Seasonal air-sea CO2 flux dynamics in Colombia’s Gorgona Marine Area during La Niña 2021–2022

Air–sea CO₂ fluxes in tropical coastal zones are strongly influenced by ENSO variability, but in situ measurements in the Eastern Tropical Pacific remain scarce. We assessed seasonal CO₂ dynamics around Gorgona Island (Panama Bight, Colombian Pacific) under La Niña 2021–2022. From November 2021 to July 2022, we conducted monthly sampling at seven stations spanning the Guapi River plume to the open ocean, measuring physical (SST, SSS, thermocline depth), chemical (TA, DIC, pH, carbonate system parameters), and biological (chlorophyll-a) variables, and estimating net CO₂ fluxes (FCO₂) with the Liss and Merlivat (1986) parameterization and atmospheric CO₂ from NOAA. La Niña featured a cool-water anomaly (−0.78 °C), enhanced precipitation (+59%) and river discharge (+44%) relative to multi-year means. The nine-month mean CO₂ flux was near neutral (−0.01049 ± 0.00014 mol C m⁻²) but strongly seasonal: six post-upwelling months showed slightly positive fluxes (0.00929 ± 0.000147 mol C m⁻²) associated with high precipitation (746.4 ± 214.7 mm), warmer SST (27.5 ± 0.4 °C), elevated pCO₂w (567 ± 97.5 µatm) and lower pH (7.869 ± 0.040), whereas three upwelling months showed slightly negative fluxes (−0.00119 ± 0.00010 mol C m⁻²) with reduced precipitation (165.8 ± 82.4 mm), cooler SST (26.5 ± 0.2 °C), lower pCO₂w (461 ± 92.8 µatm) and higher pH (7.968 ± 0.048). La Niña amplified pCO₂w variability (316–839 µatm) via vertical Ekman pumping, horizontal transport (Zonal Ekman Transport, tides), and freshwater inputs, while a persistent thermocline (10–40.1 m) restricted deep CO₂-rich waters from reaching the surface. Biological uptake further modulated outgassing, as evidenced by chlorophyll-a and ΔDIC dynamics. Overall, CO₂ fluxes were relatively low compared with other tropical estuarine and oceanic sources. These results underscore the need for sustained in situ observations in estuarine–ocean systems to refine predictive models of CO₂ fluxes under ENSO conditions.