Abstract. The incomplete spatial coverage of CO 2 partial pressure (pCO 2 ) measurements across estuary types represents a significant knowledge gap in current regional-and global-scale estimates of estuarine CO 2 emissions. Given the limited research on CO 2 dynamics in large estuaries and bay systems, as well as the sources of error in the calculation of pCO 2 (carbonic acid dissociation constants, organic alkalinity), estimates of air-sea CO 2 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 pCO 2 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 pCO 2 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 CO 2 to the atmosphere, with an area-averaged CO 2 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 pCO 2 . Whereas we used the dissociation constants of Cai and Wang (1998)