Our recent steady-state mass-balance modeling suggests that most global carbonic-acid weathering of silicate rocks occurs in the vadose zone of aquifer systems not on the surface by atmospheric CO2. That is, the weathering solute flux is nearly equal to the total global continental riverine carbon flux, signifying little atmospheric weathering by carbonic acid. This finding challenges previous carbon models that utilize silicate weathering as a control of atmospheric CO2 levels. A robust analysis utilizing global estimates of groundwater carbon concentration generated by a geospatial machine learning algorithm was coupled with recharge flux in a geographic information system environment to yield a total global groundwater carbon flux of between 0.87 and 0.96 Pg C/year to the surface environment. On discharging to the surface, the carbon is speciated between 0.01 and 0.11 Pg C/year as CaCO3; 0.35 and 0.38 Pg C/year as CO2 gas; and 0.49 and 0.51 Pg C/year as dissolved HCO3−. This total weathering carbon flux was calculated for direct ocean discharge (0.030 Pg C/year); endorheic basins (0.046 Pg C/year); cold-wet exorheic basins (0.058 Pg C/year); warm-dry exorheic basins (0.072 Pg C/year); cold-dry exorheic basins (0.115 Pg C/year); and warm-wet exorheic basins (0.448 Pg C/year).
