On the Frontline: Tracking Ocean Acidification in an Alaskan Shellfish Hatchery

The invasion of anthropogenic carbon dioxide (CO2) into the ocean is shifting the marine carbonate system such that saturation states of calcium carbonate (CaCO3) minerals are decreasing, and this is having a detrimental impact on early life stages of select shellfish species. The global, secular decrease in CaCO3 saturation states is occurring on top of a backdrop of large natural variability in coastal settings; progressively shifting the envelope of variability and leading to longer and more frequent exposure to adverse conditions. This is a great concern in the State of Alaska, a high-latitude setting vulnerable to rapid changes in the marine carbonate system, where an emerging shellfish industry plans major growth over the coming decades. Currently, the Alutiiq Pride Shellfish Hatchery (APSH) in Seward, Alaska is the only hatchery in the state, and produces many shellfish species with early life stages known to be sensitive to low CaCO3 saturation states. Here we present the first land-based OA measurements made in an Alaskan shellfish hatchery, and detail the trends in the saturation state of aragonite (Ωarag), the more soluble form of CaCO3, over a 10-month period in the APSH seawater supply. These data indicate the largest changes are on the seasonal time scale, with extended periods of sub-optimal Ωarag levels (Ωarag < 1.5) in winter and autumn associated with elevated water column respiration and short-lived runoff events, respectively. The data pinpoint a 5-month window of reprieve with favorable Ωarag conditions above the sub-optimal Ωarag threshold, which under predicted upper-bound CO2 emissions trajectories is estimated to close by 2040. To date, many species in production at APSH remain untested in their response to OA, and the data presented here establish the current conditions at APSH as well as provide a framework for hatchery-based measurements in Alaska. The current and expected conditions seen at APSH are essential to consider for this developing Alaskan industry.


Introduction
A total of nearly 155 Pg (1 Pg = 10 15 g) of carbon have been taken by the ocean due to rising atmospheric carbon dioxide (CO 2 ) levels that result largely from anthropogenic fossil fuel emissions [1]. Ocean CO 2 uptake increases the concentration of hydrogen ions (H + ; lowers seawater pH) while also decreasing carbonate ion (CO 3 2-) levels; reactions that collectively have been termed ocean acidification (OA) [2][3][4][5][6]. Decreasing CO 3 2- where O phase is the saturation state for typically either aragonite or calcite phases of CaCO 3 , Ca 2+ is the calcium concentration, and K sp is the phase specific solubility product. Of these reactions associated with OA, a decrease in O phase is a critical alteration in carbonate chemistry affecting rates of calcification and shellfish growth [7,8]. For the aragonite phase of CaCO 3 , O arag values where sub-lethal sensitivity first become evident for a number of bivalve larvae range between 1.2 and 2.0 [7]. At O arag < 1.0, thermodynamics dictate aragonite will actively dissolve in seawater [9]. An unfavorable condition for aragonite precipitation is a considerable stressor for early life stages of select shellfish species that can ultimately lead to mortality [7,8].
Due to the critical nature of this carbonate system parameter in dictating shellfish production, it is imperative for aquaculture facilities to accurately track O arag in their seawater supplies. In the State of Alaska, having this capacity is essential for an emerging shellfish industry in the midst of a shifting biogeochemical seascape [10].
Over the time scale of a year, ocean uptake of anthropogenic CO 2 results in small incremental chemical changes. On decadal time scales, these changes can shift baselines in marine carbonate chemistry with potentially significant ecosystem level impacts [11][12][13]. This is especially the case in high-latitude settings where anthropogenically-driven reductions in O arag are expected to occur ahead of changing conditions in warmer locales [14,15]. The risk associated with such wholesale biogeochemical shifts in the marine environment will potentially have farreaching effects beyond the ecosystems themselves and the dependence of coastal communities on ocean goods and services [10,15]. Aquaculture hatcheries have the capacity to be important contributors to food security in support of growing populations by ingesting seawater from adjacent coastal settings in order to produce marketable marine species. As such, hatcheries are on the frontline both in terms of being one of the first industries impacted by OA [8,11,15,16] and in their need to track and cope with the changing ocean biogeochemical conditions. It is important to note that in vulnerable coastal settings, secular changes in CO 2 are occurring in conjunction with large, in many cases understudied, carbonate system variability; demonstrated, for instance, by data collected on the United States (U.S.) west coast (Oregon) continental shelf over multiple years [17]. As shown by Harris et al. [17], large carbonate system variability does not diminish the importance of an increasing trend. Effectively, the gradual human-induced changes in O arag shift the envelope of variability such that organisms experience more adverse conditions longer and more often than they would otherwise [17][18][19].
Pacific oyster, Crassostrea gigas, is a cultured species that has been described as a "canary in the coalmine" for being an organism highly impacted by reduced O arag levels in the coastal zone [8,16,20]. Pacific oyster production failures in U.S. Pacific Northwest (PNW) shellfish hatcheries between 2005 and 2009 resulted from the ingestion of low-O arag , high-pCO 2 water, which is naturally upwelled along the coast, into shore-side hatchery facilities. As detailed by Waldbusser et al. [21], adverse O arag conditions dually impact larval Crassostrea gigas during D-hinge (prodissoconch I) shell formation by: (1) exposing calcification surfaces until completion of the D-hinge shell, and (2) increasing kinetic demands that then impacts the energy budget already limited by endogenous resources. The combination of these factors explains the sensitivity to O arag > 1 when aragonite is thermodynamically favored to precipitate. The additional anthropogenic CO 2 signal in upwelled water along the PNW is small relative to natural CO 2 levels in upwelling source waters [22], however large enough to alter the natural envelop of variability in these systems and impart a significant stress that resulted in a near collapse in a $270 Million industry for the State of Washington [23]. The hatchery production failures served as a vital catalyst in two ways [16]: (1) by unifying stakeholders and scientists in establishing robust monitoring techniques to track O arag in hatchery intake water [8], and (2) by instigating state and federal government action to fund scientific activities needed to catalog trends in coastal manifestations of OA along with their impacts on marine ecosystems and the economies they support. To this end, data are presently streaming in near real-time from a number hatchery facilities along the Pacific coast of North America [24], and guidelines have been developed to aid stakeholders in making OA related land-based measurements [25].
In this contribution, we report on land-based determinations of O arag from seawater entering an Alaskan hatchery facility in the northern Gulf of Alaska (GOA): the Alutiiq Pride Shellfish Hatchery (APSH) [26]. APSH is located in Seward, Alaska on the shore of Resurrection Bay (Fig 1) and is currently the only shellfish hatchery in the state. This is expected to change in the near future, as the Alaska Fisheries Development Foundation, Inc. currently seeks to pass an Alaska Mariculture Initiative with the aim of growing the industry to a $1 billion level in 30 years [27]. APSH is the primary producer of Blue and Red Alaskan King crab, Paralithodes platypus and Paralithodes camtschaticus, respectively. The latter crab species has been shown to be sensitive to low O arag levels [28]. APSH also produces other animals that have either previously shown sensitivity to reduced O arag (i.e. Pacific oyster [7,8]) or are believed to be vulnerable during early life stages (Panopea generosa, Crassadoma gigantea, Leukoma  staminea). Prior to the data presented here, it was unclear whether APSH had been experiencing threatening levels of O arag in their seawater supply. APSH is unique owing to its close proximity to an area of intense OA-related sampling from central Resurrection Bay to over the adjacent continental shelf [29]. Full water column measurements have been made in this area every May and September since 2008, and the GOA OA (GAKOA) mooring has been in place at the mouth of Resurrection Bay since 2011 (Fig 1) [30]. Using high-speed measurements of the seawater entering the hatchery and a compiled dataset from the adjacent coastal ocean, trends in O arag were resolved over a 10-month period and projections were made for the changing conditions at this integral site for the state of Alaska's growing shellfish industry.

Materials and Methods
From October 5, 2013 to August 5, 2014, temperature, salinity and CO 2 partial pressure (pCO 2 ) data were collected from the seawater supply flowing into APSH. This seawater supply was provided to the hatchery by the University of Alaska D. W. Hood Laboratory at the Seward Marine Center [31]. Untreated seawater was pumped to the facilities from an intake pipe extending to 75 m depth and 91 m from the laboratory into Resurrection Bay. Two seawater pumps provided flow at~200 gallons per minute (gpm), and a tangential split off the main supply line provided seawater at~1 gpm to the hydrographic and chemical sensors positioned in the preparation room of the D.W. Hood Laboratory. Temperature and salinity (calculated from conductivity) were measured with a Sea-Bird Electronics SBE 45 MicroTSG thermosalinograph [32], with salinity reported in this paper using the Practical Salinity Scale (PSS-78, dimensionless).
Seawater pCO 2 data were calculated from corrected measurements of CO 2 mixing ratio (xCO 2 ) made using a Sunburst Sensors SuperCO 2 System [33] following protocols recommended by Pierrot et al. [34] with the system theory and calculations described in detail elsewhere [35][36][37][38]. Seawater continuously flowed first through the thermosalinograph and then through two (primary and secondary) showerhead equilibrators [34,39]. The primary showerhead equilibrator supplied equilibrated carrier gas (marine air) to a non-dispersive infrared gas analyzer (LI-COR LI840A CO 2 /H 2 O) housed within the SuperCO 2 system's electronics box at a rate of~50 ml min -1 . Unaltered marine air was drawn from 0.25 inch polyethylene tubing that connected a vented water trap positioned outside of the D.W. Hood Laboratory to the secondary equilibrator. Marine air was pre-equilibrated in the secondary equilibrator, and then plumbed to the primary equilibrator as a make-up air supply to replace equilibrated carrier gas provided to the LI-COR. Pressure and temperature were continuously measured in the primary equilibrator using a Honeywell ASCX Microstructure Pressure Sensor and a Minco Fast Response RTD, respectively. Equilibrated carrier gas, three standard gases of known mixing ratio (148, 448 and 748 ppm; Scott-Marin, Inc.), and unaltered marine air were all plumbed to provide gas flow to the SuperCO 2 system's electronics box. The SuperCO 2 system was controlled using National Instruments LabVIEW software run on an HP laptop computer. The software controls data acquisition from the thermosalinograph, the pressure and temperature sensors, and the LI-COR; while also controlling Valco Instruments Co. Inc. (VICI) multi-port actuators that cycle between the gas streams plumbed to the electronics box. None of the gas streams were dried prior to analysis, and all measurements were made at 0.5 Hz. The prescribed measurement scheme controlled by the software was to supply equilibrated carrier gas from the primary equilibrator to the LI-COR continuously for 240 minutes, then cycle the actuators to consecutively allow the three standard gas streams and unaltered marine air to be measured for 90 s at 100 ml min -1 before returning to sample the carrier gas equilibrated with seawater xCO 2 . From each sequence of standard gas measurements, the final 20 s of data in the 90 s interval before the actuator changed position were used to construct calibration functions that were then interpolated in time between standard gas sequences. These functions were then used to calibrate the xCO 2 measurements of seawater equilibrated carrier gas, with an adjustment on the order of 1%. Atmospheric xCO 2 measurements made by the SuperCO 2 system were not used in this analysis, as these data do not enter in to calculations of O arag . Annual mean atmospheric pCO 2 was calculated from the well-resolved data collected on the GAKOA mooring. This annual mean atmospheric pCO 2 value was then used in anthropogenic CO 2 calculations. Corrected seawater xCO 2 was subsequently adjusted for under-pressurization in the primary showerhead equilibrator using the ratio of equilibrator to vented LI-COR cell pressure, and then converted to pCO 2 using atmospheric pressure measured by the LI-COR. The 2-s seawater pCO 2 , temperature and salinity data were quality controlled, and then bin-averaged in 5-min interval bins.
Taking advantage of the expansive and growing dataset of paired total alkalinity (TA) and salinity measurements made throughout the northern GOA coastal ocean (Fig 1; [29]), these data were compiled and used to construct a TA-salinity relationship. The compiled data were collected as part of the growing suite of carbonate system and hydrographic measurements made from oceanographic cruises that have occurred every May and September since 2008. Briefly, seawater samples were collected at specific stations and depths using a rosette equipped with 5-l Niskin bottles and a Sea-Bird Electronics 911Plus conductivity-temperature-depth (CTD) profiler. Seawater was drawn from Niskin bottles into clean 250 ml Pyrex glass reagent bottles using established sampling techniques [40], and treated with 200 μl of saturated HgCl 2 solution to halt biological activity in the samples. TA was measured in the fixed seawater samples at the University of Alaska Fairbanks (UAF) Ocean Acidification Research Center (OARC) using a VINDTA 3C [41]. Seawater certified reference materials (provided by A. G. Dickson, Scripps Institute of Oceanography) were analyzed before samples were processed to ensure measurements were accurate to within 0.1% (2 μmol kg -1 ). A robust linear regression was computed between the TA and salinity data using MathWorks MATLAB software. A small number of TA validation measurements (n = 23) were collected at APSH to verify the credibility of this linear relationship (S1 Table). 5-min binned seawater pCO 2 data, TA derived using the TA-salinity relationship, temperature and salinity were input into a customized MathWorks MATLAB version of CO2SYS [42] for carbonate system calculations. pH on the total hydrogen ion scale (pH T ), CO 3 2and O arag were computed using the equilibrium constants for the dissociation of carbonic acid from Millero [43]. Ca 2+ concentrations were calculated from salinity using the relationship from Riley and Tongudai [44]. For estimating the contribution of anthropogenic CO 2 in setting O arag variability at APSH, the mean atmospheric pCO 2 from GAKOA (398 ± 6.7 μatm) was used to calculate the sea-air difference in pCO 2 (ΔpCO 2 ). While the term "sub-lethal" has been used in the literature to describe the impact of O arag levels < 1.5 on select larval shellfish species [7,8,15], to date, there has been no evidence of die-offs directly linked to O arag conditions at APSH. Here, we consider O arag values 1.5 to be sub-optimal growth conditions for many larval species produced by APSH (e.g. Crassostrea gigas); O arag equal to 1.5 marks the threshold value below which increased stress and vulnerability is expected [7,15].

Results and Discussion
Presented here are the first O arag values from seawater entering an Alaskan shellfish hatchery, and the data indicate unique characteristics setting conditions at APSH apart from those at other hatcheries to the south along the U.S. Pacific coast within the California Current System [8]. In source waters to APSH, temperature and salinity underwent only modest changes throughout the period of observation, with temperature ranging from 5 to 11°C and salinity varying between 26 and 32 (Fig 2). pCO 2 spanned 305 to 463 μatm, and was relatively confined in comparison to the dynamic range seen in surface water measured by the GAKOA buoy at the mouth of Resurrection Bay [45]. At this nearby surface site, only limited pCO 2 data above atmospheric levels have been observed to occur between February and March [30]. The most prolonged periods of high pCO 2 levels at APSH were from January to April (Fig 2). Note that record gaps in April, May and June ranged from 2 to 20 days, and were caused by either electrical power outages in Seward or technical issues associated with the SuperCO 2 system. The TAsalinity relationship built from compiled GOA data [29] was robust (Fig 3) as indicated by a high r 2 value (0.94), low root mean square error (RMSE; 17.21 μmol kg -1 ), and close agreement with our validation samples (S1 Table). O arag , calculated using TA from the TA-salinity relationship with the pCO 2 data, showed a smaller dynamic range of relative to the measurements made within PNW hatcheries [8], with only a 0.8 unit oscillation about the 1.5 O arag threshold over the 10-month record (Fig 4). The largest abrupt changes in O arag were seen during autumn runoff events (0.5 unit), whereas during July and August the near 40-μtam diurnal signal in pCO 2 (Fig 2) amounted to only 0.1 unit daily changes in O arag (Fig 4). This small variability in O arag is greater than the at most 0.06 unit error in the O arag calculation estimated by propagating the RMSE value from the TA-salinity relationship with the error in the pCO 2 data (2.1 μatm; dominated by error in the LI-COR LI840A calibration functions). Short-term diurnal variability in the Alaskan hatchery is minimal in comparison to the extremes seen at the Whiskey Creak Hatchery (WCH) in Netarts Bay, Oregon where O arag can vary by more than a whole unit over a day [8]. O arag varies more on seasonal time scales at APSH, with highest values in July and August and low sub-optimal levels persisting from January to April. October through December is a transition period to the prolonged sub-optimal conditions in winter, which is punctuated by abrupt short-lived decreases in O arag . Assuming that May and June also are transition months to the high summer O arag levels, and that conditions in September are similar to those of August, this implies a 5-month window (May through September) of favorable O arag conditions at APSH. During these 5 months, prolonged O arag levels above the suboptimal 1.5 O arag threshold should alleviate this as a stressor for larval shellfish production, whereas during the other months of the year this may not be the case. Over the 10-month period, pH T ranged from 7.96 to 8.13 and averaged 8.03 ± 0.02. Average O arag was near the sub-optimal threshold at 1.55 ± 0.15, with levels varying from 1.2 to 2.0. pH T and O arag closely tracked each other for only 6 of the 10 months of observations (Fig 4). From January to June, pH T and O arag were tightly coupled, with O arag largely near or below the 1.5 O arag threshold. pH T and O arag became decoupled in July, August, October and November, with pH T depressed while O arag persisted above the sub-optimal O arag threshold, except during Relationship between total alkalinity (TA; μmol kg -1 ) and salinity in the northern Gulf of Alaska calculated using the data described in Fig 1. The linear fit (gray line) for these data is: TA = 48.7709*S + 606.23 μmol kg -1 (r 2 = 0.94, root mean square error = 17.21 μmol kg -1 ). This fit was calculated using the MathWorks MATLAB robust linear regression algorithm with only salinity data < 33.6 (blue dots). Measurements above this salinity range are deep samples collected over the outer continental shelf that have a steeper TA-salinity relationship (black dots). Twenty-three validation TA measurements (S1 Table) were made and are shown here as red dots. short-lived periods of abrupt decreases in O arag . Maximal error in the pH T calculations estimated by propagating pCO 2 and TA-salinity errors was 0.01 units, well below the 0.07 unit drop in pH T during the pH T − O arag . The pH T and O arag decoupling was opposite in direction to what was observed in surface water impacted by glacial melt within nearby Prince William Sound [46], and highlights the absolute requirement to track two concurrent carbonate system parameters in order to constrain the carbonate system in coastal settings [8,25,47]. In this case, the decoupling was entirely temperature driven. As eloquently described by Takahashi et al. [48], the equilibrium constants for the dissociation of carbonic acid, and the K sp-arag , are temperature and salinity dependent, with the second equilibrium constant for the dissociation of carbonic acid increasing 1.5x faster with increasing temperature than the first. This inequality in the temperature response between the two equilibrium constants drives an adjustment within the carbonate system to increase carbonic acid and CO 3 2while decreasing bicarbonate ion in order to maintain charge and mass balances. Carbonic acid equals the product of CO 2 solubility and pCO 2 , and the required increase in pCO 2 drives an increase in H + lowering pH T .
In addition to the increase in CO 3 2caused by the temperature-driven shift within the carbonate system, K sp-arag decreases with increasing temperature and the sum of these two effects raises O arag levels. For the months of July, August, October and November at APSH, utilizing pH alone as an indication of O arag would produce misleading and inaccurate results. Using output Ca 2+ , CO 3 2and K sp-arag data from CO2SYS, the O arag record was decomposed into its constituents in order to pinpoint the drivers of sub-optimal O arag levels that occurred during months outside of the window of favorable O arag conditions (Fig 4). Ca 2+ in the seawater supply to the hatchery ranged between 7.5 and 9.5 mmol kg -1 , dropping to lowest values only during brief periods when low salinity surface water was mixed vertically to the depth of the intake in October and November (Figs 2 and 4). CO 3 2levels vary over a narrower range, between 0.08 and 0.13 mmol kg -1 , and were highest during the period of sustained low pCO 2 and warm water temperatures during July and August (Figs 2 and 4). The dynamic range of CO 3 2-, 51 μmol kg -1 , was 13x the error of 3.9 μmol kg -1 calculating by propagating pCO 2 and TA-salinity errors. During winter months, Ca 2+ and K sp-arag were largely stable and near 9 mmol kg -1 and 7 (mol kg -1 ) 2 x 10 −7 , respectively. The decline in O arag over this period resulted from the steady decrease in CO 3 2to a minimum in February and March (Fig 4), consistent with the period of sustained pCO 2 above atmospheric levels (Fig 2) and indicating a respiration-driven depression in O arag . During October and November, CO 3 2decreases abruptly ( Fig   4) along with short-lived increases in pCO 2 above atmospheric levels (Fig 2). During this time, CO 3

2-
, Ca 2+ and K sp-arag all decrease precipitously during the episodic low salinity events (Figs 2 and 4). Low O arag conditions during these autumn months are triggered by runoff events and match the intensity but not the length of exposure of the respiration-driven decline in O arag seen in winter (Fig 4). In temperature-salinity space, the two differing water masses responsibly for these trends in O arag become readily apparent, and this provides a framework for avoidance of sub-optimal O arag levels (Fig 5). Temperature-salinity data collected over the 10-month period show a two-pronged distribution, with sub-optimal O arag levels seen in water masses with the lowest salinities (<29) and coldest temperatures (< 7°C). Both of these water mass types have high pCO 2 with respect to the atmosphere (Fig 2), and these temperature-salinity-pCO 2 characteristics set up a diagnostic for APSH to respond to adverse O arag conditions in the absence of real-time O arag calculations.
The 5-month window of favorable O arag conditions at APSH will shrink in the future under the continued increase in anthropogenic CO 2 emissions. Using the approach of Harris et al. [17], we can subtract the anthropogenic CO 2 signal from the APSH pCO 2 data by making the assumptions that the sea-air CO 2 disequilibrium is roughly constant as are the processes shaping the distributions of TA, salinity and temperature. In this same way, we can also add additional anthropogenic CO 2 to pinpoint when the window of favorable growth conditions may close. Target atmospheric pCO 2 levels were added to APSH ΔpCO 2 data computed using the mean atmospheric pCO 2 from GAKOA (Figs 1 and 2) with the result being an adjusted seawater pCO 2 record. O arag at the target atmospheric pCO 2 levels were calculated with CO2SYS using the adjusted seawater pCO 2 , TA, salinity and temperature data. Under a pre-industrial atmospheric pCO 2 level (pCO 2 = 280 μatm), the window for favorable O arag conditions is open the entire year (Fig 6); none of the estimated O arag levels fell below the 1.5 O arag threshold. Note that the average difference between total CO 2 calculated using both unadjusted seawater pCO 2 and adjusted seawater pCO 2 to a pre-industrial atmosphere with TA, salinity and temperature data was 52.4 μmol kg -1 . This was in close agreement with an independent estimate of the total CO 2 addition from anthropogenic CO 2 uptake on density surfaces < 26.9 kg m -3 in the North Pacific until 2006 and scaled to 2014 using a growth rate of 1.62 μmol kg -1 yr -1 (47.96 μmol kg -1 ) [49]. Atmospheric pCO 2 concentrations were then iteratively increased beyond the mean 2014 atmospheric level to estimate when the window of favorable O arag conditions may close. The atmospheric pCO 2 at which this occurred was 500 μatm (Fig 6). Under the Intergovernmental Panel on Climate Change (IPCC) representative concentration pathway (RCP) emissions scenario 8.5, this atmospheric pCO 2 threshold will be reached by approximately 2040 [50]. O arag levels at APSH in 2040 would still remain above the thermodynamic threshold for aragonite dissolution (O arag = 1; Fig 6), but clearly seawater manipulation strategies will need to be in place by this time in order to facilitate conditions favorable for larval shellfish production.

Conclusions
The window of favorable O arag conditions in source waters to APSH is gradually closing, and this biogeochemical shift is consistent across the coastal Gulf of Alaska [10,15,51]. In 3 decades, sub-optimal O arag levels will be a constant condition for this premier Alaskan hatchery (Fig 6). Manipulation strategies will need to be invoked to alter the carbonate chemistry of the seawater supply, or alternatively draw seawater from a shallower depth at the risk of experiencing a higher frequency of low O arag runoff events. Inherently this implies a decreased volume in the water column of Resurrection Bay containing O arag above the sub-optimal threshold, and Temperature (°C)-salinity diagram with contours of seawater potential density anomaly (σ t ) and Ω arag as colored dots. Note the two areas of sub-optimal Ω arag (Ω arag < 1.5; warm colors) at the lowest salinities (<29) and coldest temperatures (< 7°C). doi:10.1371/journal.pone.0130384.g005 Tracking Ocean Acidification in an Alaskan Shellfish Hatchery this has implications for natural populations of vulnerable shellfish species in the region. A limited number of studies examining biological impacts due to shifting carbonate chemistry in Alaskan coastal waters exist [28], and additional studies are essential for detailing the severity of these changes for natural and cultured assemblages. Most of the species currently in production at APSH remain untested in their response to OA. The data presented here establish the current conditions experienced by APSH, and provide a framework for hatchery-based measurements in Alaska. For Alaska's growing shellfish aquaculture industry to reach its target growth and be successful in the shifting biogeochemical climate of the coastal ocean, the implementation of robust measures for tracking O arag are key, as are strong partnerships between stakeholders and scientists. State and federal government provision of OA-related scientific research that directly supports stakeholders such as the Alaskan shellfish industry is an excellent model for simultaneously backing industry while progressing scientific initiatives.
Supporting Information S1 Table. Validation measurements (n = 23) of temperature, salinity and total alkalinity collected in the hatchery during the 10-month period of observation. (EPS) Top and bottom panels are Ω arag computed using total CO 2 (TCO 2 ) adjusted to atmospheric pCO 2 levels of 280 (blue) and 500 (red) μatm, respectively, following the approach of Harris et al. [17] by assuming sea-air CO 2 disequilibria and the processes that determine TA, temperature and salinity variability are constant in time. The vertical dashed black line in all panels is the sub-optimal 1.5 Ω arag threshold where some early life stages of marine bivalves become stressed [7,8]. An atmospheric pCO 2 of 500 is expected by 2040 if the IPCC AR5 RCP 8.5 emissions trajectory is realized. doi:10.1371/journal.pone.0130384.g006