Contrasted Chemical Weathering Rates in Cratonic Basins: The Ogooué and Mbei Rivers, Western Central Africa

Despite the absence of tectonic activity, cratonic environments are characterized by strongly variable, and in places significant, rock weathering rates. This is shown here through an exploration of the weathering rates in two inter-tropical river basins from the Atlantic Central Africa: the Ogooué and Mbei River basins, Gabon. We analyzed the elemental and strontium isotope composition of 24 water samples collected throughout these basins. Based on the determination of the major element sources we estimate that the Ogooué and Mbei rivers total dissolved solids (TDS) mainly derive from silicate chemical weathering. The chemical composition of the dissolved load and the area-normalized solute fluxes at the outlet of the Ogooué are similar to those of other West African rivers (e.g., Niger, Nyong, or Congo). However, chemical weathering rates (TZsil+ rate expressed as the release rate of the sum of cations by silicate chemical weathering) span the entire range of chemical weathering intensities hitherto recorded in worldwide cratonic environments. In the Ogooué-Mbei systems, three regions can be distinguished: (i) the Eastern sub-basins draining the Plateaux Batéké underlain by quartz-rich sandstones exhibit the lowest TZsil+ rates, (ii) the Northern sub-basins and the Mbei sub-basins, which drain the southern edge of the tectonically quiescent South Cameroon Plateau, show intermediate TZsil+ rates and (iii) the Southern sub-basins characterized by steeper slopes record the highest TZsil+ rates. In region (ii), higher DOC concentrations are associated with enrichment of elements expected to form insoluble hydrolysates in natural waters (e.g., Fe, Al, Th, REEs) suggesting enhanced transport of these elements in the colloidal phase. In region (iii), we suggest that a combination of mantle-induced dynamic uplift and lithospheric destabilization affecting the rim of the Congo Cuvette induces slow base level lowering thereby enhancing soil erosion, exhumation of fresh primary minerals, and thus weathering rates. The study points out that erosion of lateritic covers in cratonic areas can significantly enhance chemical weathering rates by bringing fresh minerals in contact with meteoric water. The heterogeneity of weathering rates amongst cratonic regions thus need to be considered for reconstructing the global, long-term carbon cycle and its control on Earth climate.

Despite the absence of tectonic activity, cratonic environments are characterized by strongly variable, and in places significant, rock weathering rates. This is shown here through an exploration of the weathering rates in two inter-tropical river basins from the Atlantic Central Africa: the Ogooué and Mbei River basins, Gabon. We analyzed the elemental and strontium isotope composition of 24 water samples collected throughout these basins. Based on the determination of the major element sources we estimate that the Ogooué and Mbei rivers total dissolved solids (TDS) mainly derive from silicate chemical weathering. The chemical composition of the dissolved load and the area-normalized solute fluxes at the outlet of the Ogooué are similar to those of other West African rivers (e.g., Niger, Nyong, or Congo). However, chemical weathering rates (TZ + sil rate expressed as the release rate of the sum of cations by silicate chemical weathering) span the entire range of chemical weathering intensities hitherto recorded in worldwide cratonic environments. In the Ogooué-Mbei systems, three regions can be distinguished: (i) the Eastern sub-basins draining the Plateaux Batéké underlain by quartz-rich sandstones exhibit the lowest TZ + sil rates, (ii) the Northern sub-basins and the Mbei sub-basins, which drain the southern edge of the tectonically quiescent South Cameroon Plateau, show intermediate TZ + sil rates and (iii) the Southern sub-basins characterized by steeper slopes record the highest TZ + sil rates. In region (ii), higher DOC concentrations are associated with enrichment of elements expected to form insoluble hydrolysates in natural waters (e.g., Fe, Al, Th, REEs) suggesting enhanced transport of these elements in the colloidal phase. In region (iii), we suggest that a combination of mantle-induced dynamic uplift and lithospheric destabilization affecting the rim of the Congo Cuvette induces slow base level lowering thereby enhancing soil erosion, exhumation of fresh primary minerals, and thus weathering rates. The study points out that erosion of lateritic covers in cratonic areas can significantly enhance chemical weathering rates by bringing fresh minerals in contact with meteoric water. The heterogeneity of weathering rates amongst cratonic regions thus need to be considered for reconstructing the global, long-term carbon cycle and its control on Earth climate.

INTRODUCTION
Over geological time scales, chemical weathering acts as a major player of the global biogeochemical cycles of elements in the Earth's Critical Zone. In particular, silicate weathering is known to consume CO 2 through mineral hydrolysis and neutralization of base cations hosted by silicate minerals, and therefore can contribute to regulating the global climate (Berner et al., 1983). Tectonic activity is thought to be a primary driver of chemical weathering (Raymo and Ruddiman, 1992;Herman et al., 2013). Orogenic uplift associated with collision between tectonic plates forms the major world mountain chains, favoring mechanical erosion, and the exposure of fresh primary mineral surfaces to meteoric water and subsequent chemical weathering. This phenomenon has long been suggested to be responsible for the gradual cooling of the Earth over the Cenozoic under the effect of global mountain uplift (Herman et al., 2013;Becker et al., 2018), although this hypothesis is strongly debated (e.g., Godderis, 2010;Willenbring and Von Blanckenburg, 2010;Von Blanckenburg et al., 2015;Norton and Schlunegger, 2017;Caves Rugenstein et al., 2019;Hilton and West, 2020;Penman et al., 2020). By contrast, tectonically quiescent cratonic areas, which represent almost 70% of the continents surface (Artyushkov et al., 2018), have long been considered to be relatively inefficient in terms of chemical weathering compared to erosive, mountainous regions (e.g., Carson and Kirkby, 1972;Stallard, 1985). Indeed, especially in the humid tropics, these low-relief settings favor the formation of deep regolith covers, chemically depleted in base cations and limiting water-bedrock interactions due to slow water percolation from the surface to the bedrock (e.g., Stallard and Edmond, 1987;Braun et al., 2005Braun et al., , 2012West, 2012;Riebe et al., 2017). While these hot and humid cratonic areas commonly exhibit low silicate weathering rates by comparison with orogenic areas (e.g., Gaillardet et al., 1999b;Moon et al., 2014) they dominate the intertropical regions surface area and, therefore, represent a significant proportion of the global delivery of dissolved matter to the oceans (Milliman and Farnsworth, 2011;Von Blanckenburg et al., 2015). In particular, according to modeling results (Goddéris et al., 2008) the net weathering budget of the intertropical cratonic areas and their role on the long term CO 2 budget may have been underestimated and need to be investigated in more detail to characterize their potential role in the Cenozoic global climate evolution.
River hydrochemical analyses are essential tools for estimating catchment-scale silicate weathering fluxes, through the quantification of the export of silicate-derived dissolved cations (Gaillardet et al., 1999b). River basins draining cratonic areas exhibit contrasted weathering rates (expressed as the drainage-area normalized flux of the sum of cations released by silicate chemical weathering, TZ + sil rate) from those recorded under boreal conditions (TZ + sil rate < 1 t km −2 yr −1 ; Millot et al., 2002;Zakharova et al., 2005Zakharova et al., , 2007Pokrovsky et al., 2015) to those measured in India in the Kavery (Pattanaik et al., 2013) and Nethravati (Gurumurthy et al., 2012) basins (TZ + sil rate > 15 t km −2 yr −1 ). Various drivers can be invoked to explain this spatial variability. For example, the Nsimi experimental watershed is located in a tectonically quiescent area of the Nyong river basin (South Cameroon Plateau) where swamp environments are widespread. There, organic-rich waters increase the mobilization and transfer of some elements generally considered as immobile (e.g., Al, Fe, Th, Zr) through colloidal transport (Oliva et al., 1999;Viers et al., 2000;Braun et al., 2005Braun et al., , 2012. Such enhancement of weathering in the presence of organic matter has also been reported for boreal Siberian rivers rich in dissolved organic matter and poor in suspended matter (Zakharova et al., 2005(Zakharova et al., , 2007Pokrovsky et al., 2015Pokrovsky et al., , 2016. Indeed, in this type of environment such environments, organo-metal complexes form, leading to the solubilization of Al, Fe, Th, and Zr and thus to the breakdown of silicate minerals (e.g., Oliva et al., 1999;Tamrat et al., 2019). The Kaveri and Nethravati basins drain the Indian craton and exhibit low annual runoff (<220 mm yr −1 ) by comparison with other tropical cratonic areas. In this region, intense monsoons can enhance weathering through strong erosional processes, exposing of felsic granulites and gneissic rocks (Pattanaik et al., 2013;Meunier et al., 2015). In this case, climate and erosion act as dominant drivers of weathering rates. Finally, in the small monolithologic basins of the Mule Hole-India tropical watershed, the presence of minor/accessory minerals (Ca-bearing minerals like epidote and apatite), the dissolution of smectite and calcite, as well as the drainage characteristics of weathering profiles have all been shown to play a key role on weathering budget (Braun et al., 2009;Violette et al., 2010). However, most of the aforementioned work has focused on a local scale (i.e., soil profiles or small watersheds) or, when dealing with weathering fluxes measured at a larger scale, the variability of weathering rates amongst cratonic environments was not considered to be part of the scope of the study.
The existence of sustained uplift due to mantle dynamics or lithospheric destabilization (Cottrell et al., 2004;Jaupart et al., 2007;Hu et al., 2018) in cratonic areas has not been considered yet as a potential driver of Earth denudation. Although much slower and occurring over larger spatial scales than mountain uplift mediated by faulting in collisional tectonic settings (Lamb and Watts, 2010;Flament et al., 2013), mantleinduced dynamic uplift or lithospheric destabilization can lead to lowering of the geomorphological base level that in turn can trigger physical erosion processes (Kusky et al., 2014) as observed for the Southern African Craton (Braun et al., 2014), North China Craton (Zhu et al., 2017) or Brazilian Shield (Rodríguez Tribaldos et al., 2017). Such slow, large-scale mantle and lithospheric dynamics, in conjunction to eustatic changes, have also been shown to control relief as well as erosion and sedimentary processes (Conrad and Husson, 2009;Guillocheau et al., 2015). Presumably, physical erosion processes sustained by mantle induced dynamic uplift or lithospheric destabilization in cratonic settings could lead to significant rock weathering through soil erosion and subsequent increased exposure of "fresh" mineral surfaces-a phenomenon we set out to address in the present study.
The Ogooué River Basin, western Central Africa, is located between the western border of the Congo cuvette and the south of the Cameroun Plateau. This intra-cratonic basin experiences a homogeneous tropical humid climate (Bogning et al., 2018(Bogning et al., , 2020 and has undergone successive and contrasted uplift phases over the Cenozoic (Guillocheau et al., 2015). The present study provides the opportunity to explore the variability in weathering fluxes and rates in the large Ogooué River Basin (drainage area of ∼ 215,000 km²) as well in the neighboring, smaller Mbei River Basin (drainage area of ∼ 1,800 km 2 ) with respect to geomorphology, tectonics, and lithology. Based on discrete hydrochemical analyses of the main Ogooué River tributaries, we assess, for the first time, the Ogooué weathering fluxes and their variability throughout the basin. These new constraints allow us to explore the main potential drivers controlling the variability in weathering rates in cratonic areas and to discuss the implications for the long-term evolution of the Earth's climate.

STUDY AREA
The equatorial Ogooué River basin covers ∼215,000 km² (location: between 3 • 9'S and 2 • 4'N and between 8 • 5'E and 14 • 3E; Figure 1). Around 85% of the basin lies within Gabon, 12% in the Republic of Congo and the remaining area in Cameroon and Equatorial Guinea. With an annual discharge of 4,750 m 3 s −1 (Bogning et al., 2018), the Ogooué is the third largest river in terms of annual discharge along the African West Coast after the Congo (∼41,000 m 3 s −1 ; Laraque et al., 2009) and the Niger (∼6,000 m 3 s −1 ) rivers (Dai and Trenberth, 2002). In addition to the Ogooué Basin, this study reports on the hydrochemistry of the smaller Mbei River (location: between 0 • 1'N and 1 • 1'N and between 10 • 7'E and 10 • 4E; Figure 1), which is a northern tributary of the Komo River discharging to the Atlantic Ocean around 150 km north of the Ogooué River outlet. The Mbei River is characterized by an annual discharge of ∼60 m 3 s −1 (ORSTOM, 1964;Njutapvoui Fokouop, 2017; data only available for the period 1964-1973) and a drainage area of ∼ 1,800 km 2 . According to the Köppen-Geiger classification, the Ogooué and Mbei basins experience a tropical savanna climate. The Ogooué Basin receives around 2,000 mm yr −1 in annual precipitation, which leads to a runoff of around 700 mm yr −1 (Mahe et al., 2013;Bogning et al., 2018;Kittel et al., 2018). The basin exhibits a bi-modal precipitation regime with wet periods from March to May and from October to December with a maximal monthly rainfall of ∼200 and ∼300 mm month −1 , respectively. The dry period extends from June to August with a minimum monthly rainfall <15 mm month −1 in July. The mean annual temperature is around 24 • C and is relatively invariant across the year. In the present work, the study area corresponds to the Ogooué Basin upstream from the Lambaréné station ( Figure 2F) and covers 206,000 km², representing 96% of the entire Ogooué Basin (Figure 1), combined to four Mbei tributaries each covering basin areas <500 km². At the Lambaréné station, the Ogooué River discharge variation follows the rainfall regime of the basin (Mahe et al., 1990). Overall, the seasonal variability in discharge (SV as quantified by the ratio between maximum and minimum monthly discharge; SV = 3.7; Figure 3) is relatively low by comparison with other intertropical rivers experiencing a monsoonal regime (e.g., Pacific Peruvian rivers: SV = 6-21; Moquet et al., 2018;Nethravati River: SV = 600-1,100;Gurumurthy et al., 2012). The elevation of the Ogooué basin, as considered in the present study, extends from 914 m.a.s.l (Lolo River sub-basin upstream) to 20 m.a.s.l. (Lambaréné station).
The Ogooué Basin can be separated into three geomorphological domains (Figure 1; Table 1). The Plateaux Batéké region (Eastern sub-basins) located close to the Congo border is composed of Cenozoic sandstone (pure quartz; Seranne et al., 2008) and is active in terms of dynamic uplift (Guillocheau et al., 2015;Weber et al., 2016). The rest of the basin is underlain by the Archean basement of the Congo craton and by remnants of Proterozoic orogenic belts. The rivers of the Northern sub-basins exclusively drain granitic rocks of the Archean craton whereas the Southern sub-basins are to a significant extent underlain by Paleoproterozoic sedimentary rocks in addition to the typical Archean cratonic rocks (gneiss and undifferencied orthogneiss; Thiéblemont et al., 2009). Some sub-basins in the Northern part (Missanga basin) and in the Southern part (Leyou basin) partly drain Paleoproterozoic volcanic and volcanic sedimentary rocks (Thiéblemont et al., 2009). The southernmost area is characterized by the presence of the Ogoulou-Offoué and Ikoy-Ikobé Faults which produce steeper relief, particularly in the Lopé sub-basin (Weber et al., 2016). The central Ogooué channel flows over Quaternary fluvial sediments laying on the Paleoproterozoic Francevilian metasediments (1,800-1,600 Ma- Seranne et al., 2008;Weber et al., 2016). The Mbei basin is essentially underlain by the Archean craton. Based on Nd radiogenic isotopes measurements, Thiéblemont et al. (2014) show that the whole Ogooué Basin is covered by a homogenous 1-3 meter-thick sandy to clayey lateritic surface cover (dated at around 3,000-2,000 years BP) named "Cover Horizon" inherited from silt-to clay-sized aeolian deposits transported from the Northern part of the Namibia desert.
The study area is essentially covered by rainforest with patchy savanna grassland and is home to a high biodiversity (e.g., Koffi et al., 2011), but has shifted following climate changes over the last 26,000 years (Kim et al., 2010;Oslisly et al., 2013). Approximately 650,000 people live in the Ogooué Basin resulting in a low population density of <3 FIGURE 1 | Study area and sampling locations. The regional map of the study area is modified from Guillocheau et al. (2015). inhabitants km −2 , mostly concentrated along the river courses (UNEP, 2010). Uranium and manganese mining activity is reported in the Southern Upper Ogooué area and future dam development may potentially affect river connectivity (Cutler et al., 2020). However, overall the contemporary anthropogenic influence on the Ogooué River system can be assumed as insignificant in terms of chemistry, water flow, and sediment fluxes.  (Bogning et al., 2018). See photo of the Ogooué River at Lambaréné City ( Figure 2F).

Sampling and in-situ Analyses
Surface water samples were collected at 24 locations in the Ogooué main channel (six samples), Ogooué's tributaries (14 samples), and in the Mbei River tributaries (four samples) in September 2017 (Figures 1, 2; Table 1).
Electrical conductivity (normalized to a temperature of 25 • C), pH and water temperature were measured in situ. We followed the protocol of Bouchez et al. (2012) for water sampling and filtration. Briefly, 10 L of surface water were sampled with a "plastic bag" at the middle of the river course, when possible, or from the bank in a zone of significant water velocity. The plastic bags were filled with a graduated, 5 L bucket and a funnel, both previously rinsed with water from the river being sampled. The bags were rinsed three times with riverine water prior to sampling. The 10 L water samples were transported in a cool box to limit their warming and were filtered 2-10 h later through PES (Millipore) 0.22-µm porosity, pre-weighted filter membranes (142 mm diameter). The filtration unit was filled three times with ∼1 L of water sample, each of which was discarded after filtration to pre-condition the filter unit and the membrane. Gloves were used to avoid contamination during the filtration process. The remaining filtered water was splitted into three aliquots: one unacidified aliquot of 30 mL for major anion analyses, one acidified (using three drops of mono-distilled nitric acid to reach a pH ∼ 2) aliquot of 1L for cation, trace element, and 87 Sr/ 86 Sr analyses, and one acidified (using HCl) aliquot of 60 mL stored in an amber glass tube for dissolved organic carbon (DOC) analysis. Alkalinity was determined on the field using the Gran method (Gran, 1950). Note that the filtered water samples (<0.22 µm) are considered to reflect primarily a "truly dissolved" fraction, but are likely to include organic and inorganic colloids, depending on the element considered. This is further discussed in section Dynamics of Organic Matter and Colloids.

Analytical Methods
All laboratory procedures and measurements were performed at the PARI (Plateau d'Analyse Haute Résolution) analytical platform of the IPGP (Institut Physique du Globe de Paris). Major cations (Ca 2+ , Na + , Mg 2+ , K + ) and major anions (Cl − , SO 2− 4 , and NO − 3 ) concentration were determined using ion chromatography (CS16 cationic colomn and AS9HC anionic column-IC5000+DIONEX THERMO). Trace elements were analyzed by Quadrupole-ICP-MS (ICP-QMS 7900 Agilent). For major element concentration, the analytical uncertainty is around 5% (95% confidence interval), and for most trace elements, the analytical uncertainty is <5%. Accuracy was checked using repeated measurements of the river water reference material SLRS-6 (St-Lawrence River water, National Research Council of Canada), with concentration measurements being mostly <10% away from the certified concentrations. DOC concentrations was quantified using a Shimadzu TOC-V CSH with a relative analytical uncertainty <5%. Major element and DOC concentrations are reported in Table 2 and trace element concentration are reported in Table 3.

Isotopic Measurements
The chemical purification of Sr (∼200 ng) was performed using extraction chromatography (Sr-SPEC resin; Eichrom) before isotope analysis according to Hajj et al. (2017) with a total blank of <0.5 ng along the Sr separation and MC-ICP-MS measurements procedures. The Sr isotopic composition was measured on a 50 ppb Sr solution within a 0.5 N HNO 3 matrix by MC-ICP-MS (Neptune, ThermoFisher Scientific). Each sample was analyzed three times. The instrumental mass fractionation was corrected for using a 86 Sr/ 88 Sr ratio of 0.1194. An average internal precision of ± 10 ppm (2σ) was obtained and the quality of the 87 Sr/ 86 Sr ratio measurements was verified by repeated analysis of the NIST SRM987 standard (n = 21; measured value: 87 Sr/ 86 Sr = 0.710320 ± 0.000015; certified range: 87 Sr/ 86 Sr = 0.71034 ± 0.00026; accepted NIST SRM987 value: 87 Sr/ 86 Sr = 0.71026 ± 0.00002; Stein et al., 1997;Swoboda et al., 2008). Uncertainty on the 87 Sr/ 86 Sr measurements are lower than 0.000043 (95% confidence interval on the three replicates; Table 2).

Discrimination of Solute Sources
To discriminate between the different river solutes sources we use a "forward method" (Garrels and MacKenzie, 1972;Moquet et al., 2011), consisting in the application of two calculation steps: atmospheric correction and quantification of the relative roles of silicate and carbonate weathering.

Atmospheric Inputs
In the studied basins, the primary source of Cl − is precipitation as evaporite rocks are virtually absent from the Ogooué and Mbei basins (Thiéblemont et al., 2009) and because anthropogenic inputs can be considered as negligible (see section Study Area).  TSS, DOC, TDS, and NICB stand for "total suspended solid", "dissolved organic carbon," "total dissolved solid," and "normalized inorganic charge balance," respectively. Uncertainty on major elements is <5%.
Frontiers in Water | www.frontiersin.org    Therefore, for all water samples, we can apply the following formula to estimate the atmospheric contribution of each solute (X) concentration ([X rain ] in µmoles l −1 ): with X= SO 2− 4 , Na + , Ca 2+ , Mg 2+ , K + , and Sr 2+ , [Cl − riv ] the total Cl − concentration in the river, and X Cl − seawater the X/Cl − molar ratio of seawater, considered here as the sole source of ions to the rain (Berner and Berner, 1987). Then the concentration of the element X corrected from the atmospheric inputs (annotated [X * ]) is calculated as: Thereafter, the " * " symbol stands for concentrations corrected from atmospheric inputs.

Contribution of Silicate and Carbonate Weathering
After correction from rainfall inputs, the dissolved load of the rivers is considered to be the result of weathering of silicate and carbonate minerals. The quantitative estimation of the sum of the cations concentrations delivered by silicate weathering ([TZ + sil ] in mg l −1 ) is: with [X sil ] the concentration of the cation X (with X = Na + , Ca 2+ , Mg 2+ and K + ) derived from silicate weathering (here in µmoles l −1 ) and M C the molar mass of the corresponding cation C (in g mol −1 ). We consider that all the K + and the Na + remaining after correction from atmospheric inputs is derived from silicate weathering only: The concentrations of Ca 2+ and Mg 2+ derived from silicate weathering can then be calculated as: As the Ogooué Basin drains silicate rocks similar to those of the Congo basin (cratonic plutonic and metamorphic lithology) to estimate [Ca sil ] and [Mg sil ] and thus [TZ + sil ] we used the value of (Ca/Na) sil = 0.35 ± 0.15 and (Mg/Na) sil = 0.24 ± 0.12 (mol/mol) determined by Négrel et al. (1993) and used for the global scale by (Gaillardet et al., 1999b). Note that these values were previously determined based on water geochemistry of rivers draining only the corresponding rock types, and not the rocks themselves, such that this method does not rely on any particular assumption regarding the congruent or incongruent character of weathering reactions.
Without alternative proton sources such as pyrite oxidation (e.g., Calmels et al., 2007;Yu et al., 2020), the main proton sources for mineral hydrolysis is atmospheric / soil CO 2 dissolution in water. Therefore, the estimation of the CO 2 consumption associated with silicate weathering (CO 2 sil in µmoles l −1 ) can be calculated following the equation: After estimating the contribution of silicate weathering to the river dissolved load, the remaining dissolved Ca 2+ and Mg 2+ is attributed to carbonate weathering, the result in terms of total concentration ([TZ + carb ] in mg l −1 ) being therefore: From the metrics mentioned above, we can finally calculate the three components of TDS concentration (mg l −1 ; Figure 5): The fluxes (F) and area-normalized fluxes (F spe ) were calculated by multiplying the concentrations of TDS, TZ + sil , and CO 2 sil by the discharge (Q) and the drainage area-normalized discharge (Q spe also named specific discharge), respectively (see the section Hydrological and Climate Data about the calculation of Q and Q spe ). For silicate weathering we also calculated the chemical denudation (D chem sil ) expressed in m Ma −1 using the concentration of solutes (Ca sil , Mg sil , K sil , Na sil , and SiO 2 in mg l −1 ), expressed as equivalent oxides (CaO, MgO, K 2 O, Na 2 O, SiO 2 in mg l −1 ) and using a rock density (d) of 2.7 g cm −3 (West et al., 2005;Bouchez and Gaillardet, 2014): 15) with M O the molar mass of O. This D chem sil parameter expresses the rate at which silicate weathering processes result in a lowering of the Earth surface.

Hydrological and Climate Data
Drainage areas and mean slopes upstream of the sampling points were extracted from the digital elevation model SRTM 90 (Shuttle Radar Topography Mission; NASA) using ArgGis 10.3 (Esri) and QGis 2.18. The lithological composition for each sub-basin was also extracted from the lithological map of Gabon using ArcGis 10.3 (Esri) (Thiéblemont et al., 2009; see Supplementary Table 1).
The daily water discharge of the Ogooué River is available only for the Lambaréné station from Bogning et al. (2018) for the period from July 2001 to August 2017 (Figure 3), meaning that the sampling campaign was performed just after the end of the precipitation record used here. Note that according to the discharge record, which does includes the period where water sampling was carried out (September 2017), the discharge pattern of the hydrological year October 2016-September 2017 is very close to the mean pattern observed for the previous 16 years. Therefore, hydrological and precipitation conditions of the year of water sampling can be considered as representative of the longer period used to constrain hydrological fluxes. We calculated an annual discharge over this period to assess an annual solute flux based on the discrete sampling performed in this study.
To calculate the annual discharge of other sampled sites we apply a statistical approach based on a regional polynomial relationship between specific discharge (Q spe ) and rainfall (P) as performed by Scherler et al. (2017) to estimate Q spe in ungauged Himalayan rivers. We first calculated the annual mean rainfall received by each studied sub-basin from the TRMM dataset (extracted for the Ogooué Basin over the period 1998-2015 according to the TRMM data-https://gpm.nasa.gov/). The calculated rainfall for these basins ranges from 1,890 to 2,692 mm.yr −1 . At the Lambaréné station, P = 2,159 mm yr −1 and Q spe is calculated according to: with Q the discharge (136 10 9 m 3 yr −1 ) and S the basin area (205.9 10 3 km²). At Lambaréné, Q spe is therefore equal to 666 mm yr −1 (Table 1).
Second, we compiled a new database (Supplementary Figure 2; Supplementary Table 3) for both specific discharge (Q spe ) and rainfall (P) for rivers of western Central Africa (Congo: Becker et al., 2018;Laraque et al., 2020; Cameroun rivers: BVET observatory database/https:// mtropics.obs-mip.fr/catalog-m-tropics and Sighomnou, 2004); Atlantic African rivers: (Lienou et al., 2008;Conway et al., 2009). We then fitted a second-order polynomial relationship to the P-Q spe data from rivers submitted to P between 1,500 mm.yr −1 and 3,000 mm yr −1 (n = 32), corresponding to the conditions measured in the Ogooué and Mbei rivers. The equation of the best fit second-order polynomial is: (17) We then applied this relationship to the P-values estimated for each sub-basin of the Ogooué and Mbei basins. We used the RMSE of the fit (161 mm yr −1 ) as a measure of the uncertainty on these Q spe estimates (see Supplementary Figure 2), and propagated this uncertainty in the solute flux calculations (see section Solute Flux Calculation). Note that for the Ogooué River at Lambaréné station, for a rainfall of 2,159 mm yr −1 , the simulated Q spe is 780 mm yr −1 while the measured Q spe was 666 mm yr −1 . The difference between simulated and measured value is therefore 113 mm yr −1 (12%) and is lower than the Q spe uncertainty considered (161 mm yr −1 ).

Solute Flux Calculation
The area-normalized fluxes (hereafter called F spe and "rates" when referring to weathering variables) of each solute was estimated by multiplying the concentration of each solute parameters values by the Q spe value estimated at each sampling sites. To compute the area-normalized fluxes for each subregion of the study area, we substracted the upstream fluxes (F) where necessary (i.e., samples number 17, 20, 26, 8, 6 along the Ogooué River and sample number 10 on the Okano River).
Our river hydrochemistry dataset features only one sampling date. Therefore, the computed solute fluxes might be affected by significant uncertainty if solute concentrations were to vary along the year. We first note that in other West African rivers like the Nyong (Viers et al., 2000), the Niger (Picouet et al., 2002) and the Congo (Laraque et al., 2009) rivers, the TDS concentration varies only slightly with discharge throughout the year. This chemostatic behavior is observed in numerous rivers in the world from small to large catchments (e.g., Godsey et al., 2009), enabling first-order estimations of annual river dissolved fluxes based on a limited sampling time resolution. Indeed, under such conditions, solute fluxes are mainly controlled by discharge variability rather than by concentration variability, thereby providing confidence to our flux estimates based on only one sampling campaign. However, a systematic bias might exist in our flux estimates since our sole sampling campaign was performed during the dry season, when the highest river solute concentrations are usually measured. Existing hydrochemical time series on the neighboring Congo River helps assessing this bias. Using only the hydrochemical data corresponding to the lowest monthly discharge (August) for the Congo River at Brazzaville (HYBAM database, http://so-hybam.org/), TDS and cations fluxes are +10% and +20%, respectively, higher than flux estimates based on monthly data covering the whole hydrological cycle. However, while informative for assessing the reliability of our flux estimates for the Ogooué and Mbei rivers, we did not consider that our data can be corrected for this bias using constraints from the Congo River, as the specifics of solute concentration-discharge relationships in rivers can be influenced by a variety of processes (Chorover et al., 2017). In particular these relationships can depend on reactive transport processes (e.g., Ameli et al., 2017;Kim et al., 2017), biological effects (e.g., Moatar et al., 2017) or, like in large rivers in hydrological conditions similar to the Ogooué river, on specific patterns of mixing between tributaries (Moquet et al., 2016;Bouchez et al., 2017). Therefore, we took the estimated bias from the Congo River as a measure of uncertainty on the Ogooué and Mbei river solute fluxes associated with the sampling strategy, i.e., ±20% and ±10% for cation and TDS fluxes, respectively.

Uncertainty Calculation
The uncertainties associated to the calculated flux values (Tables 1, 4 and related figures) take into account the propagation of the uncertainty on the major element concentration measurements (5%), a sensitivity test performed on the Ca/Na sil and Mg/Na sil ratios used in the Equations 6 and 7, the uncertainty on the discharge estimate at each sampling location (RMSE = 161 mm yr −1 ), and the uncertainty associated to our relatively loose sampling time resolution (20%). Note that the propagation of the uncertainty of Ca/Na sil and Mg/Na sil produces asymmetric values and only affects F TDS sil , F TZ + sil, F CO 2 sil , and D chem sil . In the text (sections Results and Discussion) and in the Table 4, we thus report these fluxes values as "F +x −y " where F is the central estimate, and F +x and F −y the upper and lower bound of the range of estimates, respectively. According to this method, the relative uncertainties on F TDS, F TDS * , F TDS sil , F TZ + sil , F CO 2 sil , and D chem sil range between 11 and 84% ( Table 4).

RESULTS
In order to ease the presentation of results and the discussion thereof, we divided the samples into five groups which correspond to individual basins, as well as geomorphological and lithological units (Table 1; Figures 1, 2): The Mbei tributaries, the Northern Ogooué tributaries, the Plateaux Batéké Ogooué tributaries, the Southern Ogooué tributaries and the Ogooué River main channel. Note that the Ogooué River main channel exhibit intermediate values for all parameters (pH, conductivity, solutes concentration) indicating that its composition simply results from the mixing between the composition of the upstream tributaries inputs.

Physico-Chemical Parameters
The pH of the river water samples range between 4.63 and 7.62 and the conductivity range between 5 and 78 µS cm −1 ( Table 2). Intermediate values of pH (6.62-7.10) and conductivity (16-27 µS cm −1 ) were recorded in the Ogooué mainstream. The two lowest values of pH (<5.5) and conductivity (<10 µS cm −1 ) were recorded in two tributaries draining the Plateaux Batéké (Lékoni and Passa rivers). The highest values (pH > 7.2 and conductivity > 45 µS cm −1 ) were recorded in the Southern Ogooué tributaries. The other rivers exhibit intermediate values (6.30 < pH < 7.18; 11 < conductivity < 34 µS cm −1 ). Water temperature range between 21.6 and 28.1 • C. The lowest temperature values were recorded in the Mbei tributaries (21.6-22.9 • C) whereas the Ogooué Basin samples exhibit a narrow temperature range (23.9-28.1 • C). No direct relationship between temperature and elevation is observed. The instantaneous SPM (Suspended Particulate Matter) concentration range between 6 and 70 mg l −1 ( Table 2) and its distribution does not display any clear spatial distribution.

Major Elements and Dissolved Organic Carbon
The TDS concentration (and conductivity) are variable throughout the basin and range between 11 and 75 mg l −1 ( Table 2). The highest values (55-75 mg l −1 ) were recorded in the Southern basins, the lowest values (11-23 mg l −1 ) are observed in the rivers draining the sandstone region of the Plateaux Batéké, while the other groups (Northern basins, Mbei tributaries, and Ogooué main channel) exhibit intermediate values (21-39 mg l −1 ).
The total cationic charge (TZ + , in meq l −1 ) is generally dominated by Ca 2+ , Mg 2+ , and Na + in almost equivalent proportion, while the contribution of K + to TZ + is systematically lower (7 to 12% of TZ + ). The anionic charge (TZ − , in meq l −1 ) is generally dominated by HCO − 3 (> 80% of TZ − ). Concentrations of Ca 2+ , Mg 2+ , Na + , K + , HCO − 3 , and SiO 2 are significantly correlated to the conductivity and to TDS concentration and thus followed the same spatial distribution ( Table 2). The SiO 2 contribution to TDS concentration range from 86% (Plateaux Batéké) to 31% (Southern basins) and decrease in importance as the TDS concentration increases. Interestingly, dissolved Si concentration is correlated to HCO − 3 concentration (R = 0.91; n = 24; p < 0.01; Figure 4B), showing that the alkalinity and, therefore, CO 2 consumption associated with water-rock interactions were likely due to silicate weathering in the Ogooué and Mbei basins. The concentrations of Cl − and SO 2− 4 do not follow the same spatial distribution. We did not identify any parameter controlling the SO 2− 4 concentration distribution; however, SO 2− 4 range from 3.7 to 21 µmoles l −1 which is small by comparison with the global riverine discharge-weighted average (SO 2− 4 = 108 µmoles l −1 ; Burke et al., 2018). The Mbei tributaries, which are located relatively close to the coast (∼150 km), exhibited variable Cl − concentrations (9-28 µmole l −1 ). In the Ogooué Basin, the Cl − concentrations range between 3 and 25 µmole l −1 and depend on the basin distance from the sea, as observed in the Congo Basin (Négrel et al., 1993; Figure 4A). This observation confirms that Cl − concentration was mainly controlled by the atmospheric inputs ( Figure 4A). The highest NO − 3 concentrations (5.4-10.9 µmole l −1 ) were recorded in all the Mbei tributaries and in one river draining a Northern basin (Missanga at Ndjolé), all corresponding to basins draining the southern Cameroon plateau over catchment areas smaller than 500 km². The other rivers exhibited low NO − 3 concentration (<2 µmole l −1 ).
The normalized inorganic charge balance (NICB = TZ + − TZ − TZ + + TZ − ; note that charges borne by organic matter are not taken into account in this definition of NICB) is smaller than ±10% for most samples. For 5 samples the NICB was within 11% to 23%, which reflected an excess of cationic charge relatively to the anions ( Table 2). As suggested by the weak but statistically significant correlation between DOC (see below) and NICB (R= 0.63; N = 24; p < 0.01), as previously reported for Guyana rivers (Sondag et al., 2010), and by the fact that the NICB decreased to 0% when the concentration of "inorganic" solutes increases, the on-average positive NICB is most likely due to the presence of negatively charged dissolved organic matter. Assuming a negative charge of 6 ± 0.5 µeq mg −1 of DOC , the "corrected" NICB values of less than ±10% for 23 samples and less than ±15% for all samples confirm that the deficit of negative charge relatively to cationic charge is most likely due to the presence of negatively-charged dissolved organic matter. Dissolved organic carbon concentration ranged between 0.9 and 11 mg l −1 . The highest values were recorded in the Northern basins with values ranging from 4.3 to 11 mg l −1 , while the other samples exhibited values ranging between 0.9 and 2.7 mg l −1 .

Strontium Isotope Ratios
Dissolved 87 Sr/ 86 Sr ratios range between 0.7125 and 0.7226 over the studied basins ( Table 2). According to the 87 Sr/ 86 Sr vs. Ca/Na relationship ( Figure 4C) and the end members determined by Négrel et al. (1993) for crystalline rock types in the Congo Basin, these signatures are distributed between the silicate weathering, evaporite weathering, and rain end members. As no evaporite outcrops are known in the studied area (Thiéblemont et al., 2009), the 87 Sr/ 86 Sr ratios indicate that in the region major river solutes derive mainly from silicate weathering and atmospheric inputs. The 87 Sr/ 86 Sr ratios does not display any specific pattern in terms of spatial distribution.

Discrimination of Solute Sources
According to the results of the source discrimination method explained in section Discrimination of Solute Sources, atmospheric inputs to river chemistry in the Ogooué and Mbei basins are low and contribute to <15% of the TDS for most of the rivers (Figure 5). Interestingly, according to this method, the entirety of river SO 2− 4 derive from atmospheric inputs. Again, this is consistent with the absence of known evaporite outcrops in the region. The atmospheric contribution to the river budget of other dissolved species is generally lower than 40% with a decreasing impact in the order: Mg 2+ (40 ± 21%) > K + (29 ± 16%) > Ca 2+ (27 ± 12%) > Na + (19 ± 21%) > Sr 2+ (5 ± 3%). Given the low relative input of rain to the dissolved Sr budget, no attempt was made to correct dissolved 87 Sr/ 86 Sr ratios from the rain contribution.
The Mg/Na * , Ca/Na * , and HCO 3 /Na * molar ratios of the sampled waters are consistent with the silicate end member previously defined for the Congo Basin (Négrel et al., 1993) and used at global scale by Gaillardet et al. (1999b) (Figure 6). Together with the strong relationship between HCO − 3 and SiO 2 (Figure 4B) and the overall relatively high 87 Sr/ 86 Sr ratios (Figure 4C), this observation confirms that the major solutes in the Ogooué-Mbei basins are mainly sourced from silicate weathering. More quantitatively, source discrimination calculations show that silicate weathering largely dominate the TDS export (by up to 70%; Figure 5). Carbonate weathering significantly contribute to TDS only in one Mbei tributary and one Ogooué R. tributary (Ouagna at Wagny). Such low contribution of carbonate weathering to the solute load of the study area is consistent with the lithological map, which displays only sparse carbonate outcrops in the Ogooué and Mbei basins.

Trace Elements
The concentrations of trace elements in the Ogooué and Mbei basins (Table 3; Figure 7), are generally lower than the global average . The concentrations of trace alkali and alkaline earth metals (Rb, Cs, Sr, and Ba) in the Ogooué Southern basins are higher than in the Ogooué main channel, and are the lowest in the Plateaux Batéké region. The Northern Ogooué basins and the Mbei tributaries exhibit concentrations similar to those of the Ogooué mainstream. Conversely, transition metals and REE concentrations are generally higher in the Northern basins and in the Mbei tributaries than in Southern and Plateaux Batéké tributaries. For most trace elements, Plateaux Batéké tributaries exhibit the lowest concentrations (exceptions: Al and Cd).  Based on correlation analysis with other parameters (physicochemical parameters and concentration of DOC and of major elements) across the sample set, two groups of trace elements can be distinguished. First, elements such as Be, Al, Cr, Fe, Co, Y, most REEs, Zr, Th, and U correlate positively with DOC concentration (R > 0.5; p < 0.01; Supplementary Figure 4; see for example the Fe-DOC relationship in Figure 8). In particular, the higher DOC concentration measured in the Northern Ogooué basins corresponds to higher concentration for these elements. Second, other elements such as Li, V, Cu, As, Rb, Sr, Ba, and Eu correlate positively (R > 0.5; p < 0.01) with conductivity and the concentration of most major elements, and are therefore reflective of release by rock weathering. These elements thus exhibit high concentration in  Table 1. Global river averages are from Gaillardet et al. (2014) for trace elements and from Meybeck (2003) for Na + , K + , Mg 2+ , and Ca 2+ . Elements are ordered along the y-axis from the highest to lowest values recorded in the Ogooué main channel.
the Southern basins, low concentration in the Plateaux Batéké and intermediate concentration in the other basins. We note that B concentration is strongly correlated to Cl − concentration (which might point toward a dominantly atmospheric origin of B in the Ogooué-MBei rivers), and that other elements (Ti, Mn, Ni, Nb, Cd, Cs) do not show any significant correlation with the parameters cited above. Interestingly, the elements commonly considered as weakly soluble during weathering (e.g., Al, Fe, REEs+Y) or strongly insoluble (Th, Zr) are correlated to DOC concentration in the Ogooué and Mbei basins, where they exhibit low concentration both in comparison with global rivers  Figure 7) and with other African rivers like the Nyong (Viers et al., 2000), Niger (Picouet et al., 2002), and Congo (Gaillardet et al., 1995) rivers.

Silicate Weathering Fluxes and Associated CO 2 Consumption
Using estimates of Q spe values (specific discharge; see section Hydrological and Climate Data), the TDS and silicate weathering fluxes (F TDS and F TZ + sil ) of the Ogooué and Mbei rivers and their tributaries were calculated. As carbonate weathering is a small contributor to the river dissolved load over the studied area, no attempt was made to estimate the corresponding river dissolved fluxes. At the Lambaréné station, the sampling location closest to the Ogooué outlet, the Ogooué river export a F TDS of 4.4 ± 1.1 Mt yr −1 including 2.1 ± 0.5 Mt yr −1 of F SiO 2 and 0.5 +0.2 −0.2 Mt yr −1 of F TZ + sil . The corresponding F spe are 21 ± 5.3, 10 ± 2.5 and 2.5 +0.8 −0.9 t km −2 yr −1 , respectively, and the D chem sil was 5.72 +1.46 −1.61 m Ma −1 ( Table 4). Values of the F spe TDS, F spe TZ + sil , and D chem sil parameters are particularly variable throughout the basin (Figure 9A). The highest values (F spe TDS = 29 ± 16 to 86 ± 31 t km −2 yr −1 ; F spe TZ + sil = 3.9 +2.3 −2.4 to 12.1 +4.4 −5.3 t km −2 yr −1 ; D chem sil = 7.8 +4.4 −4.5 to 22.5 +8.1 −8.6 m Ma −1 ) are recorded in the Southern basins, the lowest ones are observed in the Plateaux Batéké tributaries (F spe TDS = 18 ± 7 to 19 ± 7 t FIGURE 9 | Spatial distribution of weathering fluxes in the Ogooué Basin: (A) Color map showing the spatial distribution of specific fluxes (F spe ) of TDS, TZ + sil , SiO 2 , and associated CO 2 sil consumption F spe distribution throughout the Ogooué Basin (the darkness of the color is proportional to these rates; the extreme values are reported in the caption). Note that at this scale the basins of the Mbei tributaries are too small to be visible. (B) Relative contribution of each region in term of drainage area, discharge and fluxes of TDS, TZ + sil , SiO 2 , and CO 2 sil consumption. The uncertainties are reported in Table 4.

Silicate Weathering Rates in the Ogooué and Mbei Basins, and Comparison With Regional and Global Rates
This study presents the first TDS flux estimates for the Ogooué River. With a TDS flux F TDS = 4.4 ± 1.1 10 6 t yr −1 , the Ogooué River contributes to around 4% of the TDS flux for around 7% of the discharge of Western Africa (according to Western Africa TDS flux and discharge estimates of Milliman and Farnsworth, 2011). The Ogooué River area-normalized TDS flux (F spe TDS = 21 ± 5 t km −2 yr −1 ) is lower than the global average (30-50 t km 2 yr 1 ; Meybeck, 1976Meybeck, , 2003Milliman and Farnsworth, 2011), partly due to the fact that at the global scale, a large fraction of the river TDS flux is provided by the weathering of carbonate and evaporite rocks, which are not present in the studied area. However, the area-normalized specific fluxes of silicate-derived cationic charge and CO 2 consumption for the whole Ogooué (F spe TZ + sil = 2.5 +0.8 −0.9 t km −2 yr −1 ; F spe CO 2 sil = 115 +40 −54 10 3 moles km −2 yr −1 ) are close to the global average (F spe TZ + sil = 2.3 t km −2 yr −1 and F spe CO 2 sil = 96-117 mol km −2 yr −1 ; Gaillardet et al., 1999b;Moon et al., 2014). These numbers challenges the paradigm of cratonic areas being rather inactive in terms of silicate weathering (e.g., Gaillardet et al., 1999b;West, 2012;Moon et al., 2014).
These flux values estimated at the Ogooué outlet result from the mixing of waters derived from regions characterized by a variety of heterogeneous weathering rates. Indeed, the river TDS, TZ + sil , SiO 2 , associated CO 2 consumption F spe and the chemical silicate denudation (D chem sil ) estimated throughout the Ogooué Basin are heterogenous. These estimates cover the whole range of values yet measured in cratonic environments under all latitudes (Figure 10). As throughout the Ogooué Basin the area-normalized fluxes of TDS, TZ + sil , SiO 2 (and associated CO 2 FIGURE 10 | Silicate weathering rate (F spe TZ + sil in t km −2 yr −1 ) vs. specific discharge (Q spe ) for rivers of the Ogooué and Mbei basins, and of other rivers draining cratons at the global scale. The HYBAM database provides major elements fluxes for the Upper Negro, Branco, Xingu, and Tapajos (Moquet et al., 2016), Oyapock andMaroni (Sondag et al., 2010) and Congo (Laraque et al., 2009) for the period 2003-2019. Values of F spe TZ + sil for the HYBAM stations here correspond to the cation fluxes corrected from the atmospheric inputs, according to the method used in the present study. The remaining cations rate were considered to be only derived from silicate weathering. Error bars reflect uncertainties, the calculation of which is detailed in section Uncertainty Calculation. consumption) as well as the D chem sil are strongly correlated to one another (R > 0.87; p < 0.01), here we discuss only controls on the values of F spe TZ + sil . At the global scale, cratonic environments exhibit a large variability in silicate weathering rates (Figure 10), ranging from 0.25 t km −2 yr −1 (Slave River, Canada; West et al., 2002) to 16 t km −2 yr −1 (Kaveri River, India; Pattanaik et al., 2013). In the Ogooué Basin the lowest values obtained for the Plateaux Batéké (F spe TZ + sil < 1.3 +0.7 −0.7 t km −2 yr −1 ) are amongst the lowest silicate weathering rates measured on Earth, lower than or similar to weathering rates recorded in the Siberian craton (Zakharova et al., 2005(Zakharova et al., , 2007Pokrovsky et al., 2015) and in the Canadian shields (Millot et al., 2002). The highest values obtained for the Southern Ogooué basins (F spe TZ + sil = 3.9 +2.3 −2.4 to 12.1 +4.4 −5.3 t km −2 yr −1 ) are on the order of those measured in the Kaveri Basin (India), a region particularly active in terms of silicate weathering (Gurumurthy et al., 2012;Pattanaik et al., 2013). Interestingly, by comparison with the silicate weathering rates measured in orogenic areas, these values are in the lower range of those recorded in the Andes (F spe TZ + sil = 9-104 t km −2 yr −1 ; Moquet et al., 2011Moquet et al., , 2018, in the New Zealand Alps (F spe TZ + sil = 2-187 t km −2 yr −1 ; Moore et al., 2013) but are commensurate to, or higher than those recorded in the Himalayas (mean F spe TZ + sil = 5.78 t km −2 yr −1 , West et al., 2002), in the European Alps (F spe TZ + sil = 0.03-11 t km −2 yr −1 , Donnini et al., 2016) Négrel et al., 1993;Gaillardet et al., 1995), the Nyong (F spe TZ + sil = 0.49 t km −2 yr −1 ; Viers et al., 2000;Regard et al., 2016) or the Niger (F spe TZ + sil = 0.58 to 2.2 t km −2 yr −1 ; Picouet et al., 2002) rivers.
Among cratonic environments, no direct relationship between specific discharge and silicate weathering intensity is observed, neither at the global scale nor at the scale of the Ogooué Basin (Figure 10). The present study shows that the silicate weathering flux calculated at the outlet of the Ogooué Basin does not reflect an intrinsic property of weathering in cratonic areas (i.e., low physical erosion rates associated to low weathering rates; e.g., West, 2012), but results from the mixing of solute fluxes derived from contrasted environments in terms of weathering. In addition, our observations highlight that cratonic areas can be particularly active in terms of weathering, as in the Southern Ogooué Basin, and that the absence of tectonic activity does not necessarily imply slow weathering.

Controlling Factors of Silicate Weathering in the Ogooué Basin and Implications for the Long-Term Carbon Cycle
At the global scale, silicate weathering is controlled by a range of variables like climate, lithology, geomorphology or the presence of organic matter (e.g., Goudie and Viles, 2012) which together influence the rate of exposure of mineral surfaces to water and air, the type of water-mineral interactions, and water flowpaths.

Climate
With all things considered equal, climate is a key driver in differential weathering reaction rates, which can alter reaction temperature (e.g., Oliva et al., 2003) and water availability and potential for weathering (White and Blum, 1995;Maher and Chamberlain, 2014). Due to its equatorial position, the Ogooué Basin is characterized by a relatively constant temperature throughout the year with an average of ∼24 • C and is relatively homogenous throughout the basin. Annual rainfall estimates in the basin exhibit a significant spatial variability (from 1,890 to 2,692 mm yr −1 ) with an estimated specific discharge ranging from 512 to 1,694 ( ± 161) mm yr −1 . Interestingly, amongst the Southern Ogooué tributaries the values of F spe TZ + sil calculated significantly (p < 0.01) increase with specific discharge suggesting that within this domain rainfall distribution partly controls the F spe TZ + sil variability ( Figure 11A). This is consistent with the observation made at global scale for numerous river basins (e.g., Gaillardet et al., 1999b;West, 2012) that weathering rates increase with precipitation, and thus specific discharge. In this case, weathering rates are controlled by the amount of water available for weathering reactions. Indeed, increased water flow through regolith and rock increases the surface of contact between minerals and water (White and Blum, 1995), and the export flux of solutes when reactions occur near thermodynamic equilibrium (Maher, 2010). However, for a given Q spe value, F spe TZ + sil is highly variable throughout the basin ( Figure 11A). Therefore, at the scale of the whole studied area, climatic parameters (temperature, rainfall or runoff) alone cannot explain the significant variability in silicate weathering rates observed throughout the basin.

Lithology
Lithology is another major factor controlling the Earth surface chemical denudation (e.g., White and Blum, 1995;Hartmann et al., 2014). In the Ogooué Basin, a strong lithological contrast exists between the Plateaux Batéké in the East and the rest of the basin. Plateaux Batéké are composed of Cenozoic sandstone formations, whose mineralogy is dominated by quartz, and which are very poor in soluble elements (e.g., Thiéblemont et al., 2009), whereas the remaining part of the basin is mainly made of plutonic and metamorphic rocks. As already observed in the Congo Basin (Négrel et al., 1993;Gaillardet et al., 1995), this contrast explains the very low major cation content measured in the Ogooué tributaries draining the Plateaux Batéké.
The remaining part of the basin drains the Congo craton mainly composed of Archean and Paleoproterozoic granitic and gneissic rocks. Non-calcareous sedimentary formations (Eburnéen foreland-Francevilien D Group composed of greenish pelite with sandstone and tuffaceous intercalations; Thiéblemont et al., 2009) constitute a more important component of the underlying rocks in the Southern basins than in the Northern basins (see Supplementary Table 1). The dissolved Mg/Na * vs. Ca/Na * and HCO 3 /Na * vs. Ca/Na * ratios of the studied basins are relatively homogenous (Figure 6) and lie within the field of the granitic weathering end member previously defined for the Congo Basin (Négrel et al., 1993). This observation indicates that the river the dissolved major element composition of the Ogooué Basin derives from a relatively homogenous silicate domain. Therefore, other drivers of chemical weathering need to be invoked to explain the latitudinal gradient in chemical weathering rates observed in the Ogooué Basin (Figure 9).

Dynamics of Organic Matter and Colloids
The Northern tributaries of the Ogooué Basin drain the southern part of the Cameroon plateau and exhibits characteristics similar to those of the neighboring Nyong Basin: tectonic quiescence, low slopes, and subsequent water stagnation in swamp systems (Viers et al., 1997;Oliva et al., 1999;Braun et al., 2012). In such environments, chemical weathering contributes to around 60-70% of total (physical plus chemical) denudation (Regard et al., 2016), rapid turnover of carbon on hillside soils promote DOC exportation (Nkoue Ndondo et al., 2020), and organicrich surface waters enhance the mobilization and transfer of insoluble elements such as Al and Fe-and even of highly-charged transition metals such as Th, Ti, and Zr (Braun et al., 2005. The Northern tributaries exhibit weathering rates similar to those of rivers draining the Cameroon plateau (Nsimi brook, Nyong river, and tributaries; Figure 10) and display similar relationships between DOC and concentrations of High Field Strength Elements (Al, Fe, Th, Zr, and REEs) as those observed in these rivers (Olivié-Lauquet et al., 1999;Viers et al., 2000;Braun et al., 2012; Figure 8). In such environments, these elements are mainly transported in the stream in the colloidal phase [which is typically composed of a mixture of organic matter and Al and Fe oxy(hydroxides)], operationally included in the dissolved phase in the present study, because of the filter porosity. Such behavior has been shown to be a primary mechanism for these elements to mobilize in many surface and groundwater environments (for example, Pokrovsky et al., 2006;Pourret et al., 2007;Trostle et al., 2016;Vázquez-Ortega et al., 2016). In the Ogooué and Mbei basins, this phenomenon is exemplified by the correlation between the concentrations of these trace elements and DOC concentration (Figure 8; Supplementary Figure 4). However, these correlations can reveal either that dissolved organo-metal complexes are forming that solubilize Al, Fe, Th, and Zr (e.g., Oliva et al., 1999;Tamrat et al., 2019) and/or that these systems promote the export of colloids composed of organic matter, nanoclays, and Al and Fe oxyhydroxides (e.g., Rose et al., 1998;Olivié-Lauquet et al., 1999;Allard et al., 2004;Guinoiseau et al., 2016). In any case, we emphasize that in the neighboring small river basin of Nsimi, major dissolved cations (Ca 2+ , Mg 2+ , K + , and Na + ) and SiO 2 are typically transported in the "truly" dissolved phase, i.e., with a size < 5 kD (Viers et al., 1997). The river fluxes of these elements in the Ogooué and Mbei basins thus reflect actually dissolved weathering products, and the intermediate weathering rates calculated in this region are at least partly controlled by the content in dissolved organic matter and associated inorganic colloidal fraction. The Mbei tributaries and the Missanga River (a small Northern Ogooué tributary) exhibit the highest NO − 3 concentration (from 5.4 to 11 µmoles l −1 ) of the dataset (Table 2; Figure 11B). Some of these basins also feature the highest mean slopes of the study area (from 3.0 to 7.9%) and all drain an area <500 km². The other Northern basins exhibit low (<1.2 µmoles l −1 ) NO − 3 concentrations and low slopes <1.3%. The NO − 3 concentration values recorded in the Mbei and Missanga rivers are within the range recorded in the hillslope piezometers of the experimental Nsimi watershed in the neighboring Nyong Basin (Braun et al., 2005). There, NO − 3 concentrations decrease from 0.75 µmol −1 on hillslopes to values <0.12 µmole in downslope swampy areas and in the downstream Nyong River (Viers et al., 1997;Boeglin et al., 2003;Braun et al., 2005). This trend was attributed to the fact that mineralization of organic-rich soil horizons on hillslope favors nitrification, while denitrification occurs in the open water, swampy area (Braun et al., 2005). We suggest that similar processes explain the high NO − 3 concentrations observed in some of the Northern Ogooué basins. Indeed in these small (<500 km²), overall steep basins, hillslope soil processes favoring nitrification would dominate the catchment-scale N cycle, and thus result in high river NO − concentration and export fluxes at the outlet. The Missanga stream flows through the suburbs of Ndjolé city, such that its higher NO − 3 concentration can also be attributed to human activities. These observations show that in relatively small basins solute loads do not reflect the same processes as those of larger basins. We thus focus our following discussion of the roles of relief and erosion (section Relief and Erosion) on catchments with drainage areas > 500 km 2 , and accordingly exclude from this discussion the data from the Mbei and Missanga rivers.
Altogether, these observations suggest that the presence or absence of flat, swampy areas exert a significant control on the export fluxes of trace elements and nitrate in the Ogooué Basin.

Relief and Erosion
The Southern tributaries of the Ogooué, which exhibit the highest weathering rates of the whole basin, contain DOC concentrations similar to those of rivers draining the Plateaux Batéké, and drain rocks similar to those of the Northern basins. The main characteristic of this region amongst the sub-basins of the Ogooué Basin is the relatively steeper slopes than the Northern sub-basins and the main Ogooué channel (mean slopes = 1.5-2% for the Southern sub-basins; Table 1). It is widely acknowledged that in tectonically active areas mountain uplift triggers mechanical erosion, which in turns favors the exposure of "fresh" mineral surfaces to water and reactive gases, providing a potential explanation for the high weathering rates recorded in orogenic regions (Larsen et al., 2014) such as the Himalayas (Raymo and Ruddiman, 1992;France-Lanord et al., 2003), the Andes (Stallard, 1985;Moquet et al., 2011Moquet et al., , 2014Moquet et al., , 2018Carretier et al., 2018), or the Australian Alps (Hagedorn and Cartwright, 2009). The Ogooué Basin drains a tectonically inactive area, such that mountain uplift cannot be invoked as a control on weathering rates via increased erosion. However, Guillocheau et al. (2015) point out that the Congo Basin borders, which correspond to the Southern and Eastern parts of the Ogooué Basin, have been subjected to active uplift for the past 45 Ma. The origin of this uplift, potentially related to mantle-induced dynamic uplift or lithospheric processes (Guillocheau et al., 2015(Guillocheau et al., , 2018, has in particular promoted higher relief in the Atlantic side of the Congo Basin for the past 11.5 Ma. Conversely, the Northern Ogooué tributaries and the Mbei Basin, located in the Southern part of the South Cameroon Plateau, do not undergo such uplift (Weber et al., 2016;Guillocheau et al., 2018). Although little is known about the effects of long-term mantle or lithosphere dynamics on erosion and weathering, and excluding the lithological effect explaining the low weathering rates of the Plateaux Batéké, the North-South gradient in specific weathering fluxes over the Ogooué Basin can be attributed to a southward increase in the intensity of mantle-induced dynamic uplift, or potential to lithospheric destabilization (Cottrell et al., 2004;Jaupart et al., 2007).
As for mountain uplift in tectonically active areas, the most likely process for a positive influence of mantle-or lithosphere-induced dynamic on weathering in tectonically quiescent areas is the enhancement of erosion rates. However, whereas in mountains a tight coupling between erosion and weathering is required to maintain a finite regolith thickness over significant timescales (>10 3 -10 4 yrs), in tropical cratonic environments mantle uplift would rather sustain regolith rejuvenation (i.e., exposure of deep regolith horizons hosting relatively unweathered primary minerals to water flowpaths in the critical zone) through long-term thinning of the lateritic cover. Regardless of the exact dynamics at play, in both cases erosion remains the main driver for "fresh" material exposure to the Earth surface. Estimates of erosion rates are not available yet for the Ogooué Basin to test this hypothesis, but the steeper slopes of the Southern basins lends support to the following scenario ( Figure 11C). In the Southern Ogooué basins, higher erosion rates driven by mantle-induced dynamics or by uplift due to lithospheric instability would enhance weathering rates through accrued erosion and thinning of lateritic soils, allowing for an increase in the exposure of "fresh" mineral surfaces to reactive fluids, while in the Northern basins, the absence of uplift inhibits soil erosion and leads to the formation and the preservation of deep lateritic soils. Such scenario would imply that weathering in the Ogooué Basin operates with the "supply-limited" regime ( Figure 11C), meaning that chemical weathering rates are limited by physical erosion rates (e.g., Riebe et al., 2017).
To our knowledge, this is the first time that an influence of mantle or lithosphere dynamic (also sometimes referred to as "epeirogenic" uplift by geologists) on weathering rates is proposed for cratonic areas. This hypothesis implies that erosion could affect weathering, the global long-term carbon cycle, and climate not only through the formation of the main orogenic belts in collisional contexts (Raymo and Ruddiman, 1992) but also through the slow, large-scale dynamics of the mantle and lithosphere in cratonic areas. Cratonic environments can therefore have a significant role on the continental weathering budgets, and need to be considered carefully when evaluating the global carbon cycle (Goddéris et al., 2008;Carretier et al., 2014). As basin slope alone is an imperfect index for physical erosion, this hypothesis needs to be tested through the comparison of chemical weathering rates with erosion rates. As no sediment gauging program exists in the studied area, alternative approaches such as cosmogenic nuclides (e.g., 10 Be, Regard et al., 2016) or/and steady state calculations based on the geochemical composition of river sediments (Gaillardet et al., 1999a;Picouet et al., 2002;Louvat et al., 2008) need to be performed to better characterize the interplay between weathering and erosion in cratonic areas.

CONCLUSION
Cratonic areas located in humid tropical regions have been reported to exhibit low chemical weathering rates due to a "shielding" effect of deep, mature regolith covers. Nevertheless, thanks to their wide aerial extent on intertropical surfaces, cratonic areas represent at the global scale a significant proportion of the dissolved matter delivery to the oceans. Despite this crucial significance, assessment of chemical weathering fluxes and rates in these environments is seldom conducted, and the variability in the intra-cratonic weathering rates-especially in terms of the diversity of geomorphological setting-has generally not been considered to date.
The present study allows us to quantify the chemical weathering budget of the intertropical cratonic basins of the Mbei and the Ogooué basins, Gabon, and to explore the main drivers of weathering in this context. The chemical composition (major and trace element concentration and 87 Sr/ 86 Sr ratios) of 24 river water samples taken in September 2017 was measured. Solute source discrimination shows that atmospheric inputs account for <15% of the TDS, while silicate weathering is the main TDS source. Interestingly, the hydrochemical composition of the Ogooué River at the outlet and the whole-basin silicate weathering rate are similar to those of the world average, and similar to (if not higher than) other basins of Western Africa. The significant weathering flux of the Ogooué results from contrasted weathering regimes across the basin, which cover the range of values yet recorded in cratonic basins at the global scale.
In the Ogooué and Mbei basins, three domains, submitted to similar climate, have been identified: (i) Low weathering rates were recorded in the Plateaux Batéké (Eastern Ogooué tributaries). These low weathering rates, corresponding to the lowest values recorded in cratonic environments at the global scale, are due to a lithological effect associated with the low abundance of chemically-mobile elements in the minerals (i.e., mostly quartz) which compose the underlying Cenozoic sandstones.
(ii) Intermediate weathering rates, comparable to the Ogooué mainstream values, were recorded in the Mbei tributaries and the Northern Ogooué sub-basins. The higher DOC concentrations recorded in the Northern Ogooué sub-basins are associated with elevated concentrations of classically insoluble elements (e.g., Fe, Al, Th, Zr, REEs) suggesting a prominent control of colloidal transport in these basins. In these rivers, our data also show that active nitrification occurs on hillslopes while denitrification is promoted in swamp areas in the valleys.
(iii) The highest weathering rates were recorded in the Southern Ogooué tributaries basins. The rim of the Congo Cuvette drained by these tributaries is submitted to active mantle uplift, which enhances soil erosion, and leads to the dismantlement of soils and to increased availability of fresh mineral surfaces to reactive fluids, thereby promoting higher weathering rates. In addition to the lithological effect observed in the Plateaux Batéké and to the role of organic matter in Northern basins, the novelty of the present study is to hypothesize that mantleinduced dynamic uplift or lithospheric destabilization in cratonic areas can significantly enhance chemical weathering rates by leading to soil erosion and bringing fresh rocks in contact with meteoric water. As a corollary, the cratonic zones may also host hot-spots of weathering. These Earth surface movements of long spatial (100-1,000 km) and temporal (tens of millions years) wavelengths should be considered, along with the establishment of major orogenic belts linked to plate collision, in models of the global long-term carbon cycle and climate. The drivers of weathering in shield environments thus need deeper investigation to estimate continental weathering budgets and to constrain global-scale, long-term biogeochemical cycles.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
J-JB, JG, and JB designed the project. J-SM, J-JB, SB, and AM measured in-situ data and collected the samples. J-SM and JB performed and coordinated the chemical analyses. SB contributed to the hydroclimate data calculation. SB and J-PB contributed to the hydrological data production. EM contributed to the GIS data compilation. J-SM and JB interpreted the results with the help of J-JB and JG. J-SM wrote the manuscript with the help of JB and J-JB. JG, SC, and VR contributed to the data interpretation and revised the manuscript. M-CP contributed to the study implementation. All authors contributed to the article and approved the submitted version.

FUNDING
This study was supported by the project RALTERAC EC2CO INSU, by the International Joint Laboratory DYCOFAC (Dynamics of the forested ecosystems of Central Africa in a context of global change) and by the Programme Emergences of the City of Paris Chemical weathering of sediments in large tropical floodplains (agreement205DDEEES165). Parts of this work were supported by IPGP multidisciplinary programme PARI and by Paris-IdF region SESAME Grant No. 12015903 and by a grant overseen by the French National Research Agency (ANR) as part of the Investments d'Avenir Programme LabEx VOLTAIRE, 10-LABX-0100.