Due to the logistical complexities of implementing standard gas setups on a ship, we did not adopt the built-in δ13C-CO2 calibration program of the G2131-i CRDS system. Instead, leveraging multiple in-house standards with pre-calibrated δ13C-DIC values facilitated the correction of δ13Cmean inaccuracies. Throughout the analytical period, in-house standards were sub-sampled into 12-mL glass vials weekly and then sent to the UC Davis Stable Isotope Facility for δ13C-DIC analysis. In their approach, DIC in water was converted to headspace CO2 using phosphoric acid and analyzed using headspace equilibration technique with a Thermo Scientific GasBench II and Thermo Finnigan Delta Plus XL isotope-ratio mass spectrometer (IRMS). The δ13C-DIC values, obtained through Gasbench-IRMS method at the facility, were utilized to calibrate the CRDS measurements of δ13C-DIC.
To balance the need for frequent calibrations with the onboard sample processing efficiency, a calibration using one of the three in-house standards was conducted following analysis of every eight seawater samples. This procedure ensured each standard was assessed a minimum of three times daily. We calculated the CO2 concentration-weighted mean δ13C-CO2 (δ13Cmean) for each analysis by incorporating both raw δ13C-CO2 (δ13Craw) data and net CO2 concentration (CO2net) readings from the CRDS at every time point. The δ13Cmean values for each in-house standard, derived from itsadjacent measurements, were used in a time-based linear regression model to track the instrumental drift and estimate the value of the standard’s δ13C signal (δ13Cest) at the time of each sample measurement. This enabled the establishment of a separate three-point calibration curve (R2 > 0.999) for each measurement, incorporating the δ13Cest and the exact δ13C-DIC values of three in-house standards.
In our approach, each sample or reference material was subjected to a minimum of two and up to four consecutive measurements to achieve the preset relative standard deviation (RSD) of 0.001 for the net integration area and 0.06 for the CO2-weighted mean of δ13C-CO2. From these measurements, we selected two "valid" rounds that met our precision criteria, and the final DIC concentrations and δ13C-DIC results were always reported as an average of these two valid rounds. In addition, CRM Batch #197, #199, #201, #202, and #206 were randomly included in the sample sequence as quality checks for δ13C-DIC analysis.
Based on our previous study, the method’s uncertainty is 0.03‰ for the δ13C-DIC value (1σ). The GO-SHIP A16N cruise in 2023 analyzed 320 replicate samples from 150 CTD stations. Excluding 15 pairs of abnormal data with a δ13C-DIC difference greater than 0.2‰, the mean absolute differences was 0.06 ± 0.05‰ for δ13C-DIC (1σ, n = 305), which was within 2σ of the overall uncertainty of the method.