Dataset: Houston Galveston Bay pCO2

Field Sampling

Galveston Bay is a semi-enclosed microtidal estuary located in the northwestern Gulf of Mexico (nwGOM) (Montagna, Palmer, & Pollack, 2013). With an average water depth of 3 m and a surface area of 1554 km², Galveston Bay is the seventh largest estuary in the U.S. and the second largest estuary on the Texas coast (Bass, Torres, Irza, Proft, Sebastian, Dawson, Bedient, 2018; Morse et al., 1993; Solis & Powell, 1999). Galveston Bay receives freshwater from the Trinity River, San Jacinto River, Clear Creek, and smaller bayous and creeks, with the Trinity River providing 70% of the freshwater entering the Bay (Bass et al., 2018; Morse et al., 1993; Solis & Powell, 1999). The Bolivar Peninsula and Galveston Island separate Galveston Bay from the Gulf of Mexico (GOM), with exchange of water between the Bay and the GOM occurring through Bolivar Roads, the mouth of the Bay (Glass, Rooker, Kraus, & Holt, 2008).

Monthly cruises were conducted between October 2017 and September 2018 aboard the R/V Trident. The timing of the study allowed for examination of the factors regulating CO2 flux over the course of a year following Hurricane Harvey in late August 2017. Although the study began more than 45 days after the hurricane (the residence time of the Bay), salinity recovery of the Bay was likely still ongoing in the inner and middle sections (Du & Park, 2019; Du, Park, Dellapenna, & Clay, 2019).

During each monthly survey, a transect was run between five water sampling stations, extending northwest from the Bay mouth (Station 1) to the Five Mile Marker on the Houston Ship Channel (Station 5). One offshore cruise in the nwGOM outside Galveston Bay was conducted in October 2018. Underway pCO2 measurements were taken along a northwesterly transect from stations 1 through 5. A SUPER-CO2 System equipped with a LI-COR® LI-840A infrared gas analyzer was used to collect both water and air xCO2 after drying through a Peltier thermoelectric device. The xCO2 data, after removing residual water vapor (Honkanen et al., 2021), was converted to pCO2 at sea surface temperature assuming 100% water vapor pressure (Jiang et al., 2008). Underway seawater was taken from a steel pipe attached to the side of the research vessel, as it did not have a dedicated water intake system, and a diaphragm water pump was used to feed water to the equilibrator. In situ sea surface temperature (SST) and salinity were measured with a SeaBird Scientific SBE45® Thermosalinograph mounted parallel to the equilibrator of the SUPER-CO2 System. Prior to and following each sampling trip, the SUPER-CO2 System was calibrated using standards of known CO2 concentrations (273.3, 774.3, and 1468.7 ppm).

To calculate the pCO2 of seawater and air from measurements, the measured mole fraction of CO2 in seawater (xCO2, water) and measured equilibrator barometric pressure and xH2O were first used to calculate xCO2 in dry air (xCO2, air). This xCO2, air was then converted to pCO2 of equilibration (pCO2, eq) using measured temperature of equilibration (Teq) and water vapor pressure of equilibration, which was calculated from salinity and Teq according to methods outlined in Weiss and Price (1980). Next, SST and Teq were used to convert pCO2, eq to pCO2, water (Weiss & Price, 1980). For pCO2, air, xCO2, air was converted to pCO2, air using water vapor pressure at SST and salinity, assuming 100% humidity (Borges et al., 2004).

Meteorological Data

Three National Oceanic and Atmospheric Administration (NOAA) buoys from throughout Galveston Bay provided six-minute interval averages of continuous wind speed data (NOAA, 2022). The average wind speed for all three buoys during sampling times was calculated and applied to the timing of sampling in Galveston Bay. Prior to calculations, wind speeds were converted to a height of 10 m (u10) using the wind profile power law (Hsu, Meindl, & Gilhousen, 1994):

u1/u2 = (z1/z2)^P

where u2 is wind speed at height z2 = 10 m, u1 is the collected wind speed data at height z1, and the exponent P (0.11) around the GOM area is extracted by Hsu et al. (1994).

United States Geological Survey (USGS) streamgages for the Trinity River (gage #08066500) and San Jacinto River, east fork (SJE; gage #08070200) and west fork (SJW; gage #08068000), were used to obtain freshwater discharge (USGS, 2021). These stations were identified as the closest gages to the mouths of the rivers having complete discharge data for the period of study. Discharges of less than or equal to 45 days (residence time of the Bay) prior to flux estimates were utilized (Bass et al., 2018; Morse et al., 1993). The Texas Commission on Environmental Quality (TCEQ) performs routine water quality monitoring, and TCEQ water sampling stations were used for river endmember values from the San Jacinto (average of west fork station #11243 and east fork station #11238) and Trinity (station #10896) rivers (TCEQ, 2022). River endmember DIC was calculated from TA and pH measurements using K1 and K2 constants from Millero (1980), and pH values on the NBS scale. Seasonally weighted averages were calculated by summing the TA or DIC concentration multiplied by daily discharge values for all river measurements of that season and dividing by the sum of all discharge values for all river measurements of that season (using meteorological seasons).

Historical Data

Results from this study were compared to historical data for Galveston Bay obtained from the Surface Ocean CO2 Atlas (SOCAT) database, which provided fCO2, water and xCO2, air values, along with surface seawater salinity, temperature, and depth, with observations from 2006 and 2010 through 2016, primarily during the month of September (Bakker et al., 2016). SOCAT transects followed a similar route to our study transect, beginning near Station 4 and continuing outward into the GOM, with a side transect through the Galveston Channel, which separates Pelican Island from Galveston Island. fCO2 values were converted to pCO2 using the R package seacarb (Gattuso et al., 2022). SOCAT data were analyzed independently from the results of this study. As done previously with ship data, SOCAT xCO2, air was converted to pCO2, air by accounting for water vapor pressure based on SST and SSS, assuming 100% humidity (Borges et al., 2004).