Spatial and Temporal Heterogeneity Creates a “ BrownTide ” in Root Phenology and Nutrition

Spatial and temporal heterogeneity in plant phenology and nutrition benefits herbivores by prolonging the period in which they can forage on nutritious plants. Landscape heterogeneity can therefore enhance population performance of herbivores and may be a critically important feature of their habitat. The benefits of resource heterogeneity over space and time should extend not only to large herbivores using above-ground vegetation but also to omnivores that utilize below-ground resources. We used generalized linear models to evaluate whether spatial heterogeneity influenced temporal variation in the crude protein content of alpine sweetvetch (Hedysarum alpinum) roots in west-central Alberta, Canada, thereby potentially offering nutritional benefits to grizzly bears (Ursus arctos). We demonstrated that temporal patterns in the crude protein content of alpine sweetvetch roots were influenced by spatial heterogeneity in annual growing season temperatures and soil moisture and nutrients. Spatial heterogeneity and asynchrony in the protein content of alpine sweetvetch roots likely benefit grizzly bears by prolonging the period they can forage on high quality resources. Therefore, we have presented evidence of what we termed a “brown wave” or “brown tide” in the phenology and nutrition of a below-ground plant resource, which is analogous to the previously described “green wave” in above-ground resources.


Introduction
A central challenge in ecology lies in understanding the ways in which spatiotemporal variation in resources shape the biotic environment.While consumable above-ground net primary production limits the ability of a landscape to support herbivores [1,2], this limitation may be altered by the timing and spatial patterns of plant growth.This is especially true for spatial heterogeneity that causes temporal variation in plant phenology (e.g., budburst, flowering, and fruiting), because it provides nutritional benefits to herbivores by prolonging the period they can forage on highly nutritious immature plants [3].Thus, variable landscapes increase the ability of herbivores to access highly nutritious plants by reducing spatial and temporal autocorrelation in plant growth [4].Without the influence of spatial heterogeneity on plant phenology, plant growth would be synchronized, thereby reducing the time period during which animals could consume the most nutritious plants [5].Access to landscape heterogeneity can therefore enhance population performance of herbivores and may be a critically important feature of their habitat [4,6,7].
Many relationships between time, plant phenology, and animal foraging have been established.For example, during springtime when green vegetation is immature, plants have higher protein and lower fibre content than during the summer when plants are maturing [8,9].Migratory geese have been shown to follow this "green wave" of nutritious spring growth as they travel from temperate areas to their Arctic breeding grounds [10][11][12].This "surfing" of the green wave across space and time is also common to other migratory animals including insectivorous passerine birds [13] and large herbivores such as wildebeest (Connochaetes taurinus) [14].
Regional spatial heterogeneity in phenology can also be important to animal populations, because variation in a wide range of abiotic factors (e.g., soil moisture and nutrients, sunlight, temperature, precipitation, topography, snow accumulation) can cause plants at different locations to initiate growth at different times.For instance, in the Serengeti plains of Africa, the population viability of Thompson's gazelle (Eudorcas thomsonii) was directly influenced by access to patches of grassland that varied in phenology as a result of spatial heterogeneity in topography, soils, and rainfall [6].Likewise, habitats with greater spatial heterogeneity in the normalized difference vegetation index may have higher ungulate carrying capacities [15].Heterogeneity in elevation has also been linked to the population growth of ungulates in both North America and Europe [16], and many animals shift their distribution along altitudinal gradients in response to plant phenology.For example, several deer species, including roe deer (Capreolus capreolus) [17], red deer (Cervus elaphus) [18], sika deer (Cervus Nippon) [19], and reindeer (Rangifer tarandus) [20], migrate along elevation gradients to access newly emergent, high quality forage.
Studies of spatial and temporal herbivory patterns have necessitated that researchers focus primarily on aboveground plant growth; roots are also an important food resource for many animals including omnivorous species [21][22][23][24][25]. Yet, despite the importance of root resources to several wildlife species, relatively little is known about root phenology and nutrition [26] despite that root production accounts for 50-90% of the primary production in temperate vegetation [27,28].In west-central Alberta, Canada, roots of alpine sweetvetch (Hedysarum alpinum) are an important seasonal food for grizzly bears (Ursus arctos) [25].The nutritional quality of alpine sweetvetch roots vary temporally dependent on plant phenology; the protein content and digestibility of sweetvetch roots is highest preflowering and postseed development [29,30].This pattern is common for many short-stature perennial plants, where above-ground spring vegetation is dependent on resources that have been translocated to roots the previous fall and winter [31].Grizzly bears respond to sweetvetch phenology and nutrition by typically consuming roots during the spring and fall when they are most nutritious [25,29,32].Similar to herbivores, grizzly bears move to higher elevations and north facing slopes as spring progresses to access younger, more nutritious plants including alpine sweetvetch [29].Considering the "green wave" phenomenon observed in herbivores, the phenological and nutritional patterns associated with alpine sweetvetch root likely provide grizzly bears with a "wave" or (given the ebb and flow in nutritional characteristics) a "tide" in below-ground resources.Similar to a "green wave," this underground or "brown wave/tide" likely prolongs the availability of nutritional roots to grizzly bears.
The objective of this paper was to test the hypothesis that spatial heterogeneity in local abiotic and biotic factors influence temporal variation in the nutritional quality of alpine sweetvetch roots, thereby potentially offering nutritional benefits to species that use the roots, such as grizzly bears.In other words, our purpose was to evaluate whether a "brown tide" in the nutritional quality of alpine sweetvetch roots exited within our study region of west-central Alberta, Canada.In order to do this, we first describe general temtemporal poral trends in the crude protein content of alpine sweetvetch roots by plant phenophase.Then we evaluate the importance of spatial heterogeneity in crude protein content of roots using environmental variables in a geographic information system (GIS).Specifically, we hypothesized that the following five factors influence spatial and temporal patterns in crude protein content in alpine sweetvetch roots: (1) time (Julian day and year), (2) temperature (i.e., growing degree days which are influenced by elevation), (3) soil moisture and nutrients, (4) solar radiation (including slope and aspect), and (5) landcover category.We predict that, as well as observing general temporal trends and an elevation lag in phenology related to growing season temperature (i.e., interaction between time and growing degree days), crude protein content in alpine sweetvetch roots will be affected by local site factors including soils, solar radiation, and landcover category.Finally, we mapped spatial-temporal predictions of alpine sweetvetch root crude protein content over a 5,435 km 2 spatial scale to examine "brown tide" patterns in our study area.

Study Area.
The study area (Figure 1) is located on the eastern slopes of the Rocky Mountains and Foothills of west-central Alberta, Canada (53 • 15 N, 117 • 30 W).Average summer temperatures range from 11.9 • C in montane regions to 9.4 • C in the subalpine, while average winter temperatures are −7.8 • C and −8.9 • C, respectively [33].Annual average rainfall in the lower elevation montane area is 464 mm versus 568 mm at the higher elevation subalpine area [33].The foothills contain a diverse array of habitats, including black spruce (Picea mariana) and tamarack (Larix laricina) bogs; lodgepole pine (Pinus contorta) stands; open marshes; riparian areas; and mixed forests composed of lodgepole pine, aspen (Populous tremuloides), and white spruce (Picea glauca).Subalpine forests consist of spruce (Picea engelmannii × glauca), subalpine fir (Abies lasiocarpa), and lodgepole pine, with alpine meadows and rock outcrops occurring at higher elevations.

Sample Collections and Laboratory
Analyses.Alpine sweetvetch taproots were collected over three years (2008)(2009)(2010) from fourteen sites (Table 1).Five of the sites (CadMix, Cardinal, FicCon, FicMix, and Prospect) were located along a 90 km transect across known grizzly bear habitat.The transect was designed to observe phenological changes and growing season attributes across an elevation gradient (lower foothills, upper foothills, and subalpine habitats) by using digital time-lapse cameras as reported in Bater et al. [34].Alpine sweetvetch roots collected at camera locations were sampled just outside of observable camera images and were selected to best represent the phenology of the site.In order to capture a wider range of environmental conditions, we also collected alpine sweetvetch roots at sites in the montane and foothills subregions outside of the transect area (Drinnan, Folding, Greg, Camp, Kinky, and Watson)   and at a grizzly bear telemetry locations with evidence of sweetvetch digging (sites 1049 and 1117).All samples were frozen following collection.Nutritional analyses of crude protein content were performed in the Proximate Laboratories of the Department of Agriculture, Food, and Nutritional Sciences at the University of Alberta.Crude protein (N × 6.25) was determined via combustion analysis using a LECO TruSpec N/C Analyzer (Leco Corporation, St. Joseph, Michigan, USA).Moisture was determined by oven drying at 110 • C (a small number of samples were determined by oven drying at 105 • C).All estimates were corrected to a 100% dry matter basis.

Phenological Changes.
We grouped 112 of 117 total root samples (we could not identify the phenological stage of five samples) into six phenophases (preleaf, leaf, bud, flower, seed, dormant).For each phenological stage, mean and standard error of crude protein was estimated and an ANOVA used to test for significant differences between groups (phenophases).A histogram plot of residuals and a Shapiro-Wilk test was used to assess normality, while heterogeneity of variances was assessed by examining a residual plot and a Bartlett test.Pair-wise one-way t-tests with adjustments (Holm's) for multiple inferences were used to test for significant differences between phenophases.
Because grizzly bears in the study area often have home ranges either primarily in (1) subalpine and alpine habitats (elevation >1,700 m) or (2) foothills and montane habitats (elevation <1,700 m) [25], we also describe the range of dates of each phenophase according to these broad elevation classes.

Spatial and Temporal Heterogeneity of Crude Protein
Content in Alpine Sweetvetch Roots.To examine how spatial heterogeneity influences temporal variation in the crude protein content of alpine sweetvetch roots, we evaluated support for a series of a priori candidate models that describe both temporal and spatial heterogeneity of root crude protein content (n = 117).Spatial factors hypothesized to influence crude protein content in roots included soil moisture and nutrients, temperature, solar radiation, and landcover category (Table 2).
Soil moisture and nutrients were represented by both geological formation (GEO) [35] representing regional differences in parent material and for local soil conditions by the steady state terrain wetness index called the compound topographic index (CTI) derived from a 30 m digital elevation model (DEM) and the CTI.AML ArcInfo script from Evans [36].CTI has been positively related to soil moisture, horizon depth, percent silt, organic matter, and phosphorous [37,38] and was previously used for modelling the distribution of several grizzly bear plant foods in the study region, including alpine sweetvetch [38,39].Annual growing degree days (GDD) were used as a measure of heat accumulation (integral of temperature), which in the study area was controlled primarily by elevation (short, cool growing seasons or low GDD for the alpine and longer, warmer growing seasons or a high GDD in the montane and lower foothills).We used annual GDD with base 5 • C [40], since that is considered the minimum threshold for general plant growth and forage in Alberta [41].Global solar radiation (sum of shortwave and diffuse radiation) was again measured for each site using a 30 m DEM and the SHORTWAVC.AML and DIFFUSE.AML ArcInfo scripts from Zimmermann [42].Influence of landcover (e.g., plant community, competition for soil and light resources, influence on soil chemistry, associated organisms) were based on a remote sensing classification of landcover and a landscape canopy closure model from McDermid et al. [43].We calculated a correlation matrix (Pearson's) to ensure that model covariates were not highly correlated (r < 0.6).
Candidate models included univariate models, simple multivariate models without interactions, and multivariate models with interaction terms believed to be relevant to nutritional dynamics of roots.Candidate models were ranked for support using the Akaike information criterion with small sample size correction (AIC c ) [44].Because samples were collected over multiple years, year of collection was used to test for annual variation in nutritional quality.We used the program R [45] to fit mixed-effects generalized linear models (GLMMs) with site used as a random effect to account for multiple samples collected at a site.However, we found no difference between sites (all intercepts were statistically equivalent), which was supported by lower AIC c ranks (AIC c w i >2.0) than models created using fixed-effects models (GLMs).Site level variation was likely explained instead by site spatial and temporal variables; therefore, we only report here on fixed-effects models (GLMs) based on AIC c scores [46].The most supported model was checked for normality by examining a histogram of residuals followed by a Shapiro-Wilk test.All statistical analyses were conducted in the program R [45].
To interpret spatio-temporal patterns in the nutritional quality of alpine sweetvetch roots, we used our most supported model to predict monthly spatial patterns of crude protein content in roots from 07 May to 07 October.Model predictions were combined with an existing presence/absence model of alpine sweetvetch from Nielsen et al. [39] to limit predictions of protein content to areas where the species was predicted to occur.We did not predict protein content outside of 7 May and 7 October in order to avoid errors associated with extrapolating beyond dates that roots were collected.3).Observed phenology in the mountains lagged behind the foothills by 29 days for the preleaf stage, 7 days for the leaf stage, 26 days for the bud stage, 14 days for the flower stage, and 6 days for the seed stage.No lag, however, was observed for the onset of the dormant stage.While there was a lag in the onset of phenological stages between foothills and mountain habitats, the last observations for a particular stage in the mountains and foothills were often close to one another in dates, with the exception of plants in the seed stage lasting up to 17 days later in foothills habitats.This suggests that, despite general lags in the phenological stages of alpine sweetvetch between the two ecosystems, there was local variation potentially due to site level effects.

Spatial and Temporal Heterogeneity of Crude Protein
Content in Alpine Sweetvetch Roots.Temporal patterns in the crude protein content of alpine sweetvetch roots were influenced by spatial heterogeneity in annual growing season temperatures and soils with the most supported a priori model including factors for Julian day (nonlinear quadratic response), soil moisture and nutrients (CTI), annual growing degree days (GDD base 5 • C), and a three-way interaction between Julian day, CTI and GDD (Tables 4 and 5).This suggested that temporal patterns in crude protein content of roots were dependent not only on elevation lags (i.e., temperature) but also, local soil conditions.Specifically, during springtime roots were predicted to have higher protein content in areas with both low soil moisture and annual growing season temperatures as well as areas with high soil moisture and growing season temperatures.Crude protein content in these sites during the spring (21 May) was predicted to be as high as approximately 21%.In contrast, roots growing in areas where annual growing season temperatures were high and soil moisture was low were predicted to have lower protein values (∼16%).Areas of high soil moisture and low growing season temperatures were predicted to have very low crude protein content at approximately 4%.As spring progressed, spatial patterns in crude protein content of alpine sweetvetch roots were attenuated and equalized across the region by approximately 21 July with increases in crude protein content during late summer and fall following the mirrored response of springtime trends.Thus, during the late summer and fall, areas having high soil moisture and low growing season temperatures as well as low soil moisture and high growing season temperatures were predicted to have higher crude protein content of roots.By 21 September, crude protein content of roots was predicted to be as high as 21% for optimal sites and as low as 10% for areas of high moisture and growing season temperatures.
Figure 3 illustrates spatial-temporal predictions of crude protein content for each biweekly period from 7 May to 7 October.A lag in the decline of crude protein at high  elevations during the first half of the season is apparent.In contrast, during late summer and fall an increase in crude protein content occurred first at lower elevations followed by increases at higher elevations.On 7 May, roots with high protein content (>17% crude protein) were widely predicted for both high and low elevations.By 7 June, protein levels decreased to moderate levels (15.5-17%) for most low elevation sites, while many higher elevation mountain slopes and valleys retained high protein levels.By 7 July, the protein content of roots at most low elevations decreased to 14-15.5% crude protein content, while higher protein content root patches could still be found elevated at higher elevations.On 7 August, protein levels began to rise in low elevation river valley patches, and this pattern continued through to 7 October.During this time, the crude protein content of roots at high elevations lagged behind those at low elevations and along the bottoms of river valleys.By 7 September, low elevation valleys had moderate protein levels, while high elevation sites were still in a low protein state.By 7 October, protein levels of roots were predicted to be close to early spring (7 May) levels.

Discussion
We found that alpine sweetvetch roots had moderate protein content, which was expected given that alpine sweetvetch is a nitrogen fixing legume [47,48].Similar to Hamer and Herrero [29], crude protein content of alpine sweetvetch roots were highest early in the year, declined as phenology advanced to the flowering stage, and increased again to spring levels when dormant.The autumn rebound in crude protein content of roots suggests that a lag in phenology is the main driver of crude protein content in alpine sweetvetch roots.Despite the similarity between our results and Hamer  and Herrero [29], alpine sweetvetch roots in west-central Alberta had higher crude protein content, which may be due to differences in laboratory methods (Hamer and Herrero's [29] samples were air-dried at room temperature and were not dry matter corrected).A similar phenological pattern in crude protein content was found for yellow sweetvetch (Hedysarum sulphurescens) roots on the eastern slopes of Banff National Park and the Flathead River valley [29,30], although crude protein content of yellow sweetvetch roots was noticeably lower than that of alpine sweetvetch.Hamer and Herrero [29] also examined acid detergent fibre content of alpine sweetvetch roots as it related to phenology and found that fibre was inversely related to protein.While total dietary fibre is now considered the most appropriate measure of digestibility for bears [49], we expect that the relationship between fibre, protein, and phenology would remain the same.
We expected spatial heterogeneity in growing season temperature along an elevation gradient to influence the crude protein content of alpine sweetvetch roots.This is because high elevation sites experience later spring greenup and earlier fall frosts than low elevation sites, which results ISRN Ecology in regional differences in plant phenology and gradients of forage quality [50].Unlike Hamer and Herrero [29], a slope/aspect relationship based on our measures of solar radiation was not supported.
We found not only spatial-temporal lags due to elevation gradients but also local effects on root quality due to variations in soil moisture and nutrients.The effect of soils moisture and nutrients may be due to their influence on the timing of phenology and/or the overall protein content of alpine sweetvetch roots.For instance, the lower protein content of roots in dry low-elevation areas may be due to advanced phenology within the preleaf stage, since spring growth in wet areas may be delayed from excessive water caused by snow melt and/or precipitation.Soil moisture and nutrient composition is also an important factor for nitrogen fixation in legumes such as alpine sweetvetch, where soil moisture is often a limiting factor [51].This may be the case in wet alpine regions, where spring conditions leave those areas wet or saturated with water due to high snow packs.However, we expect that the extremely low protein value predicted for some springtime wet high-elevation root habitats may be an artefact of sampling, as fewer roots were collected in the subalpine at this time (in fact, many subalpine areas in our study region would likely be snow covered).As well, it is possible that roots with higher spring protein content growing in both dry high-elevation areas and wetter low-elevation areas may experience a combination of temperature and soil moisture conditions that are suitable to nitrogen fixation; our measure (index) of soil moisture has been correlated with soil nutrients including phosphorous [37], which is often a limiting factor in nitrogen fixation [52].Further examination of the effects of moisture and nutrients on alpine sweetvetch root quality is therefore needed.
Tracking the phenology and nutritional quality of resources over heterogeneous landscapes is an important mechanism by which animals increase diet quality, individual health, and population performance.In this paper, we provide evidence of a "brown wave" or "brown tide" in the nutritional quality alpine sweetvetch root, which suggests that the benefits of resource heterogeneity over space and time extend not only to large herbivores and thus aboveground vegetation but also to omnivores such as grizzly bears and their below-ground resources.While we demonstrate "brown tide" patterns in our study region, we have not explicitly linked this with grizzly bear foraging; however, given previous observations suggesting that grizzly bears track resource phenology along an elevation gradient [29], we expect that such a relationship does, in fact, exist.For this reason, we recommend that future research test the hypothesis that grizzly bears track the phenology and nutritional quality of alpine sweetvetch root.
This work lends to the argument that grizzly bear habitat models should more explicitly consider spatial and temporal heterogeneity in food resources [25,39].In fact, modeling has successfully predicted movement patterns for other wildlife species when animals were assumed to track variations in plant phenology within heterogeneous landscapes [14,53].While we focused on crude protein as a measure of nutritional quality, this approach can also be extended to other nutritional (e.g., carbohydrate, fibre) estimates.As well, energetic estimates (e.g., digestible protein, digestible dry matter, and net digestible energy) could also be used but would require estimates of sweetvetch abundance, density, or biomass which have not yet been developed.Such an approach may be useful in understanding the effects of environmental change on wildlife when combined with factors that regulate populations from the top down [39].

Figure 1 :
Figure 1: Study area depicting towns, alpine sweetvetch (Hedysarum alpinum) sampling sites, montane, lower foothills, upper foothills, subalpine, and alpine subregions.The study region outline indicates the extent of the area mapped for crude protein content predictions.

Figure 2 :
Figure 2: Average crude protein content (% dry matter basis) of Hedysarum alpinum roots sampled in west-central Alberta, Canada, in different phenological stages.Error bars are 1 standard error above and below the mean.Lower case letters indicate phenophases with significantly different levels of crude protein (i.e., preleaf and flower, preleaf and seed, seed and dormant, flower and dormant).

Figure 3 :
Figure 3: Nutritional landscape maps demonstrating "brown tide" patterns in west-central Alberta, Canada, by predicting the percent crude protein content of alpine sweetvetch (Hedysarum alpinum) roots from 7 May to 7 October.

Table 2 :
Hypothesized factors, variable names, and abbreviations of variables used in candidate models.

Table 3 :
Range of dates and average percent crude protein content (by dry weight) of alpine sweetvetch (Hedysarum alpinum) roots (n = 112) by phenophase for both mountain (elevation ≥1700 m) and foothills (elevation <1700 m) environments.Phenophase observations and root collection occurred from 2008-2010.There was only one observation made of the leaf phenological stage in mountain habitats.§Unique lower case lettering indicates significant differences in crude protein content.

Table 4 :
Hypothesized candidate models, model structure, Akaike Information Criterion with small sample size correction (AIC c ), difference in AIC c score from highest ranked candidate model (ΔAIC c ), and Akaike weight (w i ) of the top 10 candidate models used to predict crude protein content of alpine sweetvetch (Hedysarum alpinum) in west-central Alberta, Canada.The final model (in bold) was selected based on the highest Akaike weight (w i ) and ΔAIC c > 2.

Table 5 :
Estimated coefficients, standard error (SE), t-value, and significance levels (P) of model variables used to predict crude protein content of alpine sweetvetch (Hedysarum alpinum) roots in west-central Alberta, Canada.