Abstract. Rivers are a major source of nutrients, carbon and alkalinity to the global ocean. In this study, we firstly estimate pre-industrial riverine loads of nutrients, carbon and alkalinity based on a hierarchy of weathering and terrestrial organic matter export models, while identifying regional hotspots of the riverine exports. Secondly, we implement the riverine loads into a global ocean biogeochemical model to describe their implications for oceanic nutrient concentrations, net primary production (NPP) and air–sea CO2 fluxes globally, as well as in an analysis of coastal regions. Thirdly, we quantitatively assess the terrestrial origins and the long-term fate of riverine carbon in the ocean. We quantify annual bioavailable pre-industrial riverine loads of 3.7 Tg P, 27 Tg N, 158 Tg Si and 603 Tg C delivered to the ocean globally. We thereby identify the tropical Atlantic catchments (20 % of global C), Arctic rivers (9 % of global C) and Southeast Asian rivers (15 % of global C) as dominant suppliers of carbon for the ocean. The riverine exports lead to a simulated net global oceanic CO2 source of 231 Tg C yr−1 to the atmosphere, which is mainly caused by inorganic carbon (source of 183 Tg C yr−1) and by organic carbon (source of 128 Tg C yr−1) riverine loads. Additionally, a sink of 80 Tg C yr−1 is caused by the enhancement of the biological carbon uptake from dissolved inorganic nutrient inputs from rivers and the resulting alkalinity production. While large outgassing fluxes are simulated mostly in proximity to major river mouths, substantial outgassing fluxes can be found further offshore, most prominently in the tropical Atlantic. Furthermore, we find evidence for the interhemispheric transfer of carbon in the model; we detect a larger relative outgassing flux (49 % of global riverine-induced outgassing) in the Southern Hemisphere in comparison to the hemisphere's relative riverine inputs (33 % of global C inputs), as well as an outgassing flux of 17 Tg C yr−1 in the Southern Ocean. The addition of riverine loads in the model leads to a strong NPP increase in the tropical west Atlantic, Bay of Bengal and the East China Sea (+166 %, +377 % and +71 %, respectively). On the light-limited Arctic shelves, the NPP is not strongly sensitive to riverine loads, but the CO2 flux is strongly altered regionally due to substantial dissolved inorganic and organic carbon supplies to the region. While our study confirms that the ocean circulation remains the main driver for biogeochemical distributions in the open ocean, it reveals the necessity to consider riverine inputs for the representation of heterogeneous features in the coastal ocean and to represent riverine-induced pre-industrial carbon outgassing in the ocean. It also underlines the need to consider long-term CO2 sources from volcanic and shale oxidation fluxes in order to close the framework's atmospheric carbon budget.

The contribution of continental shelves to the marine carbon cycle is still poorly understood. Their preindustrial state is, for one, essentially unknown, which strongly limits the quantitative assessment of their anthropogenic perturbation. To date, approaches developed to investigate and quantify carbon fluxes on continental shelves have strongly simplified their physical and biogeochemical features. In this study, we enhance the global ocean biogeochemistry model HAMburg Ocean Carbon Cycle by explicitly representing riverine loads of carbon and nutrients, as well as improving the representation of organic matter dynamics in the coastal ocean. Our simulations, at a resolution of ∼0.4°, reveal a globally averaged shelf water residence time (RT) of 12–17 months, which is much shorter than the global RTs previously assumed in benchmark studies (>4 years). This shorter global RT, induced primarily through outer shelf regions with large oceanic inflows, promotes an efficient offshore transport of terrestrial and marine organic carbon (0.44 PgCyr−1) and a dissolved inorganic carbon sink from the organic cycling of carbon on the global shelf (net ecosystem productivity [NEP] equal to +0.20 PgCyr−1). In turn, this autotrophic state of continental shelves contributes to a weak global preindustrial sink of atmospheric CO2 (0.04 PgCyr−1), dominated by extensive regions with large oceanic inflows and positive NEPs, such as the Patagonian shelf, the East China Sea and the outer North Sea. The contemporary global shelf CO2 uptake of 0.15 PgCyr−1 furthermoresuggests that the anthropogenic CO2 uptake (0.11 PgCyr−1) on the global continental shelf is less efficient with respect to the open ocean.