Abstract. The incomplete spatial coverage of CO2 partial pressure (pCO2) measurements across estuary types represents a significant knowledge gap in current regional- and global-scale estimates of estuarine CO2 emissions. Given the limited research on CO2 dynamics in large estuaries and bay systems, as well as the sources of error in the calculation of pCO2 (carbonic acid dissociation constants, organic alkalinity), estimates of air–sea CO2 fluxes in estuaries are subject to large uncertainties. The Estuary and Gulf of St. Lawrence (EGSL) at the lower limit of the subarctic region in eastern Canada is the world's largest estuarine system, and is characterized by an exceptional richness in environmental diversity. It is among the world's most intensively studied estuaries, yet there are no published data on its surface-water pCO2 distribution. To fill this data gap, a comprehensive dataset was compiled from direct and indirect measurements of carbonate system parameters in the surface waters of the EGSL during the spring or summer of 2003–2016. The calculated surface-water pCO2 ranged from 435 to 765 µatm in the shallow partially mixed upper estuary, 139–578 µatm in the deep stratified lower estuary, and 207–478 µatm along the Laurentian Channel in the Gulf of St. Lawrence. Overall, at the time of sampling, the St. Lawrence Estuary served as a very weak source of CO2 to the atmosphere, with an area-averaged CO2 degassing flux of 0.98 to 2.02 mmol C m−2 d−1 (0.36 to 0.74 mol C m−2 yr−1). A preliminary analysis revealed that respiration (upper estuary), photosynthesis (lower estuary), and temperature (Gulf of St. Lawrence) controlled the spatial variability in surface-water pCO2. Whereas we used the dissociation constants of Cai and Wang (1998) to calculate estuarine pCO2, formulations recommended for best practices in open ocean environments may underestimate pCO2 at low salinities, while those of Millero (2010) may result in overestimates.