Understanding the reallocation of atmospheric anthropogenic carbon (C) in the different compartments of the Earth System is a priority in Earth Science. Global numerical modeling of the C cycle stands as one of the fundamental tools for understanding how C cycles between the atmosphere, continents, and the ocean. However, Earth System Models and other large-scale models still lack a comprehensive depiction of the role of aquatic ecosystems along the Land-to-Ocean Aquatic Continuum (LOAC) in modulating organic carbon (OC) exchanges between terrestrial ecosystems and the ocean. The capacity of aquatic ecosystems to sequester organic carbon in the sediments they accumulate (i.e., organic carbon burial (OCB)) is a fundamental process for understanding the role of the LOAC in the global C cycle. Yet, the inclusion of this process into large-scale numerical models of the C cycle is still in its early stage. Here, we review the ecosystem processes involved in OCB along the LOAC and the terminology used by different authors, OCB measurement methodologies, the structure of large-scale C models, OCB rates available in the literature, and other data sources for modeling purposes. Our goal is to pinpoint the obstacles and potential solutions for incorporating OCB along the LOAC into Earth System Models and other large-scale applications. We identify the lack of language harmonization across different scientific disciplines working with ecosystems along the LOAC as a major caveat, and suggest a controlled vocabulary about OCB to assist addressing this challenge. We have compiled an updated global data set of OBC rates across ecosystems along the LOAC (lakes, reservoirs, floodplains, and coastal ecosystems), encompassing 1163 OCB rate values corresponding to 713 individual ecosystems, and showing strong biases in its distribution across the global geography and ecosystem types. We also show that virtually no existing large-scale C model incorporates OCB along the LOAC, although several have already made first steps towards the inclusion of this process at the global scale. Finally, we analyze the challenges and potential solutions to help paving the road for integrating OCB along the LOAC in large-scale models of the C cycle, including the pressing need for a multidisciplinary perspective in OCB modeling studies that brings together researchers from the several disciplines involved in the study of the ecosystems pertaining to the LOAC.

This project has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grand agreement No. 956623. N.C. is supported by the “Ramón y Cajal” fellowship RYC-2021-033714-I, and RM and BO acknowledge funding by projects PID2020-114024GB-C32 and PID2020-114024GB-C31, all of them funded by the Ministry of Science and Innovation of Spain. Authors also acknowledge the support from the Economy and Knowledge Department of the Catalan Government through Consolidated Research Groups (ICRA-ENV 2021 SGR 01282, and FORESTREAM 2021-SGR-00530).

Peer reviewed