BACKGROUND
ecological linkages between marine and terrestrial communities are important processes structuring coastal ecosystems.
one of the most geographically broad linkages that has received increased attention in recent decades is the role of migrating salmon  to coastal watersheds  <cit> .
in addition to the contribution of salmon carcasses to primary productivity in lakes and streams  <cit> , salmon also comprise a substantial component of the diet in a diversity of marine and terrestrial taxa including pinnipeds and bears as well as an extensive history of use by first nations  <cit> .
despite the increased attention and numerous publications, there remains ambiguity and large knowledge gaps regarding the relationship between salmon and receiving ecosystems.
during an investigation from  <dig> to  <dig> at bag harbour, an intact salmon watershed on the southern reaches of haida gwaii, western canada, it became evident that foraging black bears  transfer large numbers of salmon from the stream into the riparian zone and uneaten remnants of the carcasses accumulate on the forest floor throughout the spawning period  <cit> .
these carcass remnants were widely distributed but the majority occurred within 50 m of the stream with densities averaging  <dig>  carcasses/m <dig> adjacent to the spawning zones.
multiple scavengers including flocks of corvids and dipteran larvae utilized and further dispersed these nutrients  <cit> .
subsequent investigations in other watersheds of coastal british columbia confirmed widespread transfer of salmon nutrients to the riparian zone by bears  <cit>  as well as a utilization of these nutrients by riparian vegetation  <cit> , insects  <cit>  and songbirds  <cit> .
studies in alaskan watersheds have identified both flooding and hyporheic movement as well as predator activity as a source for uploading of salmon nutrients  <cit> .
the ecological consequences of salmon-derived nutrients to riparian vegetation include increased plant growth rates  <cit> , enrichment of foliar nitrogen  <cit>  and altered plant community diversity  <cit> .
because nitrogen is often a limiting nutrient in temperate forests of the pacific northwest  <cit> , it is possible that short and/or long term differences in salmon returns to spawning rivers would co-vary with primary productivity in the riparian zone and if so, annual growth rings of these old-growth trees might contain historical information on relative abundance of salmon.
the recent identification of nitrogen isotopes in small quantities of wood from individual tree rings  <cit>  offers an additional direct assessment of the isotopic enrichment of marine-derived nutrients over time  <cit> .
we extend the investigation at bag harbour and examine the potential influences of salmon nutrients on yearly growth in sitka spruce, a dominant riparian conifer, through analyses of tree rings.
recent studies suggest that there is limited potential of tree rings for interpreting historical trends in salmon abundance due to the complexity of edaphic and climatic processes influencing tree growth, the extent of nitrogen uptake and translocation among rings as well as stable isotope fractionation  <cit> .
however, the distribution of carcasses has not been quantified in previous studies of riparian growth nor has a fine-scale  examination of both tree growth and δ15n in relation to salmon been previously undertaken.
recently, drake and naiman  <cit>  provide comparisons across multiple watersheds yet control sites were for  <dig> to 25 km distance from salmon sites with no exploration of spatial or temporal within-stand relationships to salmon abundance.
we established a sampling grid adjacent to bag harbour stream between the estuary and the headwater lake that encompassed the previously identified distribution of salmon carcasses  <cit> .
thirty-six of the largest sitka spruce  located over a range of carcass densities throughout this grid were cored and subsequently quantified for yearly radial growth.
yearly abundance of salmon in this watershed has been recorded  from  <dig> through to the present and varies by more than an order of magnitude among years  and consequently, our results and analyses are restricted to this time period.
in this study, our objectives are to: 1) characterize the spatial pattern of annual growth  and nitrogen isotopic signatures  in trees across replicated gradients of salmon carcass densities, 2) identify the temporal signatures  of radial growth and nitrogen isotope signatures in relation to known yearly salmon escapement and 3) assess the predictors of both tree ring growth  and δ15n levels using an information theoretics approach that relies on general linear mixed models .
potential predictors of growth and δ15n values in annual tree rings included environmental , spatial  and total nitrogen as variables.
our replicated transects perpendicular to the stream channel across sharp gradients in carcasses and throughout the length of the stream along a broad range of salmon spawning densities allowed multiple paired comparisons of growth and nitrogen isotope signatures of trees.
although not directly exploring the use of tree ring data as a proxy for historical salmon abundance in this study, our results are useful for evaluating tree ring growth and δ15n values as potential proxies.
methods
study area
bag harbour watershed , located in southeastern haida gwaii, british columbia, canada , is largely composed of old-growth forest, mainly sitka spruce , western hemlock  and western redcedar .
alder  are uncommon in the riparian zone of this watershed.
one salmon-bearing stream, which ranges from  <dig> to 20 m in width and usually less than  <dig>  m depth, is fed by a small lake, located ~ <dig>  km upstream from the estuary.
the majority of salmon spawn in gravel beds located between 100 m and 800 m upstream, with limited spawning occurring above 800 m. chum  make up the majority  of spawning salmon, although smaller numbers of coho  and pink  salmon are also present .
black bears  are the dominant vector of salmon from the stream to the forest floor and a diversity of species utilize the abandoned carcasses  <cit> .
field methods
in fall  <dig>  a riparian sampling grid was established across a previously documented gradient of high to zero salmon carcass densities and extending in a 100 m band along the stream from the estuary to the lake  <cit>  .
at each of  <dig> sites along the length of the stream, separated by approximately 100 m, the largest sitka spruce closest to the stream  was chosen for coring.
at  <dig> of these sites, one or two additional of the largest spruce were cored along a perpendicular transect into the forest, the first near the outer edge of the distribution of salmon carcasses  and the second, when present, approximately twice the distance from the stream and typically comprising the most distant spruce identified from the stream.
spruce were less abundant above the spawning reaches of the stream and when present were close to the stream channel.
we used the presence or absence of salmon carcasses as a general proxy for the spatial input of salmon nutrients into the forest.
we refer to trees located at sites without evidence of carcasses  as ’control or reference’ trees.
these trees reflect low to trace salmon-nutrient loading.
tree core sampling and growth index
each spruce was cored at breast height using a  <dig>  cm × 40 cm increment borer, measured for diameter at breast height  and distance from the stream and distance upstream from the estuary.
presence or absence of salmon carcass remnants around the base of the tree  was recorded.
cores were oven dried at 60oc for at least four weeks, mounted on plastic sheaths, sanded with fine-grade sandpaper and digitally scanned.
images of cores were used to individually measure tree ring widths to the nearest  <dig>  mm using coorecorder  with problematic rings visually assessed under a dissecting microscope.
tree ring series were crosschecked for accuracy using default settings in cofecha  <cit> .
individual tree ring series were then fitted to a single smoothing spline with a 50% frequency cutoff of 100 years, a procedure that removes age-related growth trends in individual series, using the program arstan  <cit> .
the use of a 100-year spline as opposed to the default 32-year spline is that mid to high frequency variations are retained in trends less than 100 years and long-term growth related trends are removed .
we refer to this metric as the arstan-modified growth index or growth index in subsequent analyses.
stable isotopes
stable isotope samples were collected by excising individual tree rings from dried cores using a scalpel blade.
when rings were very small  and there was insufficient mass for analyses, samples from two or more sequential years were combined.
we analyzed a relatively complete ring series  on  <dig> of the trees comprising six paired comparisons  at different positions  along the length of the stream channel  and two additional control trees  at  <dig> and 700 m above the upper reaches of the spawning zones.
in addition, we analysed  <dig> rings from each of  <dig> additional trees distributed from the estuary to the headwater lake  comprising two time periods .
these two periods were chosen as part of a coastwide assessment of nitrogen isotope signatures in conifer tree rings and characterize the earliest period  where salmon abundance is known for most watersheds and a recent period  that largely excludes use of more recent sapwood rings.
tree ring samples were ground for five minutes using a wig-l-bug amalgamator.
samples  were packaged in tin capsules and analyzed for total nitrogen and isotopes of nitrogen at the university of california davis stable isotope facility using a pdz europa anca-gsl elemental analyzer interfaced to a pdz europa 20- <dig> isotope ratio mass spectrophotometer .
nitrogen isotope results are reported in δ notation  relative to the international standard of atmospheric n <dig>  in parts per thousand.
isotopic abundances are calculated by:
δ15n=rsample/rstandard
where r equals the ratio of 15 n/14 n stable isotopes.
reproducibility of δ15n values on replicated wood samples averaged +/-  <dig> ‰.
environmental data
temporal data included salmon abundance  data for bag harbour  modified using a three-year running average from the previous three years.
this three-year average with a one-year lag  <cit>  was used as a conservative temporal estimate of potential salmon nutrient input.
precipitation and temperature data were obtained from nearby sandspit for the majority of years, with data from tlell or alternatively masset used when sandspit data were not available .
total precipitation and mean maximum monthly temperature data were generated from september in the previous growth year to october of the growth year for any given tree ring, as these periods more closely reflect the biological year  <cit> .
generalized linear mixed models
we assessed the influence of salmon abundance and other potential predictors on tree growth using the arstan-modified growth index and δ15n isotopic enrichment separately, using an information theoretics approach.
predictors were selected due to their availability but most importantly, their possible influence on tree growth and δ15n values.
for both model sets, response variables were first examined with all potential predictor variables to determine the nature of the potential relationship.
predictor variables were also examined for collinearity.
statistical analyses were performed using spss v. <dig> .
predictors of the growth index were assessed using general linear mixed models  for the period 1948- <dig>  individual trees  were specified as random factor subjects with year as the repeated measure.
covariance was specified as ar <dig> in order to account for the autoregressive nature of the data.
growth index, the response variable, was transformed closer to normality using a box-cox transformation and salmon abundance, precipitation and temperature were log transformed.
thirty-two biologically feasible models were constructed using the predictor variables salmon abundance, precipitation, temperature and the interaction with distance upstream  and distance into the forest .
spatial predictors were only included as interaction terms due to the consequence of arstan growth index, which standardizes tree series, thus removing the main spatial effects of growth, along with age-related trends.
the null model  was also included.
models were ranked using aikaike’s information criteria for finite sample sizes .
where more than one model was ranked with an aic value of > <dig>  model averaging of the confidence set of candidate models was performed on parameter estimates and variances.
the confidence set of candidate models was based on those with a cumulative sum of >95% aicc weight.
predictors of δ15n isotopic enrichment for individual tree rings were also assessed using multiple general linear mixed models  for the period 1948- <dig> .
tree ring isotopic signatures obtained in years more recent than  <dig> were excluded due to elevated total nitrogen and variable δ15n isotopic values in the biologically active sapwood .
individual trees were specified as random factor subjects with year as the repeated measure.
covariance was specified as autoregressive  in order to account for the autoregressive nature of the data.
distance into the forest and tree dbh demonstrated an unacceptably high level of collinearity.
as there was a strong relationship between δ15n and tree dbh, we adjusted δ15n isotopic enrichment values based on a dbh correction function generated using a linear regression of dbh and the δ15n values of all trees:
logy′ij=logyij−vjlogxi−logx
where y’ij is the adjusted value of δ15n for year i, yij is the original δ15n value, vj is the regression slope, xi is the annual radial growth of year i and x the δ15n signature for the year  <dig>  we present data for both unadjusted and dbh-adjusted isotopic signatures but only use the latter for analysis using glmms.
dbh-adjusted δ15n isotopic signature, the response variable, was transformed closer to normality using a box-cox transformation, growth index was log transformed and distance upstream  was converted to a categorical value  due to a non-linear relationship with the response variable.
twenty-three biologically feasible models were constructed using the predictor variables salmon abundance , growth index, percent nitrogen, categorical distance upstream , distance into the forest  and interaction terms.
interaction terms involving spatial predictors and the growth index were not included.
the null model  was also included.
models were ranked using aikaike’s information criteria for finite sample sizes .
RESULTS
spatial trends
in the  <dig> transects between estuary and the headwater lake, trees in the carcass zone were larger  than adjacent control trees  but did not differ in average age .
the largest average annual growths  were found in trees adjacent to the major spawning reaches of the stream but average growth did not vary with respect to distance into the forest .
among the  <dig> trees with stable isotope data, δ15n signatures  ranged from - <dig>  to  <dig>  among trees, with the highest values occurring with carcasses  than in the control trees .
data were similar for 1950- <dig> .
we examined the interaction between tree size  and isotope values and found a significant positive slope  with a similar trend for both carcass and non-carcass trees although only significant in the former .
percent nitrogen was marginally but non-significantly higher for carcass trees  than control trees  .
temporal trends
we have relatively complete isotope data  for  <dig> trees.
within most of these , the unadjusted δ15n signatures varied approximately 5‰ over the 53 years with some trees  showing a depletion historically and others  exhibiting higher signatures in the past.
five of the six trees found further from the stream had lower signatures  relative to trees located closer to the stream, with a single exception , but in this case, salmon carcasses were present in both sites.
the highest signatures of any tree , reaching 12‰, occurred in the major bear-foraging area.
these signatures were ~10‰ higher than the corresponding control tree  50 m into the riparian zone but this difference was greatly reduced historically as demonstrated by the increased isotopic signatures in the control tree during earlier periods.
when we corrected isotope data for size of tree in each year, the temporal trends suggest the control tree increased by 10‰ achieving similar isotopic signatures to bh <dig> during the earliest sampling period in  <dig>  smaller increases  of size-adjusted values occurred historically in three additional carcass trees  and in three control trees .
the lowest unadjusted and adjusted signatures occurred in the two control trees  above the upper extent of the spawning zone and both of these showed a 3-5‰ increase back to  <dig>  we evaluated the relationship with salmon abundance which fluctuated from ~ <dig> in the late 1990s to ~ <dig>  in the late 1940′s.
five of the  <dig> trees showed a non-significant negative association while nine showed a positive relationship between size-adjusted δ15n and salmon abundance of which four  were statistically significant .
temporal trends for glmm variables
to provide a visual context for variables used in the glmms, we characterize the mean values for standardized growth index, δ15n, percent nitrogen, salmon abundance, temperature and total precipitation .
standardized growth was highly variable, with clear peaks in mean growth over the study period .
adjusted δ15n was also highly variable with decline in mean values towards the present .
percent nitrogen showed low variance among years apart from the recent decades which showed a striking increase in percentages concordant with the outer active sapwood .
salmon escapement fluctuates over time with two large pulses occurring in the late 1940′s and early 1960′s  and a general decline with time from the 1960′s to the recent.
abundances in bag harbour prior to the mid-1940′s are not currently known but are anticipated to have been higher  <cit> .
the environmental variables, mean temperature and total precipitation, were calculated for the biological year.
on average, temperatures generally increase with time  and bag harbour and the surrounding area receives between  <dig> and 2100 mm of precipitation  for the biological year .
predictors of growth using glmms
due to a lack of discrimination for the highest ranked model using aicc, the top two ranked models that best predicted the growth index were averaged and the reported cumulative weight  generated using both model terms .
leading models include the temporally dynamic predictor variable salmon abundance, temperature and the interaction terms salmon*temperature and salmon* distance into the forest.
summing parameter weights across the candidate set of best-supported models demonstrated the high relative importance of salmon, temperature and the interaction between salmon and temperature; spatial interaction terms with salmon were lesser contributors to the top model set .
although parameters estimates were small, particularly for the influence of temperature , we note that the growth index is a heavily modified metric and both means and variances are standardized between tree series.
salmon, the interaction between salmon abundance and temperature had positive parameter estimates, as did the interaction between salmon abundance and the distance into the forest.
the interactions between  <cit>  salmon abundance and temperature and  <cit>  salmon abundance and distance into the forest were similar, with a larger influence of salmon abundance on growth at higher temperatures and similarly, a larger effect of salmon abundance with increasing distance into the forest .
model predictors are annual salmon abundance , mean maximum temperature  for the biological year , distance from stream , distance upstream , total precipitation for the biological year  and derived annual tree ring width .
models were ranked using aicc values.
for predictors of the growth index, the top two models were averaged.
* interaction effect.
* interaction effect.
*interaction effect.
predictors of δ15n using glmms
in contrast to the predictors of growth, the two spatial variables and their interaction best predicted δ15n isotopic enrichment in individual tree rings .
summing parameter weights across the top three best-supported models  demonstrated the high importance of the spatial predictors of distance into the forest, distance upstream and their interaction with far lower support for the salmon and growth .
distance into the forest, distance upstream and their interaction negatively all had negative parameter estimates .
the interaction between distance upstream and into the forest is interpreted as having a lesser importance of distance into the forest with distance upstream .
the two spatial predictors  were significant .
although four of the  <dig> trees yielded a significant positive covariation between isotopic signatures and salmon abundance, this effect is accentuated on trees occurring at low carcass densities.
notably, our model-based inference did not identify salmon abundance as an important predictor.
we also ran the glmm using the full temporal complement  in order to include the last decade of high nitrogen levels in the sapwood and found the identical top-ranked models relative to the original model-based selection .
* interaction effect.
discussion and
CONCLUSIONS
the spatial patterns of carcass distribution, with declining carcasses with distances upstream and into the forest, delineated our study area and provided the spatial framework for our analysis.
of the multiple predators and scavengers observed at bag harbour, haida gwaii, during the five separate salmon spawning runs , black bears were the dominant species responsible for the movement of salmon from the stream into riparian zones  <cit> .
in  <dig>  eight bears uploaded  <dig> salmon  to the riparian zone  <cit>  and consumed approximately one-half of each carcass leaving 330 g of carcass remnants/m <dig> or 10 g n/m <dig> <cit> .
these estimates are similar to those documented in alaska, where gende et al.
<cit>  estimated that bears contributed 5 mg n/m <dig> to 10 m strip adjacent to a salmon stream.
in addition, brown bears contributed approximately 5 mg n/m <dig> in a 500 m region adjacent to an alaskan stream, of which 97% was added via urine  <cit> .
in addition to nitrogen provided by salmon carcasses, this input was estimated to account for 10-25% of the forest nitrogen budget  <cit> .
despite the quantitative estimates for quantity and spatial extent of carcasses, we note that abundance and spatial extent of the carcass distribution at bag harbour probably comprises only a fraction of historical levels; in the late 1940s and early 1960′s, the size of the spawning run was approximately seven times the abundance in  <dig>  further, the bear trails that occurred beyond the high carcass zone and above the upper reaches of the spawning area  <cit>  broaden the potential contribution of urine-derived nitrogen at much greater distances than examined.
adjacent to the main spawning area, where salmon carcass densities were highest, we also found high annual radial tree growth.
however, high growth was also found further into the riparian zones where carcasses were much less abundant but bear trails and assumed urine deposition were prevalent.
other studies on the annual growth of alaska sitka spruce  that grew adjacent to salmon streams found that these trees grew at more than three times the rate of trees at reference sites above waterfalls of unknown distance from the salmon spawning areas  <cit> .
this elevated productivity of trees led to further examinations of the temporal relationship between salmon abundance and tree growth at the watershed level.
after an initial investigation  <cit> , an examination of nine sites detected positive relationships between tree ring growth and salmon abundance at five locations, which allowed for the historical reconstruction of salmon runs at five locations  <cit> .
similar to drake and naiman  <cit> , our results generated using indexed growth clearly suggest that salmon abundance has a positive influence on tree growth.
however, instead of using reference sites located at a distance  , we used within-watershed trees across a sharp ecological gradient from high to very low or zero salmon-nutrient loading.
relying on model-based inference, we found that the abundance of salmon and temperature are major predictors of the tree growth index, with all other predictors consisting of salmon abundance and its interaction with temperature and distance into the forest.
interestingly, the interaction between salmon and distance into the forest was positive.
initially considered counterintuitive, this positive relationship may be the result of elevated nitrogen levels in areas with higher salmon loading where, as a consequence, trees may not always be nitrogen-limited.
their growth, although high, may not correlate as closely with salmon as trees that are more nitrogen-limited and experience lower levels of salmon-derived nitrogen.
considered initially as controls, the growth response of trees at larger distances from the spawning gravel beds  to salmon abundance provides evidence that the spatial effects of salmon on tree growth extend at least 90 m into the riparian zone and likely further.
temperature and the interaction between temperature and salmon abundance were also highly ranked as predictors, suggesting that differencing techniques for the removal of environmental signals from a chronology  may, at the minimum, remove some of the influences of salmon and subsequently alter reconstructions.
patterns of δ15n enrichment tracked what is known of the spatial transfer of salmon into the forest, with higher levels of enrichment found adjacent to the spawning area and declining with distance into the forest and above the upper reaches of the spawning zone.
not surprisingly, the use of glmm confirmed our visual assessment; spatial variables were found to be the strongest predictors of δ15n at bag harbour.
this finding is similar, although on a much finer spatial resolution, to hocking and reimchen  <cit> , who found that chum salmon abundance was a leading predictor of δ15n variation in plants and invertebrates at the watershed level.
similarly, hilderbrand et al.
<cit>  found that not only δ15n signatures of white spruce decline with distance from salmon-spawning streams, but that the activity of radio-collared brown bears were highly correlated with this same gradient of depletion.
the mechanism of the influence of salmon on tree growth is generally considered to be related to the input of salmon-derived nitrogen to nitrogen-limited trees.
in addition to growth being positively influenced by salmon abundance, δ15n derived from salmon and subsequently laid down in annual tree rings may provide an additional proxy for salmon abundance.
an experimental study that monitored the uptake of δ15n by western redcedar trees  found that δ15n uptake in the fall was distributed to leaves and stems in the spring, which substantiates our use of a one year lag effect  <cit>  but also furthers the idea that salmon-derived δ15n could be linked to salmon abundance, similar to growth.
studies of natural δ15n levels in tree rings have suggested that they may reflect patterns of atmospheric nitrogen deposition  and fertilization studies have demonstrated that tree ring δ15n changes may reflect the historical n regime in a forest  <cit> .
kiernan and johnson  <cit>  found relatively weak support for the idea that δ15n levels in tree rings can be linked to salmon abundance and suggested that nitrogen dynamics in the adjacent stream, soil and within-tree may have acted to obscure their ability to discriminate salmon-derived and autochthonous sources of n.
major factors that remain poorly understood, and which potentially confound the use of δ15n levels as a predictor of salmon abundance in this study and others, include denitrification, nitrogen isotope fractionation, translocation of n isotopes across tree rings  and local soil and stream nitrogen dynamics  <cit> .
in particular, denitrification, where 14 n is preferentially lost in the system and higher proportions of 15 n remain, could be particularly problematic if it varied with spatial patterns of salmon loading .
similar δ15n declines with distance away from the salmon spawning area  have been previously attributed to salmon-derived nutrients but these conclusions were strengthened by a lack of similar declines on nearby non-salmon bearing streams  <cit> .
while we have no prior expectation that denitrification varied spatially over the study site we cannot discriminate denitrification from salmon loading using δ15n if the two are similar.
the close spatial coupling between the distribution of salmon carcasses at bag harbour and the δ15n signatures of trees, particularly the steep gradients in signatures  over short distances , suggest that tree ring δ15n signatures in this watershed may be proxies for historical nutrient uploading in this watershed.
although there are a number of feasible alternative explanations, we make several possible inferences if these δ15n patterns are linked to salmon.
in the riparian zone adjacent to the current major spawning and bear feeding areas, there has been general temporal continuity of the nutrient uploading during the last five decades.
yet, in the upper reaches of the watershed, between the present spawning zone and the headwater lake, and in two trees further downstream in the low carcass zone, where there is no current salmon transfer, δ15n signatures were historically equivalent to those observed in the major feeding areas, but have become depleted by approximately 5‰ by  <dig>  suggesting a historically wider distribution of salmon-derived nutrients.
consistent with this suggestion, there was also a depletion in yearly isotopic signatures, relative to historic periods, and which are more similar to the yearly variability in trees in present high carcass densities.
there was a limited isotopic signature of the two major pulses of salmon recorded in the late 1940s and early 1960′s apart from two of the low-carcass density trees .
bears abandon a greater number of carcasses and consume a smaller proportion of each carcass when salmon escapement is relatively high  <cit>  yet it is likely that numbers of bears in a watershed are determined on average annual salmon escapement, rather than yearly fluctuations.
if so, this would decouple any immediate covariation and reduce the chances of detecting year to year fluctuations, thus limiting the application of isotopic proxies for salmon abundance to a longer time frame.
tree ring information has been used as a proxy for constructing historical salmon abundances at multiple watersheds with riparian zone stand chronologies and nearby salmon-free controls  <cit>  and at the regional level in climate-sensitive series  <cit> , both with some success.
similar techniques have been employed to examine historical seabird activity extending from the input of guano to forested nesting colonies  <cit> .
our study, although not completing salmon reconstructions and offering only a relatively short time-scale analysis, provides insight into the use of growth indexes and tree ring δ15n as proxies for historical salmon abundance.
we observed that δ15n signatures exhibit a relatively consistent decrease from large to small trees both in high and low carcass density zones allowing calibration of δ15n in cores that extend from the cambium towards the center.
without a dbh-correction, δ15n values tend to remain somewhat stable over time.
yet, because we detected a strongly positive relationship between tree dbh and δ15n and subsequently used this relationship as a calibration factor, δ15n values demonstrate a decline over time.
we also found elevated isotopic signatures in the late 1940′s and 1950′s of several trees that are currently in low carcass zones, suggesting elevated salmon uploading at that time.
this size factor has consequences for potential use and improvement of salmon abundance reconstructions , although further investigation is required.
our study, which reaffirms the importance of salmon to riparian productivity and contributes new insight for the complex, spatiotemporal influences of salmon abundance, temperature and tree location relative to salmon upon conifer growth and δ15n levels.
both of these variables provide strong spatial indicators.
the use of δ15n for temporal reconstructive purposes remains challenging.
we speculate that the appropriate choice of trees in riparian zones across ecological gradients of salmon carcass densities and bear activity and the isotopic correction for yearly changes in tree size may allow improved assessment of tree ring growth as a historical proxy for salmon abundance and distribution during previous centuries in watersheds where ancient trees still persist.
abbreviations
aicc: aikaike’s information criteria for finite sample sizes; bh: bag harbour; dbh: diameter at breast height; glmm: general linear mixed model.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
tr conceived of, carried out fieldwork and statistical analyses for the study.
cf completed tree ring measurements, data collection and glmm statistical analyses.
both authors participated in drafting the manuscript and have read and approved the final manuscript.

 ecological linkages between marine and terrestrial communities are important processes structuring coastal ecosystems. one of the most geographically broad linkages that has received increased attention in recent decades is the role of migrating salmon to coastal watersheds. in addition to the contribution of salmon carcasses to primary productivity in lakes and streams, salmon also comprise a substantial component of the diet in a diversity of marine and terrestrial taxa including pinnipeds and bears as well as an extensive history of use by first nations. despite the increased attention and numerous publications, there remains ambiguity and large knowledge gaps regarding the relationship between salmon and receiving ecosystems. during an investigation from to at bag harbour, an intact salmon watershed on the southern reaches of haida gwaii, western canada, it became evident that foraging black bears transfer large numbers of salmon from the stream into the riparian zone and uneaten remnants of the carcasses accumulate on the forest floor throughout the spawning period. these carcass remnants were widely distributed but the majority occurred within 50 m of the stream with densities averaging carcasses/m adjacent to the spawning zones. multiple scavengers including flocks of corvids and dipteran larvae utilized and further dispersed these nutrients. subsequent investigations in other watersheds of coastal british columbia confirmed widespread transfer of salmon nutrients to the riparian zone by bears as well as a utilization of these nutrients by riparian vegetation, insects and songbirds. studies in alaskan watersheds have identified both flooding and hyporheic movement as well as predator activity as a source for uploading of salmon nutrients. the ecological consequences of salmon-derived nutrients to riparian vegetation include increased plant growth rates, enrichment of foliar nitrogen and altered plant community diversity. because nitrogen is often a limiting nutrient in temperate forests of the pacific northwest, it is possible that short and/or long term differences in salmon returns to spawning rivers would co-vary with primary productivity in the riparian zone and if so, annual growth rings of these old-growth trees might contain historical information on relative abundance of salmon. the recent identification of nitrogen isotopes in small quantities of wood from individual tree rings offers an additional direct assessment of the isotopic enrichment of marine-derived nutrients over time. we extend the investigation at bag harbour and examine the potential influences of salmon nutrients on yearly growth in sitka spruce, a dominant riparian conifer, through analyses of tree rings. recent studies suggest that there is limited potential of tree rings for interpreting historical trends in salmon abundance due to the complexity of edaphic and climatic processes influencing tree growth, the extent of nitrogen uptake and translocation among rings as well as stable isotope fractionation. however, the distribution of carcasses has not been quantified in previous studies of riparian growth nor has a fine-scale examination of both tree growth and δ15n in relation to salmon been previously undertaken. recently, drake and naiman provide comparisons across multiple watersheds yet control sites were for to 25 km distance from salmon sites with no exploration of spatial or temporal within-stand relationships to salmon abundance. we established a sampling grid adjacent to bag harbour stream between the estuary and the headwater lake that encompassed the previously identified distribution of salmon carcasses. thirty-six of the largest sitka spruce located over a range of carcass densities throughout this grid were cored and subsequently quantified for yearly radial growth. yearly abundance of salmon in this watershed has been recorded from through to the present and varies by more than an order of magnitude among years and consequently, our results and analyses are restricted to this time period. in this study, our objectives are to: 1) characterize the spatial pattern of annual growth and nitrogen isotopic signatures in trees across replicated gradients of salmon carcass densities, 2) identify the temporal signatures of radial growth and nitrogen isotope signatures in relation to known yearly salmon escapement and 3) assess the predictors of both tree ring growth and δ15n levels using an information theoretics approach that relies on general linear mixed models. potential predictors of growth and δ15n values in annual tree rings included environmental, spatial and total nitrogen as variables. our replicated transects perpendicular to the stream channel across sharp gradients in carcasses and throughout the length of the stream along a broad range of salmon spawning densities allowed multiple paired comparisons of growth and nitrogen isotope signatures of trees. although not directly exploring the use of tree ring data as a proxy for historical salmon abundance in this study, our results are useful for evaluating tree ring growth and δ15n values as potential proxies. methods study area bag harbour watershed, located in southeastern haida gwaii, british columbia, canada, is largely composed of old-growth forest, mainly sitka spruce, western hemlock and western redcedar. alder are uncommon in the riparian zone of this watershed. one salmon-bearing stream, which ranges from to 20 m in width and usually less than  m depth, is fed by a small lake, located ~  km upstream from the estuary. the majority of salmon spawn in gravel beds located between 100 m and 800 m upstream, with limited spawning occurring above 800 m. chum make up the majority of spawning salmon, although smaller numbers of coho and pink salmon are also present. black bears are the dominant vector of salmon from the stream to the forest floor and a diversity of species utilize the abandoned carcasses. field methods in fall a riparian sampling grid was established across a previously documented gradient of high to zero salmon carcass densities and extending in a 100 m band along the stream from the estuary to the lake. at each of sites along the length of the stream, separated by approximately 100 m, the largest sitka spruce closest to the stream was chosen for coring. at of these sites, one or two additional of the largest spruce were cored along a perpendicular transect into the forest, the first near the outer edge of the distribution of salmon carcasses and the second, when present, approximately twice the distance from the stream and typically comprising the most distant spruce identified from the stream. spruce were less abundant above the spawning reaches of the stream and when present were close to the stream channel. we used the presence or absence of salmon carcasses as a general proxy for the spatial input of salmon nutrients into the forest. we refer to trees located at sites without evidence of carcasses as ’control or reference’ trees. these trees reflect low to trace salmon-nutrient loading. tree core sampling and growth index each spruce was cored at breast height using a  cm × 40 cm increment borer, measured for diameter at breast height and distance from the stream and distance upstream from the estuary. presence or absence of salmon carcass remnants around the base of the tree was recorded. cores were oven dried at 60oc for at least four weeks, mounted on plastic sheaths, sanded with fine-grade sandpaper and digitally scanned. images of cores were used to individually measure tree ring widths to the nearest  mm using coorecorder with problematic rings visually assessed under a dissecting microscope. tree ring series were crosschecked for accuracy using default settings in cofecha. individual tree ring series were then fitted to a single smoothing spline with a 50% frequency cutoff of 100 years, a procedure that removes age-related growth trends in individual series, using the program arstan. the use of a 100-year spline as opposed to the default 32-year spline is that mid to high frequency variations are retained in trends less than 100 years and long-term growth related trends are removed. we refer to this metric as the arstan-modified growth index or growth index in subsequent analyses. stable isotopes stable isotope samples were collected by excising individual tree rings from dried cores using a scalpel blade. when rings were very small and there was insufficient mass for analyses, samples from two or more sequential years were combined. we analyzed a relatively complete ring series on of the trees comprising six paired comparisons at different positions along the length of the stream channel and two additional control trees at and 700 m above the upper reaches of the spawning zones. in addition, we analysed rings from each of additional trees distributed from the estuary to the headwater lake comprising two time periods. these two periods were chosen as part of a coastwide assessment of nitrogen isotope signatures in conifer tree rings and characterize the earliest period where salmon abundance is known for most watersheds and a recent period that largely excludes use of more recent sapwood rings. tree ring samples were ground for five minutes using a wig-l-bug amalgamator. samples were packaged in tin capsules and analyzed for total nitrogen and isotopes of nitrogen at the university of california davis stable isotope facility using a pdz europa anca-gsl elemental analyzer interfaced to a pdz europa 20- isotope ratio mass spectrophotometer. nitrogen isotope results are reported in δ notation relative to the international standard of atmospheric n in parts per thousand. isotopic abundances are calculated by: δ15n=rsample/rstandard where r equals the ratio of 15 n/14 n stable isotopes. reproducibility of δ15n values on replicated wood samples averaged +/- ‰. environmental data temporal data included salmon abundance data for bag harbour modified using a three-year running average from the previous three years. this three-year average with a one-year lag was used as a conservative temporal estimate of potential salmon nutrient input. precipitation and temperature data were obtained from nearby sandspit for the majority of years, with data from tlell or alternatively masset used when sandspit data were not available. total precipitation and mean maximum monthly temperature data were generated from september in the previous growth year to october of the growth year for any given tree ring, as these periods more closely reflect the biological year. generalized linear mixed models we assessed the influence of salmon abundance and other potential predictors on tree growth using the arstan-modified growth index and δ15n isotopic enrichment separately, using an information theoretics approach. predictors were selected due to their availability but most importantly, their possible influence on tree growth and δ15n values. for both model sets, response variables were first examined with all potential predictor variables to determine the nature of the potential relationship. predictor variables were also examined for collinearity. statistical analyses were performed using spss v.. predictors of the growth index were assessed using general linear mixed models for the period 1948- individual trees were specified as random factor subjects with year as the repeated measure. covariance was specified as ar in order to account for the autoregressive nature of the data. growth index, the response variable, was transformed closer to normality using a box-cox transformation and salmon abundance, precipitation and temperature were log transformed. thirty-two biologically feasible models were constructed using the predictor variables salmon abundance, precipitation, temperature and the interaction with distance upstream and distance into the forest. spatial predictors were only included as interaction terms due to the consequence of arstan growth index, which standardizes tree series, thus removing the main spatial effects of growth, along with age-related trends. the null model was also included. models were ranked using aikaike’s information criteria for finite sample sizes. where more than one model was ranked with an aic value of > model averaging of the confidence set of candidate models was performed on parameter estimates and variances. the confidence set of candidate models was based on those with a cumulative sum of >95% aicc weight. predictors of δ15n isotopic enrichment for individual tree rings were also assessed using multiple general linear mixed models for the period 1948-. tree ring isotopic signatures obtained in years more recent than were excluded due to elevated total nitrogen and variable δ15n isotopic values in the biologically active sapwood. individual trees were specified as random factor subjects with year as the repeated measure. covariance was specified as autoregressive in order to account for the autoregressive nature of the data. distance into the forest and tree dbh demonstrated an unacceptably high level of collinearity. as there was a strong relationship between δ15n and tree dbh, we adjusted δ15n isotopic enrichment values based on a dbh correction function generated using a linear regression of dbh and the δ15n values of all trees: logy′ij=logyij−vjlogxi−logx where y’ij is the adjusted value of δ15n for year i, yij is the original δ15n value, vj is the regression slope, xi is the annual radial growth of year i and x the δ15n signature for the year we present data for both unadjusted and dbh-adjusted isotopic signatures but only use the latter for analysis using glmms. dbh-adjusted δ15n isotopic signature, the response variable, was transformed closer to normality using a box-cox transformation, growth index was log transformed and distance upstream was converted to a categorical value due to a non-linear relationship with the response variable. twenty-three biologically feasible models were constructed using the predictor variables salmon abundance, growth index, percent nitrogen, categorical distance upstream, distance into the forest and interaction terms. interaction terms involving spatial predictors and the growth index were not included. the null model was also included. models were ranked using aikaike’s information criteria for finite sample sizes.  spatial trends in the transects between estuary and the headwater lake, trees in the carcass zone were larger than adjacent control trees but did not differ in average age. the largest average annual growths were found in trees adjacent to the major spawning reaches of the stream but average growth did not vary with respect to distance into the forest. among the trees with stable isotope data, δ15n signatures ranged from - to among trees, with the highest values occurring with carcasses than in the control trees. data were similar for 1950-. we examined the interaction between tree size and isotope values and found a significant positive slope with a similar trend for both carcass and non-carcass trees although only significant in the former. percent nitrogen was marginally but non-significantly higher for carcass trees than control trees. temporal trends we have relatively complete isotope data for trees. within most of these, the unadjusted δ15n signatures varied approximately 5‰ over the 53 years with some trees showing a depletion historically and others exhibiting higher signatures in the past. five of the six trees found further from the stream had lower signatures relative to trees located closer to the stream, with a single exception, but in this case, salmon carcasses were present in both sites. the highest signatures of any tree, reaching 12‰, occurred in the major bear-foraging area. these signatures were ~10‰ higher than the corresponding control tree 50 m into the riparian zone but this difference was greatly reduced historically as demonstrated by the increased isotopic signatures in the control tree during earlier periods. when we corrected isotope data for size of tree in each year, the temporal trends suggest the control tree increased by 10‰ achieving similar isotopic signatures to bh during the earliest sampling period in smaller increases of size-adjusted values occurred historically in three additional carcass trees and in three control trees. the lowest unadjusted and adjusted signatures occurred in the two control trees above the upper extent of the spawning zone and both of these showed a 3-5‰ increase back to we evaluated the relationship with salmon abundance which fluctuated from ~ in the late 1990s to ~ in the late 1940′s. five of the trees showed a non-significant negative association while nine showed a positive relationship between size-adjusted δ15n and salmon abundance of which four were statistically significant. temporal trends for glmm variables to provide a visual context for variables used in the glmms, we characterize the mean values for standardized growth index, δ15n, percent nitrogen, salmon abundance, temperature and total precipitation. standardized growth was highly variable, with clear peaks in mean growth over the study period. adjusted δ15n was also highly variable with decline in mean values towards the present. percent nitrogen showed low variance among years apart from the recent decades which showed a striking increase in percentages concordant with the outer active sapwood. salmon escapement fluctuates over time with two large pulses occurring in the late 1940′s and early 1960′s and a general decline with time from the 1960′s to the recent. abundances in bag harbour prior to the mid-1940′s are not currently known but are anticipated to have been higher. the environmental variables, mean temperature and total precipitation, were calculated for the biological year. on average, temperatures generally increase with time and bag harbour and the surrounding area receives between and 2100 mm of precipitation for the biological year. predictors of growth using glmms due to a lack of discrimination for the highest ranked model using aicc, the top two ranked models that best predicted the growth index were averaged and the reported cumulative weight generated using both model terms. leading models include the temporally dynamic predictor variable salmon abundance, temperature and the interaction terms salmon*temperature and salmon* distance into the forest. summing parameter weights across the candidate set of best-supported models demonstrated the high relative importance of salmon, temperature and the interaction between salmon and temperature; spatial interaction terms with salmon were lesser contributors to the top model set. although parameters estimates were small, particularly for the influence of temperature, we note that the growth index is a heavily modified metric and both means and variances are standardized between tree series. salmon, the interaction between salmon abundance and temperature had positive parameter estimates, as did the interaction between salmon abundance and the distance into the forest. the interactions between salmon abundance and temperature and salmon abundance and distance into the forest were similar, with a larger influence of salmon abundance on growth at higher temperatures and similarly, a larger effect of salmon abundance with increasing distance into the forest. model predictors are annual salmon abundance, mean maximum temperature for the biological year, distance from stream, distance upstream, total precipitation for the biological year and derived annual tree ring width. models were ranked using aicc values. for predictors of the growth index, the top two models were averaged. * interaction effect. * interaction effect. *interaction effect. predictors of δ15n using glmms in contrast to the predictors of growth, the two spatial variables and their interaction best predicted δ15n isotopic enrichment in individual tree rings. summing parameter weights across the top three best-supported models demonstrated the high importance of the spatial predictors of distance into the forest, distance upstream and their interaction with far lower support for the salmon and growth. distance into the forest, distance upstream and their interaction negatively all had negative parameter estimates. the interaction between distance upstream and into the forest is interpreted as having a lesser importance of distance into the forest with distance upstream. the two spatial predictors were significant. although four of the trees yielded a significant positive covariation between isotopic signatures and salmon abundance, this effect is accentuated on trees occurring at low carcass densities. notably, our model-based inference did not identify salmon abundance as an important predictor. we also ran the glmm using the full temporal complement in order to include the last decade of high nitrogen levels in the sapwood and found the identical top-ranked models relative to the original model-based selection. * interaction effect. discussion and  the spatial patterns of carcass distribution, with declining carcasses with distances upstream and into the forest, delineated our study area and provided the spatial framework for our analysis. of the multiple predators and scavengers observed at bag harbour, haida gwaii, during the five separate salmon spawning runs, black bears were the dominant species responsible for the movement of salmon from the stream into riparian zones. in eight bears uploaded salmon to the riparian zone and consumed approximately one-half of each carcass leaving 330 g of carcass remnants/m or 10 g n/m. these estimates are similar to those documented in alaska, where gende et al. estimated that bears contributed 5 mg n/m to 10 m strip adjacent to a salmon stream. in addition, brown bears contributed approximately 5 mg n/m in a 500 m region adjacent to an alaskan stream, of which 97% was added via urine. in addition to nitrogen provided by salmon carcasses, this input was estimated to account for 10-25% of the forest nitrogen budget. despite the quantitative estimates for quantity and spatial extent of carcasses, we note that abundance and spatial extent of the carcass distribution at bag harbour probably comprises only a fraction of historical levels; in the late 1940s and early 1960′s, the size of the spawning run was approximately seven times the abundance in further, the bear trails that occurred beyond the high carcass zone and above the upper reaches of the spawning area broaden the potential contribution of urine-derived nitrogen at much greater distances than examined. adjacent to the main spawning area, where salmon carcass densities were highest, we also found high annual radial tree growth. however, high growth was also found further into the riparian zones where carcasses were much less abundant but bear trails and assumed urine deposition were prevalent. other studies on the annual growth of alaska sitka spruce that grew adjacent to salmon streams found that these trees grew at more than three times the rate of trees at reference sites above waterfalls of unknown distance from the salmon spawning areas. this elevated productivity of trees led to further examinations of the temporal relationship between salmon abundance and tree growth at the watershed level. after an initial investigation, an examination of nine sites detected positive relationships between tree ring growth and salmon abundance at five locations, which allowed for the historical reconstruction of salmon runs at five locations. similar to drake and naiman, our results generated using indexed growth clearly suggest that salmon abundance has a positive influence on tree growth. however, instead of using reference sites located at a distance, we used within-watershed trees across a sharp ecological gradient from high to very low or zero salmon-nutrient loading. relying on model-based inference, we found that the abundance of salmon and temperature are major predictors of the tree growth index, with all other predictors consisting of salmon abundance and its interaction with temperature and distance into the forest. interestingly, the interaction between salmon and distance into the forest was positive. initially considered counterintuitive, this positive relationship may be the result of elevated nitrogen levels in areas with higher salmon loading where, as a consequence, trees may not always be nitrogen-limited. their growth, although high, may not correlate as closely with salmon as trees that are more nitrogen-limited and experience lower levels of salmon-derived nitrogen. considered initially as controls, the growth response of trees at larger distances from the spawning gravel beds to salmon abundance provides evidence that the spatial effects of salmon on tree growth extend at least 90 m into the riparian zone and likely further. temperature and the interaction between temperature and salmon abundance were also highly ranked as predictors, suggesting that differencing techniques for the removal of environmental signals from a chronology may, at the minimum, remove some of the influences of salmon and subsequently alter reconstructions. patterns of δ15n enrichment tracked what is known of the spatial transfer of salmon into the forest, with higher levels of enrichment found adjacent to the spawning area and declining with distance into the forest and above the upper reaches of the spawning zone. not surprisingly, the use of glmm confirmed our visual assessment; spatial variables were found to be the strongest predictors of δ15n at bag harbour. this finding is similar, although on a much finer spatial resolution, to hocking and reimchen, who found that chum salmon abundance was a leading predictor of δ15n variation in plants and invertebrates at the watershed level. similarly, hilderbrand et al. found that not only δ15n signatures of white spruce decline with distance from salmon-spawning streams, but that the activity of radio-collared brown bears were highly correlated with this same gradient of depletion. the mechanism of the influence of salmon on tree growth is generally considered to be related to the input of salmon-derived nitrogen to nitrogen-limited trees. in addition to growth being positively influenced by salmon abundance, δ15n derived from salmon and subsequently laid down in annual tree rings may provide an additional proxy for salmon abundance. an experimental study that monitored the uptake of δ15n by western redcedar trees found that δ15n uptake in the fall was distributed to leaves and stems in the spring, which substantiates our use of a one year lag effect but also furthers the idea that salmon-derived δ15n could be linked to salmon abundance, similar to growth. studies of natural δ15n levels in tree rings have suggested that they may reflect patterns of atmospheric nitrogen deposition and fertilization studies have demonstrated that tree ring δ15n changes may reflect the historical n regime in a forest. kiernan and johnson found relatively weak support for the idea that δ15n levels in tree rings can be linked to salmon abundance and suggested that nitrogen dynamics in the adjacent stream, soil and within-tree may have acted to obscure their ability to discriminate salmon-derived and autochthonous sources of n. major factors that remain poorly understood, and which potentially confound the use of δ15n levels as a predictor of salmon abundance in this study and others, include denitrification, nitrogen isotope fractionation, translocation of n isotopes across tree rings and local soil and stream nitrogen dynamics. in particular, denitrification, where 14 n is preferentially lost in the system and higher proportions of 15 n remain, could be particularly problematic if it varied with spatial patterns of salmon loading. similar δ15n declines with distance away from the salmon spawning area have been previously attributed to salmon-derived nutrients but these conclusions were strengthened by a lack of similar declines on nearby non-salmon bearing streams. while we have no prior expectation that denitrification varied spatially over the study site we cannot discriminate denitrification from salmon loading using δ15n if the two are similar. the close spatial coupling between the distribution of salmon carcasses at bag harbour and the δ15n signatures of trees, particularly the steep gradients in signatures over short distances, suggest that tree ring δ15n signatures in this watershed may be proxies for historical nutrient uploading in this watershed. although there are a number of feasible alternative explanations, we make several possible inferences if these δ15n patterns are linked to salmon. in the riparian zone adjacent to the current major spawning and bear feeding areas, there has been general temporal continuity of the nutrient uploading during the last five decades. yet, in the upper reaches of the watershed, between the present spawning zone and the headwater lake, and in two trees further downstream in the low carcass zone, where there is no current salmon transfer, δ15n signatures were historically equivalent to those observed in the major feeding areas, but have become depleted by approximately 5‰ by suggesting a historically wider distribution of salmon-derived nutrients. consistent with this suggestion, there was also a depletion in yearly isotopic signatures, relative to historic periods, and which are more similar to the yearly variability in trees in present high carcass densities. there was a limited isotopic signature of the two major pulses of salmon recorded in the late 1940s and early 1960′s apart from two of the low-carcass density trees. bears abandon a greater number of carcasses and consume a smaller proportion of each carcass when salmon escapement is relatively high yet it is likely that numbers of bears in a watershed are determined on average annual salmon escapement, rather than yearly fluctuations. if so, this would decouple any immediate covariation and reduce the chances of detecting year to year fluctuations, thus limiting the application of isotopic proxies for salmon abundance to a longer time frame. tree ring information has been used as a proxy for constructing historical salmon abundances at multiple watersheds with riparian zone stand chronologies and nearby salmon-free controls and at the regional level in climate-sensitive series, both with some success. similar techniques have been employed to examine historical seabird activity extending from the input of guano to forested nesting colonies. our study, although not completing salmon reconstructions and offering only a relatively short time-scale analysis, provides insight into the use of growth indexes and tree ring δ15n as proxies for historical salmon abundance. we observed that δ15n signatures exhibit a relatively consistent decrease from large to small trees both in high and low carcass density zones allowing calibration of δ15n in cores that extend from the cambium towards the center. without a dbh-correction, δ15n values tend to remain somewhat stable over time. yet, because we detected a strongly positive relationship between tree dbh and δ15n and subsequently used this relationship as a calibration factor, δ15n values demonstrate a decline over time. we also found elevated isotopic signatures in the late 1940′s and 1950′s of several trees that are currently in low carcass zones, suggesting elevated salmon uploading at that time. this size factor has consequences for potential use and improvement of salmon abundance reconstructions, although further investigation is required. our study, which reaffirms the importance of salmon to riparian productivity and contributes new insight for the complex, spatiotemporal influences of salmon abundance, temperature and tree location relative to salmon upon conifer growth and δ15n levels. both of these variables provide strong spatial indicators. the use of δ15n for temporal reconstructive purposes remains challenging. we speculate that the appropriate choice of trees in riparian zones across ecological gradients of salmon carcass densities and bear activity and the isotopic correction for yearly changes in tree size may allow improved assessment of tree ring growth as a historical proxy for salmon abundance and distribution during previous centuries in watersheds where ancient trees still persist. abbreviations aicc: aikaike’s information criteria for finite sample sizes; bh: bag harbour; dbh: diameter at breast height; glmm: general linear mixed model. competing interests the authors declare that they have no competing interests. authors’ contributions tr conceived of, carried out fieldwork and statistical analyses for the study. cf completed tree ring measurements, data collection and glmm statistical analyses. both authors participated in drafting the manuscript and have read and approved the final manuscript. 
the marine-terrestrial transfer of salmon provides a substantial pulse of nutrients to receiving ecosystems along the pacific coast of north america and has been shown to enhance productivity and isotopic signatures of conifers and other riparian vegetation. in a small salmon-bearing watershed in haida gwaii, canada, the transfer and distributional pattern of salmon carcasses into the riparian zone by black bears provided a spatial basis for investigating the influence of salmon on sitka spruce tree ring growth and nitrogen isotopic signatures across a gradient of salmon carcass densities in relation to salmon escapement. however patterns observed in individual trees, particularly those at increased distance into the forest, suggest positive relationships with historical salmon abundance. tree diameter demonstrated a positive relationship with δ15n signatures for trees with and without salmon carcass influence. using an information theoretics approach with general linear mixed models, we show that salmon abundance, mean annual temperature and the interaction terms salmon abundance*temperature and salmon abundance*distance into the forest best predict tree growth. although salmon abundance has been previously linked to annual conifer growth and δ15n levels, our approach demonstrates the need to incorporate additional predictors including tree size and opens up the prospect of their dual use as historical proxies for salmon abundance. annual growth was found to be highest in the high salmon carcass zone and δ15n signatures closely tracked the known distribution of salmon carcasses at distances into the forest and upstream. in addition, spatial variables and their interaction are the strongest predictors of δ15n signatures.
 the marine-terrestrial transfer of salmon provides a substantial pulse of nutrients to receiving ecosystems along the pacific coast of north america and has been shown to enhance productivity and isotopic signatures of conifers and other riparian vegetation. an explicitly spatial, within-watershed investigation of the influence of salmon on conifers has never been previously investigated. in a small salmon-bearing watershed in haida gwaii, canada, the transfer and distributional pattern of salmon carcasses into the riparian zone by black bears provided a spatial basis for investigating the influence of salmon on sitka spruce tree ring growth and nitrogen isotopic signatures across a gradient of salmon carcass densities in relation to salmon escapement.  annual growth was found to be highest in the high salmon carcass zone and δ15n signatures closely tracked the known distribution of salmon carcasses at distances into the forest and upstream. tree diameter demonstrated a positive relationship with δ15n signatures for trees with and without salmon carcass influence. using an information theoretics approach with general linear mixed models, we show that salmon abundance, mean annual temperature and the interaction terms salmon abundance*temperature and salmon abundance*distance into the forest best predict tree growth. in addition, spatial variables and their interaction are the strongest predictors of δ15n signatures. however patterns observed in individual trees, particularly those at increased distance into the forest, suggest positive relationships with historical salmon abundance.  using a replicated spatial sampling design across a sharp gradient in salmon nutrient loading, our study provides clear evidence that the temporal pattern in an allochthonous nutrient source and an interaction with temperature and spatial location influences conifer growth. although salmon abundance has been previously linked to annual conifer growth and δ15n levels, our approach demonstrates the need to incorporate additional predictors including tree size and opens up the prospect of their dual use as historical proxies for salmon abundance. spatial subsidysalmonnitrogensitka spruceconiferriparianstable isotope 