Thank you very much for your reply. If I understand correctly, GWB should have a known bulk HCO3- in the fluid or fugacity to calculate the corresponding fugacity or dissolved carbonates, respectively. The software cannot calculate carbon distribution between gas and liquid phases based on a total carbon equal to a sum of CO2 in the gas phase and aqueous carbonate species, right?