Comment on “Soil carbon fraction responses to grazing intensity and texture in a semiarid grassland” (Graham et al., 2026)

Comment on “Soil carbon fraction responses to grazing intensity and texture in a semiarid grassland” (Graham et al., 2026) | Authorea try { document.documentElement.classList.add('js'); } catch (e) { } var _gaq = _gaq || []; _gaq.push(['_setAccount', 'G-8VDV14Y67G']); _gaq.push(['_trackPageview']); (function() { var ga = document.createElement('script'); ga.type = 'text/javascript'; ga.async = true; ga.src = ('https:' == document.location.protocol ? 'https://ssl' : 'http://www') + '.google-analytics.com/ga.js'; var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s); })(); Skip to main content Preprints Collections Wiley Open Research IET Open Research Ecological Society of Japan All Collections About About Authorea FAQs Contact Us Quick Search anywhere Search for preprint articles, keywords, etc. Search Search ADVANCED SEARCH SCROLL This is a preprint and has not been peer reviewed. Data may be preliminary. 27 February 2026 V1 Latest version Share on Comment on “Soil carbon fraction responses to grazing intensity and texture in a semiarid grassland” (Graham et al., 2026) Authors : Tanzeel Muzaffar 0009-0007-2138-5391 [email protected] and Luna Zayter Authors Info & Affiliations https://doi.org/10.22541/au.177222627.76243306/v1 95 views 68 downloads Contents Abstract Information & Authors Metrics & Citations View Options References Figures Tables Media Share Abstract This Letter to the Editor provides a collegial commentary on Graham et al. (2026), focusing on statistical modeling structure and soil carbon stock calculation methods relevant to grazing–carbon research. Comment on “Soil carbon fraction responses to grazing intensity and texture in a semiarid grassland” (Graham et al., 2026) To the Editor: Graham et al. , 2026 examined how over 80 years of cattle grazing at different intensities affected soil organic carbon (SOC) fractions (particulate OM – POM-C vs. mineral-associated OM – MAOM-C) in a semiarid mixed-grass prairie. They sampled 0–30 cm soils in six pastures (high, medium, low grazing) and two texture classes (clay loam vs. silty clay), then used Bayesian mixed-effects models (R brms ) to analyze carbon in 0–7.5, 7.5–15, and 15–30 cm layers. They report that grazing intensity (GI) did not significantly change total MAOM-C stocks, but that surface-layer POM-C was significantly higher under medium and high grazing in clay-loam soils (with no grazing effect in silty clay). Overall POM:MAOM ratios were higher in coarser-textured soils and under heavier grazing, implying greater accrual of labile C but also greater potential loss. Notably, POM-C and MAOM-C were only modestly coupled (r ≈ 0.3) across samples, suggesting that increases in POM did not translate into proportional increases in the more stable MAOM pools (Becker et al. , 2022). We applaud the focus on carbon fractions and Bayesian modeling, but wish to raise several technical points for consideration. First, regarding the statistical approach: the Bayesian brms framework is well suited for unbalanced designs, greater transparency in prior specification and posterior summaries would further strengthen interpretability. With only six pastures and a complex design, the choice of priors and model structure can materially affect results (McElreath, 2020). We suggest re-running a single hierarchical model with depth as a random effect (or as a factor with a common error term) to capture overall depth trends. A multilevel model could tighten uncertainty bounds and avoid multiple comparison issues. We also urge reporting of full posterior summaries (means, credible intervals) rather than just “significance”. For instance, stating that GI “did not affect” MAOM-C would be clearer if phrased in terms of overlapping 95% credible intervals or posterior probability (e.g. P ≈ 0.50). In short, sensitivity checks on priors and a hierarchical depth model (e.g. carbon stock ~ GI * Texture * Depth + (1 | Pasture) ) would strengthen confidence in the inferences. Second, the calculation of SOC stocks merits scrutiny. Graham et al. computed stocks via SOC% × bulk density × fixed depth. However, bulk density can vary with grazing and texture; fixed-depth samples therefore compare different soil masses. Note that a plausible reduction in BD (e.g. 1.5→1.1 g·cm⁻³) can cause a ~17% underestimation of SOC change under fixed-depth sampling (Fowler et al. , 2023). In other words, non-uniform BD makes fixed-depth stocks non-comparable. To avoid this bias, we recommend recalculating SOC on an equivalent soil mass (ESM) basis (Wendt and Hauser, 2013). An ESM approach would ensure all treatments reflect the same soil mass, as demonstrated by (Fowler et al. , 2023). Applying an ESM correction here might alter the reported grazing effects. For example, heavy grazing could reduce BD, meaning that the observed POM-C gain might partly reflect sampling deeper, lighter soil. Recomputing C stocks on equal-mass samples would clarify whether GI truly increased total C in clay loam soils or simply in the same sampled volume. Third, the biological interpretation of POM versus MAOM merits further nuance. The authors correctly note that POM:MAOM ratios increased with grazing and texture, inferring that POM gains are vulnerable. Indeed, MAOM is widely viewed as the more stable fraction and POM as more labile (Villarino et al. , 2023). They explicitly state “MAOM is generally considered to be a more stable fraction, whereas POM is more readily decomposable”. Thus, the higher POM:MAOM under heavy grazing suggests a transient C pool. Given the 80-year timeframe, one might question whether the system is at equilibrium: has new POM fully converted to MAOM, or could current POM still be turning over? We note that long-term grazing exclosure studies often find little change in total SOC. In fact, Derner et al. reported that even after or its turnover, and grazing effects did not interact with texture (Derner, Augustine and Frank, 2019). Fourth, we recommend explicitly testing the texture × grazing interaction. The results suggest GI effects only in clay loam, none in silty clay. Was a formal GI*Texture term included? If not, a unified model with the interaction would confirm whether the difference is statistically robust. Such an interaction might arise if, for example, silty clays already had high baseline POM and BD differences. Additionally, rather than modelling each layer separately, a multivariate model across all depths could capture covariance between depths. In summary, a single model with GI, Texture, Depth as predictors might simplify interpretation and highlight any significant interaction effects. In conclusion, these points do not call the data itself into question but suggest that some results might change under alternative analyses. For example, if an ESM correction reduces the reported clay-loam POM increase, then grazing’s carbon benefit may be smaller than stated. Likewise, a hierarchical model might reveal GI effects at other depths or tighten credibility of non-effects. Ultimately, addressing these methodological issues will clarify how confidently we can use heavy grazing to sequester grassland carbon. We appreciate the value of Graham et al. , 2026 long-term dataset and believe that refining the analysis as suggested will enhance its utility for grazing management. We would welcome further dialogue or data sharing to jointly explore these questions. Sincerely, Tanzeel Muzaffar, VIT Bhopal University, kothrikalan, Sehore 466114, India [email protected] Luna Zayter, American University of Beirut, Beirut 1107 2020, Lebanon [email protected] References Becker, A.E. et al. (2022) “Surface-soil carbon stocks greater under well-managed grazed pasture than row crops,” Soil Science Society of America Journal , 86(3), pp. 758–768.Derner, J.D., Augustine, D.J. and Frank, D.A. (2019) “Does Grazing Matter for Soil Organic Carbon Sequestration in the Western North American Great Plains?,” Ecosystems , 22(5), pp. 1088–1094. Available at: https://doi.org/10.1007/s10021-018-0324-3.Fowler, A.F. et al. (2023) “A simple soil mass correction for a more accurate determination of soil carbon stock changes,” Scientific Reports , 13(1), p. 2242. Available at: https://doi.org/10.1038/s41598-023-29289-2.Graham, C. et al. (2026) “Soil carbon fraction responses to grazing intensity and texture in a semiarid grassland,” Soil Science Society of America Journal , 90(1), p. e70184. Available at: https://doi.org/10.1002/saj2.70184.McElreath, R. (2020) Statistical Rethinking: A Bayesian Course with Examples in R and Stan . 2nd ed. Chapman and Hall/CRC. Available at: https://doi.org/10.1201/9780429029608.Villarino, S.H. et al. (2023) “A large nitrogen supply from the stable mineral-associated soil organic matter fraction,” Biology and Fertility of Soils , 59(7), pp. 833–841. Available at: https://doi.org/10.1007/s00374-023-01755-z.Wendt, J.W. and Hauser, S. (2013) “An equivalent soil mass procedure for monitoring soil organic carbon in multiple soil layers,” European Journal of Soil Science , 64(1), pp. 58–65. Available at: https://doi.org/10.1111/ejss.12002. Information & Authors Information Version history V1 Version 1 27 February 2026 Copyright This work is licensed under a Creative Commons Attribution 4.0 International License Keywords bayesian modelling equivalent soil mass grazing intensity semiarid grasslands soil organic carbon Authors Affiliations Tanzeel Muzaffar 0009-0007-2138-5391 [email protected] VIT Bhopal University View all articles by this author Luna Zayter American University of Beirut Faculty of Health Sciences View all articles by this author Metrics & Citations Metrics Article Usage 95 views 68 downloads .FvxKWukQNSOunydq8rnd { width: 100px; } Citations Download citation Tanzeel Muzaffar, Luna Zayter. Comment on “Soil carbon fraction responses to grazing intensity and texture in a semiarid grassland” (Graham et al., 2026). Authorea . 27 February 2026. DOI: https://doi.org/10.22541/au.177222627.76243306/v1 If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download. For more information or tips please see 'Downloading to a citation manager' in the Help menu . Format Please select one from the list RIS (ProCite, Reference Manager) EndNote BibTex Medlars RefWorks Direct import Tips for downloading citations document.getElementById('citMgrHelpLink').addEventListener('click', function() { popupHelp(this.href); return false; }); $(".js__slcInclude").on("change", function(e){ if ($(this).val() == 'refworks') $('#direct').prop("checked", false); $('#direct').prop("disabled", ($(this).val() == 'refworks')); }); View Options View options PDF View PDF Figures Tables Media Share Share Share article link Copy Link Copied! Copying failed. Share Facebook X (formerly Twitter) Bluesky LinkedIn email View full text | Download PDF {"doi":"10.22541/au.177222627.76243306/v1","type":"Article"} Now Reading: Share Figures Tables Close figure viewer Back to article Figure title goes here Change zoom level Go to figure location within the article Download figure Toggle share panel Toggle share panel Share Toggle information panel Toggle information panel Go to previous graphic Go to next graphic Go to previous table Go to next table All figures All tables View all material View all material xrefBack.goTo xrefBack.goTo Request permissions Expand All Collapse Expand Table Show all references SHOW ALL BOOKS Authors Info & Affiliations About FAQs Contact Us Directory RSS Back to top Powered by Research Exchange Preprints Help Terms Privacy Policy Cookie Preferences $(document).ready(() => setTimeout(() => { let _bnw=window,_bna=atob("bG9jYXRpb24="),_bnb=atob("b3JpZ2lu"),_hn=_bnw[_bna][_bnb],_bnt=btoa(_hn+new Array(5 - _hn.length % 4).join(" ")); $.get("/resource/lodash?t="+_bnt); },4000)); (function(){function c(){var b=a.contentDocument||a.contentWindow.document;if(b){var d=b.createElement('script');d.innerHTML="window.__CF$cv$params={r:'9fe45ee14c1a52ad',t:'MTc3OTIwNzAxNA=='};var a=document.createElement('script');a.src='/cdn-cgi/challenge-platform/scripts/jsd/main.js';document.getElementsByTagName('head')[0].appendChild(a);";b.getElementsByTagName('head')[0].appendChild(d)}}if(document.body){var a=document.createElement('iframe');a.height=1;a.width=1;a.style.position='absolute';a.style.top=0;a.style.left=0;a.style.border='none';a.style.visibility='hidden';document.body.appendChild(a);if('loading'!==document.readyState)c();else if(window.addEventListener)document.addEventListener('DOMContentLoaded',c);else{var e=document.onreadystatechange||function(){};document.onreadystatechange=function(b){e(b);'loading'!==document.readyState&&(document.onreadystatechange=e,c())}}}})();