Carbonate system in the Cabo Frio upwelling

Oct 11, 2023

Scientific Reports volume 13, Article number: 5292 (2023) Cite this article

301 Accesses

6 Altmetric

Metrics details

The quantitative assessment of the carbonate system represents one of the biggest challenges toward the "Sustainable Development Goals" defined by the United Nations in 2015. In this sense, the present study investigated the Spatio-temporal dynamics of the carbonate system and the effects of the El Niño and La Niña phenomena over the Cabo Frio upwelling area. The physical characterization of the site was carried out through data on wind speed and sea surface temperature. Water samples were also collected during the oceanographic cruise onboard the Diadorim R/V (Research Vessel). From these samples, the parameters of absolute and practical salinity, density, pH, total alkalinity, carbonate, calcite, aragonite, bicarbonate dissolved inorganic carbon, carbon dioxide, partial pressure of carbon, calcium, and total boron were obtained. The highest average concentration of bicarbonate in S1 (2018 µmol/kg) seems to contribute to the dissolved inorganic carbon values (2203 µmol/kg). The values of calcite saturation state, aragonite saturation state, and carbonate were higher on the surface of each station (calcite saturation state = 4.80–5.48; aragonite saturation state = 3.10–3.63, and carbonate = 189–216 µmol/kg). The mean values of pH were similar in the day/night samples (7.96/7.97). The whole carbonate system was calculated through thermodynamic modeling with the Marine Chemical Analysis (AQM) program loaded with the results of the following parameters: temperature, salinity, total alkalinity, and pH parameters. This manuscript presents original data on the carbonate system and the "acidification" process influenced by the Cabo Frio upwelling, which directly depends on the El Niño and La Niña phenomena oscillations in the sea surface temperature.

Carbon dioxide (CO2) sources, transport mechanisms, and transformations are essential in oceanography field studies1,2. Inorganic CO2 can exhibit significant spatial and temporal variability within the same water mass since the ocean content is dependent on processes such as the atmospheric exchange through sea surface and organic matter degradation (both autochthonous and allochthonous derived)3.

The reduction in seawater pH caused by the increase in CO2 in this compartment can lead to a decrease in marine carbonate (Reaction 1), a process also known as ocean acidification (OA)4. The coastal ocean waters are naturally subjected to daily, seasonal, and even annual pH variations amplified by OA5. Seawater's pH oscillation affects the carbonate system speciation by reducing the amounts of \({\text{CO}}_{3}^{2-}\) while increasing CO2 and \({\text{HCO}}_{3}^{-}\) contents, interfering in the natural processes of photosynthesis and calcification of marine organisms, thus creating negative ecological, social, and economic impacts6.

The reduction in the amount of \({\text{CO}}_{3}^{2-}\) available in the water will reduce the ocean's capacity to remove the CO2 released into the atmosphere by human activities. The absorption of H+ and CO2 by \({\text{CO}}_{3}^{2-}\) reduces the capacity of shallow waters to retain CO2. Several authors have linked the calcium saturation state (Ω) with reducing the calcification capacity of organisms linked to a drop in \({\text{CO}}_{3}^{2-}\) availability5,7,8. As \({\text{CO}}_{3}^{2-}\) concentration decreases in the seawater (Reaction 1), there is a reduction in the carbonate saturation state (Ω) (Eq. 1). The Ω has been implied in the reduction of calcification of marine organisms. The sea animals that present carbonate structures such as skeletons, shells, and spines are the most affected by OA9,10.

The speciation and quantification of the carbonate system are considered challenges for the Sustainable Development Goals (SDG) defined in 2015 by the United Nations11 for the next nine years. One of these challenges is establishing an analytical protocol and implementing a monitoring program for OA. Data obtained in studies of OA (pH, TA, [\({\text{HCO}}_{3}^{-}\)], [\({\text{CO}}_{3}^{2-}\)], [CO2]aq, ρCO2, Ωcalc, Ωarag) are also essential to validate regional and global models of CO2 flow between ocean and atmosphere interface. OA studies have presented little progress in this direction1,12,13,14,15. The main difficulties related to the implementation of an OA Monitoring Program in coastal and oceanic waters are (1) the absence of a database of the carbonate system, (2) the lack of a unified protocol for the determination of pH and total alkalinity (TA), (3) the non-disclosure of the polynomial precision of the carbonate system data, and (4) lack of an integrated and open access data repository on CO2 fluxes in coastal and ocean waters16.

An OA monitoring network requires keeping records of the principal chemical parameters (i.e., pH and total alkalinity), which allows, for example, the determination of the aragonite saturation state (ΩArag) and a complete description of the carbonate system14. The suitable parameters for this purpose are defined through the balance of the reactions (2–6) that occur when CO2 dissolves in seawater2,17 as follows:

Climatic phenomena such as the El Niño Southern Oscillation can affect the intensity of upwelling by increasing the SST18. This is conditioned to the power of events that affect the tropical and subtropical cyclone and anticyclone systems, changing the intensity of the winds nearby the upwelling region and leading to an increase in the SST, decreasing the strength of the upwelling. The influence of these phenomena in the Cabo Frio upwelling region has not yet been fully understood, nor have their implications for the carbonate saturation state. This experiment is part of an extensive study on the feasibility of implementing a unified protocol in a monitoring program for the acidification of coastal and offshore areas of the Brazilian ocean, such as (1) upwelling sites, (2) coral reefs, (3) oceanic islands, and (4) oil platforms. The present investigation originally elucidates the space–time dynamics of the carbonate system in the Cabo Frio resurgence and assesses La Niña and El Niño in the resurgence phenomenon.

Located on the coast of Arraial do Cabo (Rio de Janeiro, Brazil, Fig. 1), the study site occurred between latitudes 22° 58′–23° 06′ South and longitudes 42° 10′–42° West. The Brazilian coast presents a NE-SW orientation in this region, causing the NE wind to blow parallel to the coast and the shallow coastal waters to move towards the open sea due to Ekman's transport19. Cabo Frio coastal upwelling phenomenon brings the deepest, coldest, and richest nutrient waters20.

Study area. Sampling stations are represented by dots and respective numbers (S1, S2, S3, and S4). This map was generated with the ArcMap software v. 10.8.2 (https://www.esri.com/en-us/arcgis/about-arcgis/overview), by Senez-Mello, TM.

Campos Basin is directly under the influence of the South Atlantic Subtropical Anticyclone (SASA)21, the main feature of the atmospheric circulation over the South Atlantic Ocean that affects the Brazilian weather and climate22. This feature is responsible for thermodynamic stability and low-intensity northeast winds that predominate in the southeast region of Brazil21. Teleconnection patterns (wave propagation in the upper atmosphere) can change the intensity and location of the South Atlantic subtropical high at any time of year, and changes in the positioning of the SASA generates a significant change in the wind pattern23.

Sun et al.24 analyzed data from 1979 to 2015 to compare the austral summer mean location of the SASA with the Southern Annular Mode (SAM)25 and extended multivariate El Niño–Southern Oscillation (ENSO) index (MEI)26. The authors found that when the SASA shifts poleward, the SAM was in a positive phase of La Niña, and when it shifts equatorward, the SAM was in a negative phase during El Niño. The wind speeds observed across Brazil can also describe this shift in the SASA pattern.

Three water masses occur at Arraial do Cabo from the shore to 500 m deep. The Tropical Water (TW) mass is found on the outer side of the continental shelf in the first 200 m and is characterized by a temperature range between 27.37 and 28.26 °C and a salinity range between 36.44 and 37.55 psu1,27,28. The South Atlantic Central Water (SACW) mass, ranging from 142 to 567 m in depth, is characterized by temperatures between 13.33 and 15.59 °C and salinity between 35.41 and 35.78 psu1,20,27,29,30,31. The SACW is a water mass rich in nutrients that reach the surface through the coastal Cabo Frio upwelling phenomenon20,32. The Coastal Water (CW) mass is found in the inner part of the continental shelf33, resulting from the mixture of continental waters, TW, and SACW. The CW is characterized by temperatures above 23 °C and salinity below 34 psu34. The primary water flux in the region is the Brazil Current System (BC). This current acts until approximately 500 m deep, carrying the TW and the SACW towards the south31.

The climate in the study area is classified as hot-semiarid35, characterized by intense regional evaporation and reduced rainfall if compared to adjacent areas, yielding an arid climate36.

Wind velocity and surface sea temperature (SST) data from the Cabo Frio (RJ) region, organized in a temporal series, were used to characterize the marine upwelling phenomenon. The daily averages of the wind velocity represented the NE (0°–45° direction) and SW (180°–270° direction) quadrants from September 2006 to December 2016. The data series of SST daily averages referred to the period between 1994 and 2016 (National Institute of Meteorology—www.inmet.gov.br). The data series of SST weekly averages referred to the period between 1994 and 2016 (data from the Admiral Paulo Moreira Marine Research Institute—IEAPM).

The daily average wind velocity from the NE and SW quadrants was tabulated with the respective SST (Sea Surface Temperature) averages from 2006 to 2016. In cases where the data of average wind velocity did not have corresponding SST data in the temporal series, the SST was estimated using the Kriging interpolation method. The estimate was obtained from the interpolation of the SST data in a regular grid, where: (X) corresponded to the time in days, (Y) to the time in years, and (Z) to the daily and annual SST37. A univariate statistic of residual values of SST resulting from the interpolation model was used to evaluate the standard error of the data, and to indicate the function that best fits the input data. The Q–Q (quantile–quantile) relation was used to analyze the wind velocity, NE, and SW, compared to its respective SST38. The data were compared to those of the El Niño and La Niña periods checking the influence of these events on wind velocity and SST39. Trend curves were also obtained relating the average wind velocity with the SST. Pearson's R2 coefficient was used to assess data correlation.

The Diadorim R/V from the Brazilian Navy collected the water samples on January 20, 2016, using a Northeast/Southwest oriented transect. The first station (S1) was located between Cabo Frio Island and the continent, and the furthermost station (S4) was located 14 km from the coast, at the 100 m isobath (Table 1, Fig. 1).

The sampling campaign involved a two-scale analysis. The spatial scale was held perpendicular to the coastline (S1–S4), Fig. 1. The ship was anchored in station 12 for 12 h for the temporal scale. Both campaigns were performed on the same day. In the spatial approach, water was sampled from the surface (~ 3 m) between 6:30 and 10:00 h (UTC) using a pump adapted to a hose without forming bubbles. The other water samples below the surface were collected through a Niskin bottle of 10 L for the middle (half of the total depth) and bottom (~ 5 m above the seafloor).

In the temporal approach, in station 2, the surface water samples were collected hourly (12 samples total), between 13:00 to 01:00 h (UTC), while middle and bottom water samples were collected in alternate hours: 13, 15, 17, 21, 23, and 1 h.

Water measurements, such as temperature, depth, and dissolved oxygen, were performed with a CTD vessel (Midas Valeport).

The modeling and the procedural calibrations described in this manuscript were performed with the Marine Chemical Analysis (AQM) program1,2,16,17. The AQM is a package of thermodynamic equations, executed via MS Excel, which can predict the complex composition of the marine carbonate system. This package is based on measurements that can be relatively inexpensively (pH, temperature, and alkalinity), reducing the overall costs of ocean acidification monitoring programs. The AQM program is available upon request to the corresponding author's email.

Non-parametric Kruskal Wallis test was chosen for comparisons between groups. All statistical tests were performed using the Statistica 7.0 software (TIBCO) with a significance level set at p < 0.05.

The analytical method was based on international procedures for studies involving the chemistry of inorganic carbon dioxide in marine waters15,40,41 with the necessary adaptations.

For the determination of TA, water samples were collected and filtered in a Nalgene filtration system through GF/F filters before being transferred to BOD type flasks (300 mL, Kimble brand) and immediately analyzed40.

The potentiometric determination was conducted with duplicate samples in an open thermostated glass cell, where 3 mL (to obtain v1) and 10 mL (to obtain v2) of HCl 0.1 M were added to each 100 mL sample42. The method consists in determining the slope of the line by obtaining two points for the function of Gran (F): F (1) defined by v1 and F (2) defined by v242. A Thermo Scientific Orion Star potentiometer coupled to the Orion glass reference electrode cell, model 8102BNUWP was used for potentiometric determinations. The pH electrode was calibrated daily with "Tris" buffer (0.04 m) for sample readings (maximum 12 samples per day). Due to the reduced number of samples per day, the short period of the oceanographic cruise, and the constant working conditions (electricity source, solutions, and equipment), we chose to verify the electrode performance at the beginning and the end of the oceanographic cruise. The electrode's percent efficiency ranged between 99.49 and 99.54% concerning the theoretical Nernst value (59 mV). More details are available in hydrogen potential (pH).

The analytical precision and accuracy were calculated from five replicates of the reference material (Dickson–CRM, for oceanic CO2 measurements, batch 104)43, which obtained a 95% recovery rate from the expected value (Table 2). The calculated TA was obtained by the AQM program through the equation: TA (µmol/kg) = 660 + 47.6S, defined by Hunter44 for waters of the Atlantic and Pacific oceans by the GEOSECS Program. The normalized total alkalinity (NTA) was obtained by the AQM program using the equation: NTA (µmol/kg) = TA (µmol/kg) × 35/Salinity (g/kg), where 35 was assumed to be the representative salinity of the water masses.

The total pH of the water samples collected during the cruise was determined in the "wet laboratory" as follows: pHT (= − log([H+] + [(HSO4−]/co), where co is the thermodynamic concentration (1 mol/kg-soln).

The internal solution of the combined pH electrode was filled up with 0.7 m NaCl to reduce the potential liquid junction. The electrode's electromotive force (emf) was related to the molar concentration of the proton [H+], as shown in Eq. (2).

where: E° is the standard electrode potential, determined by titrating a 0.7 m NaCl solution with 0.179 M HCl45. The pHT (total scale) values were measured immediately after each collection at a constant temperature of 25 ℃ in a thermostatic cell connected to a microprocessed thermostatic bath with external circulation (Qimis) to avoid temperature bias46. The determinations were made by the Thermo Scientific Orion Star potentiometer coupled to the Orion glass reference electrode with a 0.7 m NaCl outer chamber filling solution, model 8102BNUWP. The analytical slope for the electrode was within ± 0.13 mV (theoretical Nernst value at 25 °C). The electrode was calibrated with a "Tris" buffer (0.04 m) prepared in the laboratory47, where pH values were assigned by spectrophotometry (m-cresol method)16,40,48. The "Tris" buffer allows the accuracy of 0.001 pH units units47,49. Subsequently, using the AQM program, the pH results were corrected for the temperature recorded at the sampling moment (pHt = pH25 + A + Bt + Ct2)50.

The determination of Ca and TB was conducted using a MIP OES (microwave-induced plasma optical emission spectrometer, 4200 MP-AES, Agilent brand). The external analytical curves were made with monoelementary standards (1000 mg/L, VHG®) with concentrations ranging from 0.1 to 10 mg/L in an ultrapure water matrix. A matrix influence test was conducted in which it was found that both the boron and calcium signals did not show any significant difference between the ultrapure water matrices and the 500 mg/L NaCl solution. The calculated calcium and boric acid were also obtained using the equations: [Ca2+]T = 2.938 × 10−4 × S15, and [B]T = 0.000416 × (S/35)51. The analyzed versus calculated values of Ca and TB in water samples showed a relative error (RE%) from 0.2 to 8% (Table 3).

All parameters from the inorganic CO2 system (CO2, \({\text{CO}}_{3}^{2-}\), \({\text{HCO}}_{3}^{-}\), DIC, ρCO2, ΩCalc, and ΩArag) were calculated using the carbonate system dissociation constant K52 defined as follows:

\({lnk}_{B}^{*}\)53

\({lnk}_{Si}^{*}\)54

\({lnk}_{1}^{*}\)(H3PO4)55

\({lnk}_{2}^{*}\) (\({H}_{2}{PO}_{4}^{-}\))55

\({lnk}_{3}^{*}\) (\({HPO}_{4}^{2-}\))55

\({lnk}_{2}^{*}\) (\({CO}_{3}^{2-}\))56

The aqueous concentrations (CO2(aq)) and the partial pressure (ρCO2) were calculated from the variables of temperature, salinity, pH, and TA and by using the thermodynamic and stoichiometric constant K (\({pk}_{1}^{o}\), \({pk}_{2}^{o}\), \({pk}_{1}^{*}\), and \({pk}_{2}^{*}\))57,58. The AQM was also used in this phase, aiding the calculations.

The CO2 flux equation between oceans and the atmosphere is defined between aqueous CO2 and saturated CO2, defined as follows (Eq. 3):

Aqueous CO2 and saturated CO2 are components that characterize the "balance" of the flux equation (Eq. 3) determined in this study using the AQM program. Wind velocity has a significant effect on the gas transfer equation. The relationship between gas exchange and wind velocity can have non-linear effects on calculating the gas transfer velocities for particular wind velocity measurements, depending on the sampling design and the wind velocity1. The most critical parameter in the gases transfer velocity equation (Eq. 4) in terms of the function with wind velocity is based on the gaseous exchange coefficient (kT in cm/h):

where: u is the wind velocity module at 10 m from the surface in m/s, Sc is the Schmidt number of CO2 in seawater59,60, and 660 is the Sc value in seawater at 20 °C. The Schmidt number is defined as the water kinematic viscosity divided by a gas diffusion coefficient defined as follows (Eq. 5) by a polynomial:

The kT calculation (Eq. 4) considers that the wind velocity (u) has a fundamental quadratic dependence on the CO2 flux calculation. Typically, the wind velocity sampling for the kT calculation takes the climatological average into account at 10 m above the water surface61. However, such a procedure can include an error related to the average wind incompatibility and the CO2 sampling in situ. In addition, reading the wind speed on the sea surface (10 m) can produce an inaccurate result due to the orographic average of the wind induced by the vegetation. Thus, the closer the wind measurement is performed to the water surface, the more significant the effect of wind on Sc will be2.

The physical–chemical data obtained at the four stations are represented in Figs. 2, 3, and 4. The lowest temperatures were recorded at the bottom of each station, varying from 13.9 °C (S2) to 14.7 °C (S1), showing the presence of SACW (colder and less saline waters). In station 1, the water emerged, showing lower temperatures and salinities, like the bottom samples from the other stations. On the surface of station 1, the temperature was 18.2 °C, and the salinity was 35.59 psu, where the temperature and salinity seem to characterize SACW2. The SACW was observed in S2, S3, and S4 at intermediate depths, where the temperature ranged from 15.4 to 16.0 °C and salinity between 35.55 and 35.74 psu. The presence of SACW was observed at the bottom of all stations. On the surface of S3–S4, water temperatures ranged from 22.2 to 25.7 °C, and salinities were greater than 36 psu. When increasing the distance from the coast (from S2–5 km), the temperature and salinity increased on the surface, characterizing the TW, see Figs. 2 and 3. Salinity did not present significant variations between stations and depths (Fig. 3). The highest values were recorded at each station surface, varying from 35.59 (S1) to 36.37 psu (S2), and the lowest values were recorded at the bottom, varying from 35.39 (S2) and 35.41 psu (S4).

Temperature of the longitudinal samples from the surface, middle, and bottom of the water column.

Salinity of the longitudinal samples from the surface, middle, and bottom of the water column.

Dissolved oxygen in the longitudinal samples from the surface, middle, and bottom of the water column.

All stations showed the highest DO concentrations (193–216 µmol/L) in superficial samples and the lowest concentrations (180–196 µmol/L) in bottom-most samples, Fig. 4. The highest DO contents were obtained in S2, S3, and S4, between 10 and 20 m depth, where respiration may predominate over productivity, making this water mass richer in organic and inorganic carbon and nutrients20. Thus, when SACW comes to the surface, it brings water rich in CO2 and nutrients and O2 depleted. In shallow water, the increase in oxygen solubility is favored by photosynthesis and gas exchange with the atmosphere (Fig. 4).

It was observed through a daily data series (from September 2006 to December 2016) of the study area that the NE and SW wind speeds increased from the year 2010. Between 2010 and 2015, this increase was simultaneous with the periods in which SST showed a downward trend favoring higher NE and SW wind speeds and a greater intensity of Cabo Frio's resurgence from 2010 to 2015. Between 2012 and 2013, the El Niño and La Niña events were observed, which favored the most considerable variations in SST given by the significant deviations from the mean, verified by the lowest SST concerning the other periods of the time series (Figs. 5 and 6).

Wind velocity at quadrants SW and NE (blue line) with the respective variation of sea surface temperature (SST) at the Cabo Frio/RJ upwelling area (circle with a vertical bar in red), showing the average temperature and the standard deviation between 09/2006 and 12/2016, according to mild (Weak), and moderate to strong (Mod. to Str.) El Niño and La Niña events.

Temperature Grid of the Sea Surface Temperature (SST) at the Cabo Frio upwelling obtained by Kriging the SST data (1995–2016). Blue spots (cold): longer-term and intensity of the upwelling phenomenon. Red blots (hot): longer-term and intensity of temperature anomalies. Ellipse: SST in the period correspondent to the sampling campaign. Hachured rectangle: Str—increased wind magnitude period (from 2010 to 2016); Weak—absence of the events of a mild El Niño (from 2012 to 2013); and Mod. To Str—moderate to strong winds (from 2006 to 2010).

In contrast, before 2010, moderate to strong El Niño and La Niña events affected the Cabo Frio resurgence phenomenon by reducing its intensity, increasing SST, and decreasing wind speed (Figs. 5 and 6). Elias62 had already observed a similar pattern in the Cabo Frio resurgence, but on a larger time scale as the El Niño periods, considering stronger winds responsible for changes in the resurgence process pattern, generating reflexes in the resurgence process local. The author also found that the Cabo Frio resurgence phenomenon is related to the frequency of the El Niño and La Niña events, causing changes in SST and the wind dynamics throughout the seasons. Thus, affecting the nutrient transport and dispersion from the sea bottom to the surface in the upwelling vicinity.

During the water sample collection period in January 2016, according to Figs. 5 and 6, sea surface temperatures were rising (on average from 20 to 22 °C) due to the strong El Niño that occurred in 2015. Also, a decrease in the average wind speed from 10 to 8 m/s was observed (Fig. 5). This decrease in wind speed corroborated with the increase in temperature of the water samples collected in the oceanographic campaign of January 2016 (22.7ºC average). This average temperature is even lower when compared to the average found between the years 2006 and 2010 (Figs. 5 and 6).

When the El Niño and La Niña events are scarce or absent, it is observed that the SST decreases and the wind speeds up, increasing the intensity of the resurgence of Cabo Frio (Figs. 5 and 6). It was also observed that the duration of the more moderate events of El Niño and La Niña decreased the intensity of the resurgence, causing an increase in SST and a decrease in the speed of NE and SW winds (Figs. 5 and 6). This relationship between SST and wind speeds can be confirmed by comparing wind speed quantiles with SST in a Q–Q plot dispersion graph. This statistical test showed that the variables (wind speed x SST) have similar behavior in both wind directions (R2NE = 0.99 and R2SW = 0.94). In addition, the linear regression showed that the lower the SST, the greater the wind speed (Fig. 7).

Quantile–Quantile plot (QQ-plot) displaying the SST and Wind Velocity data distribution from 2006 to 2016 in the Cabo Frio upwelling area. Blue dots: SW wind. Green dots: NE wind. Red lines: best-fit linear equation. Greenline: NE wind intensity trend curve; blue line: SW wind intensity trend curve, both fitted to 95% prediction bands.

The SST data and El Niño and La Niña presented in Fig. 6, which pointed to an influence on the wind pattern in the Cabo Frio upwelling area (Fig. 5), may be related to the movement of lunar nutation. The most evident effect is the water imprisonment in a nodal cycle or Saros in each hemisphere due to the alteration of the gravity center between the Moon and the Earth63.

This 18.6-year cycle coincides from 1993 to 2011, with a minimum lunar orbit declination (− 18.35°) in 1998 and the maximum (− 28.65°) in 2006. The minimum lunar declination occurred again in the following Saros cycle in 2016 (U.S. Naval Ephemeris Observatory). This observation coincided with the trapping of the warmest waters in the Southern Hemisphere, which may influence the decrease in the frequency of the Cabo Frio resurgence phenomenon and increase the effects of El Niño and La Niña on the resurgence by decreasing the winds. From 2010 onwards, the upwelling frequency became more intense, with a more significant decrease in SST causing the highest observed wind speeds due to the diminishing effects of El Niño and La Niña on the upwelling. Exceptionally, concerning the other years of the time series, from 2010 onwards, there was a more significant variation in SST (on average of − 3 °C) and wind speed (on average of + 4 m/s) (Fig. 5). Similar effects of these recent changes in upwelling systems have already been described64 to understand the aragonite saturation state in monsoon areas in Indonesia65.

The Fig. 8 showed that pH did not have a pattern in the profiling of each station or between stations. The average value between stations ranged from 7.93 to 7.96. On the other hand, TA concentrations showed the highest surface values in all stations (2587–2614 µmol/kg), with the highest concentration in S3 (2614 µmol/kg). [CO32−] seems to contribute significantly (R = 0.95; p < 0.05) to the TA values from the seasons.

pH and Alkalinity of the longitudinal samples from the surface, middle, and bottom water column.

The highest average concentration of \({\text{HCO}}_{3}^{-}\) in S1 (2018 µmol/kg), seems to contribute strongly (R = 0.92; p < 0.05) to the DIC values (2203 µmol/kg). Part of the [\({\text{HCO}}_{3}^{-}\)] in S1 can be explained by the increase in ρCO2 generated by the decomposition of organic matter of natural origin present in the resurgence of water66,67 (reactions 7–9).

SACW also contributes as a source of \({\text{HCO}}_{3}^{-}\), in S1 through the upwelling phenomenon. The values of ΩCalcite, ΩAragonite, and \({\text{CO}}_{3}^{2-}\) were higher on the surface of each station (ΩCalcite = 4.80–5.48; ΩAragonite = 3.10–3.63, and \({\text{CO}}_{3}^{2-}\) = 189–216 µmol/kg), Fig. 9. The mixture of water masses favored by the bathymetry seems to influence the high values of the carbonate system in surface waters. SACW, rich in nutrients compared to other water masses in the region66, provides a photosynthetic activity that affects the carbonate system equilibrium constants2,17, favoring the increase of the saturation state, of the concentrations of the carbonate and the TA. Equilibrium reactions 4 and 5 elucidate this process.

Calcite saturation state (ΩCalcite) and carbonate concentration (µmol/kg) of the longitudinal samples from the water column's surface, middle, and bottom.

Recent research68,69,70 suggested that Ω from seawater does not control the rate of calcification (calcifying fluid). That is, we cannot merely link the availability of carbonate (affected by the decrease in pH values) to the rate of calcification. Until now, no \({\text{CO}}_{3}^{2-}\) transporter has been found in calcifying organisms (for example, coccolithophorids). Conversely, there is ample evidence of \({\text{HCO}}_{3}^{-}\) transporters. The bicarbonate within the calcifying fluid provides the formation of calcium carbonate (Reaction 6).

Reactions involving the formation of calcium carbonate are dependent on the electrochemical gradient (\({H}^{+}\)) between the marine environment and the organisms' tissue and suitable Ω (higher) in the cytoplasmic fluid. Several gaps still need to be clarified concerning the calcification of marine organisms6,71,72. The pH is a masterful measure for its participation in various chemical equilibrium reactions (proton concentration), through which accurate data on the speciation and quantification of the carbonate system can be obtained at reduced costs1,13,73.

The events of El Niño and La Niña can strongly affect the parameters of the carbonate system since these parameters are affected by SACW. The direct influence is the bicarbonate transport, while the indirect influence is the nutrient transport by the SACW stimulating the photosynthetic activity of the surface waters (reactions 4 and 5), removing the CO2 from the water. It is essential to highlight that in the sampling period (Jan/2016), there was an extreme El Niño event, affecting the upwelling phenomenon and favoring an increase in SST and a decrease in the speed of NE and SW winds. The opposite scenario occurred before 2010 when the SST averages were higher, and the wind speeds were lower than in the following years (Figs. 5, 6, and 7).

The values of the carbonate system discussed above showed that the studied area has a low corrosive effect, unlike the other upwelled waters generally presenting pH values below 7.80 and reduced carbonate saturation state (< 1.0)1. The low corrosive effect was probably due to the water masses mixture (SACW and TW) and the low-intensity upwelling by decreased NE wind speeds 5 days before the collection (Figs. 5 and 6), characterizing the effects of the intense El Niño on the Cabo Frio upwelling system.

In temporal sampling (S2), the effect of sunlight on the carbonate system parameters over 12 h was not elucidated (Figs. 10 and 11). The mean values of pH, TA, DIC, ρCO2, ΩCa, and ΩAr were similar in the day/night samples, respectively: pH = 7.96/7.97; TA = 2.626/2.632 µmol/kg; DIC = 2.231/2.215 µmol/kg; ρCO2 = 523/518 µmol/kg; ΩCa = 4.6/4.6, and ΩAr = 3.0/3.0. In this time scale scenario, the physical variables (temperature and salinity) were determinants in the dynamics of the carbonate system's parameters concerning the biological processes (respiration and photosynthesis). It was also evident, through Fig. 10, that there was a mixture between the TW layers (temperature above 18 °C and salinity 36–37.37 psu) and SACW (temperature below 21 °C and salinity between 34.29 and 36.19 psu) over the sampling period.

Temporal sampling. Temperature (°C), salinity (psu), aragonite, and calcite saturation state. X-axis: hour–day. Y-axis: coordinates of stations. Z-axis: Depth-m. Data sampled from 1:00 p.m. to 1:00 a.m. on the 20th and 21st of January 2016.

Temporal sampling. Total Alkalinity (µmol/kg), bicarbonate (µmol/kg), Dissolved Inorganic Carbon (µmol/kg), Potential Hydrogenionic (pH). Dissolved Oxygen (µmol/kg) and Carbon Dioxide Partial Pressure (atm). X-axis: hour–day. Y-axis: coordinates of stations. Z-axis: Depth-m. Data sampled from 1:00 p.m. to 1:00 a.m. on the 20th and 21st of January 2016.

The averages of the carbonate system parameters between the surface (TW) and bottom (SACW) layers were distinct, respectively: ΩCa = 5.2 and 3.9; ΩAr = 3.4 and 2.6; TA = 2642 and 2624 µmol/kg; \({\text{HCO}}_{3}^{-}\) = 2026 and 2217 µmol/kg; DIC = 2182 and 2284 µmol/kg; pH = 7.99 and 7.94; OD = 218 and 181 µmol/kg, and ρCO2 = 483 and 566 µmol/kg. The higher DIC values in SACW were due to the higher values of ρCO2 and HCO3− (2026–2217 µmol/kg), already discussed in the carbonate system section. Reaction 1 explains the higher \({\text{HCO}}_{3}^{-}\) concentrations in this water mass. The raised ρCO2 found in SACW was caused by the decomposition of autochthonous organic matter in upwelled waters66,67 (reactions 7–9). Higher concentrations of \({\text{CO}}_{3}^{2-}\) on the surface (189–202 µmol/kg) seem to contribute to the higher values of the calcite and aragonite saturation state (Reaction 6) affected by the resurgence over 12 h. The higher values of the carbonate saturation state on the surface (TW) are favored by the photosynthetic activity that increases the concentration of carbonate ions (see Reaction 1).

The model used in Figs. 10 and 11 to spatialize the temperature and the carbonate system over the collection time (X) and depth (Y) clearly shows a mixture between TW and SACW water masses. Reduced upwelling intensity (20 to January 21, 2016) suggested the influence of the NE wind during the day. The lower intensity of the resurgence in the collection period from January 20 to January 21, 2016 (Fig. 6) suggests that the effects of the local NE wind blowing intensely towards the continent (onshore) during the day by a more significant temperature gradient between ocean and continent, can also promote the entrapment or stacking of surface water masses along the coast. On the other hand, with increasing water temperature, the winds blow towards the ocean (offshore) at night, accentuating the resurgence effect of SACW and its approach to the coast. This resurgence process during the night can be observed by the lower values of temperature and omega (Fig. 10). Furthermore, the higher TA, DIC, and [\({\text{HCO}}_{3}^{-}\)] values and pH decrease during the night endorse the effect of the upwelled SACW (Fig. 11).

The carbonate system of the present study, when compared to other study areas, showed a lower corrosive effect of SACW, probably due to the phenomena of El Niño and La Ninã over the Cabo Frio upwelling decreasing its intensity (see item 3.1.2). Studies1 of the carbonate system carried out in the South Atlantic (Ocean Basin Sergipe/Alagoas, Brazil) showed a more corrosive water from SACW (250 m deep), with the following values: pH = 7.74 mol/kg-sun; ρCO2 = 953 µmol/kg; ΩCa = 1.9; OD = 146 µmol/kg, and ΩAr = 1.2. Another study65 pointed to the effect of upwelling events on the aragonite saturation state, showing lower values (2.97–3.44) when compared to the area not affected by upwelling (4.45–3.57). An earlier study5 also described the same corrosive effect of upwelled acidified water on the continental shelf. All these authors suggested the corrosive effect of upwelled waters (pH < 7.75 and ΩAr < 1.0).

Recently, coastal areas have been attracting significant attention due to CO2 fluxes74,75,76,77. These researchers generated these CO2 fluxes as they were motivated by attempts to understand the effect CO2 has on ocean acidification, and coastal areas' role in CO2 sequestration. Irrespective of their motivations, the values generated by these CO2 fluxes are inaccurate in terms of the interface between water and air2,61,78. Indeed, fundamental physical K, such as the Schmidt number and the gaseous exchange coefficient (KT) components of the flux equation, are limited when considered in field experiments once they are empirically generated in controlled laboratory conditions. Such K do not consider other factors that interface with the CO2 flux, such as: water/air turbulence interface, air bubbles, surfactant substances, and precipitation. What is the associated error for estimating the flow with these physical factors that are not considered and the assumed laboratory conditions? The average flow of CO2 from water to the atmosphere in this work was 0.12 mmol/m2/day, which is relatively low when compared to the Caravelas estuary (5–1.377 mmol/m2/day)16, and average flow of CO2 found in upwelled waters off the coast of Chile (1.6–2.19 mmol/m2/dia)79.

This study presented new preliminary data on the dynamics of "acidification" in the Cabo Frio upwelling governed by the El Niño and La Niña phenomena, which control SST. Furthermore, it demonstrated the importance of implementing seasonal studies with protocols involving more sensitive analytical techniques with greater precision and reduced cost to generate more data and promote knowledge about the carbonate system along the Brazilian coast1.

The El Niño Southern Oscillation might affect the upwelling intensity by raising the SST. This climate phenomenon influences the carbonate system since these parameters are affected by SACW. Calcium carbonate seems responsible for the increased values of TA in all studied stations. On the other hand, in the S1 station, DIC was influenced by bicarbonate concentrations from upwelled waters. In temporal sampling (S2), the effect of sunlight on the carbonate system parameters over 12 h was not observed.

Lastly, the results of our study emphasize the importance of ENOS phenomena, and the nodal cycle should be considered in studies of acidification of the oceans in resurgence areas. Interdisciplinary studies with the implementation of a specific protocol for temporal and seasonal scales are necessary to understand the effects on biota and support climate change models.

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Ramos e Silva, C. A. et al. Inventory of water masses and carbonate system from Brazilian's northeast coast: Monitoring ocean acidification. PLoS ONE 17, e0271875 (2022).

Article CAS PubMed PubMed Central Google Scholar

Ramos e Silva, C. A., Dávalos, P. B., Da Silva, M. P. & de Miranda, L. B. Variability and transport of inorganic carbon dioxide in a tropical estuary. J. Oceanogr. Mar. Res. 05, 155 (2017).

Google Scholar

Chester, R. Marine Geochemistry (Academic Division of Unwin Hyman Ltd, 1990). https://doi.org/10.1007/978-94-010-9488-7.

Book Google Scholar

Cai, W. J. et al. Acidification of subsurface coastal waters enhanced by eutrophication. Nat. Geosci. 4, 766–770 (2011).

Article ADS CAS Google Scholar

Feely, R. A., Sabine, C. L., Hernandez-Ayon, J. M., Ianson, D. & Hales, B. Evidence for upwelling of corrosive ‘acidified’ water onto the continental shelf. Science (80-). 320, 1490–1492 (2008).

Article ADS CAS Google Scholar

Fassbender, A. J., Sabine, C. L. & Feifel, K. M. Consideration of coastal carbonate chemistry in understanding biological calcification. Geophys. Res. Lett. https://doi.org/10.1002/2016GL068860 (2016).

Article Google Scholar

Smith, H. E. K. et al. Predominance of heavily calcified coccolithophores at low CaCO3 saturation during winter in the Bay of Biscay. Proc. Natl. Acad. Sci. USA. 109, 8845–8849 (2012).

Article ADS CAS PubMed PubMed Central Google Scholar

Takahashi, T. The fate of industrial carbon dioxide. Science (80-). 305, 352–353 (2004).

Article CAS Google Scholar

Hofmann, G. E. et al. The effect of Ocean acidification on calcifying organisms in marine ecosystems: An organism-to-ecosystem perspective. Annu. Rev. Ecol. Evol. Syst. https://doi.org/10.1146/annurev.ecolsys.110308.120227 (2010).

Article Google Scholar

Watson, S. A. Giant clams and rising CO2: Light may ameliorate effects of ocean acidification on a solar-powered animal. PLoS ONE https://doi.org/10.1371/journal.pone.0128405 (2015).

Article PubMed PubMed Central Google Scholar

United Nations. United Nations Summit on Sustainable Development 2015: Informal summary. in 70th Session of the General Assembly 1–10 (2015).

Garcia, H. E. et al. Data management strategy to improve global use of ocean acidification data and information. Oceanography 28, 226–228 (2015).

Google Scholar

Marion, G. M. et al. pH of seawater. Mar. Chem. 126, 89–96 (2011).

Article CAS Google Scholar

McLaughlin, K. et al. Core principles of the California current acidification network: Linking chemistry, physics, and ecological effects. Oceanography 28, 160–169 (2015).

Article Google Scholar

Millero, F. J. Carbonate constants for estuarine waters. Mar. Freshw. Res. 61, 139 (2010).

Article ADS CAS Google Scholar

Ramos e Silva, C. A. et al. Acidificação dos oceanos em um sopro: Prática educacional para construção de conhecimento das mudanças globais. Rev. Experiências em Ensino Ciências 12, 1–25 (2017).

ADS Google Scholar

Ramos e Silva, C. A., Sternberg, L. S. L., Dávalos, P. B. & de Souza, F. E. S. The impact of organic and intensive farming on the tropical estuary. Ocean Coast. Manag. 141, 55–64 (2017).

Article Google Scholar

Lübbecke, J. F., Burls, N. J., Reason, C. J. C. & Mcphaden, M. J. Variability in the South Atlantic anticyclone and the Atlantic Niño mode. J. Clim. https://doi.org/10.1175/JCLI-D-14-00202.1 (2014).

Article Google Scholar

Silva, L. C. C. Estudo Das Massas D’água e da Circulação Geostrófica na Região Sudeste da Bacia Do Brasil. COPPE/UFRJ (Universidade Federal do Rio de Janeiro, 2006).

Valentin, J. L. Analyse des paramètres hydrobiologiques dans la remontée de Cabo Frio (Brésil). Mar. Biol. https://doi.org/10.1007/BF00392407 (1984).

Article Google Scholar

Dereczynski, C. P. & Menezes, W. F. Meteorologia da bacia de campos. in Meteorologia e Oceanografia 1–54 (Elsevier, 2015). https://doi.org/10.1016/B978-85-352-6208-7.50008-8.

Reboita, M. S., Ambrizzi, T., Silva, B. A., Pinheiro, R. F. & da Rocha, R. P. The South Atlantic subtropical anticyclone: Present and future climate. Front. Earth Sci. 7, 1–15 (2019).

Article ADS Google Scholar

Gilliland, J. M. & Keim, B. D. Position of the South Atlantic anticyclone and its impact on surface conditions across Brazil. J. Appl. Meteorol. Climatol. 57, 535–553 (2018).

Article ADS Google Scholar

Sun, Y. et al. OCO-2 advances photosynthesis observation from space via solar-induced chlorophyll fluorescence. Science. 358, 1–6 (2017).

Article Google Scholar

Marshall, G. J. Trends in the Southern Annular Mode from observations and reanalyses. J. Clim. 16, 4134–4143 (2003).

2.0.CO;2" data-track-action="article reference" href="https://doi.org/10.1175%2F1520-0442%282003%29016%3C4134%3ATITSAM%3E2.0.CO%3B2" aria-label="Article reference 25" data-doi="10.1175/1520-0442(2003)0162.0.CO;2">Article ADS Google Scholar

Wolter, K. & El Timlin, M. S. Niño/Southern Oscillation behaviour since 1871 as diagnosed in an extended multivariate ENSO index (MEI.ext). Int. J. Climatol. 31, 1074–1087 (2011).

Article Google Scholar

Da Silveira, I. C. A., Schmidt, A. C. K., Campos, E. J. D., De Godoi, S. S. & Ikeda, Y. A corrente do Brasil ao largo da costa leste brasileira. Braz. J. Oceanogr. 48, 171–183 (2000).

Article Google Scholar

Emílsson, I. The shelf and coastal waters off southern Brazil. Bol. do Inst. Ocean. 11, 101–112 (1961).

Article Google Scholar

Almeida, R. A. F. F., Nobre, P., Haarsma, R. J. & Campos, E. J. D. D. Negative ocean-atmosphere feedback in the South Atlantic Convergence Zone. Geophys. Res. Lett. 34, L18809 (2007).

Article ADS Google Scholar

da Silva, A. C., Araújo, M. & Bourlès, B. Variação sazonal da estrutura de massas de água na plataforma continental do Amazonas e área oceânica adjacente. Rev. Bras. Geofísica 23, 145–157 (2005).

Article Google Scholar

Stramma, L. & England, M. On the water masses and mean circulation of the South Atlantic Ocean. J. Geophys. Res. Ocean. 104, 20863–20883 (1999).

Article ADS Google Scholar

Campos, E. J. D., Velhote, D. & da Silveira, I. C. A. Shelf break upwelling driven by Brazil Current Cyclonic Meanders. Geophys. Res. Lett. 27, 751–754 (2000).

Article ADS Google Scholar

Kampel, M. Estimativa da produção primária e biomassa fitoplanctônica através de sensoriamento remoto da cor do oceano e dados in situ na costa sudeste brasileira. Situ (2003).

Castro, F. de & Miranda, L. de. Physical oceanography of the western Atlantic continental shelf located between 4 N and 34 S. in The Sea vol. 11 209–251 (Wiley, 1998).

Barbiéri, E. B. Ritmo climático e extração de sal em Cabo Frio. Rev. Bras. Geogr. 37, 23–109 (1975).

Google Scholar

Barbiéri, E. B. Cabo Frio e Iguaba Grande, dois microclimas distintos a um curto intervalo espacial. In: Lacerda, L.D. de; Araújo, D.S.D. de; Cerqueira, R. & Turcq, B. (orgs.). Restingas: Origem, Estrutura, Processos. An. do Simpósio sobre Restingas Bras. CEUFF, Niterói. p. 413–425. 477 (1984).

Isaaks, E. H. & Srivastava, R. M. Spatial continuity measures for probabilistic and deterministic geostatistics. Math. Geol. https://doi.org/10.1007/BF00892982 (1988).

Article MathSciNet MATH Google Scholar

Gnanadesikan, R. Methods for statistical data analysis of multivariate observations. J. Am. Stat. Assoc. 92, 1656 (1997).

Article Google Scholar

NOAA. GeoPlatform. NOAA Climate. https://noaa.maps.arcgis.com/home/group.html?id=b4e31886df6740b68f35953cfe15079f#overview (2019). Accessed 2019.

DOE. Handbook of methods for the analysis of the various parameters of the carbon dioxide system in sea water; version 2. ORNL/CDIAC-74 (Eds Dickson, A G & Goyet, C.) 2 (1994).

Pratt, K. W. Measurement of pHT values of Tris buffers in artificial seawater at varying mole ratios of Tris: Tris·HCl. Mar. Chem. 162, 89–95 (2014).

Article CAS Google Scholar

Van Den Berg, C. M. G. & Rogers, H. Determination of alkalinities of estuarine waters by a two-point potentiometric titration. Mar. Chem. 20, 219–226 (1987).

Article Google Scholar

Dickson, A. G. Certificate of analysis. Reference material for oceanic CO2 measurements Batch 104 (Bottled on June 11, 2010). Mar. Phys. Lab. 1–3 (2010).

Hunter, K. A. The temperature dependence of pH in surface seawater. Deep. Res. Part I Oceanogr. Res. Pap. https://doi.org/10.1016/S0967-0637(98)00047-8 (1998).

Article Google Scholar

Ramos e Silva, C. A., Liu, X. & Millero, F. J. Solubility of siderite (FeCO3) in NaCl solutions. J. Solution Chem. 31, 97–108 (2002).

Article Google Scholar

Perez, F. F. & Fraga, F. A precise and rapid analytical procedure for alkalinity determination. Mar. Chem. 21, 169–182 (1987).

Article CAS Google Scholar

Millero, F. J. et al. The use of buffers to measure the pH of seawater. Mar. Chem. https://doi.org/10.1016/0304-4203(93)90199-X (1993).

Article Google Scholar

Clayton, T. D. & Byrne, R. H. Spectrophotometric seawater pH measurements: total hydrogen ion concentration scale calibration of m-cresol purple and at-sea results. Deep. Res. Part I. https://doi.org/10.1016/0967-0637(93)90048-8 (1993).

Article Google Scholar

Millero, F. J. The pH of estuarine waters. Limnol. Oceanogr. 31, 839–847 (1986).

Article ADS CAS Google Scholar

Millero, F. J. Thermodynamics of the carbon dioxide system in the oceans. Science (80-). 59, 661–677 (1995).

CAS Google Scholar

Millero, F. J. Chemical Oceanography 2nd edn. (CRC Press, Taylor & Francis, 1996).

Google Scholar

Lueker, T. J., Dickson, A. G. & Keeling, C. D. Ocean pCO2 calculated from dissolved inorganic carbon, alkalinity, and equations for K1 and K2: Validation based on laboratory measurements of CO2 in gas and seawater at equilibrium. Mar. Chem. 70, 105–119 (2000).

Article CAS Google Scholar

Dickson, A. G. Thermodynamics of the dissociation of boric acid in synthetic seawater from 273.15 to 318.15 K. Deep Sea Res. Part A. Oceanogr. Res. Pap. 37, 755–766 (1990).

Article ADS CAS Google Scholar

Millero, F. J. The physical chemistry of natural waters. Pure Appl. Chem. 57, 1015–1024 (1985).

Article CAS Google Scholar

Yao, W. & Millero, F. J. The chemistry of the anoxic waters in the Framvaren Fjord. Norway. Aquat. Geochem. 1, 53–88 (1995).

Article CAS Google Scholar

Roy, R. N. et al. The dissociation constants of carbonic acid in seawater at salinities 5 to 45 and temperatures 0 to 45 °C. Mar. Chem. 44, 249–267 (1993).

Article CAS Google Scholar

Millero, F. J. Chemical Oceanography 3rd edn, Vol. 3 (CRC Press, Taylor & Francis, 2006).

Google Scholar

Millero, F. J., Graham, T. B., Huang, F., Bustos-Serrano, H. & Pierrot, D. Dissociation constants of carbonic acid in seawater as a function of salinity and temperature. Mar. Chem. 100, 80–94 (2006).

Article CAS Google Scholar

Jähne, B. et al. On the parameters influencing air–water gas exchange. J. Geophys. Res. 92, 1937 (1987).

Article ADS Google Scholar

Jähne, B., Heinz, G. & Dietrich, W. Measurement of the diffusion coefficients of sparingly soluble gases in water. J. Geophys. Res. 92, 10767 (1987).

Article ADS Google Scholar

Wanninkhof, R. Relationship between wind speed and gas exchange over the ocean. J. Geophys. Res. 97, 7373 (1992).

Article ADS Google Scholar

Elias, L. M. C. Variabilidade Interanual da Ressurgência de Cabo Frio. UFRJ (Universidade Federal do Rio de Janeiro, 2009).

Yasuda, I. Impact of the astronomical lunar 18.6-y tidal cycle on El-Niño and Southern Oscillation. Sci. Rep. 8, 15206 (2018).

Article ADS PubMed PubMed Central Google Scholar

Huang, J. et al. Integrating hierarchical bioavailability and population distribution into potential eco-risk assessment of heavy metals in road dust: A case study in Xiandao District, Changsha City, China. Sci. Total Environ. 541, 969–976 (2016).

Article ADS CAS PubMed Google Scholar

Xue, L. et al. Aragonite saturation state in a monsoonal upwelling system off Java, Indonesia. J. Mar. Syst. 153, 10–17 (2016).

Article Google Scholar

Coelho-Souza, S. A., López, M. S., Guimarães, J. R. D., Coutinho, R. & Candella, R. N. Biophysical interactions in the cabo frio upwelling system, Southeastern Brazil. Braz. J. Oceanogr. 60, 353–365 (2012).

Article Google Scholar

Guenther, M. et al. Plankton trophic structure and particulate organic carbon production during a coastal downwelling-upwelling cycle. Mar. Ecol. Prog. Ser. https://doi.org/10.3354/meps07458 (2008).

Article Google Scholar

Allemand, D., Tambutté, É., Zoccola, D. & Tambutté, S. Coral calcification, cells to reefs. in Coral Reefs: An Ecosystem in Transition 119–150 (Springer Netherlands, 2011). https://doi.org/10.1007/978-94-007-0114-4_9.

Cyronak, T., Schulz, K. G. & Jokiel, P. L. The Omega myth: What really drives lower calcification rates in an acidifying ocean. ICES J. Mar. Sci. 73, 558–562 (2016).

Article Google Scholar

Taylor, A. R., Brownlee, C. & Wheeler, G. L. Proton channels in algae: Reasons to be excited. Trends Plant Sci. 17, 675–684 (2012).

Article CAS PubMed Google Scholar

Hurd, C. L. et al. Ocean acidification as a multiple driver: How interactions between changing seawater carbonate parameters affect marine life. Mar. Freshw. Res. 71, 263 (2019).

Article Google Scholar

Waldbusser, G. G., Hales, B. & Haley, B. A. Calcium carbonate saturation state: On myths and this or that stories. ICES J. Mar. Sci. 73, 563–568 (2016).

Article Google Scholar

Bockmon, E. E. & Dickson, A. G. An inter-laboratory comparison assessing the quality of seawater carbon dioxide measurements. Mar. Chem. 171, 36–43 (2015).

Article CAS Google Scholar

Troxler, T. G. et al. Component-specific dynamics of riverine mangrove CO2 efflux in the Florida coastal Everglades. Agric. For. Meteorol. 213, 273–282 (2015).

Article ADS Google Scholar

Bates, N. R. Seasonal variability of the effect of coral reefs on seawater CO2 and air-sea CO2Exchange. Limnol. Oceanogr. 47, 43–52 (2002).

Article ADS CAS Google Scholar

Zhai, W., Dai, M., Cai, W.-J., Wang, Y. & Hong, H. The partial pressure of carbon dioxide and air–sea fluxes in the northern South China Sea in spring, summer and autumn. Mar. Chem. 96, 87–97 (2005).

Article CAS Google Scholar

Zhai, W. & Dai, M. On the seasonal variation of air–sea CO2 fluxes in the outer Changjiang (Yangtze River) Estuary, East China. Sea. Mar. Chem. 117, 2–10 (2009).

Article CAS Google Scholar

Borges, A. V. et al. Variability of the gas transfer velocity of CO2 in a macrotidal estuary (the Scheldt). Estuaries 27, 593–603 (2004).

Article CAS Google Scholar

Lefèvre, N. Observations of pCO2 in the coastal upwelling off Chile: Spatial and temporal extrapolation using satellite data. J. Geophys. Res. 107, 3055 (2002).

Article ADS Google Scholar

Download references

The authors declare that no funds, grants, or other support were received during the preparation of this manuscript. All authors certify that they have no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript.

Postgraduate Program in Dynamics of Oceans and Earth, Federal Fluminense University, Niteroi, 24230-971, Brazil

Carlos Augusto Ramos e Silva, Livia Viana de Godoy Fernandes, Thaise Machado Senez Mello & Nicole Silva Caliman Monteiro

Department of Marine Biology, Federal Fluminense University, Niteroi, 59075-970, Brazil

Carlos Augusto Ramos e Silva

Center for Study of Water, Biomass, and Oil (NAB), Federal Fluminense University, Niterói, 24210-330, Brazil

Carlos Augusto Ramos e Silva, Humberto Marotta, Anderson Araújo Rocha & Raimundo Nonato Damasceno

Agricultural School of Jundiaí, Federal University of Rio Grande do Norte, Macaíba, 59280-000, Brazil

Flavo Elano Soares de Souza

Sedimentary and Environmental Processes Laboratory (LAPSA), Federal Fluminense University, Niterói, 24210-346, Brazil

Humberto Marotta

Admiral Paulo Moreira Marine Research Institute, Arraial do Cabo, 28930-000, Brazil

Flavio da Costa Fernandes, Ricardo Coutinho & Lohengrin Dias de Almeida Fernandes

Ecosystems and Global Change Laboratory (LEMG), International Laboratory of Global Change (LINCGlobal), Niterói, Brazil

Ludmila Caetano dos Santos

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

C.A.R.e.S.—Project administration, conceptualization, methodology, field survey, data acquisition, investigation, writing (original draft and review), supervision, and validation. L.V.G.F.—Data collection, formal analysis, writing (original draft). F.E.S.S.—Formal analysis, writing (original draft, review, and editing). H.M.—Field survey, sample preparation, data acquisition, and writing (original draft). F.C.F.—Field survey, sample preparation, data acquisition, and writing (original draft). T.S.S.-M.—Data curation, visualization, writing (original draft, review, editing, and final validation), submitting agent. N.S.C.M.—Formal analysis, writing (original draft). A.A.R.—Formal analysis, writing (original draft). R.C.—Data acquisition, formal analysis, conceptualization, supervision, and validation. L.D.A.F.—Data acquisition, formal analysis, and investigation. R.N.D.—Investigation, supervision, validation, writing (review and editing). L.C.S.—Data acquisition, formal analysis, and writing (original draft). All authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Correspondence to Carlos Augusto Ramos e Silva.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Silva, C.A.R.e., de Godoy Fernandes, L.V., de Souza, F.E.S. et al. Carbonate system in the Cabo Frio upwelling. Sci Rep 13, 5292 (2023). https://doi.org/10.1038/s41598-023-31479-x

Download citation

Received: 24 November 2022

Accepted: 13 March 2023

Published: 31 March 2023

DOI: https://doi.org/10.1038/s41598-023-31479-x

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.