Spatial changes in soil stable isotopic composition in response to carrion decomposition

Decomposition provides a critical mechanism for returning nutrients to the surrounding environment. In terrestrial systems, animal carcass, or carrion, decomposition results in a cascade of biogeochemical changes. Soil microbial communities are stimulated, resulting in transformations of carbon (C) and nitrogen (N) sourced from the decaying carrion soft tissues, changes to soil pH, electrical conductivity, and oxygen availability as microbial communities release CO2 and mineralize organic N. While many of the rapid changes to soil biogeochemistry observed during carrion decomposition return to background or starting conditions shortly after soft tissues are degraded, some biogeochemical parameters, particularly bulk soil stable δ15N isotopic composition, have the potential to exhibit prolonged perturbations, extending for several years. The goal of this study was to evaluate the lateral and vertical changes to soil stable isotopic composition 1 year after carrion decomposition in a forest ecosystem. Lateral transects extending 140 cm from three decomposition “hotspots” were sampled at 20 cm intervals, and subsurface cores were collected beneath each hotspot to a depth of 50 cm. Bulk soil stable isotopic composition (δ15N and δ13C) indicated that 1 year after complete soft tissue removal and decay, soils were significantly 15N enriched by 7.5± 1.0 ‰ compared to control soils up to 60 cm from the hotspot center, and enrichment extended to a depth of 10 cm. Hotspot soils also contained 10 % more N compared to control soils, indicating that decomposition perturbs N pools. Our results demonstrate that carrion decomposition has the potential to result in long-term changes to soil biogeochemistry, up to at least 1 year after soft tissue degradation, and to contribute to bulk soil stable isotopic composition.


Introduction
Nutrient hotspots form from the introduction of carbon (C) and nitrogen (N)-rich compounds into an ecosystem, resulting in elevated reaction rates compared to surrounding regions (McClain et al., 2003).For terrestrial and aquatic systems, hotspots may be sourced from fallen trees (Lodge et al., 2016), annual deposition of deciduous leaves (Vidon et al., 2010), animal scat (Erskine et al., 1998;van der Waal et al., 2011), or animal carcasses (Parmenter and Lamarra, 1991;Carter et al., 2007;Wheeler et al., 2014;Wheeler and Kavanagh, 2017).Hotspots sourced from animal carcasses, also referred to as carrion hotspots, significantly alter surface and belowground soil physiochemistry and plant communities in terrestrial ecosystems (Carter et al., 2007;Keenan et al., 2018b).These alterations can have significant long-term impacts; for example, large animal carcasses had measurable effects on a prairie ecosystem for at least 5 years (Towne, 2000), and for a decade or more in the Arctic (Danell et al., 2002).In addition to providing a critical source of C and N, carrion hotspots are important sources of ecosystem heterogeneity (Towne, 2000;Bump et al., 2009b) and promote biodiversity (Barton et al., 2013).
Carrion decomposition occurs in a series of physical (Payne, 1965) and biogeochemical (Keenan et al., 2018b) stages.The breakdown and release of animal tissues provides S. W. Keenan et al.: Spatial changes in soil stable isotopic composition a labile source of nutrients for insect and vertebrate scavengers as well as soil microfauna and microbiota (i.e., nematodes, bacteria, and fungi).Studies evaluating the consequences of carrion decay on soil biogeochemistry have monitored decomposition on a range of timescales, from days (Metcalf et al., 2013;Macdonald et al., 2014;Keenan et al., 2018b;Szelecz et al., 2018) to years (Towne, 2000;Bump et al., 2009a;Keenan et al., 2018a) and in different climatic and geographic settings, including temperate forests (Melis et al., 2007;Cobaugh et al., 2015;Keenan et al., 2018b) and Australian rangeland (Macdonald et al., 2014), as well as under controlled laboratory settings (Carter et al., 2008(Carter et al., , 2010)).Some of the key changes that occur in soils following the deposition and decomposition of carrion include changes to pH (both increases and decreases), increased electrical conductivity, decreased oxygen availability, increased gas fluxes (CO 2 , CH 4 , N 2 O, H 2 S), elevated rates of microbially driven C and N cycling, and increased dissolved compounds available to microbes (NH + 4 , NO − 3 , Ca 2+ , and SO 2− 4 ) (Melis et al., 2007;Aitkenhead-Peterson et al., 2012;Keenan et al., 2018b).
Many of the rapid, pulsed perturbations to soil C and N pools observed at carrion hotspots, such as elevated microbial respiration rates (measured as CO 2 release) and changes to soil pH, return to background biogeochemical conditions during the skeletal stage of decomposition, when soft tissues have been largely or completely degraded by insect and vertebrate scavengers (Cobaugh et al., 2015;Keenan et al., 2018b).However, certain biogeochemical measures, including soil stable δ 15 N composition, have been observed to remain enriched in soils collected at carrion hotspots compared to background soils for a protracted period of time, up to several years (Bump et al., 2009a;Wheeler and Kavanagh, 2017).Soil stable isotopic composition integrates all biogeochemical activity within the soil as well as inputs from plant or animal matter.In contrast with δ 15 N enrichment, no changes in soil δ 13 C composition have been observed in surface soils of decomposition hotspots (Wheeler and Kavanagh, 2017;Keenan et al., 2018b).A variety of studies have demonstrated the potential for natural abundances of 15 N to be used as a tracer of ecological processes, including N input from animals (urea and feces) in N-limited and isolated ecosystems (Erskine et al., 1998) and input of marine taxa (salmon carcasses) to terrestrial and riparian areas (Kline et al., 1990;Koyama et al., 2005).
While 15 N enrichment due to carrion decomposition has been demonstrated in previous work, these studies were limited to surface soils (maximum sampling depth of 10 cm) from the center of the hotspots (Bump et al., 2009a;Wheeler and Kavanagh, 2017).This has left a gap in our understanding of the spatial extent of carcass enrichment, which is ultimately necessary for quantifying ecosystem impacts of these decomposition inputs.Given the potential for natural abundance 15 N to serve as a long-term tracer of decomposition processes, the goal of this study was to evaluate spatial changes in stable 15 N enrichment at a carrion hotspot 1 year post-decay.In particular, the lateral and vertical extent of stable isotope changes as a result of enhanced biogeochemical reactions in a hotspot is largely unknown.Soils beneath and adjacent to former carrion hotspots (up to ∼ 40 cm, the extent of visible fluid migration) were expected to remain 15 N enriched 1 year after decay.Additionally, isotopic enrichment was expected to persist to at least 10 cm depth, the maximum depth examined in previous studies.

Study area and sample collection
The study site was a mixed deciduous forest in east Tennessee (36 • 0 1.0 N, 84 • 13 1.6 W, ∼ 330 m elevation).Soils were part of the Fullerton-Pailo complex and characterized as Typic Paleudults (Soil Survey Staff, 2018).The A horizon extended to approximately 20 cm depth.Five ∼ 23 kg nuisance North American beaver (Castor canadensis) carcasses were placed frozen within scavenger prevention enclosures (1.19 × 0.74 × 0.81 m) and allowed to decay naturally, starting 31 July 2016.As part of a separate study, approximately 75 g of surface soil (0-5 cm depth) was collected a total of five times during decay beneath each animal (Keenan et al., 2018b).For this study, soils were collected on 8 August 2017, 1 year after decomposition, and after bones had been removed from the site.Soils were taken from surface transects as well as from cores obtained at depth below three carrion hotspots.Approximately 30 g of soil from the top 0-5 cm were collected using a 3 cm diameter auger within the hotspot -an elliptical area 40 to 80 cm in diameter of visibly discolored soil (Figs. 1, 2).Surface samples were additionally collected along a linear transect radiating perpendicular from the longest axis of the elliptical hotspot at 20 cm intervals up to 140 cm (Fig. 1).Within the hotspots, soils were cored to a depth of 50 cm using a 10 cm diameter auger; cores were partitioned into depth intervals of 0-5, 5-10, 10-15, 15-20, 20-30, 30-40, and 40-50 cm depth (Fig. 1).
All soil samples were homogenized to a uniform consistency in the field by hand (changing nitrile gloves between samples), removing any rocks, roots, leaves, or vegetation larger than 2 mm.Samples were transported to the lab and processed immediately.Aliquots were oven-dried in triplicate at 105 • C for 48 h to determine gravimetric moisture (Tables 1, S1 in the Supplement).Once dried, subsamples were powdered in an agate mortar and pestle and stored in 1.5 mL tubes until subsequent isotopic analyses.Samples (∼ 25 mg for surface soils; ∼ 50 mg for soils below 20 cm depth) were transferred into 5×9 mm tin capsules (Costech).Isotopic analyses were conducted at Washington University in St. Louis.Samples, standards, and blanks were loaded into a Costech Zero Blank autosampler and combusted in Biogeosciences, 16, 3929-3939, 2019    Soils collected from the center (0-5 cm depth) of all five hotspots and a composite control sample (pooled soil collected from five locations ∼ 3 m from each hotspot) were also analyzed for microbial respiration rates (as evolved CO 2 ), ammonium, nitrate, nitrification potential, pH, electrical conductivity, dissolved organic C and N, and protein content, building from a previous study at the site and following the methods described by Keenan et al. (2018b).In brief, headspace CO 2 was measured immediately after placing and sealing soil into 60 mL serum vials, and also after 24 h (LI-820, Licor Inc.).Vacuum-filtered (1 µm; Ahlstrom, glass microfiber) soil extracts (10 g soil: 40 mL 0.5 M K 2 SO 4 ) were collected after shaking for 4 h at 150 rpm at room temperature and were frozen at −20 • C until the subsequent colorimetric analysis of ammonium and nitrate (Rhine et al., 1998;Doane and Horwath, 2003).Aliquots were oxidized with a persulfate solution to measure dissolved organic carbon (DOC) as evolved CO 2 and dissolved organic nitrogen (DON) colorimetrically as nitrate (Doyle et al., 2004).Nitrification potential was determined colorimetrically using a modified chlorate block method optimized for microplates (Belser and Mays, 1980;Keeney and Nelson, 1982;Hart, 1994).Soil pH and electrical conductivity were measured from a soil slurry (3 g soil: 6 mL deionized water) using a handheld multi-parameter meter (Orion A329, Thermo Scientific).Protein content was determined using the Bradford assay (Wright and Upadhyaya, 1996;Redmile-Gordon et al., 2013).Because the goal of this study was to focus specifically on stable isotopes as long-term tracers in carrion hotspots, only surface soils from the five remnant hotspot centers were processed for full physiochemistry.

Stable isotope analyses
The contribution of carcass-derived nitrogen to bulk soil stable isotopic composition in surface transects was determined using a linear two-member isotope mixing model (Wheeler and Kavanagh, 2017;Keenan et al., 2018b), using bulk control soil δ 15 N composition (0.1 ‰) as one end-member and beaver decomposition fluid (10.2 ‰) as the other.Decomposition fluid is the by-product of microbial and autolytic processes acting on a carcass after animal death.Fluids consist of amino acids, dead and live microbial cells, urea, water, and lipids, and they represent one of the primary mechanisms for return of host's tissues to the surrounding environment.Decomposition fluid isotopic composition was previously determined, using fluids collected from three decomposing beavers left on a shallow plastic tray to intercept fluids (Keenan et al., 2018b).The linear equation for the isotope mixing model (Wheeler and Kavanagh, 2017) was defined as follows: where CDN is the carcass-derived N (%), TEM is the average δ 15 N of soil from the treatment condition (sampling interval along the surface transects), SEM is the end-member control soil stable isotopic composition (0.1 ‰), and FEM is the end-member isotopic composition of decomposition fluids (10.2 ‰).The contribution of control soil-derived δ 15 N to measured treatment conditions was calculated by subtracting CDN (%) from 100 %.
To track changes in δ 15 N between surface soil and soil collected at depth, 15 N values were calculated by subtracting the δ 15 N value of soil at each depth from values obtained at the surface of the hotspot and control sampling locations.Negative 15 N values indicate that surface soils are 15 N enriched compared to soils at depth (Martinelli et al., 1999).

Statistical analyses
Data were analyzed using SigmaPlot to test for significant differences between treatments and controls.For both surface and depth transects, data from the three transects were treated as replicates for subsequent statistical analyses.Significance (p < 0.05) was determined based on one-way ANOVA analyses with Holm-Sidak post hoc testing.Significant differences between control and hotspot soils were determined using paired t tests at each sampling depth or transect interval using R (R Core Team, 2018, version 3.5.0).

Surface soil biogeochemical changes during decomposition
During carrion decomposition, fluids sourced from the carcass were released into the surrounding environment (Fig. 2).The pulse of nutrient-rich fluids resulted in significant longterm changes to surrounding soil physiochemistry following fluid degradation by soil microbial communities (Table 1 and Table S1).In particular, after 1 year of decay, soil pH was significantly lower than control, initial, and pre-decay soils (p < 0.001; F = 59.317).In addition, bulk soil δ 15 N remained significantly enriched compared to control and starting soil isotopic composition (p < 0.001; F = 27.948).
Other physicochemical parameters, including conductivity, microbial respiration, DOC, DON, ammonium, nitrate, and nitrification potential all returned to background conditions after 1 year.With the exception of the 1 year samples, data were previously published in Keenan et al. (2018b) and are included here for comparison.

Lateral changes in stable isotopic composition
Soils were significantly 15 N enriched within the visible carrion hotspot (mean soil composition 7.5 ± 1.0 ‰) and up to 60 cm from the hotspot center (2.2 ± 0.5 ‰) compared to composite control soils (0.1 ‰) (Fig. 3, Table S2) (paired t test, p = 0.016).Soil δ 15 N values gradually declined, reaching background abundance values around 80 cm from the hotspot center.In contrast, there were no significant differences between control and hotspot soil δ 13 C and no differences as a function of distance from the hotspot center (one-way ANOVA, p = 0.464; F = 1.004).C/N ratios were lower within the hotspot and exhibited a gradual and significant increase with increasing distance from the hotspot center (Fig. 3).However, there were no significant differences between hotspot and control soil C/N.The influence of carrion decomposition on soil stable δ 15 N isotopic composition decreases with increasing distance from the hotspot (Figs.3a, S1 in the Supplement).Soil C/N composition follows a linear trend, increasing by 0.07 per centimeter from the hotspot center (Fig. S1).Based on linear two-member isotope mixing models, carcass-derived fluids exhibit a linear decrease in contribution to soil isotopic composition with increasing distance from the hotspot.Carcasses contribute to soil stable δ 15 N isotopic composition up to 60 cm from the hotspot center (Fig. 4), an area that was beyond the decomposition island and was not visibly discolored (Fig. 2).

Vertical changes in stable isotopic composition
Soil collected at depth beneath the three mortality hotspots was significantly 15 N enriched compared to control soils up to 10 cm depth (Fig. 5).Surface hotspot soils were also enriched at 30 cm depth compared to the control.Control soils became more 15 N enriched with increasing depth.There was no significant difference between control and hotspot soil δ 13 C and C/N values, and both exhibited the same trends with depth.Soils exhibited 13 C enrichment with increasing depth and a decline in C/N ratios.
Control soils exhibited a strong positive linear relationship between the negative of the natural log of bulk soil % N and stable isotopic composition, reflecting decreasing N and C availability with increasing depth.Decomposition results in a shift in the observed hotspot soil N isotopic discrimination (D, or the slope of the linear regressions) (Natelhoffer and Fry, 1988), leading to a breakdown of a strong linear relationship compared to control soils (Fig. 6).D does not change for C isotopes in control or hotspot soils.
In general, soils exhibit a trend of increasing 15 N (the difference between soil δ 15 N value at a specific depth and δ 15 N at the surface) with depth, reflecting 15 N enrichment in deep forest soils (Martinelli et al., 1999).Hotspot soils exhibited lower 15 N values compared to control soils, indicating little change in 15 N enrichment with depth (Table 2).Control soils displayed increasing 15 N with depth, a pattern globally observed in forest soils (Martinelli et al., 1999) (Table 2).Combined with D (Fig. 6), hotspot and control vertical profiles have distinct N sources and exhibit distinct N pools with depth.
Soils collected from 0-5 cm exhibited a significant contribution from decomposition fluids (Fig. 5).As depth increases beyond 10 cm, there is no change to the proportional contribution of decomposition fluid to bulk stable isotopic composition.By 30-40 cm and 40-50 cm depths, hotspot soil δ 15 N and δ 13 C compositions are similar to control soils (Fig. 5), indicating limited, if any, input from decomposition fluids.

Discussion
Soils associated with carrion decomposition hotspots retained biogeochemical markers of vertebrate decay at least 1 year after soft tissue degradation.Within the hotspots, soils remained 15 N enriched compared to control locations.At this specific location and soil type, decomposing animals have the potential to exert long-term changes, here at least 1 year, on surface and subsurface soil stable isotopic composition.
The beaver carcasses used in this study, which were between 20 and 25 kg in mass, resulted in measurable changes to soil biogeochemistry down to 10 cm depth and up to 60 cm away from the hotspot center, beyond the area that was visibly discolored from decomposition fluids.The total % N measured in soils can be used to approximate the contribution of beaver N to soil.During active decomposition, hotspot soils contained 36 % more N compared to control soils (0.362 % N vs. 0.267 %).After 1 year, hotspot soils still contained 10 % more N than control soils (0.285 % N vs. 0.260 %), reflecting a loss of ∼ 28 % of the beaver-derived N in 1 year.Based on prior calculations, each beaver introduced approximately 0.47-0.79kg N and 2.7-4.6 kg C to the soil (Keenan et al., 2018b).After 1 year, approximately 0.05-0.08kg N derived from the carcass remains within the soil.The linear two-member isotope mixing model indicates that the contribution of carcass-derived fluids to soil iso-topic composition decreases with increasing distance from the hotspot (Fig. 4).However, the observed isotopic enrichment likely reflects a contribution of a 15 N-enriched input (carcass fluid) as well as subsequent diagenetic fractionation driven by soil microorganisms.It is likely that the elevated N availability provided by a decomposing carcass would result in additional nitrification and denitrification, which would increase the δ 15 N independent of source mixing.Nitrification potential rates were elevated during earlier stages of decomposition but were not different from control soils after 1 year (Table 1).The two-member mixing model assumes a simple mixing of soil and carcass-derived N, but diagenetic fractionation is also likely involved in the observed δ 15 N patterns along the surface transects.
The contribution of carcass-derived N to soils at depth as well as laterally is influenced by a variety of physical and climatic variables.Here, decay occurred during the summer in east Tennessee, with an average high of 32.2 • C for the month of August.Carcasses were exposed to measurable precipitation 6 out of 10 d after placement, preventing soft tissues from significant desiccation and supporting abundant blowfly larvae and other insect activity.Blowfly larvae migrating away from the carcasses on the surface and within the soil to pupate likely provided an important physical mechanism to distribute beaver-enriched N to surrounding soils.Blowfly larvae can move up to 10 m away from the carcass, and typically extend down into the soil up to 10 cm depth, depending on the soil substrate properties (Gomes et al., 2006).As blowfly larvae disperse, they have the potential to physically transport decomposition fluids acquired internally or externally during feeding, release excrement during migration, or die, leaving their tissues to degrade.An estimated 66 % of pupae that disperse to pupate die en route (Putman, 1977).Rainfall may have also contributed to the downward movement of decomposition fluids.
The temporal persistence of isotopic enrichment hotspots is currently unknown but is likely to be ecosystem, carrion type, and carrion mass-specific.A larger carcass would be expected to result in greater lateral and vertical dispersal of carrion-derived fluids, as well as greater changes to ecosystem processes, because of the greater volume of decomposing soft tissue (Baruzzi et al., 2018).In addition, larger carcasses may host a larger and longer-lived insect community (Parmenter and MacMahon, 2009), including blowfly larvae, which may further nutrient dispersal and may impact a larger area.Bump et al. (2009b) observed elevated foliar δ 15 N values in plants growing on sites impacted by deer carcass (∼ 56 kg) decomposition at least 2.5 years after decay in a temperate hardwood forest, suggesting a long-lived hotspot signature.In some ecosystems, such as the Arctic tundra, isotopic enrichment is likely to persist for even longer based on perturbations to C and N surrounding muskox after 5 to 10 years of decay (Danell et al., 2002).
Increasing δ 13 C and δ 15 N values with depth in soils has previously been observed in a variety of soil types and cli-matic conditions (Natelhoffer and Fry, 1988;Martinelli et al., 1999;Billings and Richter, 2006).Changes to δ 13 C with depth are due to progressive cycling of C through microbial biomass (Liang et al., 2017), where selective preservation and biochemical fractionation together lead to 13 C-enriched organic C in soil (Natelhoffer and Fry, 1988;Billings and Richter, 2006).While we observed a similar increase in δ 13 C with depth, we did not see a significant change in 13 C as a result of carcass enrichment.Wheeler and Kavanagh (2017) similarly did not observe a change in soil δ 13 C following carrion decomposition, which may be due to the degradation of carcass-derived C (and uptake by microbes and blowfly larvae) combined with elevated background C in soil compared to N.
Increasing δ 15 N values with depth reflects two broad biochemical processes leading to fractionation, both likely driven by microbial activities.First, the preferential excretion of 15 N-depleted compounds during catabolism and anabolism leaves the residual microbial cells and soil 15 N enriched.Second, kinetic fractionation associated with gaseous N loss is also known to result in enrichment, depending on the microbial communities present and N mineralization rates (Evans, 2001;Robinson, 2001;Liang et al., 2017).Over time, as soil profiles develop, accretion of 15 N-enriched microbial cells, particularly fungi, leads to isotopic enrichment at depth (Billings and Richter, 2006).In contrast, plant and leaf litter are the dominant contributors to N pools in surface soils in most temperate forest ecosystems (Vidon et al., 2010), resulting in surface soils that are isotopically depleted compared to the soils at depth.Decomposition hotspots, however, disrupt the expected pattern (Fig. 5), causing surface enrichment and likely leave a lasting impact on soil stable isotopic composition.
For systems at or near steady state conditions, the difference in isotopic enrichment between soils at depth and the surface ( 15 N) provides a way to compare soils from different geographic and climatic locations (Martinelli et al., 1999), and was used here to compare hotspot soils and those collected at control locations. 15N values observed in the control depth profile are within the expected range observed in temperate forests worldwide (2.7 ‰ to 9.1 ‰) (Table 2).However, as a consequence of carrion inputs and decay, 15 N values are more similar to those observed in tropical forest ecosystems (1.1 ‰ to 4.3 ‰).In tropical systems, lower 15 N values are thought to reflect more open N cycling with elevated N losses (nitrification, nitrate leaching, and ammonia volatilization) under conditions of elevated total N inputs (Martinelli et al., 1999).Whether our observed changes in 15 N are due to elevated N cycling rates, disequilibrium effects across the soil profile due to changing N inputs from a system dominated by atmospheric dry and wet deposition of nitrate and ammonium to one with carrionsourced N, or both, is not known.
The observed stable isotopic discrimination (D or the slope of the linear regressions) did not differ for control and hotspot soil δ 13 C (Fig. 6), suggesting that C cycling and pools in soils 1 year after carrion decay are not altered.In contrast, D values for δ 15 N were different between control and hotspot soils, which reflects a loss of the linear relationship and indicates distinct N sources for the two soil profiles.This also emphasizes that decaying carrion provide an important and potentially distinct N pool for soil ecosystems that have the potential to mask natural (background) processes that control soil profile 15 N gradients with depth.In addition, differences in D values between the two soils suggest that there may be less discrimination occurring within hotspot soils compared to control soils, likely due to the rapid input of an isotopically enriched N pool (Evans, 2001).
Hotspot soils received the input of beaver-derived fluids (10.2 ± 0.4 ‰) (Keenan et al., 2018b) as well as soft and hard tissues (1.0 ‰ to 4.0 ‰ for beaver bone collagen from Minnesota; Fox-Dobbs et al., 2007).Stable isotopic composition of surface soils strongly suggests that decomposition fluids are a significant contributing source to bulk soil stable isotopic composition up to 60 cm from the hotspot center (Fig. 4), and δ 15 N-enriched values in the soil profile at depth also suggest some contributions up to 10 cm depth (Fig. 5).Beyond 10 cm depth, control soils and hotspot soils are indistinguishable, suggesting that decomposition fluids do not significantly influence soil stable isotopic composition.Rather, natural δ 15 N enrichment due to soil accretion processes can explain the observed soil stable isotopic composition.

Conclusions
The decay of ∼ 23 kg North American beaver carcasses resulted in rapid (within days) and long-lived (up to 1 year) 15 N enrichment in forest soils up to 10 cm depth and ∼ 60 cm distal.Observed 15 N enrichment at depth and laterally is likely due to a combination of physical movement of fluids during decomposition and the transport of fluids by insects, particularly blowfly larvae.In this system, rainfall during decomposition may have also acted as a physical transport mechanism.While likely to be significantly influenced by carcass size, climate, and soil type, decomposition has the potential to exert long-lived influences on soil stable isotopic composition.

Figure 1 .
Figure 1.Schematic cross section view of the locations of soil samples (stars) collected from each of three carrion decomposition sites.Dashed line represents the hotspot -the area of visibly discolored soil.Soils collected at depth extended to the B horizon.The visibly discolored area of soil due to carrion hotspot formation extended approximately 35-40 cm from the hotspot center along the surface and to a few centimeters depth.

Figure 2 .
Figure 2. View of a beaver after placement (a) and during advanced decay (b) to demonstrate the lateral migration of carcass-derived fluids during decay.Both photos are from the same animal, and (b) were taken during advanced decay (8 August 2016).Visible extent of fluid migration is outlined in the white dashed line.

Figure 3 .
Figure 3. Lateral changes in soil stable (a) δ 15 N and (b) δ 13 C isotopic composition and (c) C/N ratios extending from carrion hotspot centers.Soil was visibly discolored 35-40 cm from the center (here, 0 cm distance).Letters indicate soil samples taken at discrete distances from hotspot center that were not significantly different based on an ANOVA and Holm-Sidak post hoc test.Asterisks denote the significant differences between control and hotspot soils (t test) at each respective distance.The dashed line represents control surface soil (0-5 cm) composition.

Figure 4 .
Figure 4. Results of linear two-member isotope mixing distinguishing the contributions of soil and carcass fluid to bulk soil stable isotopic composition.

Figure 5 .
Figure 5. Stable isotopic composition and C/N ratios for soils beneath carrion hotspots (closed circles) and at a control location (stars).Letters indicate hotspot soil samples as a function of depth that were not significantly different based on post hoc testing, and asterisks indicate significant differences between control and hotspot soils at each depth interval (both based on one-way ANOVA, p < 0.05).

Figure 6 .
Figure 6.Bulk soil stable isotopic composition and corresponding negative natural log % N (a) and % C (b) for hotspot and control soils with depth.Linear regressions were fit to hotspot and control datasets.

Table 1 .
Keenan et al. (2018b)hemical data during 1 year of decomposition.Letters indicate hotspot soil samples within each measured dataset (i.e., pH) that were not significantly different based on one-way ANOVA (p < 0.05) and Holm-Sidak post hoc testing.Asterisks indicate significant differences between control and treatment soils.N.M. indicates parameters were not measured.Control samples were homogenized into a single representative sample and do not have standard deviations.Data, except for 1 year post decay, were previously published inKeenan et al. (2018b).

Table 2 .
Differences in soil δ 15 N at depth and δ 15 N in surface soils for hotspot and control depth profiles.