Elephant megacarcasses increase local nutrient pools in African savanna soils and plants

Reed, Courtney G.; Budny, Michelle L.; du Toit, Johan T.; Helcoski, Ryan; Schimel, Joshua P.; Smit, Izak P. J.; Strydom, Tercia; Tallian, Aimee; Thompson, Dave I.; van Coller, Helga; Lemoine, Nathan P.; Burkepile, Deron E.

doi:https://doi.org/10.5194/egusphere-2024-1514

Preprints

https://doi.org/10.5194/egusphere-2024-1514

Preprints

30 May 2024

| 30 May 2024

Elephant megacarcasses increase local nutrient pools in African savanna soils and plants

Courtney G. Reed, Michelle L. Budny, Johan T. du Toit, Ryan Helcoski, Joshua P. Schimel, Izak P. J. Smit, Tercia Strydom, Aimee Tallian, Dave I. Thompson, Helga van Coller, Nathan P. Lemoine, and Deron E. Burkepile

Abstract. African elephants (Loxodonta africana) are the largest extant terrestrial mammals, with bodies containing enormous quantities of nutrients. Yet we know little about how these nutrients move through the ecosystem after an elephant dies. Here, we investigated the initial effects (1–26 months post-death) of elephant megacarcasses on savanna soil and plant nutrient pools in Kruger National Park, South Africa. We hypothesized that: (H1) elephant megacarcass decomposition would release nutrients into soil, resulting in higher concentrations of soil nitrogen (N), phosphorus (P), and micronutrients near the center of carcass sites; (H2) carbon (C) inputs to the soil would stimulate microbial activity, resulting in increased soil respiration potential near the center of carcass sites; and (H3) carcass-derived nutrients would move from soil into plants, resulting in higher foliar nutrient concentrations near the center of carcass sites. To test our hypotheses, we identified 10 elephant carcass sites split evenly between nutrient-poor granitic and nutrient-rich basaltic soils. At each site, we ran transects in the four cardinal directions from the center of the gravesite, collecting soil and grass (Urochloa mosambicensis) samples at 0, 2.5, 5, 10, and 15 m. We then analyzed samples for CNP and micronutrient concentrations and quantified soil microbial respiration potential. We found that concentrations of soil nitrate, ammonium, ¹⁵N, P, sodium, and potassium were elevated closer to the center of carcass sites (H1). Microbial respiration potentials were positively correlated with soil organic C, and both respiration and organic C decreased with distance from the carcass (H2). Finally, we found evidence that plants were readily absorbing carcass-derived nutrients from the soil, with foliar %N, ¹⁵N, iron, potassium, and sodium significantly elevated closer to the center of carcass sites (H3). Together, these results indicate that elephant megacarcasses release ecologically consequential pulses of nutrients into the soil, which then move into above-ground nutrient pools in plants. These localized nutrient pulses may drive spatiotemporal heterogeneity in plant diversity, herbivore behavior, and ecosystem processes.

Received: 21 May 2024 – Discussion started: 30 May 2024

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this preprint. The responsibility to include appropriate place names lies with the authors.

Download & links

Preprint (PDF, 1643 KB)

Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
Preprint (1643 KB)

Supplement (1461 KB)

Download & links

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Journal article(s) based on this preprint

25 Mar 2025

Elephant megacarcasses increase local nutrient pools in African savanna soils and plants

Biogeosciences, 22, 1583–1596, https://doi.org/10.5194/bg-22-1583-2025,https://doi.org/10.5194/bg-22-1583-2025, 2025

Short summary

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2024-1514', Shawn J. Leroux, 17 Jun 2024

Summary
Animals impact elemental cycling in many direct and indirect ways. Evidence from several biomes demonstrates that even after death, animal carcasses can change the biogeochemistry of ecosystems and these impacts can be long lasting. Most studies of carcass impacts on ecosystems, however, are done on small to medium (1kg to 200kg) sized animals. In this contribution, the authors investigate the effects of elephant megacarcasses on the biogeochemistry of soils and plants. The authors report significant effects of elephant carcasses on components of soil and plant elemental cycling and they discuss how these effects may be important components of spatiotemporal heterogeneity in ecosystems.
General comments
1) Overall, I found the writing good. The authors have crafted a nice narrative that makes a compelling case that megacarcasses can be important parts of ecosystems and therefore we need to learn more about the impacts of these carcasses on ecosystems.
2) I have a few questions about the analysis. The effective sample size is 10. Obviously, it is hard to find carcasses (I would have great difficulty in finding 10 fresh moose carcasses in my system!) but the authors are trying to squeeze a lot of information out of very few data points. I have the following specific questions about the analysis:
i) While I like the transect approach, the design may have been stronger if the authors had random transects (ie, transects with no known carcass) like Risch et al. work. This would strengthen inference.
ii) lines 180-183. The author’s approach to checking for normality of response data does not seem sound to me. The assumption of normality (for linear models) is normality in light of the model, i.e., investigating the normality of residuals is a more common approach to this. Either way, it is often better to avoid transforming the data and generalized linear models do allow for a lot of flexibility to fit different error distributions. For example, the gamma family in glm is very flexible and can handle log-normal data sets. Did the authors try different families of error distributions before transforming their data?
iii) lines 187-189. How many data points did the authors have per estimated parameter in the most complex model here?
iv) line 194. This is fine but I think Burham & Anderson would say that any model within deltaAIC of 2 of the null model should not be considered to be supported. In several cases, the authors interpret top models that are ranked above the null but within deltaAIC of 2 of the null as supported (e.g., lines 217-218, 218-221).
v) what R^2 are the authors reporting? In the captions of Tables S1 and S2 (thank you for providing full AIC and coefficient tables), the authors state “R^2 is the proportion of variance explained by a model”. This is unclear. These are mixed models, and the most common approach is to report the marginal R^2 and conditional R^2. Is the R^2 in these tables one of those or another pseudo R^2? This is critical for many reasons but most importantly, given the small sample size and large number of mixed-models, I would expect at least one of the models to not converge. There are many indicators when a mixed-model does not converge and one of the best is when the marginal R^2 = conditional R^2. Without having both of these pieces of information, the reader is unable to adequately assess the fit of the models. Other indicators of models not converging are coefficients estimates or errors that are very large or very small (i.e., 0 – see next comment).
vi) I am confused by the magnitude of Table S2 sodium and iron coefficients and/or the scale of reported on the y-axis of Figure 5 for these. The iron coefficients in Table S2 seem small relative to the Figure 5c ? Or am I misreading things?
3) the reporting of results could be improved. I recommend, the authors report: top ranked models (AIC + measure of independent fit like R^2). Then report effect size or relationships (coefficients). I found key statistics to be missing throughout. Statements like “Phosphate concentrations were greater in granitic soils…” would be more informative if they included the coefficient + error in parenthesis. Coefficients can be reported for the top-ranked model or from model averaged results when there are several competing models.
4) in section 3.2 I think the reader may be more interested in coefficients and confidence intervals around those relationships than p-values that are currently reported.
Specific comments:
5) I found the use of three different terms that mean similar things (nutrient flows, ecosystem processes, nutrient availability) in the introductory sentence confusing. I recommend the authors replace “nutrient availability” with “ecosystem processes” or “nutrient flows”. Surely living animals (not just carcasses) influence nutrient availability (which is just a part of a continual nutrient cycle).
6) line 83. I believe there is no “e” at the end of the citation Risch et al.
7) lines 96-111. How do these elephants die? As someone with no experience with megacarcasses, I would appreciate some insight on the causes of death. Most large herbivore deaths in my empirical systems are from predation which I assume is not the case for elephants.
8) really excellent job with clear hypotheses and nice work carrying forward these hypotheses throughout the ms – really makes the job easier for the reader.
9) lines 132-133 Why 10cm deep core? Is that mineral soil only?
10) the discussion is well done – concise and touches on all hypotheses.
11) Figure 1 is an outstanding visual!
12) in figures 2-5 I recommend the authors consider reminding the reader of the sampling resolution because the jitter of points makes it impossible to see what distances were measured below 5m.

Citation: https://doi.org/10.5194/egusphere-2024-1514-RC1
- AC1: 'Reply on RC1', Courtney Reed, 21 Jun 2024
  
  Thank you for this thoughtful review! We will work on incorporating your feedback and will post a full response once we have done so.
  
  Citation: https://doi.org/10.5194/egusphere-2024-1514-AC1
- AC2: 'Reply on RC1', Courtney Reed, 04 Sep 2024
  
  Thank you for your feedback, which was very helpful in improving the manuscript. Our responses are attached, with reviewer comments in black and author responses in blue.
  
  Citation: https://doi.org/10.5194/egusphere-2024-1514-AC2
RC2:
'Comment on egusphere-2024-1514', Sarah Keenan, 27 Jun 2024
Summary:
Reed and coauthors present a well-written study on the impacts of megacarcasses (elephants) to soil biogeochemistry after up to 2 years of decomposition. The authors examined 10 carcass hotspots with 5 carcasses each on two different soil types. They quantified soil major and trace element chemistry as well as plants associated with the hotspots to determine if carcasses influenced soil N and P chemistry and if those elements were subsequently enriched in vegetation. The current version of this manuscript does not adequately describe the methods in enough detail to make the work reproduceable. Additionally, the handling of the data for statistical analyses is strange and non-standard. The discussion needs to be re-written to better emphasize the importance of the work (as framed in the introduction). I think this work has potential to be an important contribution, but there needs to be some major revisions.

General comments:
The importance of this study in adding to our knowledge about nutrient transfer at carrion hotspots is not emphasized clearly in the discussion. The introduction frames how megacarcasses may be “functionally different than smaller carcasses” but never returns to this aspect in the discussion, which is really where this work could add to our knowledge. Adding more to the discussion would help address this issue and would make the impact of the work clearer.

Parts of the results belong in the discussion, and I’ve tried to highlight those below in specific comments.

The methods need significantly more specific details, highlighted in specific comments. Additionally, there were no control soil or plant samples examined here. Please describe in the methods why there were no controls.

The handling of the data for statistical analyses is non-standard and not clearly justified. If data were non-normally distributed (it seems like some datasets were and some were not), why not just use a non-parametric statistical test rather than log-transforming the data? It is a bit strange to log-transform some data but not all. The approach of adding 0.001 to zero values is also not correct (described below in specific comments).

The presentation of elemental data for soil and plant composition is non-standardized throughout. Some data (i.e., iron) are presented as mg/kg (is this soil dry weight?), while others are presented as % (Ca% of what?) in the same figure (figure 5 for example). Other data are presented as mg/L (figure 2). Part of this confusion is from the missing details in the methods that clearly explain how these data were generated. In several of the figures there is a statement about back-transformed data, which is also confusing.

I can appreciate that finding carcasses that have decomposed for the same amount of time is challenging, but 1 month to 26 months is a huge range of time (at least from what we know from not megacarcasses). The biogeochemical processes occurring at a carcass decaying after 1 month postmortem is very different than a carcass that has been decaying for 26 months (from smaller carcasses). It would be useful to see some of the data, particularly ammonium, plotted as a function of postmortem interval (months) even if that is not a variable that could be included in statistical analyses because of the small sample size. It would also be helpful to see if the postmortem interval for the 10 carcasses is evenly distributed between the two soil types or if one has more fresh carcasses and the other has older carcasses, that could help with interpretation of the results.

I think it may be useful to add some photos to supplemental information (or even the main text) showing what the carcasses/sites looked like (maybe representative images from a fresher carcass and one that is older).

Specific comments:
Lines 99-100: There should be more details provided on the soil type and what makes the granitic soils “nutrient poor” compared to soils developed from a basalt protolith. Because soil type becomes an important part of this study, the details of the soil types need to be expanded in the introduction.

Line 132: Include a citation or discuss why soil samples were collected to a depth of 10 cm rather than the upper 5 cm. For decomposition studies, typically the upper 5 cm is examined, not the upper 10 cm.

Line 145: More details are needed beyond “measurements of soil ion concentrations”. What instrumentation was used? What specific extraction protocol was followed? I’m assuming deionized water was used (1:2 soil to deionized water?), but those details are not provided. How long were samples mixed (shaking platform?), what speed, etc.

Line 148: “mass spectrometry”—elaborate on what this means with respect to instrumentation used to analyze cations. Here and throughout the methods, please also include what standards were used for the different analysis types.

Lines 146-150: Clarify if these analyses were conducted on the water extracts.

Line 152: Were stable isotope analyses conduct on oven-dried soil? 10 g is an exceptionally large amount of soil—how much was actually analyzed with EA-IRMS? Were samples powdered prior to combustion?

Line 154 (and throughout with respect to stable nitrogen isotope results): The authors refer to “¹⁵N” measurements, but surely this should be presented as the ratio of 15/14N and in delta notation? In the methods here there also needs to be more description of the standard, the materials used for linearity, and the analytical precision of the instrument.

Line 175: More details on the ICP-MS are needed, including standards, detection limits, etc. Additionally, were these samples digested in nitric acid? Water? How long were they microwaved?

Line 182: Adding some random number to each variable is not a standard way to handle data that are zero in your dataset (or if it is, there is no citation here and I am not familiar with that approach). Typically for geochemical data (like what was generated with ICP-MS), you can replace zero values with ½ the detection limit to remove non-zero data. There are other more technical ways to deal with zero values from a statistical standpoint, but the ½ the detection limit is the easiest and has the longest history of use. Please justify the use of your approach or re-run the analyses following a standard method for handling non-zero data in a geochemical dataset.

Line 255: The part of the sentence that reads “…we found evidence that N from carcasses had moved from soils into plants” does not belong in the results section. This is interpretation and should be moved to the discussion.

Lines 256-258: Similar comment as above where the content of this sentence is interpretation and should be moved to the discussion.

Lines 295-297: I’m not quite sure I understand the logic presented here. First, soil microbial biomass was not measured. The respiration potential (through production of CO2) was measured, but heterotrophic activity (which is how respiration can be interpreted) consumes oxygen. I think the phrasing here needs to be re-worked to not imply that the soil respiration (and the communities producing CO2) are not necessarily the same that are driving nitrification.

Lines 304-305: There are prior studies that demonstrate the impact of increased organic C during decomposition on soil microbial processes that should be cited here (see studies by DeBruyn and colleagues)

Lines 310-311: I’m not sure that this is phrased correctly. Phosphorus (predominantly as phosphate) is considered immobile in soil partly because of low solubility because it is often sorbed with Ca, Fe, Al or organics, and the release of P is tightly controlled by soil (or fluid) pH. N does not face the same sorts of sorption immobilization constraints. I think if you rephrased it to clarify that P and N are held within different reservoirs within soils that make them behave differently (and add some citations), that would help.

Lines 324-328: As mentioned above, because the composition of the two soil types were not included, this part of the discussion is not supported by the results. It’s unclear if the authors here are trying to say that the basaltic soils contain more nutrients after being impacted by decomposition or if the native state of the soils (background conditions) are more nutrient rich. I think if the introduction described the background chemistry of the two soil types this would be better supported. Additionally, basalt and granite contain different types of minerals and therefore additional sources of elements like P. I don’t know what the specific mineralogy is of these two rock types in KNP, but it might be worth exploring. In particular, the presence of apatite (Ca-P bearing mineral) in the granite might also be contributing to elevated P measured in the granite soils.

Figure 5 (and others): I’m a little bit confused by the figure caption. I think it would be better to present this as selected elements plotted as a function of distance. These are not the results of generalized linear mixed models, but the selection of how to present the data were informed by the models.
Citation: https://doi.org/10.5194/egusphere-2024-1514-RC2
- AC3: 'Reply on RC2', Courtney Reed, 04 Sep 2024
  
  Thank you for your feedback, which was very helpful in improving the manuscript. Our responses are attached, with reviewer comments in black and author responses in blue.
  
  Citation: https://doi.org/10.5194/egusphere-2024-1514-AC3
RC3:
'Comment on egusphere-2024-1514', Anonymous Referee #3, 05 Jul 2024

This manuscript presents data on the influence of elephant carcasses on nutrient availability in South African savanna soils. It would be a surprise if a decaying elephant did not increase nutrient concentrations in the proximity of the carcass, but there are some interesting differences among nutrients in terms of the distance over which the effects extend. There are parallels with other nutrient hotspots in tropical ecosystems, including glades in African savannas (e.g. Augustine 2003) and leafcutter ant nests in tropical forests (e.g. Hudson et al. 2009) - it would be worth introducing these into the discussion for comparison. I have several questions about methodology and results that should be addressed before this manuscript could be acceptable for publication.
Line 38 – these are not graves. A grave is an excavation for burial.
Line 77 – what about the amounts of cations in an elephant?
Line 89 – cations are not micronutrients. Is there direct evidence for widespread (or any) cation limitation of growth in savanna ecosystems?
Line 99 – please include classifications for the granitic and basaltic soils in one of the internationally recognized systems. In Soil Taxonomy these are presumably Inceptisols / Alfisols and Vertisols, respectively?
Line 138 – freezing soil has implications for subsequent measurements of extractable nutrients (e.g. Turner and Romero 2009). This should be mentioned here. Given the apparently very high values for some measurements (see below) I suspect that pretreatment had a major impact on results.
Line 145/146 – please explain the difference between phosphate and plant-available P. As written, it appears they were both measured in the water extracts. Most plant-available P tests are not conducted in water (e.g. Olsen, Mehlich, Bray, etc).
Line 158 – was there any inorganic C in the samples? Savanna Vertisols developed in basalt can have considerable carbonate concentrations, albeit often in subsoil. Soil pH values would help indicate this possibility – how did carcasses affect soil pH?
Line 161 – was moisture standardized prior to the incubations?
Line 182 – an alternative is to set values to ½ detection limit.
Line 188 – I understand that there were insufficient carcasses to allow inclusion of carcass age in models. However, major differences would be expected between carcasses aged 1 month vs 2.5 years. Is there any way to provide an indication of the magnitude of the age effect? How would distance effects look if young carcasses were excluded, for example?
Line 228 – It is perhaps not surprising that P concentrations showed little variation with distance, given that P was measured in water extracts (i.e. the extraction is recovering a relatively small pool of soluble P).
Line 230 – how is plant-available P defined here?
Line 286 – Soil extractable nutrients should be expressed on the basis of dry soil, not volume.
Line 286 - these very high extractable nitrogen concentrations are presumably in part a consequence of soils being frozen prior to analysis. Another factor is time between sampling and freezing - or storage prior to freezing. Please provide a statement about sample treatment prior to analysis (time from sampling to freezing, storage conditions during this time, etc, as relevant).
Line 308 – the high available P and tissue N:P ratios (see below) indicate that there is no P limitation here. This might limit the likely influence of carcasses on foliar P, as found here (i.e. foliar P is not a strong indicator of the extent to which carcasses change P availability in general).
Line 322 – water-extractable P represents a tiny proportion of the total soil P, so reliance on this procedure probably limits the possibility of detecting change in soil P.
Line 323- 330 – this text largely repeats results.
Line 340 - missing here is a discussion of the ecological consequences of the findings. What are the implications for plant and microbial ecology in savanna ecosystems?
Figure 2B, C – these values should be presented in mg/kg soil.
Figure 2B – these are extremely high nitrate concentrations, even out to 15 m. For example, 100 mg/L is equivalent to 200 mg/kg based on a 1:2 soil to water extraction ratio. Extractions done quickly after sampling and in 2 M KCl are in the range of 1-5 mg/kg. This seems to be a clear indication of storage effects.
Figure 2B, C – are these values as NO3/NH4 or on an N basis?
Figure 2D – please express stable isotope ratios as δ¹⁵N. This may be how the results are presented, but this is not clear from the units.
Figure 2F – these are very high available P concentrations for a natural ecosystem, although there is no mention of the method used.
Figure 4 – it looks like foliar N:P ratios are around 4 (2% N, 0.5%P) – these very low values that suggest strong N limitation. This is incompatible with the very high nitrate values presented in Figure 2. This further indicates storage problems with N measurements.
Augustine, D. J. (2003). Long-term, livestock-mediated redistribution of nitrogen and phosphorus in an East African savanna. Journal of Applied Ecology, 40(1), 137-149.
Hudson et al. (2009). Temporal patterns of nutrient availability around nests of leaf-cutting ants (Attaolombica) in secondary moist tropical forest. Soil Biology and Biochemistry, 41(6), 1088-1093.
Turner, B. L., & Romero, T. E. (2009). Short-term changes in extractable inorganic nutrients during storage of tropical rain forest soils. Soil Science Society of America Journal, 73, 1972-1979.

Citation: https://doi.org/10.5194/egusphere-2024-1514-RC3
- AC4: 'Reply on RC3', Courtney Reed, 04 Sep 2024
  
  Thank you for your feedback, which was very helpful in improving the manuscript. Our responses are attached, with reviewer comments in black and author responses in blue.
  
  Citation: https://doi.org/10.5194/egusphere-2024-1514-AC4

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2024-1514', Shawn J. Leroux, 17 Jun 2024

Summary
Animals impact elemental cycling in many direct and indirect ways. Evidence from several biomes demonstrates that even after death, animal carcasses can change the biogeochemistry of ecosystems and these impacts can be long lasting. Most studies of carcass impacts on ecosystems, however, are done on small to medium (1kg to 200kg) sized animals. In this contribution, the authors investigate the effects of elephant megacarcasses on the biogeochemistry of soils and plants. The authors report significant effects of elephant carcasses on components of soil and plant elemental cycling and they discuss how these effects may be important components of spatiotemporal heterogeneity in ecosystems.
General comments
1) Overall, I found the writing good. The authors have crafted a nice narrative that makes a compelling case that megacarcasses can be important parts of ecosystems and therefore we need to learn more about the impacts of these carcasses on ecosystems.
2) I have a few questions about the analysis. The effective sample size is 10. Obviously, it is hard to find carcasses (I would have great difficulty in finding 10 fresh moose carcasses in my system!) but the authors are trying to squeeze a lot of information out of very few data points. I have the following specific questions about the analysis:
i) While I like the transect approach, the design may have been stronger if the authors had random transects (ie, transects with no known carcass) like Risch et al. work. This would strengthen inference.
ii) lines 180-183. The author’s approach to checking for normality of response data does not seem sound to me. The assumption of normality (for linear models) is normality in light of the model, i.e., investigating the normality of residuals is a more common approach to this. Either way, it is often better to avoid transforming the data and generalized linear models do allow for a lot of flexibility to fit different error distributions. For example, the gamma family in glm is very flexible and can handle log-normal data sets. Did the authors try different families of error distributions before transforming their data?
iii) lines 187-189. How many data points did the authors have per estimated parameter in the most complex model here?
iv) line 194. This is fine but I think Burham & Anderson would say that any model within deltaAIC of 2 of the null model should not be considered to be supported. In several cases, the authors interpret top models that are ranked above the null but within deltaAIC of 2 of the null as supported (e.g., lines 217-218, 218-221).
v) what R^2 are the authors reporting? In the captions of Tables S1 and S2 (thank you for providing full AIC and coefficient tables), the authors state “R^2 is the proportion of variance explained by a model”. This is unclear. These are mixed models, and the most common approach is to report the marginal R^2 and conditional R^2. Is the R^2 in these tables one of those or another pseudo R^2? This is critical for many reasons but most importantly, given the small sample size and large number of mixed-models, I would expect at least one of the models to not converge. There are many indicators when a mixed-model does not converge and one of the best is when the marginal R^2 = conditional R^2. Without having both of these pieces of information, the reader is unable to adequately assess the fit of the models. Other indicators of models not converging are coefficients estimates or errors that are very large or very small (i.e., 0 – see next comment).
vi) I am confused by the magnitude of Table S2 sodium and iron coefficients and/or the scale of reported on the y-axis of Figure 5 for these. The iron coefficients in Table S2 seem small relative to the Figure 5c ? Or am I misreading things?
3) the reporting of results could be improved. I recommend, the authors report: top ranked models (AIC + measure of independent fit like R^2). Then report effect size or relationships (coefficients). I found key statistics to be missing throughout. Statements like “Phosphate concentrations were greater in granitic soils…” would be more informative if they included the coefficient + error in parenthesis. Coefficients can be reported for the top-ranked model or from model averaged results when there are several competing models.
4) in section 3.2 I think the reader may be more interested in coefficients and confidence intervals around those relationships than p-values that are currently reported.
Specific comments:
5) I found the use of three different terms that mean similar things (nutrient flows, ecosystem processes, nutrient availability) in the introductory sentence confusing. I recommend the authors replace “nutrient availability” with “ecosystem processes” or “nutrient flows”. Surely living animals (not just carcasses) influence nutrient availability (which is just a part of a continual nutrient cycle).
6) line 83. I believe there is no “e” at the end of the citation Risch et al.
7) lines 96-111. How do these elephants die? As someone with no experience with megacarcasses, I would appreciate some insight on the causes of death. Most large herbivore deaths in my empirical systems are from predation which I assume is not the case for elephants.
8) really excellent job with clear hypotheses and nice work carrying forward these hypotheses throughout the ms – really makes the job easier for the reader.
9) lines 132-133 Why 10cm deep core? Is that mineral soil only?
10) the discussion is well done – concise and touches on all hypotheses.
11) Figure 1 is an outstanding visual!
12) in figures 2-5 I recommend the authors consider reminding the reader of the sampling resolution because the jitter of points makes it impossible to see what distances were measured below 5m.

Citation: https://doi.org/10.5194/egusphere-2024-1514-RC1
- AC1: 'Reply on RC1', Courtney Reed, 21 Jun 2024
  
  Thank you for this thoughtful review! We will work on incorporating your feedback and will post a full response once we have done so.
  
  Citation: https://doi.org/10.5194/egusphere-2024-1514-AC1
- AC2: 'Reply on RC1', Courtney Reed, 04 Sep 2024
  
  Thank you for your feedback, which was very helpful in improving the manuscript. Our responses are attached, with reviewer comments in black and author responses in blue.
  
  Citation: https://doi.org/10.5194/egusphere-2024-1514-AC2
RC2:
'Comment on egusphere-2024-1514', Sarah Keenan, 27 Jun 2024
Summary:
Reed and coauthors present a well-written study on the impacts of megacarcasses (elephants) to soil biogeochemistry after up to 2 years of decomposition. The authors examined 10 carcass hotspots with 5 carcasses each on two different soil types. They quantified soil major and trace element chemistry as well as plants associated with the hotspots to determine if carcasses influenced soil N and P chemistry and if those elements were subsequently enriched in vegetation. The current version of this manuscript does not adequately describe the methods in enough detail to make the work reproduceable. Additionally, the handling of the data for statistical analyses is strange and non-standard. The discussion needs to be re-written to better emphasize the importance of the work (as framed in the introduction). I think this work has potential to be an important contribution, but there needs to be some major revisions.

General comments:
The importance of this study in adding to our knowledge about nutrient transfer at carrion hotspots is not emphasized clearly in the discussion. The introduction frames how megacarcasses may be “functionally different than smaller carcasses” but never returns to this aspect in the discussion, which is really where this work could add to our knowledge. Adding more to the discussion would help address this issue and would make the impact of the work clearer.

Parts of the results belong in the discussion, and I’ve tried to highlight those below in specific comments.

The methods need significantly more specific details, highlighted in specific comments. Additionally, there were no control soil or plant samples examined here. Please describe in the methods why there were no controls.

The handling of the data for statistical analyses is non-standard and not clearly justified. If data were non-normally distributed (it seems like some datasets were and some were not), why not just use a non-parametric statistical test rather than log-transforming the data? It is a bit strange to log-transform some data but not all. The approach of adding 0.001 to zero values is also not correct (described below in specific comments).

The presentation of elemental data for soil and plant composition is non-standardized throughout. Some data (i.e., iron) are presented as mg/kg (is this soil dry weight?), while others are presented as % (Ca% of what?) in the same figure (figure 5 for example). Other data are presented as mg/L (figure 2). Part of this confusion is from the missing details in the methods that clearly explain how these data were generated. In several of the figures there is a statement about back-transformed data, which is also confusing.

I can appreciate that finding carcasses that have decomposed for the same amount of time is challenging, but 1 month to 26 months is a huge range of time (at least from what we know from not megacarcasses). The biogeochemical processes occurring at a carcass decaying after 1 month postmortem is very different than a carcass that has been decaying for 26 months (from smaller carcasses). It would be useful to see some of the data, particularly ammonium, plotted as a function of postmortem interval (months) even if that is not a variable that could be included in statistical analyses because of the small sample size. It would also be helpful to see if the postmortem interval for the 10 carcasses is evenly distributed between the two soil types or if one has more fresh carcasses and the other has older carcasses, that could help with interpretation of the results.

I think it may be useful to add some photos to supplemental information (or even the main text) showing what the carcasses/sites looked like (maybe representative images from a fresher carcass and one that is older).

Specific comments:
Lines 99-100: There should be more details provided on the soil type and what makes the granitic soils “nutrient poor” compared to soils developed from a basalt protolith. Because soil type becomes an important part of this study, the details of the soil types need to be expanded in the introduction.

Line 132: Include a citation or discuss why soil samples were collected to a depth of 10 cm rather than the upper 5 cm. For decomposition studies, typically the upper 5 cm is examined, not the upper 10 cm.

Line 145: More details are needed beyond “measurements of soil ion concentrations”. What instrumentation was used? What specific extraction protocol was followed? I’m assuming deionized water was used (1:2 soil to deionized water?), but those details are not provided. How long were samples mixed (shaking platform?), what speed, etc.

Line 148: “mass spectrometry”—elaborate on what this means with respect to instrumentation used to analyze cations. Here and throughout the methods, please also include what standards were used for the different analysis types.

Lines 146-150: Clarify if these analyses were conducted on the water extracts.

Line 152: Were stable isotope analyses conduct on oven-dried soil? 10 g is an exceptionally large amount of soil—how much was actually analyzed with EA-IRMS? Were samples powdered prior to combustion?

Line 154 (and throughout with respect to stable nitrogen isotope results): The authors refer to “¹⁵N” measurements, but surely this should be presented as the ratio of 15/14N and in delta notation? In the methods here there also needs to be more description of the standard, the materials used for linearity, and the analytical precision of the instrument.

Line 175: More details on the ICP-MS are needed, including standards, detection limits, etc. Additionally, were these samples digested in nitric acid? Water? How long were they microwaved?

Line 182: Adding some random number to each variable is not a standard way to handle data that are zero in your dataset (or if it is, there is no citation here and I am not familiar with that approach). Typically for geochemical data (like what was generated with ICP-MS), you can replace zero values with ½ the detection limit to remove non-zero data. There are other more technical ways to deal with zero values from a statistical standpoint, but the ½ the detection limit is the easiest and has the longest history of use. Please justify the use of your approach or re-run the analyses following a standard method for handling non-zero data in a geochemical dataset.

Line 255: The part of the sentence that reads “…we found evidence that N from carcasses had moved from soils into plants” does not belong in the results section. This is interpretation and should be moved to the discussion.

Lines 256-258: Similar comment as above where the content of this sentence is interpretation and should be moved to the discussion.

Lines 295-297: I’m not quite sure I understand the logic presented here. First, soil microbial biomass was not measured. The respiration potential (through production of CO2) was measured, but heterotrophic activity (which is how respiration can be interpreted) consumes oxygen. I think the phrasing here needs to be re-worked to not imply that the soil respiration (and the communities producing CO2) are not necessarily the same that are driving nitrification.

Lines 304-305: There are prior studies that demonstrate the impact of increased organic C during decomposition on soil microbial processes that should be cited here (see studies by DeBruyn and colleagues)

Lines 310-311: I’m not sure that this is phrased correctly. Phosphorus (predominantly as phosphate) is considered immobile in soil partly because of low solubility because it is often sorbed with Ca, Fe, Al or organics, and the release of P is tightly controlled by soil (or fluid) pH. N does not face the same sorts of sorption immobilization constraints. I think if you rephrased it to clarify that P and N are held within different reservoirs within soils that make them behave differently (and add some citations), that would help.

Lines 324-328: As mentioned above, because the composition of the two soil types were not included, this part of the discussion is not supported by the results. It’s unclear if the authors here are trying to say that the basaltic soils contain more nutrients after being impacted by decomposition or if the native state of the soils (background conditions) are more nutrient rich. I think if the introduction described the background chemistry of the two soil types this would be better supported. Additionally, basalt and granite contain different types of minerals and therefore additional sources of elements like P. I don’t know what the specific mineralogy is of these two rock types in KNP, but it might be worth exploring. In particular, the presence of apatite (Ca-P bearing mineral) in the granite might also be contributing to elevated P measured in the granite soils.

Figure 5 (and others): I’m a little bit confused by the figure caption. I think it would be better to present this as selected elements plotted as a function of distance. These are not the results of generalized linear mixed models, but the selection of how to present the data were informed by the models.
Citation: https://doi.org/10.5194/egusphere-2024-1514-RC2
- AC3: 'Reply on RC2', Courtney Reed, 04 Sep 2024
  
  Thank you for your feedback, which was very helpful in improving the manuscript. Our responses are attached, with reviewer comments in black and author responses in blue.
  
  Citation: https://doi.org/10.5194/egusphere-2024-1514-AC3
RC3:
'Comment on egusphere-2024-1514', Anonymous Referee #3, 05 Jul 2024

This manuscript presents data on the influence of elephant carcasses on nutrient availability in South African savanna soils. It would be a surprise if a decaying elephant did not increase nutrient concentrations in the proximity of the carcass, but there are some interesting differences among nutrients in terms of the distance over which the effects extend. There are parallels with other nutrient hotspots in tropical ecosystems, including glades in African savannas (e.g. Augustine 2003) and leafcutter ant nests in tropical forests (e.g. Hudson et al. 2009) - it would be worth introducing these into the discussion for comparison. I have several questions about methodology and results that should be addressed before this manuscript could be acceptable for publication.
Line 38 – these are not graves. A grave is an excavation for burial.
Line 77 – what about the amounts of cations in an elephant?
Line 89 – cations are not micronutrients. Is there direct evidence for widespread (or any) cation limitation of growth in savanna ecosystems?
Line 99 – please include classifications for the granitic and basaltic soils in one of the internationally recognized systems. In Soil Taxonomy these are presumably Inceptisols / Alfisols and Vertisols, respectively?
Line 138 – freezing soil has implications for subsequent measurements of extractable nutrients (e.g. Turner and Romero 2009). This should be mentioned here. Given the apparently very high values for some measurements (see below) I suspect that pretreatment had a major impact on results.
Line 145/146 – please explain the difference between phosphate and plant-available P. As written, it appears they were both measured in the water extracts. Most plant-available P tests are not conducted in water (e.g. Olsen, Mehlich, Bray, etc).
Line 158 – was there any inorganic C in the samples? Savanna Vertisols developed in basalt can have considerable carbonate concentrations, albeit often in subsoil. Soil pH values would help indicate this possibility – how did carcasses affect soil pH?
Line 161 – was moisture standardized prior to the incubations?
Line 182 – an alternative is to set values to ½ detection limit.
Line 188 – I understand that there were insufficient carcasses to allow inclusion of carcass age in models. However, major differences would be expected between carcasses aged 1 month vs 2.5 years. Is there any way to provide an indication of the magnitude of the age effect? How would distance effects look if young carcasses were excluded, for example?
Line 228 – It is perhaps not surprising that P concentrations showed little variation with distance, given that P was measured in water extracts (i.e. the extraction is recovering a relatively small pool of soluble P).
Line 230 – how is plant-available P defined here?
Line 286 – Soil extractable nutrients should be expressed on the basis of dry soil, not volume.
Line 286 - these very high extractable nitrogen concentrations are presumably in part a consequence of soils being frozen prior to analysis. Another factor is time between sampling and freezing - or storage prior to freezing. Please provide a statement about sample treatment prior to analysis (time from sampling to freezing, storage conditions during this time, etc, as relevant).
Line 308 – the high available P and tissue N:P ratios (see below) indicate that there is no P limitation here. This might limit the likely influence of carcasses on foliar P, as found here (i.e. foliar P is not a strong indicator of the extent to which carcasses change P availability in general).
Line 322 – water-extractable P represents a tiny proportion of the total soil P, so reliance on this procedure probably limits the possibility of detecting change in soil P.
Line 323- 330 – this text largely repeats results.
Line 340 - missing here is a discussion of the ecological consequences of the findings. What are the implications for plant and microbial ecology in savanna ecosystems?
Figure 2B, C – these values should be presented in mg/kg soil.
Figure 2B – these are extremely high nitrate concentrations, even out to 15 m. For example, 100 mg/L is equivalent to 200 mg/kg based on a 1:2 soil to water extraction ratio. Extractions done quickly after sampling and in 2 M KCl are in the range of 1-5 mg/kg. This seems to be a clear indication of storage effects.
Figure 2B, C – are these values as NO3/NH4 or on an N basis?
Figure 2D – please express stable isotope ratios as δ¹⁵N. This may be how the results are presented, but this is not clear from the units.
Figure 2F – these are very high available P concentrations for a natural ecosystem, although there is no mention of the method used.
Figure 4 – it looks like foliar N:P ratios are around 4 (2% N, 0.5%P) – these very low values that suggest strong N limitation. This is incompatible with the very high nitrate values presented in Figure 2. This further indicates storage problems with N measurements.
Augustine, D. J. (2003). Long-term, livestock-mediated redistribution of nitrogen and phosphorus in an East African savanna. Journal of Applied Ecology, 40(1), 137-149.
Hudson et al. (2009). Temporal patterns of nutrient availability around nests of leaf-cutting ants (Attaolombica) in secondary moist tropical forest. Soil Biology and Biochemistry, 41(6), 1088-1093.
Turner, B. L., & Romero, T. E. (2009). Short-term changes in extractable inorganic nutrients during storage of tropical rain forest soils. Soil Science Society of America Journal, 73, 1972-1979.

Citation: https://doi.org/10.5194/egusphere-2024-1514-RC3
- AC4: 'Reply on RC3', Courtney Reed, 04 Sep 2024
  
  Thank you for your feedback, which was very helpful in improving the manuscript. Our responses are attached, with reviewer comments in black and author responses in blue.
  
  Citation: https://doi.org/10.5194/egusphere-2024-1514-AC4

Peer review completion

AR: Author's response | RR: Referee report | ED: Editor decision | EF: Editorial file upload

ED: Reconsider after major revisions (22 Sep 2024) by Erika Buscardo

AR by Courtney Reed on behalf of the Authors (25 Oct 2024) Author's response Author's tracked changes Manuscript

ED: Referee Nomination & Report Request started (25 Oct 2024) by Erika Buscardo

RR by Anonymous Referee #3 (27 Oct 2024)

Suggestions for revision or reasons for rejection

The authors provide careful revisions to my comments, which are very much appreciated. I have a couple of remaining comments.

On the available N, the authors make a convincing case that freezing didn't affect the results greatly. However, this leaves the potential influence of storage time between sampling and freezing, which can lead to considerable changes in nitrate and ammonium concentrations. How long between sampling soil and freezing it? This is not stated in the revised description. Comparison with literature values is only helpful if those data were collected with minimal storage time, but I expect that almost all studies mentioned involve some storage time before extraction, freezing, etc. Often samples are stored for many days (fridge etc) before extraction, leading to 100+ pp concentrations of available N. This is a well-known artifact that is still not widely recognized in ecological studies. Studies that extract immediately - e.g. using in-field KCl extraction or rapidly upon return to the laboratory - tend to find very low extractable N concentrations. Anyway, as the authors mention, the samples were all treated in the same way, so the patterns with distance are presumably realistic. However, the stark contrast between the apparent strong N limitation based on foliar data and extremely high N availability based on KCl extractions is a conundrum. If someone can do in-field KCl extractions sometime at Kruger (especially seasonally?) it would be very helpful in interpreting the nutrient situation there. If such studies have already been done, it would be useful to compare them to the results presented here.

In response to my question about evidence for Ca limitation in savannas the authors state that evidence does exist. They provide citations and they have amended their manuscript to include this. However, one citation is a review from 1986 that provides no evidence for cation limitation in natural plant communities anywhere (in fact, they state that "There is virtually no direct evidence that native terrestrial ecosystems are nutrient limited"), another citation is a classic paper on uplift of nutrients by plants but does not provide direct evidence of cation limitation. A third citation concludes that productivity in woodland savanna is co-limited by N and P. Another review deals with N and P and also does not provide any evidence of cation limitation. The final paper deals with potassium dynamics in cropping systems where K is an issue because N and P are supplied as fertilizer. So, I ask again - is there any evidence for cation limitation of plant growth in natural communities on African savannas? Appropriate evidence would be an increase in productivity following addition of calcium or potassium to a savanna plant community.

Please note that orders in Soil Taxonomy are capitalized (Alfisols, Andisols, etc).

Hide

RR by Shawn J. Leroux (04 Nov 2024)

RR by Sarah Keenan (21 Nov 2024)

Suggestions for revision or reasons for rejection

General comments/major concerns:
The authors addressed some of the concerns raised in the initial review, but there are still areas that need some significant clarification and correction related to the conversion of concentrations measured in soil chemistry with the conversion to mg/kg (dry weight?). The authors describe (line 175) that water soluble nitrate and phosphate were extracted using 100 ml soil and 200 ml deionized water. I am unfamiliar with measuring soil in ml rather than mass, and I am concerned that this is not a standardized approach and does not provide a meaningful measure of those ions. This extends to the ammonium extraction—if this was also conducted as a 1:2 ratio (100 ml soil to 200 ml water), that is also a non-standardized approach. Is this a typo and this should be 100 mg of soil? This confusion extends further to the calculations for converting mg/L, which as described, that math would end up with the units of mL/mL of soil, which is not a unit I’ve seen before with soil chemistry. The concentrations should be based on a gram dry weight basis, not volume of soil. There should be some correction for dry mass of the soil.

Further, despite describing this conversion process to present ammonium concentrations in soil in a standardized way (mg/some gram dry weight basis), the authors later (lines 370-372) describe ammonium in soil as mg/L. The figures of ammonium show data as mg/kg. The values presented (assuming there was a typo and the data are actually mg/kg gram dry weight basis) are also significantly smaller than other studies with smaller carcasses (e.g., Quagiotto et al., 2019), which is the opposite of the discussed results and conclusions stated. Here, the maximum concentration of ammonium reported (based on your figure 2C) is ~150 mg/kg. The max concentration of ammonium in Quagiotto et al. (2019) is ~800 mg/kg. The opposite is true of nitrate, which your results suggest is 500-1500 mg/kg whereas in Quagiotto et al. it is ~2 mg/kg. This is exceptionally high nitrate, even compared to some of the prior data generated from Kruger and included in the author’s responses to reviewer #3. In sum, it is really difficult to figure out how the data were generated and then handled to get to the results presented for the soil geochemistry.

There is something strange going on with the soil nitrogen data (excluding the %N and nitrogen isotopes) that does not make sense from the perspective of what should be expected (high ammonium, low nitrate), what is typical of ecosystems (mentioned by reviewer #3 in the first review), and compared to other studies at Kruger.

Additionally, the phosphate seems really low, which is also maybe a conversion issue?

The concerns related to the soil geochemical data are significant and have impacts on the discussion and conclusions presented. Until these data issues are resolved, aspects of the discussion are not supported by the current results.

Additionally, one suggested change is to replace the words “post-death” with either “post-mortem” or “postmortem.” While “post-death” is not technically incorrect, the convention and terminology used more broadly is “postmortem.”

Specific comments:
- Lines 172-174: It may be better to cite an EPA or USGS protocol rather than these papers.
- Line 280: I think it’s great to evaluate change in soil chemistry as a function of carcass time since death, but could the authors explain why an exponential decay would be expected?
- Line 182: The authors state, “To determine whether soil anions were distinct…” but then list cations. At the end of this same sentence, you don’t measure cations or anions “using microwave assisted digestion.”
- Line 331: Can you clarify what is meant by the term “variance.” Is this change in mean?
- Line 366: The added text, “but much greater in magnitude in these much larger carcasses” is not supported by the data presented
- Lines 367-372: As discussed above, why is mg/L being used here. This part of the discussion is not supported by the results presented.
- Line 379: The way this sentence reads suggests that nitrate is sourced from heterotrophic microbes. Consider revising this sentence to clarify that the carbon story relates to heterotrophic microbes, but the nitrate does not necessarily inform heterotrophic microbial activity.
- Line 386-387: The units here for respiration rates are not presented on a gram dry weight basis—is this just a typo?
- Lines 442-444: Because of some of the confusing regarding the data as presented, this sentence is not supported clearly by the data. The discussion should have some more specific comparisons to smaller carcasses and their impacts on ecosystems to support this sentence (and whole paragraph).
- Line 457: I would be careful with listing “nitrate” here as being sourced from the carcasses. Ammonium is sourced from the carcasses, but nitrate forms as microbes transform the ammonium into nitrate through nitrification. The increase in nitrate only occurs later in decomposition once the system returns to being oxygenated (not examined here). The carcass itself is not releasing/introducing nitrate as the nitrogen pool. I would delete the reference to nitrate here.
- Line 462: While I agree that this type of work could help us understand the impacts of megaherbivore declines on ecosystems, I would remove the term “Anthropocene.” This is not an officially recognized term (it was rejected this year by the geological community). I would just replace that word with “at present” or “in modern ecosystems” or something similar.
- In all figures, use appropriate superscript or subscript on axis labels

Hide

ED: Reconsider after major revisions (26 Nov 2024) by Erika Buscardo

AR by Courtney Reed on behalf of the Authors (17 Dec 2024) Author's response Author's tracked changes Manuscript

ED: Referee Nomination & Report Request started (18 Dec 2024) by Erika Buscardo

ED: Publish subject to technical corrections (12 Jan 2025) by Erika Buscardo

AR by Courtney Reed on behalf of the Authors (17 Jan 2025) Manuscript

Journal article(s) based on this preprint

25 Mar 2025

Elephant megacarcasses increase local nutrient pools in African savanna soils and plants

Biogeosciences, 22, 1583–1596, https://doi.org/10.5194/bg-22-1583-2025,https://doi.org/10.5194/bg-22-1583-2025, 2025

Short summary

Supplement

https://doi.org/10.5194/egusphere-2024-1514-supplement

Viewed

Total article views: 906 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	Supplement	BibTeX	EndNote
605	165	136	906	51	32	36

HTML: 605
PDF: 165
XML: 136
Total: 906
Supplement: 51
BibTeX: 32
EndNote: 36

Views and downloads (calculated since 30 May 2024)

Month	HTML	PDF	XML	Total
May 2024	14	3	4	21
Jun 2024	293	77	13	383
Jul 2024	71	22	7	100
Aug 2024	33	12	1	46
Sep 2024	28	9	99	136
Oct 2024	14	7	2	23
Nov 2024	18	7	1	26
Dec 2024	12	10	0	22
Jan 2025	42	11	8	61
Feb 2025	68	4	1	73
Mar 2025	12	3	0	15

Cumulative views and downloads (calculated since 30 May 2024)

Month	HTML	PDF	XML	Total
May 2024	14	3	4	21
Jun 2024	293	77	13	383
Jul 2024	71	22	7	100
Aug 2024	33	12	1	46
Sep 2024	28	9	99	136
Oct 2024	14	7	2	23
Nov 2024	18	7	1	26
Dec 2024	12	10	0	22
Jan 2025	42	11	8	61
Feb 2025	68	4	1	73
Mar 2025	12	3	0	15

Viewed (geographical distribution)

Total article views: 967 (including HTML, PDF, and XML) Thereof 967 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 25 Mar 2025

Download

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Preprint (1643 KB)
Metadata XML

Short summary

We seek to understand the ecological legacies of elephants after they die. We sampled elephant carcasses in South Africa and found that they release nutrients into soil, which then enter plants and are available for consumption by other herbivores. This research reveals a key way that these elephants contribute to nutrient cycling on the savanna after death. It also highlights an important process that may be lost on savannas in areas where elephant populations are in decline.


Total:	0
HTML:	0
PDF:	0
XML:	0