Lateral, Vertical, and Temporal Variability of Seawater Carbonate Chemistry at Hog Reef, Bermuda

Spatial and temporal carbonate chemistry variability on coral reefs is influenced by a combination of seawater hydrodynamics, geomorphology, and biogeochemical processes, though their relative influence varies by site. It is often assumed that the water column above most reefs is well-mixed with small to no gradients outside of the benthic boundary layer. However, few studies to date have explored the processes and properties controlling these multi-dimensional gradients. Here, we investigated the lateral, vertical, and temporal variability of seawater carbonate chemistry on a Bermudan rim reef using a combination of spatial seawater chemistry surveys and autonomous in situ sensors. Instruments were deployed at Hog Reef measuring current flow, seawater temperature, salinity, pHT, pCO2, dissolved oxygen (DO), and total alkalinity (TA) on the benthos, and temperature, salinity, DO, and pCO2 at the surface. Water samples from spatial surveys were collected from surface and bottom depths at 13 stations covering ∼3 km2 across 4 days. High frequency temporal variability in carbonate chemistry was driven by a combination of diel light and mixed semi-diurnal tidal cycles on the reef. Daytime gradients in DO between the surface and the benthos suggested significant water column production contributing to distinct diel trends in pHT, pCO2, and DO, but not TA. We hypothesize these differences reflect the differential effect of biogeochemical processes important in both the water column and benthos (organic carbon production/respiration) vs. processes mainly occurring on the benthos (calcium carbonate production/dissolution). Locally at Hog Reef, the relative magnitude of the diel variability of organic carbon production/respiration was 1.4–4.6 times larger than that of calcium carbonate production/dissolution, though estimates of net organic carbon production and calcification based on inshore-offshore chemical gradients revealed net heterotrophy (−118 ± 51 mmol m–2 day–1) and net calcification (150 ± 37 mmol CaCO3 m–2 day–1). These results reflect the important roles of time and space in assessing reef biogeochemical processes. The spatial variability in carbonate chemistry parameters was larger laterally than vertically and was generally observed in conjunction with depth gradients, but varied between sampling events, depending on time of day and modifications due to current flow.


INTRODUCTION
Coral reefs are some of the most ecologically diverse and economically valuable ecosystems in the world, providing ecosystem services such as coastline protection, fisheries provisioning, tourism revenue, nutrient cycling, and habitat creation that amount to billions of dollars annually (e.g., Spurgeon, 1992;Moberg and Folke, 1999;Brander et al., 2007;Costanza et al., 2014). However, anthropogenic stressors on local and global scales, including temperature rise (e.g., Hughes et al., 2003), overfishing (e.g., Jackson et al., 2001), hypoxia (Nelson and Altieri, 2019), and ocean acidification (e.g., Hoegh-Guldberg et al., 2007) are threatening reef function around the world. Coral reefs may be particularly susceptible to the effects of ocean acidification due to the large presence of marine calcifiers whose calcium carbonate (CaCO 3 ) structures serve as the reef foundation (Kleypas et al., 1999;Hoegh-Guldberg et al., 2007;Kleypas and Yates, 2009). Increasingly acidified conditions resulting from changes in seawater carbonate chemistry are predicted to drive declines in calcification by organisms such as corals (e.g., Langdon et al., 2000; for a metanalysis see Chan and Connolly, 2013) and increases in bioerosion and CaCO 3 dissolution rates (e.g., Andersson and Gledhill, 2013;Schönberg et al., 2017) which could result in many reefs potentially shifting to net dissolving instead of net calcifying by the end of the century (Hoegh-Guldberg et al., 2007;Silverman et al., 2009;Eyre et al., 2018).
An individual reef 's response to acidified conditions, however, may be strongly dependent on the biogeochemistry of the particular system . Locally, carbonate chemistry is in part controlled by the metabolic processes occurring on the benthos such as calcification, CaCO 3 dissolution, primary production, and respiration Kleypas et al., 2011). The net sum of total calcification and CaCO 3 dissolution is referred to as Net Community Calcification (NCC); whereas, the net sum of primary production and total respiration is referred to as Net Community Production (NCP). Changes in NCC affect the balance of both total alkalinity (TA) and dissolved inorganic carbon (DIC) in a ratio of 2:1, whereas changes in NCP in reef environments mainly affect DIC (Zeebe and Wolf-Gladrow, 2001;Andersson and Gledhill, 2013). Consequently, positive NCC results in a decrease in seawater pH while positive NCP results in an increase in pH, and vice versa for negative NCC and NCP, respectively (Zeebe and Wolf-Gladrow, 2001). Concurrent with these metabolic processes, physical factors such as the reef geomorphology, hydrodynamics, and residence time also influence local variability in carbonate chemistry, as these physical properties control the flow trajectory and contact time between seawater and the benthos Zhang et al., 2012;Falter et al., 2013;Lowe and Falter, 2015;Page et al., 2018). It has been hypothesized that these metabolic and physical processes together may modulate local seawater chemistry to the extent they can either alleviate or exacerbate acidification over reefs compared to the open ocean (Bates et al., 2010;Anthony et al., 2011Anthony et al., , 2013Andersson et al., 2014;Cyronak et al., 2014). However, how these biological and physical drivers interact and contribute to local seawater carbonate chemistry variability in different coral reef habitats and across vertical, lateral, and temporal scales has not been fully characterized Takeshita et al., 2018).
Many previous studies have demonstrated carbonate chemistry variability as a function of time at the surface or at the bottom in numerous coral reef locations (e.g., Gray et al., 2012;Sutton et al., 2014;DeCarlo et al., 2017) with fewer studies characterizing the spatial and lateral variability across different habitats (Suzuki and Kawahata, 2003;Watanabe et al., 2006;Andersson et al., 2014). In terms of vertical variability, most studies assume that the water column above a shallow reef is well-mixed with small to no gradients above the benthic boundary layer (e.g., Falter et al., 2013;Page et al., 2018). However, given that flow regime, residence time, and benthic productivity have the potential to form strong gradients in water chemistry from upstream to downstream Takeshita et al., 2018), there may also be gradients between surface and bottom that are not often quantified (Watanabe et al., 2006;Andersson et al., 2007). Spatial gradients in any direction also change through time, on hourly to seasonal (e.g., Gagliano et al., 2010;Gray et al., 2012;Guadayol et al., 2014) and annual to decadal time scales (e.g., Pelejero et al., 2005;Wei et al., 2015). This variability has been characterized in Bermuda over shorter timescales (days) (e.g., Takeshita et al., 2018) and longer timescales (seasons to years) (e.g., Bates et al., 2010;Andersson et al., 2014;Yeakel et al., 2015;Bates, 2017). However, there is a scarcity of studies examining vertical, lateral, and temporal variations in reef carbonate chemistry simultaneously.
Compared to low latitude, tropical coral reefs, the high latitude reefs in Bermuda experience significant seasonal and inter-annual variations in seawater chemistry as a result of variations in seawater temperature and local reef metabolism Courtney et al., 2017), as well as larger scale ocean processes linked to the North Atlantic Oscillation (NAO) (Yeakel et al., 2015). In addition, this ocean region has been identified as an area with one of the fastest changing chemical environments under ocean acidification . Furthermore, it has been suggested that the coral community structure of Bermuda, which is made up of coral taxa with relatively stress-tolerant life histories, may represent one potential stable state of coral reefs in a future high CO 2 world , making it an important area to study. Here, the main objective was to compare the spatial variability of seawater chemistry laterally across a reef, at surface and at depth, and through time, in order to capture variability in four dimensions. We also estimated the relative contributions of biogeochemical processes to variability in DIC and TA locally and relative to the open ocean endmember during the study period. Understanding the variability of local seawater chemistry, both spatially and temporally, as well as the drivers of this variability, will be important for predicting the scope and magnitude of ocean acidification impacts on coral reefs.

Study Site
The study was conducted at Hog Reef (32.457 • N, 64.835 • W), a rim-reef on the northern edge of the Bermuda carbonate platform (Figure 1). Depth varies across Hog Reef, ranging from ∼4 to 25 m with a mean (± 1 SD) of 10.3 ± 3.3 m, generally increasing offshore but with high structural complexity and variability (Courtney et al., 2016;Takeshita et al., 2018 ; Figures 1, 2). Mean (± 1 SD) benthic cover at Hog Reef has been previously described as 35 ± 3% macroalgae, 27 ± 5% hard coral, 20 ± 4% turf algae, 16 ± 2% soft coral, and < 5% sand, rock, rubble, and coralline algae (Marine Environmental Program, 2007;Courtney et al., 2016). Residence time of the water on this rim reef has been estimated to be about 1-4 days using various tracer methods (Morris et al., 1977;Venti et al., 2012). Tides in this area can be characterized as mixed semi-diurnal, with a mean tidal range of 0.90 ± 0.11 m. Tide information (mean sea level, MSL; m) data were accessed from the NOAA Tides and Currents Database for the nearest tide station (St. Georges Island, Bermuda Station 1 ). Bathymetry data was accessed from Sutherland et al. (2014) and the satellite image of Hog Reef is courtesy of the Digital Globe Foundation.

Autonomous Sensors
In order to capture high frequency variability and diel cycles on the reef, a suite of autonomous sensors was deployed at Hog Reef (Figure 1d and Table 1). To characterize water flow, an Aquadopp acoustic doppler current profiler (ADCP, 1 MHz; Nortek) was deployed at the bottom (∼7 m water depth) from September 1-7, 2017. The Aquadopp measured current speeds and directions in 0.5 m depth bins from 0.4 m above the profiler to 5.4 m above the profiler. However, due to interference in the first two cells by nearby biota and reef structures, only cells from 1.4 to 5.4 m above the profiler were analyzed. Averages of 1 Hz data over 1 min were stored at 10 min intervals.
Seawater pH T (defined on the total H + scale), pCO 2 , and TA were measured with a suite of instruments from Sunburst Sensors between September 1-22, 2017, including a Submersible Autonomous Moored Instrument for pH (SAMI-pH), two SAMI-CO 2 that also recorded temperature and dissolved O 2 with Aanderaa 4531 Oxygen Optodes, a CTD (Microcat, Sea-Bird Scientific), and the recently developed SAMI for Alkalinity (SAMI-alk), which is one of the first in situ instruments for autonomous measurements of total alkalinity (Spaulding et al., 2014). SAMI-pH accuracy and precision are ± 0.005 and ± 0.001, respectively, based on repeat analyses of tris buffer. SAMI-CO 2 accuracy and precision are ± 5 and ± 1 µatm, respectively, based on comparisons with an infrared CO 2 analyzer. The Aanderaa Oxygen Optode accuracy is reported as < 1 µM. Accuracy for the CTD temperature and conductivity is ± 0.002 • C and ± 0.0003 S m −1 , respectively. Spaulding et al. (2014) found an accuracy and precision for SAMI-alk of −2 ± 13 µmol kg −1 based on comparison with discrete samples collected during a deployment. All instruments collected measurements every 15 minutes, except 1 https://tidesandcurrents.noaa.gov/stationhome.html?id=2695540#info the SAMI-alk which collected measurements every 2 h (with one measurement omitted per day due to in situ calibrations at varying times, i.e., 11 measurements per day). Initial offsets of sensor data were corrected based on water samples collected at the time of deployment and no other offsets were applied. Initial SAMI-CO 2 calibrations are referenced to an infrared method calibrated with NOAA certified gases. See Table 1 for a summary of the deployment specifics.
In addition to the Sunburst Sensors instruments, time series data for air and seawater pCO 2 were obtained from a NOAA PMEL MAPCO2 buoy moored on Hog Reef (32.457 • N, 64.835 • W). The MAPCO2 buoy measures the mole fraction of CO 2 (xCO 2 ) in seawater and air equilibrated with seawater every 3 h using a LI-COR LI-820 CO 2 infrared gas analyzer calibrated with a reference gas that meets World Meteorological Organization standards (Sutton et al., , 2019. Accuracy and precision of these xCO 2 measurements are estimated to be < ± 2 and < ± 0.7 µmol mol −1 , respectively . The buoy is also equipped with a Sensirion SHT71 sensor to measure relative humidity and temperature . Sea surface temperature and salinity were also continuously recorded by a Sea-Bird Electronics 16 plus V2, with accuracies of ± 0.005 • C and ± 0.0005 S/m for temperature and conductivity, respectively. Air and seawater pCO 2 were calculated from these parameters using equations from Zeebe and Wolf-Gladrow (2001), with uncertainty estimated to be < 2.0 µatm .

Spatial Surveys of Carbonate Chemistry
Autonomous moored sensor measurements were coupled with discrete water samples collected at 13 stations covering an area of ∼3 km 2 on Hog Reef (Figure 1). A total of 4 spatial surveys were completed via boat, each over a span of ∼2 h, with afternoon surveys (∼13:00-15:00) on September 1 and September 4 and morning surveys (∼10:00-12:00) on September 4 and September 7. Spatial surveys consisted of three overlapping transects, with the PMEL MAPCO2 buoy as a common center point, and spanned the width of the rim reef (Figures 1b,c). Using a 5 liter Niskin bottle, water samples were collected at the surface (0.5-1.0 m) and bottom (4-20 m depending on the station's local water depth) at all stations, with one intermediate sample collected at the deepest stations (Figure 1d). In situ temperature and salinity were measured with a YSI 556 Handheld Multiparameter Instrument (accuracy: temperature ± 0.15 • C, salinity ± 1%). Water samples for DIC and TA (N = 114) were collected according to best practices in 250 ml Pyrex glass sample bottles and immediately fixed with 100 µl of a saturated HgCl 2 solution (Dickson et al., 2007).
DIC was analyzed using an Automated Infra-Red Inorganic Carbon Analyzer (AIRICA, Marianda) with a Li-COR 7000 as detector. TA was analyzed via open-cell potentiometric acid titration system developed by the Dickson Lab at Scripps Institution of Oceanography (Dickson et al., 2007). Some samples (n = 29) were reanalyzed using a small volume, opencell Metrohm Titrando automatic titrator system to validate results due to occasional occurrences of non-ideal titration curves. Accuracy (± precision) of TA and DIC measurements were calculated as the mean offsets (± SD) from Certified Reference Materials (CRM) provided by the Dickson Lab and were 1.3 ± 3.5 µmol kg −1 (n = 57) for DIC, −0.46 ± 1.5 µmol kg −1 (n = 24) for the Dickson TA system, and 1.76 ± 3.0 µmol kg −1 (n = 21) for the Metrohm TA system. CO2SYS for Excel (Lewis and Wallace, 1998) and MATLAB (van Heuven et al., 2011), were used to calculate additional carbonate chemistry parameters using the first and second dissociation constants of carbonic acid (K 1 and K 2 ) from Mehrbach et al. (1973) refit by Dickson and Millero (1987), the dissociation constants of bisulfate (K HSO4 ) from Dickson (1990) and total boron by Uppström (1974), and total scale pH (pH T ).

Current Data
Current profiler data was collected in earth coordinates and corrected for the magnetic declination of the earth at the deployment site to account for differences between magnetic and true north. The magnetic declination at Hog Reef (32.457 • N 64.835 • W) on 1 September 2017 was 14.74 • W ± 0.33 • of true north 2 . All data are reported in compass coordinates where 0 • is north and 90 • is east. A principal component analysis was then performed on the corrected, depth-averaged u (east) and v (north) current velocities for both the entire duration of the deployment, as well as for individual days using the princax function in MATLAB (developed by Rich Signell, United States Geological Survey). The angle of maximum variance (θ) and the major and minor axes of a principal ellipse were computed to identify the angle containing the most variance. Mean depthaveraged flow velocity was also calculated and denoted on plots with a quiver arrow emanating from the instrument deployment location, with angle and length denoting current direction and speed, respectively.

Variability and Frequency Analyses
Variability for each parameter recorded by the autonomous sensors was calculated as the mean daily range (maximumminimum) (from 00:00 h to 23:45 h, or as close to these times as sampling frequency would allow) ± one standard deviation. Variability in measured parameters from the spatial FIGURE 2 | Depth-averaged current velocity (a) over the entire deployment (September 1-7) and for individual days corresponding to dates of the water chemistry spatial surveys (b) September 1, (c) September 2, (d) September 4, and (e) September 7 overlaid on the bathymetry of the deployment site. White stars correspond to individual velocity measurements over that time period and the arrow beginning at the black dot corresponds to the mean velocity vector. The ellipse represents the variability (1 SD) of the velocity on that day (i.e., the velocity variability (1 SD) traces out an ellipse starting at the black dot pointing toward the ellipse-for example, over the full deployment period, the mean flow is approximately off-reef to the north but varies strongly in the east/west direction, with minimal variability southwards). The origin for the velocity plots (individual measurements, ellipse, and mean velocity vector) is at the mooring represented by the black circle on each plot and x and y-axis scales of 0.1 m/s are included in light gray. Variables measured are: xCO 2 , mole fraction of carbon dioxide; pCO 2 , partial pressure of carbon dioxide; T, temperature; DO, dissolved oxygen; S, salinity; pH T , total scale pH; TA, total alkalinity; current speed; current direction; and P, pressure.
surveys was calculated as the mean range of a parameter ± 1 standard deviation on one sampling date over the surface or bottom (for lateral variability) across all sampling dates. Vertical variability was calculated as the mean of the absolute differences in a parameter between surface and bottom across the survey area for one sampling date ± 1 standard deviation, across all sampling dates. Frequency analyses were performed on all parameters measured by the autonomous sensors in order to determine the dominant frequencies of variability using the plomb function in MATLAB to create Lomb-Scargle periodograms (Figure 3; raw data are shown in Supplementary Figure 1). A high-pass filter (highpass; pass band frequency of 1 day −1 ) was used to remove additional low frequency signals and distill the diel frequency signal in temperature at surface and bottom. A Savitzky-Golay filter (sgolayfilt; order 3 and frame length 7) and then a high-pass filter (highpass; pass band frequency of 1 day −1 ) was used to distill the frequency signals in total alkalinity at 7 m as well as calculated relative NCC (r NCC-TA ). These filtered total alkalinity, r NCC-TA , and temperature data were only used to create the Lomb-Scargle periodograms (Figure 3, right column) which are plotted next to the raw, unfiltered time series data (Figure 3, main column; for a comparison, see Supplementary Figures 4-6).

Biogeochemical Drivers of Variability
Instrument time series data were used to approximate the relative contribution from biogeochemical processes to the observed variability in DIC and TA at Hog Reef. A time series of DIC at 7 m was calculated from the time series of TA (from the SAMI-alk at 7 m) and pCO 2 (from the SAMI-CO 2 at 7 m) FIGURE 3 | (Left) Temporal variability in environmental conditions and seawater chemistry on Hog Reef for September 2017 at surface (lighter colors) and bottom (darker colors). Temperature, salinity, dissolved oxygen, TA, pH T , and pCO 2 were measured with autonomous moored sensors. Measurements of bottom temperature (at 7 m depth) from multiple sensors were essentially indistinguishable, thus only one was plotted. Measured TA and calculated pH T and pCO 2 from discrete TA and DIC seawater samples were plotted as circles ("SIO" samples were analyzed in the Andersson Lab at Scripps Institution of Oceanography and "UM" samples were analyzed in the DeGrandpre Lab at the University of Montana Missoula). Mean sea level data were accessed from the NOAA Tides and Currents Database (St. George's Island). Relative NCC and NCP (r NCC-TA , r NCP-DO , and r NCP-DIC ) were estimated from sensor-derived water chemistry data using a mean anomaly (as described in the text). Gray bars indicate nighttime hours. (Right) Lomb-Scargle periodograms of autonomous moored sensor parameters as well as calculated r NCC-TA , r NCP-DO , and r NCP-DIC . Figure 3), interpolated to a frequency of every 2 h from 1 September 2017 at 15:00 to 22 September 2017 at 11:00 in order to match the lower sampling frequency of the SAMI-alk. Note that DIC calculated from pH T and TA also gave similar DIC values. TA, DIC, and DO were normalized (nTA, nDIC, nDO) to mean salinity from the interpolated time series (36.59 PSU).

(Supplementary
For a well-mixed reef habitat, NCC is routinely calculated based on Eulerian measurements of TA according FIGURE 4 | Lateral spatiotemporal variability at the surface in temperature, DIC, TA, and pH T (calculated from measured DIC and TA) on September 1, 2, 4, and 7. Plots also contain bathymetry contours (black) and the mean depth-averaged velocity vector averaged over the sampling period (red).
to (e.g., Silverman et al., 2007;Langdon et al., 2010;Koweek et al., 2015a): (1) where nTA refers to the difference between nTA on the reef and offshore, τ is the mean residence time, ρ is the seawater density, and h is the mean height of the water column. Furthermore, NCP can be calculated based on measurements of DIC and TA or DO with corrections for gas exchange: where nDIC and nDO refer to the differences between the reef and offshore, and F CO2,gas and F O2,gas represent sea-air CO 2 and O 2 gas exchange, respectively. Relative to biogeochemical processes on coral reefs with relatively short residence times (hours to days), sea-air CO 2 gas exchange is often small and therefore assumed to be negligible while O 2 flux is more strongly affected (e.g., Yeakel et al., 2015). Here, we used two different approaches to assess the biogeochemical processes underlying the observed chemical condition and variability at Hog Reef. In one approach, NCC and NCP were calculated according to Eqs 1 and 2 relative to the offshore values of TA and DIC measured at the Bermuda Atlantic Time-series Station (BATS), adopting the best estimates of the mean depth (10.3 ± 3.3 m) and residence time at Hog Reef (2.5 ± 0.4 days; Venti et al., 2012). However, because neither τ nor h were directly measured in this study and even small errors in these parameters can generate large errors in NCC and NCP (Courtney and Andersson, 2019), these calculations are admittedly uncertain. A Monte Carlo analysis was used to propagate the potential errors of τ and h through iterative calculations of NCC and NCP using randomly generated numbers (n = 10,000) within the range of uncertainties for each parameter. The resulting mean and range were used as approximations of the mean ± uncertainty in NCC and NCP. Furthermore, it is important to recognize that these estimates of NCC and NCP based on inshore-offshore chemical gradients represent integrated metabolic signals over multiple tidal cycles and habitats beyond the Hog Reef study site because currents here are not unidirectional and seawater is moving back and forth across this site multiple times before escaping (Venti et al., 2012;Takeshita et al., 2018).
In a second approach, the local biogeochemical variability was assessed based on relative changes in NCC and NCP, rather than actual rates. In the absence of detailed τ or h measurements, the relative contributions of NCC and NCP can be assessed by their ratio, termed R OI by Suzuki (1998), which represents the ratio of organic to inorganic carbon production (Suzuki, 1998;Suzuki and Kawahata, 2003): However, in the present case, to highlight the diel variability and the relative magnitude of the numerator and denominator of Eqs 4 and 5 at Hog Reef, the temporal variabilities of these terms were assessed independently relative to the local mean over the sampling period rather than the offshore conditions. We refer to these terms as r NCC-TA , r NCP-DIC , and r NCP-DO : It is important to note that r NCC-TA , r NCP-DIC , and r NCP-DO are not metabolic rates, nor do they reflect the net trophic or calcification status of the water at Hog Reef. Instead these terms simply highlight the approximate relative magnitude of variability in TA and DIC attributed to NCC and NCP at Hog Reef and its surrounding habitats on diel timescales. The spatial representation of these terms roughly corresponds to the hypothetical reef area covered by a parcel of seawater over one tidal cycle. This has previously been referred to as the "effective reef footprint" which translates to the spatial and temporal extent of the benthic communities' influence on the overlying seawater chemistry (Courtney et al., 2016). The temporal variability in r NCC-TA , r NCP -DIC , and r NCP-DO was calculated as the mean daily range (from 00:00 h to 00:00 h, or as close to these times as sampling frequency would allow) ± 1 standard deviation.

Fine Scale Temporal and Vertical Variability
Vertical profiles of current speeds and directions as well as pressure recorded at Hog Reef clearly showed the mixed semidiurnal tidal cycle, with a mean current speed of 0.08 ± 0.04 m s −1 (Supplementary Figure 2). Tide station data revealed a mean tidal range of 0.90 ± 0.11 m during the deployment (Figure 3 Consequently, water flowed predominantly offshore at Hog Reef (Figure 2a), with significant variance in the east/west (approximately along-reef) direction resulting in an elliptical benthic footprint with a semi-major axis of 1.2 km and semiminor axis of 0.8 km based on the hypothetical distance traveled by a parcel of water over 6 hours. Additional principal component analyses of the depth-averaged velocities for individual days corresponding to our spatial survey sampling dates revealed that θ varied between these days, changing from 79 • on September 1 (Figure 2b), 103 • on September 2 (Figure 2c), 86 • on September 4 (Figure 2d), and 120 • on September 7 (Figure 2e). The direction of the mean quiver for each day revealed that mean flow varied from along-reef and slightly inshore on September 1 (Figure 2b), to off-reef to the northeast on September 2 (Figure 2c), to slightly off-reef to the northwest on September 4 (Figure 2d), and off-reef approximately north on September 7 (Figure 2e) with differing degrees of variance. High frequency measurements by autonomous sensors at Hog Reef showed distinct diel variability in temperature, DO, pH T , and pCO 2 , but less so with respect to TA, which was confirmed by a frequency analysis of all data (Figure 3, right column). A secondary local maximum and minimum of some of these parameters were also apparent on most days (Figure 3), with secondary peaks in the power spectra at ∼2 day −1 , reflecting the influence from the mixed semidiurnal tide (Figure 3). This was particularly prominent for the TA data. Surface and bottom (7 m) temperatures generally agreed, with the exception of September 3-4 with a maximum vertical difference of 0.84 • C (Figure 3). Mean temperature (± 1 SD) during the deployment was 28.38 ± 0.40 • C at the surface and 28.36 ± 0.38 • C at the bottom. The mean diel temperature range (± 1 SD) in surface and bottom waters over the month-long deployment was 0.38 ± 0.15 and 0.40 ± 0.11 • C, respectively (Figure 3). Salinity remained relatively constant throughout the deployment with a mean (± 1 SD) of 36.64 ± 0.06 PSU at the surface and 36.59 ± 0.04 PSU at the bottom. We note that there was drift between the salinity measurements at the surface and 7 m beginning midway through the deployment resulting in higher mean salinity at the surface compared to the bottom. Calculations of density at the surface and bottom produced an unstable water column during the latter half of the deployment, with a mean density (± 1 SD) over the deployment of 1023.6 ± 0.2 kg m −3 at the surface and 1023.5 ± 0.1 kg m −3 at the bottom, suggesting that the drift was a sensor error and not a real signal. Given a lack of post-deployment cross-calibration between instruments, we cannot confirm which sensor experienced the drift, though the drift may be an artifact of potential biofouling due to the absence of a pump for the CTD deployed at the bottom.
DO concentrations at the surface during peak production hours were higher than at the bottom, though the mean (± 1 SD) concentrations at surface and bottom over the deployment were similar (226 ± 22 and 220 ± 18 µmol kg −1 , respectively). The mean diel range (± 1 SD) in DO, however, appeared qualitatively different at the surface (78 ± 29 µmol kg −1 ) and bottom (54 ± 16 µmol kg −1 ) (Figure 3). Furthermore, mean TA and pH T (± 1 SD) at the bottom were 2331 ± 16 µmol kg −1 and 7.953 ± 0.025, with mean diel ranges (± 1 SD) of 26 ± 13 µmol kg −1 and 0.07 ± 0.02, respectively (Figure 3; no autonomous measurements were made at the surface of these parameters, Table 1). Changes in pH T positively tracked changes in DO (Spearman's rank correlation coefficient, ρ = 0.61, p < 0.001) which also mirrored changes in pCO 2 over the deployment (pH T and pCO 2 : ρ = −0.99, p < 0.001; pCO 2 , and DO: ρ = −0.66, p < 0.001). Mean seawater pCO 2 (± 1 SD) measured by the SAMI-CO 2 instruments was 518 ± 35 µatm at the surface and 516 ± 30 µatm at the bottom, with mean diel ranges (± 1 SD) of 70 ± 25 and 76 ± 28 µatm, respectively (Figure 3), the differences of which are within the instruments' accuracy. Mean buoy seawater pCO 2 measured at the surface, however, was slightly lower (508 ± 30 µatm; mean ± 1 SD) than that measured by the SAMI-CO 2 at the surface, with the difference between these values outside of the buoy instrument's accuracy. However, buoy pCO 2 had a similar mean diel range (± 1 SD) of 70 ± 25 µatm compared to the SAMI-CO 2 instruments. Buoy pCO 2 measured in air was relatively constant (387 ± 3 µatm; mean ± 1 SD) with a small diel signal (5 ± 6 µatm; mean range ± 1 SD), which however, was disproportionally influenced by two anomalous measurements on September 11 and 16th which were outside the range of expected accuracy and precision.
Overall, parameters that were measured (TA) or calculated (pH T and pCO 2 ) from discrete bottle samples collected at the instrument site during spatial surveys showed good agreement with instrument data for the majority of samples, though there were disagreements on a few occasions (Figure 3). The degraded precision compared to the fundamental replicate reproducibility is likely due to spatial and temporal mismatches between the sensors and discrete sampling. Bottle samples were collected as close to the sensors as physically possible from a boat, but small differences in sampling depth as well as water movement and/or potential boundary layers around the instruments may have led to these observed disagreements. In addition, there is also some error introduced to measurements of pCO 2 and pH T calculated from DIC and TA that may contribute in part to these differences (Millero, 2007).

Lateral Spatiotemporal Variability
Spatial surveys over Hog Reef revealed noticeable variability in most parameters across the 3 km 2 survey area. Some lateral gradients appeared to follow inshore to offshore depth gradients, but the structure of these gradients varied between sampling dates and with changes in the direction of bulk water flow as measured by the current meter (Figures 4, 5). Mean surface temperature (± 1 SD) across the four surveys was 28.38 ± 0.15 • C at the surface and 28.28 ± 0.14 • C at the bottom whereas mean salinity (± 1 SD) was the same at both surface and bottom (36.6 ± 0.0 PSU) at the precision of the YSI. The mean ranges (± 1 SD) of these parameters were 0.41 ± 0.14 • C and 0.1 ± 0.0 PSU across the surface, and 0.38 ± 0.12 • C and 0.1 ± 0.0 PSU across the bottom across the four surveys, respectively (Figures 4, 5). Temperature was lower on the northwestern side of the survey area over the deeper reef drop off, and higher in the southeast above the shallow reef lagoon area (Figures 4, 5) as would be expected due to depthdependent heating.
DIC and TA were lower over the shallower, inshore areas of the reef, and higher over the deeper drop off (Figures 4, 5).

Vertical Spatiotemporal Variability
Vertical variability was small for all parameters of the spatial survey area and with no clear spatial trend (Figure 6). The mean temperature and salinity differences (± 1 SD) between surface and bottom were 0.11 ± 0.07 • C and 0.0 ± 0.0 PSU, respectively, across the four surveys, which are less than the precision of the YSI 556 Handheld Multiparameter Instrument measurements themselves (Figure 6). Mean absolute differences (± 1 SD) in carbonate chemistry parameters between surface and bottom were 5 ± 2 µmol kg −1 for DIC, 6 ± 3 µmol FIGURE 5 | Lateral spatiotemporal variability at the bottom in temperature, DIC, TA, and pH T (calculated from measured DIC and TA) on September 1, 2, 4, and 7. Plots also contain bathymetry contours (black) and the mean depth-averaged velocity vector averaged over the sampling period (red). kg −1 for TA, and 0.01 ± 0.01 for pH T across the four surveys ( Figure 6) and were close to the precision of the respective measurements.

Biogeochemical Drivers of Carbonate Chemistry Variability
On the local scale, r NCP-DIC , r NCP-DO , and r NCC-TA were calculated relative to a nDIC mean of 2043 µmol kg −1 , nDO mean of 220 µmol kg −1 , and nTA mean of 2331 µmol kg −1 . The variability of r NCP-DO and r NCP-DIC followed a diel cycle, with maximum values in the late afternoon and minimum values in the early morning (as confirmed by the frequency analysis; Figure 3). However, the mean diel range (± 1 SD) in r NCP-DO (51 ± 15 µmol kg −1 ) was larger than the mean diel range for r NCP-DIC (32 ± 11 µmol kg −1 ) (Figure 3). The diel variability in r NCP-DO and r NCP-DIC also corresponded to the diel cycle of pH T measured independently at Hog Reef, i.e., peaks in r NCP-DIC and r NCP-DO corresponded to peaks in pH T (Figure 3).
Diel variability in r NCC-TA was smaller than both r NCP-DO and r NCP-DIC (13 ± 6 µmol kg −1 ; mean range ± 1 SD) and r NCC-TA did not show any strong correlation to either r NCP-DO or r NCP-DIC with a correlation coefficient (r) < 0.1 in both cases. The local diel variabilities in r NCP-DIC and r NCC-TA were also clearly evident in biogeochemical assessments based on inshore-offshore chemical gradients (Figures 7B,C). Offshore TA and DIC at the surface at BATS 3 in September 2017 was 2396 and 2046 µmol kg −1 , respectively, which combined with the mean depth and residence time of the Hog Reef area, revealed NCP of −118 ± 51 mmol m −2 day −1 (mean ± 1 SD) and NCC of 150 ± 37 mmol CaCO 3 m −2 day −1 . Consequently, these calculations indicate that the integrated metabolic signals encompassing multiple tidal cycles and potentially multiple habitats were net heterotrophic and net calcifying, with the diel variabilities characterized on FIGURE 6 | Vertical spatiotemporal variability of the difference between surface and bottom temperature, DIC, TA, and pH T (calculated from measured DIC and TA; bottom values subtracted from surface values) on September 1, 2, 4, and 7. Plots also contain bathymetry contours (black) and the mean depth-averaged velocity vector averaged over the sampling period.
Hog Reef superimposed on these larger scale net metabolic signals (Figure 7).

Fine Scale Temporal and Vertical Variability
The results of this study suggest that variability in the parameters measured at Hog Reef was mainly diel, and was likely influenced by a combination of light, benthic and water column production, calcification and dissolution, and local tidal flow. In particular, temperature, DO, pH T , and pCO 2 exhibited distinct local extremes in the afternoon and after midnight each day. This dominant frequency of variability of ∼ 1 day −1 was linked to the daily light cycle and its effects on seawater temperature and primary production (Figure 3). Seawater pCO 2 variability was inverse to that of DO and pH T , reflecting CO 2 uptake by primary production during the day and release via respiration at night. These diel patterns are ubiquitous for shallow marine systems and have been documented for many reef systems in Bermuda (Bates et al., 2001(Bates et al., , 2010Sutton et al., 2014;Page et al., 2017;Takeshita et al., 2018) and around the globe (e.g., Hofmann et al., 2011;Gray et al., 2012;Guadayol et al., 2014;Koweek et al., 2015a,b;DeCarlo et al., 2017;Page et al., 2018;Cyronak et al., 2020).
In contrast to the autonomous measurements of temperature, DO, pH T , and pCO 2 , which revealed distinct diel trends, autonomous measurements of TA revealed higher frequency temporal variability and a less distinct diel trend (Figure 3). This may be related to the fact that the relative magnitude of the diel variability in calcification and CaCO 3 dissolution (r NCC ) was small compared to that of photosynthesis and respiration (r NCP-DO and r NCP-DIC ) at Hog Reef, making changes in seawater TA less pronounced compared to changes in DIC, DO, pH T , and pCO 2 (Figure 3 and Supplementary Figure 3; Takeshita et al., 2018). Furthermore, frequency analyses revealed that TA varied strongly at frequencies of both 1 day −1 (diel light cycle) as well as 2 day −1 (relating to the tidal cycle) (Figure 3, right). This suggests that flow dynamics had a strong influence on the temporal variability in TA, which depending on the flow trajectory and history, would reveal a chemical signature that did not necessarily track the diel light cycle, and led to the prominence of variability at a higher frequency (Figure 3). We hypothesize that the reason other parameters (i.e., DO, pH T , and seawater pCO 2 ) revealed a strong diel trend that was not muddled by the tidal cycle is because photosynthesis and respiration were also quantitatively important in the water column, as opposed to just the benthos, whereas calcification and CaCO 3 dissolution mainly occurred on the benthos. Consequently, changes in the seawater TA appeared to mainly reflect biogeochemical processes that predominantly occurred in the benthos (calcification and dissolution) with the chemical signal modified by hydrodynamics, while changes in DIC, DO, pH T , and pCO 2 were also strongly influenced by water column photosynthesis and respiration. This hypothesis is supported by the fact that there were notable differences in the magnitude of DO variability at the surface and the bottom throughout the deployment period. While DO minima at night were similar for surface and bottom, daytime maxima were much higher in the surface, suggesting that there was higher total net primary productivity measured at the surface compared to the benthos during the day. Furthermore, the largest differences between surface and bottom DO occurred during neap tides, potentially indicating amplified gradients during times of reduced water flow and mixing (e.g., Falter et al., 2013;Lowe and Falter, 2015). Long et al. (2019) reported that water column production at Hog Reef accounted for more than half (58%) of the total oxygen flux, compared to just 39% from the benthos, during an overlapping study interval with our observations (Long et al., 2019), lending further support to the hypothesized contribution from water column productivity compared to the benthos.
The observed gradient in DO between surface and bottom was at times partly corroborated by independent pCO 2 measurements by the SAMI-CO 2 at the bottom and the MAPCO2 buoy at the surface. Surface pCO 2 measured by the buoy was consistently lower than pCO 2 measured at 7 m by the SAMI-CO 2 , as would be expected given the DO gradient between surface and 7 m. However, because of inconsistencies between measurements of pCO 2 by the MAPCO2 and SAMI-CO 2 at the surface, it is difficult to make any strong conclusions based on these observations. It should also be noted that variations in seawater temperature and NCC between day (net calcification) and night (decreased calcification and increased CaCO 3 dissolution) influence seawater pCO 2 opposite to changes resulting from the day-night cycle in NCP (Cyronak et al., 2020), resulting in a dampened pCO 2 diel amplitude. Consequently, pCO 2 measurements are less sensitive relative to the DO measurements as a measure of net productivity differences between the surface and bottom. On coral reefs in general, the benthos is assumed to be the main contributor to coral reef organic carbon production with only minor influence from the water column. Few studies have actually quantified the relative contributions of each (e.g., Long et al., 2019), but based on the present observations and the Long et al. (2019) study, it is clear that this needs to be assessed in more detail in future studies. The relative contributions from the water column and benthos to NCP are likely to change as a function of depth so that the contribution from the benthos becomes increasingly important as the depth decreases. This may in fact explain the observations from Cyronak et al. (2020) that coral reef pH T variability can be largely predicted based on depth and mean pH T conditions. Compared to our observations at Hog Reef, other studies that have characterized TA autonomously on coral reef environments in Kāne'ohe Bay, Hawai'i (Spaulding et al., 2014, same sensor as this study; Briggs et al., 2020, solid-state sensor) and Heron Island, Australia (McMahon et al., 2018, VINDTA 3c) observed more distinct diel variability and a greater daily amplitude in TA that closely tracked the diel variability in other carbonate chemistry parameters (i.e., pH T and pCO 2 ). These study locations were characterized by much shallower depths than Hog Reef, and thus, higher benthic biomass to water volume ratio, resulting in greater influence by benthic processes on water chemistry and consequently larger changes in seawater TA. The higher biomass to water volume ratio was also evident from greater daily changes in seawater pH T and pCO 2 at these locations (Drupp et al., 2013;McMahon et al., 2018;Page et al., 2018). It should be noted that making high quality autonomous measurements of TA on coral reefs is an important and active area of research since it provides an additional master variable to constrain the carbonic acid system. Furthermore, because calcification and CaCO 3 dissolution are central processes on coral reefs, and strongly influenced by global environmental change (Hoegh-Guldberg et al., 2007;Hughes et al., 2018), these processes are critical to monitor. Direct measurements of TA (as opposed to calculated TA) offer the most accurate and precise approach to estimate NCC. Thus, high frequency measurements of TA are important for understanding the drivers and dynamics of calcification and dissolution processes on coral reefs and could reveal temporal trends, relationships, and the effects of episodic events (e.g., storms) otherwise not detected from low frequency measurements (Spaulding et al., 2014;McMahon et al., 2018;Briggs et al., 2020).

Lateral and Vertical Spatiotemporal Variability
Lateral spatial gradients in carbonate chemistry across Hog Reef appeared to be influenced by the bathymetry and the general direction of water flow across or along the reef (Figures 5, 6). The observed gradients in DIC, TA, and pH T generally decreased from offshore to inshore as a result of the uptake of TA and DIC via calcification and primary production, with an exacerbated signal in shallower areas where smaller volumes of water were being modified by the reef community (such as in Falter et al., 2013;Takeshita et al., 2018). The current data and spatial survey data together highlight the influence of hydrodynamics on the structure of the observed gradients (Figures 5, 6). On most days, the mean direction of depth-averaged flow was oriented in line with the irregularities in the structure of the inshore-offshore gradients observed (e.g., September 2, 4, 7 in Figures 5, 6). It should also be noted that September 4 was an extremely calm and sunny day, and thus, favorable for generating strong chemical gradients. While current speed and direction do not fully account for the gradients observed, the present results highlight the need to combine spatial water chemistry surveys with complementary hydrodynamic data, in order to better interpret observed chemical gradients. Furthermore, given the much larger gradients observed laterally compared to vertically, we find that the general assumption of a well-mixed water column over coral reefs is mostly a reasonable assumption at the spatial scale and depth of the current study, but the presence of stratification and/or disparate contributions from the benthos and water column biogeochemical processes, may invalidate this assumption. Furthermore, this assumption varies from place to place and needs to be given site-specific considerations.
Lateral spatial gradients in seawater carbonate chemistry have previously been characterized in Bermuda (e.g., Andersson et al., 2014;Venti et al., 2014;Takeshita et al., 2018) and in many other reef systems (e.g., Watanabe et al., 2006;Shamberger et al., 2011;Albright et al., 2015;Courtney et al., 2018;Page et al., 2018). Combined with mean residence time and depth, these gradients provide the information required to calculate rates of NCC and NCP, but accurately characterizing residence time is a non-trivial task associated with large uncertainties (Lowe et al., 2009;Venti et al., 2012;Teneva et al., 2013;Courtney and Andersson, 2019). The present data highlight the potential challenges in quantifying local rates of NCC and NCP based on spatial chemical gradients, as the water on the reef flows in different directions depending on the time of day and depth. In addition, the history of any given water parcel flowing over different habitats depends on the tide, current speeds and directions, and benthic composition. Compared to Hog Reef, there are clearly other locations that are more suitable for this task where even-depth and quasiunidirectional flow or no flow (e.g., slack tide) may assist in quantifying rates of NCC and NCP. Other approaches based on gradients in the benthic boundary layer (e.g., benthic gradient flux or eddy covariance: Long et al., 2013;Takeshita et al., 2016Takeshita et al., , 2018 may also provide a way to accurately characterize rates at locations such as Hog Reef but are limited to a small benthic footprint (∼100 m 2 ). Although the objective of this study was not to quantify metabolic rates, the different approaches explored here to assess biogeochemical drivers (local temporal variability, inshore-offshore gradients, and TA-DIC relationships) revealed the importance of the temporal and spatial scales considered, and how they might be related (Figures 7, 8).

Biogeochemical Drivers of Carbonate Chemistry Variability
Calculations of r NCP-DIC , r NCP-DO , and r NCC-TA using a mean anomaly highlighted the variability of localized processes occurring on Hog Reef and adjacent habitats (i.e., habitat-specific biogeochemical processes and the respective mixing between offshore, rim reef, and lagoon seawater). On average, these terms roughly represented a benthic reef area on the order of 3 km 2 , but there were times when the footprint may have been larger (Figure 2). The magnitude of the local, short-term variability in organic carbon processes (r NCP-DIC and r NCP-DO ) was greater than that of inorganic carbon processes (r NCC-TA ) occurring at Hog Reef during the study period with the diel range of both r NCP-DO and r NCP-DIC exceeding that of r NCC-TA (4.6 ± 2.4 µmol kg −1 for r NCP-DO : r NCC-TA and 2.8 ± 1.4 µmol kg −1 for r NCP-DIC : r NCC-TA ; mean ratio ± SD) (Figure 3). This was in agreement with a previous study at the same site that measured benthic NCC and NCP based on chemical gradients in the benthic boundary layer (Takeshita et al., 2018). Furthermore, the temporal variability in r NCP-DIC agreed well with r NCP-DO calculated from changes in DO, though the amplitude of the diel variability in r NCP-DO was larger than that of r NCP-DIC . This disagreement may be in part due to deviation from the assumed photosynthetic quotient of 1 (Smith and Marsh, 1973;Takeshita et al., 2018;Bolden et al., 2019) with a Type II linear regression between r NCP-DO and r NCP-DIC yielding a slope of 1.8 (p < 0.001), but could also arise from other factors (e.g., sea-air gas exchange was not accounted for and DIC was not measured directly but calculated from pCO 2 and TA; Supplementary Figure 3).
Direct estimates of NCP and NCC based on the chemical gradients between Hog Reef and offshore (BATS), showed substantial net heterotrophy and net calcification ( Figure 7C). These estimates reveal the cumulative influence on seawater chemistry from a reef area potentially larger than Hog Reef and for seawater residence time extending multiple tidal cycles, with the influence from the local daily metabolism superimposed on this larger scale reef metabolic signal. The generality of these results (i.e., heterotrophy and calcification) is probably robust, but the absolute values need to be viewed with caution since the flow trajectory, mean depth, and residence time are associated with substantial uncertainty (Courtney and Andersson, 2019). Nonetheless, the estimated mean rates of NCC and NCP are comparable to previous estimates from Bermuda in September (NCC: ∼50-200 mmol m −2 day −1 , NCP: ∼ 0-−250 mmol m −2 day −1 ; (Bates et al., 2010;Yeakel et al., 2015;Courtney et al., 2016Courtney et al., , 2017Takeshita et al., 2018), although several of these estimates used similar methodology and are thus not quite independent comparisons (Yeakel et al., 2015;Courtney et al., 2016Courtney et al., , 2017, and Takeshita et al. (2018) only measured NCP and NCC over 3 days.
The importance of scale, illuminated by the different approaches used to assess the importance of biogeochemical processes on Hog Reef, was also evident from graphical assessment and correlation analysis of TA-DIC data. In contrast to many other studies that have characterized TA and DIC on coral reefs over temporal or spatial scales (e.g., Cyronak et al., 2018 and references therein), the current study did not observe a strong correlation between TA and DIC across space (Figure 8), which was also the case in Takeshita et al. (2018). The lack of a clear TA-DIC relationship could be explained by both the local dominance of r NCP over r NCC as well as the differences between surface and benthic production (see previous discussion). However, when placed in context with TA and DIC data collected at other reef sites in Bermuda, covering inshore waters from Bailey's Bay to the rim reef at Hog Reef, the data from this study fit well along a linear relationship of nTA and nDIC (normalized to a salinity of 36.7, the mean salinity across FIGURE 7 | (A) r NCC-TA (blue) and r NCP-DIC (black) calculated relative to the local means; (B) r NCC-TA (blue) and r NCP-DIC (black) calculated relative to the offshore value at the BATS Station; (C) Rates of NCC (blue) and NCP (black) (mean ± shaded uncertainty) calculated relative to the offshore value at the BATS Station. data from all studies) with a slope close to 1 (slope = 0.97, R 2 = 0.86; Figure 8). This implies that NCC and NCP are closely balanced at this larger spatial and temporal scale, and that relatively small changes in mean pH T occur as a function of time and space (Suzuki and Kawahata, 2003;Andersson and Gledhill, 2013;Andersson et al., 2014;Cyronak et al., 2018). However, this does not mean that pH T is invariable; for a given habitat, seawater pH T may vary largely on diel timescales depending on the relative importance of the organic to inorganic carbon cycle Takeshita et al., 2018). These observations are directly related to the functional scale (biological, temporal, and spatial) of the study and has been discussed extensively in Cyronak et al. (2018) and Takeshita et al. (2018).

CONCLUDING REMARKS
The results of this study highlight that seawater carbonate chemistry on coral reefs varies considerably over time and space, with variability linked to diel cycles and local hydrodynamics. By measuring variability in carbonate chemistry and other parameters simultaneously over time and at multiple spatial scales, we showed that physical and biological processes interact on a reef on different scales to result in the net variability in carbonate chemistry and environmental parameters. The general assumption of a well-mixed water column with small vertical chemical gradients (compared to lateral or temporal gradients) over shallow coral reefs may be a reasonable assumption for many reefs, but is site dependent and needs to be assessed on a case-by-case basis. Notably, the present study showed a decoupling between the variability of TA and other carbonate chemistry parameters, which could be attributed to different contributions from benthic and water column biogeochemical processes as well as hydrodynamics. This observation deserves additional attention in other reef locations characterized by different depth, flow, and community compositions. In general, because of their relatively high latitude (32 • N), Bermudan reefs experience higher seasonal variability in many properties compared to more tropical reefs. It has been predicted that future changes in ocean pH may increase the variability, and thus the range, of extreme pCO 2 values that reefs are exposed to due to reduced seawater buffering capacity (Jury et al., 2013;Shaw et al., 2013). However, whether naturally more variable reef environments extend protections to corals and reef organisms under a regime of changing ocean pH is unclear (Dufault et al., 2012;Comeau et al., 2014;Camp et al., 2016;Rivest et al., 2017;Kapsenberg and Cyronak, 2019). Characterizing the presentday variability of carbonate chemistry at multiple scales, as well as the interactions between hydrodynamics, geomorphology, and biogeochemistry, on reefs such as those in Bermuda will FIGURE 8 | Salinity normalized (36.7 PSU) TA and DIC for water sample data from sites across the Bermuda carbonate platform: Crescent Reef (triangles), Bailey's Bay (circles), and various expeditions to Hog Reef (squares) ranging from 2010 to 2017 with a Type II linear regression for all data (black) and pH contours. Only data from July to October of each year were used from Courtney et al. (2017;yellow). All data from Takeshita et al. (2016;green) were collected in July 2015. All bottle sample data from the present study were collected in September 2017 (blue).
help us better predict how reefs will change as a result of future acidification.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: BCO-DMO Dataset #783568 https://www.bco-dmo.org/dataset/783568 and BCO-DMO Project #73787 https://www.bco-dmo.org/project/737878.

AUTHOR CONTRIBUTIONS
AA, TC, HP, and MD designed the study. TC, HP, AA, and CB deployed the instruments and collected the data. AP performed the laboratory work and analyzed the data with inputs from AA, TC, and SG. AP wrote the first draft of the manuscript with input from AA. All authors contributed to subsequent drafts of the manuscript.

ACKNOWLEDGMENTS
We thank Dr. Matthew Long and Jennie Rheuban (Woods Hole Oceanographic Institution) for assistance during sensor deployment and recovery and Dr. Tyler Cyronak for processing and providing imagery of Hog Reef.