We show here that increased variability of temperature and pH synergistically negatively affects the energetics of intertidal zone crabs. Under future climate scenarios, coastal ecosystems are projected to have increased extremes of low tide-associated thermal stress and ocean acidification-associated low pH, the individual or interactive effects of which have yet to be determined. To characterize energetic consequences of exposure to increased variability of pH and temperature, we exposed porcelain crabs, Petrolisthes cinctipes, to conditions that simulated current and future intertidal zone thermal and pH environments. During the daily low tide, specimens were exposed to no, moderate or extreme heating, and during the daily high tide experienced no, moderate or extreme acidification. Respiration rate and cardiac thermal limits were assessed following 2.5 weeks of acclimation. Thermal variation had a larger overall effect than pH variation, though there was an interactive effect between the two environmental drivers. Under the most extreme temperature and pH combination, respiration rate decreased while heat tolerance increased, indicating a smaller overall aerobic energy budget (i.e. a reduced O2 consumption rate) of which a larger portion is devoted to basal maintenance (i.e. greater thermal tolerance indicating induction of the cellular stress response). These results suggest the potential for negative long-term ecological consequences for intertidal ectotherms exposed to increased extremes in pH and temperature due to reduced energy for behavior and reproduction.

Increased atmospheric CO2 from anthropogenic activity is warming air and ocean temperatures (climate change, ocean warming), while also decreasing ocean pH and carbonate ion saturation (ocean acidification) (Solomon et al., 2007). Increased air temperatures from climate change are projected to include increased frequency, intensity and duration of heat waves (Solomon et al., 2007). Ocean warming will have additional global consequences including increased energy of oceanic storms, reduced ecosystem services from coastal marine habitats such as coral reefs, and shifts in species distributions and abundance (Solomon et al., 2007). Ocean acidification may also reduce calcification in planktonic organisms (Bednaršek et al., 2014; Orr et al., 2005) though the greatest concern could be for coastal benthic calcifying organisms living in the intertidal zone (Feely et al., 2008; Kelly et al., 2013; Pespeni et al., 2013). Intertidal zone habitats already exhibit considerable variability in pH (Hofmann et al., 2011) and temperature (Stillman and Tagmount, 2009), but ocean acidification, ocean warming, climate change and hypoxia are likely to result in greater intertidal zone extremes in temperature and pH in the future (Doney et al., 2012; Melzner et al., 2013; Shaw et al., 2013; Solomon et al., 2007), which could pose a significant physiological stress on intertidal zone organisms (Waldbusser and Salisbury, 2014).

Physiological stress associated with increased environmental variability is likely to have repercussions for population distribution, growth rates, extinction risk and fitness (Boyce et al., 2006). There is evidence that pH and thermal variability can have significant effects on physiological and biochemical characteristics of marine invertebrates and fish (Alenius and Munguia, 2012; Feldmeth et al., 1974; Lowe and Heath, 1969; Otto, 1973; Podrabsky and Somero, 2004). Additionally, exposure to thermal variability can cause shifts in gene expression patterns that set limits for physiological function (Podrabsky and Somero, 2004), and exposure to increased pH variability can negatively affect performance (Alenius and Munguia, 2012). Because coastal environments are predicted to experience increases in the frequency and the magnitude of temperature and pH extremes as a result of increased atmospheric CO2 (Doney et al., 2012; Shaw et al., 2013), there is a pressing need to determine how organisms within these ecosystems will respond to such variability.

Inferences regarding the impacts of ocean acidification and ocean warming on coastal marine organisms have chiefly been made from studies that expose animals to constant levels of temperature or pH (Daufresne et al., 2009; Whiteley, 2011) or the constant combination of the two environmental drivers (Arnberg et al., 2012; Matoo et al., 2013). However, during low tide, air temperatures in the intertidal zone fluctuate dramatically over short periods of time (e.g. 20°C in 6 h) (Stillman and Somero, 1996) and reach extremes when low tides coincide with hot days (Finke et al., 2007; Helmuth, 2002; Stillman and Somero, 2000). Marine intertidal organisms tolerate these natural temperature fluctuations, yet it is known that they are limited in their scope to tolerate future warming (Stillman, 2003). Additionally, pH routinely fluctuates by ≥0.5 pH units between day and night, and by ≥1 pH unit seasonally on temperate rocky shores (Wootton et al., 2008), pH on coral reefs varies by up to 0.5 pH units during day–night cycles (Birkeland et al., 2008), and pH changes rapidly by up to 0.5 pH units during upwelling of CO2-rich water on the California coast (Feely et al., 2008; Hofmann et al., 2011). While many important inferences regarding organismal responses to ocean warming and ocean acidification have been made in studies that use constant pH and temperature, those inferences may not be ecologically relevant for intertidal zone organisms, for whom daily and tidal variability in temperature and pH are the norm (Waldbusser and Salisbury, 2014). Such studies must be refined to better reflect both the natural variation in these environments and the predicted changes in near-shore habitats.

Because temperature is arguably the most important factor influencing the physiology and distribution of ectothermic organisms (Angilletta, 2009; Cossins and Bowler, 1987) and as crustaceans are efficient acid–base regulators (Melzner et al., 2009), we hypothesized that exposure to daily variation in temperature would elicit a larger physiological response than exposure to daily variation in pH, but that there would be an interactive effect of the two environmental drivers because of the effect of temperature on acid–base regulation processes (e.g. ion pumping). Levels of experimental variability were determined by present-day (Hofmann et al., 2011; Stillman and Tagmount, 2009) and future temperature and pH scenarios for coastal regions following models of the Intergovernmental Panel on Climate Change (IPCC) and coastal oceanographic studies (Feely et al., 2008; Waldbusser and Salisbury, 2014).

The porcelain crab Petrolisthes cinctipes (Randall) was used as a model intertidal organism to identify the physiological consequences of increased environmental variability on an organism that already experiences considerable variability in nature. We used respiration rate as a proxy for the total energy available for biological processes, and changes in upper thermal tolerance to identify shifts in the cellular stress response (CSR), a critical component of basal maintenance (Sokolova, 2013). We observed a reduced metabolic rate accompanied by increased thermal tolerance under future temperature and pH variability, indicating a potential decline in ecological performance of these crabs due to energetic constraints.

Exposure to the combination of increasing extremes of pH and temperature reduced respiration rates by up to 25% (ANOVA, F2,57=3.17, P<0.05, Fig. 1A). Exposure to increasing levels of pH variability in isolation did not affect respiration rates (ANOVA, F2,55=0.08, P=0.96, Fig. 1A), while exposure to increasing temperature extremes in isolation reduced respiration rates by up to 17%, though that reduction was not statistically significant (ANOVA, F2,58=1.69, P=0.17, Fig. 1A).

Exposure to the combination of the greatest extremes of pH and temperature variation had the strongest effect on thermal tolerance limits, increasing the critical thermal maximum (CTmax) by nearly 2°C above controls (ANOVA, F2,56=10.9, P<0.001, Fig. 1B). Thermal variability in isolation increased CTmax (ANOVA, F2,55=3.16, P<0.05, Fig. 1B,) by up to ~1°C at 25°C and no further increases in CTmax were observed at 30°C exposure (Fig. 1B). Increasing pH variability in isolation also had a positive effect on CTmax (ANOVA, F2,59=3.15, P<0.05, Fig. 1B) by up to 1°C between pH 8.1 and 7.15 (Fig. 1B).

The striking, interactive physiological effect of temperature and pH variability is apparent in the plots of respiration rate versus thermal tolerance (Fig. 2). Increased pH variability in isolation had a positive effect on the relationship between these two response variables, whereas increased temperature variability in isolation had the opposite effect (Fig. 2). In combination, these effects did not offset, but instead, the combined effect of pH and temperature variability was distinct and greater at both levels of variation (present day and projected future) (Fig. 2).

We observed a synergistic interaction of pH and temperature variability on decreasing energy production (respiration rate) and an additive relationship between the two on increasing basal maintenance (i.e. survival, via the CSR). Reduced constant pH alone has been shown to reduce the metabolism of embryonic and larval porcelain crabs (Carter et al., 2013; Ceballos-Osuna et al., 2013) and larval urchins (Dorey et al., 2012), resulting in reduced scope for growth. Yet there are complexities in metabolic responses whereby not all individuals reduce metabolism to the same degree, if at all (Carter et al., 2013; Ceballos-Osuna et al., 2013), suggesting the potential for adaptive responses to ocean acidification exists. Meta-analyses across many taxa indicate reduced calcification and survival in response to ocean acidification (Kroeker et al., 2013). In combination with increased constant temperature, reduced pH has been shown to increase the metabolism of coastal marine organisms such as oysters (Matoo et al., 2013), clams (Matoo et al., 2013), shrimp (Arnberg et al., 2012) and coral (Cumbo et al., 2013). Our results indicate that the addition of environmentally relevant variability of pH and temperature could alter the trajectory of metabolic responses from elevation to metabolic depression. Metabolic depression may be an energy-saving mechanism for these crabs (Sokolova, 2013) or a method to reduce oxidative damage (Teranishi and Stillman, 2007) during times of extreme environmental variability. Reduced aerobic metabolism can also be an indicator of reduced ecological performance (Sokolova et al., 2012).

Decapod crustaceans possess a well-developed ability to actively reduce hemolymph pH during acidification stress (Melzner et al., 2009). In addition to active acid–base regulation, some species of porcelain crabs in the genus Petrolisthes are capable of passively reducing hemolymph pH during times of aerial stress by means of leg membranes that aid in passive respiration (Stillman and Somero, 1996; Stillman and Somero, 2000) and acid–base regulation (Vargas et al., 2010). pH extremes in isolation did not alter metabolic rates of adult P. cinctipes in the present study, while early life history stages of P. cinctipes (lacking the leg membrane structures) showed negative effects on energetics in response to lowered pH (Carter et al., 2013; Ceballos-Osuna et al., 2013), indicating that leg membrane structures may represent an exaptation for withstanding pH stress in the intertidal zone, providing an energetically inexpensive mechanism for acid–base regulation.

Multi-stressor studies of pH and temperature have illustrated the complex interactions between these variables, which may depend on species, life history stage and geographic range studied (Harvey et al., 2013). In corals, the combined effects of decreased pH and increased warming (increased overall stress) increased rates of microbioerosion of skeletal structures (Reyes-Nivia et al., 2013). Intertidal zone organisms may be able to withstand increased warming and acidification, but with performance effects that could lead to reduced ecological resilience. The exposure to the combination of constant low pH and constant elevated temperature has been shown to reduce growth in intertidal snails (Melatunan et al., 2013) and activity in shrimp (Dissanayake and Ishimatsu, 2011). These kinds of organismal tradeoffs for stress tolerance would be expected to occur for P. cinctipes if exposure to increased stress reduced the energy available for key processes such as growth or activity.

Mechanisms that underlie metabolic suppression during exposure to increased environmental variability are unknown, but may involve adjustments in cellular machinery that change mitochondrial abundance or capacity, enzyme activity, ion and gas transport, and membrane composition, as these are central to generalized responses to altered environments (Hochachka and Somero, 2002; Sokolova, 2013). The CSR is a generalized response to cellular perturbation, whereby separate stressors induce a common stress response pathway involving protein stabilization, DNA repair and free radical scavenging (Kültz, 2005). Petrolisthes cinctipes is known to respond to heat stress by upregulation of molecular chaperone genes (i.e. CSR) and downregulation of mitochondrial respiration genes (Teranishi and Stillman, 2007), and that response is known to shift with habitat temperature variability (Stillman and Tagmount, 2009). Following the initial stressor (e.g. thermal stress), the products of stress response pathways may remain, protecting against a subsequent, yet different, stressor (e.g. pH stress), resulting in cross-tolerance (Kültz, 2005). Acidification has been shown to induce CSR heat shock proteins hsp70 and hsp90 in gill tissue of zebrafish, demonstrating that cellular responses to acidification share signaling pathways and cross-talk with cellular responses to thermal stress (Sinclair et al., 2013; Tiedke et al., 2013). Even without clear mechanisms, the observed decreases in metabolism combined with increased upper temperature tolerance (i.e. basal maintenance/survival) suggest that limited energy may remain for discretionary metabolic demands [growth, reproduction and/or behavior (e.g. Sokolova, 2013)] and may result in negative ecological consequences due to increased variability in temperature and pH.

Our current findings suggest that acclimatization capacity may be beneficial in the short term (increased tolerance for extreme temperatures), yet detrimental in the long term (reduced overall energy). The potential for evolutionary adaptive responses to interacting stressors is unknown, but will depend on evolutionary capacity of the pathways and tolerance mechanisms involved (Sunday et al., 2014). Assessing the standing genetic diversity of genes involved in those pathways across populations from distinct thermal habitats [e.g. hotspots (Helmuth et al., 2011)] or pH gradients [e.g. upwelling centers (Feely et al., 2008)] could reveal the evolutionary potential that these crabs have to adapt to future change. The role of standing genetic diversity and local adaptation in responding to a single stressor, ocean acidification, has been demonstrated in the purple sea urchin, Strongylocentrotus purpuratus (Pespeni et al., 2013), red abalone, Haliotis rufescens (De Wit and Palumbi, 2013), and the intertidal zone carnivorous snail Concholepas concholepas (Lardies et al., 2014). Here, we highlight the need to understand the potential to adapt to multiple interacting stressors.

In conclusion, by highlighting the importance of environmentally realistic experiments to characterize organismal responses to ocean acidification and ocean warming, we have shown that the coastal impacts of ocean acidification and ocean warming will likely result in a physiological response that will reduce overall energy budgets and have ecological consequences as a result of reductions in energy available for behavior, growth and reproduction.

Porcelain crab collection and maintenance

Male, adult P. cinctipes were collected from Fort Ross, CA, USA (38.5143°N, 123.2438°W) during low tides from September 2012 to January 2013 for a total of five collection dates. Specimens were immediately transported to the Romberg Tiburon Center for Environmental Studies (Tiburon, CA, USA; transport time ~90 min), where they were placed in recirculating, 20 μm filtered, UV-sterilized seawater (12±0.5°C, salinity 33.3±0.2, ambient pH 8.09±0.08) for 24 h. Animals were not fed during this 24 h period. On the second day, specimens were acclimated in one of nine combinations (see below) of present day and future variability in pH (CO2 acidified water) and low-tide air temperature.

Seawater was obtained during high tide from a 15 m deep pipe in San Francisco Bay, settled in a stagnant tank for ≥1 week to remove sediment, filtered to 20 μm, treated with activated carbon, UV sterilized, and foam fractionated with a protein skimmer before being added to the recirculating system. In the recirculating system, water flowed continuously past a UV sterilizer before being delivered to any tank, and outflow from any tank was pumped past 100 μm and 20 μm mechanical filters, activated carbon, and a UV sterilizer before being returned to the reservoir. The total volume of the recirculating system was ~10,000 l and ~4000 l of new seawater was exchanged before each acclimation experiment.

Experimental aquaria

Acclimation experiments were performed in acrylic cylinders (one crab per cylinder, 7×15 cm diameter×height) that had 5.3 cm of reticulated foam filter (Aquatic Ecosystems, PF7) as a base to prevent upwelling of water into the cylinder (Fig. 3). Cylinders were placed in insulated and lidded acrylic tanks (Aqua Logic, USA), where they were arranged in a 5×8 cylinder grid. Each tank had a standpipe to set an upper water height that submerged the crabs by 5 cm, as well as a lower standpipe to set a lower water level that exposed each cylinder to air, but left enough water to regulate the temperature of the air space in the tanks (see below). Each tank received a continuous flow of recirculated seawater. While crabs were immersed, the water in each tank was maintained at 11±0.5°C. All animals were fed every other day with equal amounts of Shellfish Diet© (Reed Mariculture Inc., www.reedmariculture.com) containing a mixture of 30% Isochrysis sp., 20% Pavlov sp., 20% Thalossiosira weissflogii and 30% Tetraselmis sp.

Daytime emersion and temperature ramping

Petrolisthes cinctipes of Northern California (Fort Ross) usually experience a diurnal tidal cycle with low tides coinciding with dawn and dusk, while P. cinctipes that are ~450+ miles north (Cape Arago, OR, USA) usually experience low tides during peak sunlight hours (see Stillman and Tagmount, 2009). Crabs perch in moist areas under rocks during low tides, usually gregariously. We chose to expose the specimens in our study to air during a ‘low tide’ at 12:00 h each day (Fig. 3). Solenoid valves (Hayward PVC/CPVC) controlled by a timer opened to enable tanks to drain over 25 min, so that by 12:30 h each day the crabs were exposed to air. During low tide, air temperatures were regulated within each sealed tank by vigorously bubbling the water below the cylinders to thermally equilibrate the air and water (Fig. 3).

To increase air temperature during low tide, immersion heaters in the water, regulated by temperature controllers whose temperature sensors were in the air, heated the air at 2.6 or 3.7°C h−1 for 5 h (maximum air temperature of 25 and 30°C, respectively, Fig. 4B). Water was chilled back to 11±0.5°C before the crabs were resubmerged. At 18:30 h, the solenoid valves closed to fill the tanks, submerging the crabs and ending the low tide period (Fig. 4B). Crabs in the tank that had no heating experienced 11.0±0.5°C air temperatures during aerial exposure (Fig. 4B).

Night-time pH variability

To alter the pH each crab experienced, water was pumped from reservoirs containing ambient or CO2-acidified water. Water was delivered through vinyl tubing with 0.94 l h−1 irrigation drippers (Netafim 71) that were threaded through the lids of the tanks and suspended above each cylinder (Fig. 4A). The irrigation drippers ensured equal flow rates to each cylinder. Pumps delivered water for 7 h (03:00 h to 10:00 h) from the ambient pH water tank (~pH 8.09±0.08), or for 5 h (04:00 h to 09:00 h) from the low pH water tank (pH 7.59±0.08 or 7.14±0.07, depending on experimental treatment, Fig. 4A). Water in all cylinders was automatically re-equilibrated to ambient conditions through air exposure and through vigorous aeration of the circulating water by air pumps and protein skimmers by 10:30 h; thus, all crabs experienced the same pH and temperature conditions prior to the period of emersion (Fig. 4). pH was monitored in randomly chosen cylinders at 1 min intervals throughout each experiment using an Omega pH electrode attached to an Omega pH/H-SD1 pH data logger. The electrode was moved to a different cylinder each day. The probe was calibrated every 5 days with NBS buffers 4, 7 and 10 (Fisher, SB101, SB110, SB115). In addition to logging pH with the probe, pH was also determined spectrophotometrically biweekly during the experiment (see below).

Water was acidified to pH 7.6 and 7.15 using compressed CO2 gas (Praxair). A pH electrode (Fisher Scientific, Accumet pH Submersible Combination Electrode no. 13-620-Ap56) and pH controller (Professional, pH-ORP 2002C controller) activated a solenoid valve on the CO2 tank regulator (Widdicombe and Needham, 2007). The pH electrode was placed in an insulated and covered ~600 l CO2 mixing tank containing water held at 12±1°C. When the pH was above the setpoint of the pH controller, CO2 gas was introduced into the mixing tanks using a wooden air diffuser. To prevent background respiration from altering the PCO2, the water was continuously circulated through a 0.2 μm filter. Electrodes were calibrated once per week during experiments using a two-point calibration with NBS buffers 7 and 10 (Fisher, SB110, SB115). The CO2 mixing tank was additionally bubbled with ambient air for oxygenation. Each day, the mixing tank was filled with recirculated seawater between 10:00 h and 12:00 h, CO2 acidified, and pumped into the cylinders 16–18 h later (see above). Salinity did not change during the short period that the water was being acidified in the tank.

Water chemistry

Seawater samples to determine pH were collected bi-weekly from random cylinders between 08:45 h and 09:00 h when pumps were delivering acidified or ambient water to the cylinders. Seawater pH was monitored daily with electrodes and measured spectrophotometrically using m-Cresol Purple sodium salt dye (PharmaSpec UV-1700, Shimadzu, Columbia, MD, USA; Sigma-Aldrich no. 211761, St Louis, MO, USA) following the modified protocol of DOE (DOE, 1994). Seawater samples for alkalinity were collected at the same times as pH. Total alkalinity by mass (μmol kg−1 seawater) was calculated through a linear Gran procedure using a Metrohm dosimat 765 pH-Meter (Metrohm Herisau Switzerland) (Gran, 1952). Nitrogen species (nitrate, nitrite and ammonium) were analyzed by colorimetric assay using test kits (Salifert, American Pharmaceuticals), and were consistently undetectable above 0 mg l−1.

Additional parameters of the seawater carbonate chemistry including [HCO3−], [CO32−], PCO2, Ωcalcite and Ωaragonite (where Ω is the saturation state) were calculated using SeaCarb 2.4.8 (Lavigne et al., 2013) in the statistical computing program R v 3.0.1 (R Development Core Team, 2013). Seawater salinity, temperature, pH and alkalinity were used as input parameters (flag 8) using the total pH scale (pHT). Total phosphate and total silicate were assumed to be 0 mol kg−1, though measurements were not taken. Total alkalinity was high relative to open ocean water for our study in all pH treatments (3133±133 μequiv kg−1 at pH 8.1, 2702±360 μequiv kg−1 at pH 7.6 and 2873±235 μequiv kg−1 at pH 7.15, Table 1), reflecting the influence of terriginous input in our coastal/estuarine water source (Duarte et al., 2013). Calcite and aragonite saturation states did not fall below 1 at pH 8.1 and 7.6, though aragonite did fall to 0.6 at pH 7.15 (Table 1).

Oxygen consumption rate measurements

Oxygen consumption rate measurements were conducted after 17 days of exposure to the experimental treatments, using closed respirometry vials (N=17–20 crabs per treatment). Metabolic rates were measured at 18°C, the maximal average body temperature of P. cinctipes in nature (Stillman and Tagmount, 2009). Crabs were collected from the acclimation tanks at 10:00 h, following the low pH exposure, and were immediately placed in 70 ml glass vials (one crab per vial, Fisher Scientific) containing 18°C, 0.2 μm filtered seawater at ambient PCO2 (salinity 33, pH 8.05) and a 4×4 cm piece of mesh to stabilize the crabs and minimize disturbance during measurement. Vials were sealed with a plastic screw lid containing a gas-tight rubber gasket. A magnetic impeller built into the vial lid help to ensure the water in the vial was well mixed and prevented the formation of oxygen concentration gradients. Planar optode spots (diameter 0.5 cm; PreSens Precision Sensing GmbH) were glued to the inside of each vial and calibrated using O2-saturated and O2-free sterile seawater prior to each sampling experiment. A Fibox T3 oxygen meter (PreSens Precision Sensing GmbH) was used to measure the oxygen saturation in each vial using PST6v541 software (PreSens Precision Sensing GmbH). Once sealed, the vials were placed in a dark, recirculating water bath at 18°C for a 10 min recovery period. Oxygen was measured every 10 min for 40 min in a custom-milled aluminium block that maintained temperature and aligned the O2 probe and sensor spot. Oxygen saturation of the water never went below 75–100% O2 saturation. A ‘blank’ vial containing no crab was run along-side each set of animals. After respirometry, each crab was towel dried and fresh mass was measured to the nearest 0.1 g. Slopes were calculated by the decrease of oxygen in the vial per hour and normalized to fresh mass to get respiration rates in units of μmol O2 min−1 g−1.

Upper thermal tolerance limits of cardiac performance (CTmax)

After 18 days, thermal tolerance limits were indexed by cardiac CTmax (Stillman, 2003). Specimens were collected from the acclimation tanks at 10:00 h immediately following low pH exposure. A small piece of cork (1.5×1.9 cm) was glued to the carapace of each crab using cyanoacrylate glue (superglue) as a handle for immobilizing each crab. Pinholes through the carapace were made laterally on either side of the heart and 25 μm diameter impedance electrodes (ceramic coated copper wire) were inserted into each hole and glued into place with cyanoacrylate glue. Specimens were immersed in an aerated, temperature-controlled seawater bath at 12°C. Impedance electrodes were connected to a PowerLab 16sp data acquisition (ADInstruments) unit and crab heart rates were recorded continuously using the Chart v. 5 recording software (ADInstruments). Surgery never resulted in mortality and specimens were given a 30 min recovery period at a constant temperature of 12°C in which heart rates decreased to a baseline rate (see supplementary material Fig. S1) before ramping began. Temperature was increased at 0.1°C min−1 for 5 h, a rate determined to be environmentally realistic for a thermally extreme day (Stillman and Tagmount, 2009). Heart rate was monitored for a total of 5 h and was expressed as beats min−1 or transformed to the natural logarithm of beats min−1 for Arrhenius plots. An Arrhenius break temperature (ABT) represents the temperature during a thermal ramp at which heart rate dramatically changes (typically from a steady increase with temperature to an abrupt decrease). ABTs were determined in Microsoft Excel by fitting linear regressions to ascending and decreasing heart rate data (pre- and post-ABT). The temperature at which these regressions intersected was considered the ABT (Stillman, 2003). ABTs were cross-checked using the segmented package in the statistical computing program R (v 3.0.1), which generated ABTs identical to those generated in Microsoft Excel. Data shown were calculated from Microsoft Excel.

Statistical analyses

All statistical analyses were carried out using the R statistical program (v 3.0.1). Metabolic rate and CTmax data were analyzed by two-way ANOVA (with nightly pH spike and daily temperature spike as independent or combined fixed factors depending on best-fit modeling). Planned comparisons for specific pairwise analysis were done with one-way ANOVA. All residual errors were tested for normality using QQ normality plots, and data were log transformed to meet assumptions of normality or non-parametric Kruskal–Wallace analysis was applied. Tukey HSD post hoc comparisons were made to identify significant differences between treatments. Values reported in the text are means ± 1 s.e.m.

We would like to thank Michael Garland and Lars Tomanek for logistical support, Kristine Okimura for assistance with water chemistry analysis, Tessa Page, Eric Armstrong, Rachel Diner and Melanie Schiffer for assistance with specimen collection, and George Somero and Anne Todgham for their helpful comments on drafts of the manuscript.

Funding

This material is based upon work supported by the National Science Foundation [grant no. 1041225 to J.H.S.].

Alenius
B.
,
Munguia
P.
(
2012
).
Effects of pH variability on the intertidal isopod, Paradella dianae
.
Mar. Freshw. Behav. Physiol.
45
,
245
-
259
.
Angilletta
M. J.
(
2009
).
Thermal Adaptation: A Theoretical and Empirical Synthesis
.
Oxford, UK
:
Oxford University Press
.
Arnberg
M.
,
Calosi
P.
,
Spicer
J. I.
,
Tandberg
A. H. S.
,
Nilsen
M.
,
Westerlund
S.
,
Bechmann
R. K.
(
2012
).
Elevated temperature elicits greater effects than decreased pH on the development, feeding and metabolism of northern shrimp (Pandalus borealis) larvae
.
Mar. Biol.
160
,
2037
-
2048
.
Bednaršek
N.
,
Feely
R. A.
,
Reum
J. C. P.
,
Peterson
B.
,
Menkel
J.
,
Alin
S. R.
,
Hales
B.
(
2014
).
Limacina helicina shell dissolution as an indicator of declining habitat suitability owing to ocean acidification in the California Current Ecosystem
.
Proc. Biol. Sci.
281
,
20140123
.
Birkeland
C.
,
Craig
P.
,
Fenner
D.
,
Smith
L.
,
Kiene
W. E.
,
Riegl
B.
(
2008
).
Geologic setting and ecological functioning of coral reefs in American Samoa
. In
Coral Reefs of the USA
(ed.
Riegl
B.
,
Dodge
R. E.
), pp.
741
-
765
.
New York, NY
:
Springer
.
Boyce
M. S.
,
Haridas
C. V.
,
Lee
C. T.
The Nceas Stochastic Demography Working Group
(
2006
).
Demography in an increasingly variable world
.
Trends Ecol. Evol.
21
,
141
-
148
.
Carter
H. A.
,
Ceballos-Osuna
L.
,
Miller
N. A.
,
Stillman
J. H.
(
2013
).
Impact of ocean acidification on metabolism and energetics during early life stages of the intertidal porcelain crab Petrolisthes cinctipes
.
J. Exp. Biol.
216
,
1412
-
1422
.
Ceballos-Osuna
L.
,
Carter
H. A.
,
Miller
N. A.
,
Stillman
J. H.
(
2013
).
Effects of ocean acidification on early life-history stages of the intertidal porcelain crab Petrolisthes cinctipes
.
J. Exp. Biol.
216
,
1405
-
1411
.
Cossins
A.
,
Bowler
K.
(
1987
).
Temperature Biology of Animals
.
New York, NY
:
Chapman and Hall
.
Cumbo
V. R.
,
Fan
T. Y.
,
Edmunds
P. J.
(
2013
).
Effects of exposure duration on the response of Pocillopora damicornis larvae to elevated temperature and high pCO2
.
J. Exp. Mar. Biol. Ecol.
439
,
100
-
107
.
Daufresne
M.
,
Lengfellner
K.
,
Sommer
U.
(
2009
).
Global warming benefits the small in aquatic ecosystems
.
Proc. Natl. Acad. Sci. USA
106
,
12788
-
12793
.
De Wit
P.
,
Palumbi
S. R.
(
2013
).
Transcriptome-wide polymorphisms of red abalone (Haliotis rufescens) reveal patterns of gene flow and local adaptation
.
Mol. Ecol.
22
,
2884
-
2897
.
Dissanayake
A.
,
Ishimatsu
A.
(
2011
).
Synergistic effects of elevated CO2 and temperature on the metabolic scope and activity in a shallow-water coastal decapod (Metapenaeus joyneri; Crustacea: Penaeidae)
.
ICES J. Mar. Sci.
68
,
1147
-
1154
.
DOE
(
1994
).
Handbook of Methods for the Analysis of the Various Parameters of the Carbon Dioxide System in Sea Water, Version 2
(ed.
Dickson
A. G.
,
Goyet
C.
),
ORNL/CDIAC-74
.
Oak Ridge, TN
:
Carbon Dioxide Information and Analysis Center
.
Doney
S. C.
,
Ruckelshaus
M.
,
Duffy
J. E.
,
Barry
J. P.
,
Chan
F.
,
English
C. A.
,
Galindo
H. M.
,
Grebmeier
J. M.
,
Hollowed
A. B.
,
Knowlton
N.
, et al. 
. (
2012
).
Climate change impacts on marine ecosystems
.
Ann. Rev. Mar. Sci.
4
,
11
-
37
.
Dorey
N.
,
Melzner
F.
,
Martin
S.
,
Oberhänsli
F.
,
Teyssié
J.-L.
,
Bustamante
P.
,
Gattuso
J.-P.
,
Lacoue-Labarthe
T.
(
2012
).
Ocean acidification and temperature rise: effects on calcification during early development of the cuttlefish Sepia officinalis
.
Mar. Biol.
160
,
2007
-
2022
.
Duarte
C. M.
,
Hendriks
I. E.
,
Moore
T. S.
,
Olsen
Y. S.
,
Steckbauer
A.
,
Ramajo
L.
,
Carstensen
J.
,
Trotter
J. A.
,
McCulloch
M.
(
2013
).
Is ocean acidification an open-ocean syndrome? Understanding anthropogenic impacts on seawater pH
.
Estuaries Coasts
36
,
221
-
236
.
Feely
R. A.
,
Sabine
C. L.
,
Hernandez-Ayon
J. M.
,
Ianson
D.
,
Hales
B.
(
2008
).
Evidence for upwelling of corrosive “acidified” water onto the continental shelf
.
Science
320
,
1490
-
1492
.
Feldmeth
R. C.
,
Stone
E. A.
,
Brown
J. H.
(
1974
).
An increased scope for thermal tolerance upon acclimating pupfish (Cyprinodon) to cycling temperatures
.
J. Comp. Physiol.
89
,
39
-
44
.
Finke
G.
,
Navarrete
S.
,
Bozinovic
F.
(
2007
).
Tidal regimes of temperate coasts and their influences on aerial exposure for intertidal organisms
.
Mar. Ecol. Prog. Ser.
343
,
57
-
62
.
Gran
G.
(
1952
).
Determination of the equivalence point in potentiometric titrations. Part ll
.
The Analyst
77
,
661
-
671
.
Harvey
B. P.
,
Gwynn-Jones
D.
,
Moore
P. J.
(
2013
).
Meta-analysis reveals complex marine biological responses to the interactive effects of ocean acidification and warming
.
Ecol. Evol.
3
,
1016
-
1030
.
Helmuth
B.
(
2002
).
How do we measure the environment? Linking intertidal thermal physiology and ecology through biophysics
.
Integr. Comp. Biol.
42
,
837
-
845
.
Helmuth
B.
,
Yamane
L.
,
Lalwani
S.
,
Matzelle
A.
,
Tockstein
A.
,
Gao
N.
(
2011
).
Hidden signals of climate change in intertidal ecosystems: what (not) to expect when you are expecting
.
J. Exp. Mar. Biol. Ecol.
400
,
191
-
199
.
Hochachka
P. W.
,
Somero
G. N.
(
2002
).
Biochemical Adaptation: Mechanism and Process in Physiological Evolution
.
New York, NY
:
Oxford University Press
.
Hofmann
G. E.
,
Smith
J. E.
,
Johnson
K. S.
,
Send
U.
,
Levin
L. A.
,
Micheli
F.
,
Paytan
A.
,
Price
N. N.
,
Peterson
B.
,
Takeshita
Y.
, et al. 
. (
2011
).
High-frequency dynamics of ocean pH: a multi-ecosystem comparison
.
PLoS ONE
6
,
e28983
.
Kelly
M. W.
,
Padilla-Gamiño
J. L.
,
Hofmann
G. E.
(
2013
).
Natural variation and the capacity to adapt to ocean acidification in the keystone sea urchin Strongylocentrotus purpuratus
.
Glob. Chang. Biol.
19
,
2536
-
2546
.
Kroeker
K. J.
,
Kordas
R. L.
,
Crim
R.
,
Hendriks
I. E.
,
Ramajo
L.
,
Singh
G. S.
,
Duarte
C. M.
,
Gattuso
J.-P.
(
2013
).
Impacts of ocean acidification on marine organisms: quantifying sensitivities and interaction with warming
.
Glob. Chang. Biol.
19
,
1884
-
1896
.
Kültz
D.
(
2005
).
Molecular and evolutionary basis of the cellular stress response
.
Annu. Rev. Physiol.
67
,
225
-
257
.
Lardies
M. A.
,
Arias
M. B.
,
Poupin
M. J.
,
Manríquez
P. H.
,
Torres
R.
,
Vargas
C. A.
,
Navarro
J. M.
,
Lagos
N. A.
(
2014
).
Differential response to ocean acidification in physiological traits of Concholepas concholepas populations
.
J. Sea Res.
90
,
127
-
134
.
Lavigne
A. H.
,
Epitalon
J.
,
Gattuso
J.
(
2013
).
seacarb: Seawater Carbonate Chemistry with R. R package version 2.4.8
. Available at: http://cran.r-project.org/package=seacarb.
Lowe
C.
,
Heath
W.
(
1969
).
Behavioral and physiological responses to temperature in the desert pupfish Cyprinodon macularius
.
Physiol. Zool.
42
,
53
-
59
.
Matoo
O. B.
,
Ivanina
A. V.
,
Ullstad
C.
,
Beniash
E.
,
Sokolova
I. M.
(
2013
).
Interactive effects of elevated temperature and CO2 levels on metabolism and oxidative stress in two common marine bivalves (Crassostrea virginica and Mercenaria mercenaria)
.
Comp. Biochem. Physiol.
164A
,
545
-
553
.
Melatunan
S.
,
Calosi
P.
,
Rundle
S.
,
Widdicombe
S.
,
Moody
A.
(
2013
).
Effects of ocean acidification and elevated temperature on shell plasticity and its energetic basis in an intertidal gastropod
.
Mar. Ecol. Prog. Ser.
472
,
155
-
168
.
Melzner
F.
,
Gutowska
M. a.
,
Langenbuch
M.
,
Dupont
S.
,
Lucassen
M.
,
Thorndyke
M. C.
,
Bleich
M.
,
Pörtner
H.-O.
(
2009
).
Physiological basis for high CO2 tolerance in marine ectothermic animals: pre-adaptation through lifestyle and ontogeny?
Biogeosciences
6
,
2313
-
2331
.
Melzner
F.
,
Thomsen
J.
,
Koeve
W.
,
Oschlies
A.
,
Gutowska
M. a.
,
Bange
H. W.
,
Hansen
H. P.
,
Körtzinger
A.
(
2013
).
Future ocean acidification will be amplified by hypoxia in coastal habitats
.
Mar. Biol.
160
,
1875
-
1888
.
Orr
J. C.
,
Fabry
V. J.
,
Aumont
O.
,
Bopp
L.
,
Doney
S. C.
,
Feely
R. A.
,
Gnanadesikan
A.
,
Gruber
N.
,
Ishida
A.
,
Joos
F.
, et al. 
. (
2005
).
Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms
.
Nature
437
,
681
-
686
.
Otto
G.
(
1973
).
Temperature tolerance of the mosquitofish, Gambusia affinis (Baird and Girard)
.
J. Fish Biol.
5
,
575
-
585
.
Pespeni
M. H.
,
Sanford
E.
,
Gaylord
B.
,
Hill
T. M.
,
Hosfelt
J. D.
,
Jaris
H. K.
,
LaVigne
M.
,
Lenz
E. A.
,
Russell
A. D.
,
Young
M. K.
, et al. 
. (
2013
).
Evolutionary change during experimental ocean acidification
.
Proc. Natl. Acad. Sci. USA
110
,
6937
-
6942
.
Podrabsky
J. E.
,
Somero
G. N.
(
2004
).
Changes in gene expression associated with acclimation to constant temperatures and fluctuating daily temperatures in an annual killifish Austrofundulus limnaeus
.
J. Exp. Biol.
207
,
2237
-
2254
.
Reyes-Nivia
C.
,
Diaz-Pulido
G.
,
Kline
D.
,
Guldberg
O.-H.
,
Dove
S.
(
2013
).
Ocean acidification and warming scenarios increase microbioerosion of coral skeletons
.
Glob. Chang. Biol.
19
,
1919
-
1929
.
Shaw
E. C.
,
Munday
P. L.
,
McNeil
B. I.
(
2013
).
The role of CO2 variability and exposure time for biological impacts of ocean acidification
.
Geophys. Res. Lett.
40
,
4685
-
4688
.
Sinclair
B. J.
,
Ferguson
L. V.
,
Salehipour-Shirazi
G.
,
Macmillan
H.
(
2013
).
Cross-tolerance and cross-talk in the cold: relating low temperatures to desiccation and immune stress in insects
.
Integr. Comp. Biol.
53
,
545
-
556
.
Sokolova
I. M.
(
2013
).
Energy-limited tolerance to stress as a conceptual framework to integrate the effects of multiple stressors
.
Integr. Comp. Biol.
53
,
597
-
608
.
Sokolova
I. M.
,
Frederich
M.
,
Bagwe
R.
,
Lannig
G.
,
Sukhotin
A. A.
(
2012
).
Energy homeostasis as an integrative tool for assessing limits of environmental stress tolerance in aquatic invertebrates
.
Mar. Environ. Res.
79
,
1
-
15
.
Solomon
S.
,
Qin
D.
,
Manning
M.
,
Chen
Z.
,
Marquis
M.
,
Averyt
K. B.
,
Tignor
M.
,
Miller
H. L.
(
2007
).
Climate Change: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the intergovernmental Panel on Climate Change
.
Cambridge; New York, NY
:
Cambridge University Press for IPCC
.
Stillman
J. H.
(
2003
).
Acclimation capacity underlies susceptibility to climate change
.
Science
301
,
65
.
Stillman
J.
,
Somero
G.
(
1996
).
Adaptation to temperature stress and aerial exposure in congeneric species of intertidal porcelain crabs (genus Petrolisthes): correlation of physiology, biochemistry and morphology with vertical distribution
.
J. Exp. Biol.
199
,
1845
-
1855
.
Stillman
J. H.
,
Somero
G. N.
(
2000
).
A comparative analysis of the upper thermal tolerance limits of eastern Pacific porcelain crabs, genus Petrolisthes: influences of latitude, vertical zonation, acclimation, and phylogeny
.
Physiol. Biochem. Zool.
73
,
200
-
208
.
Stillman
J. H.
,
Tagmount
A.
(
2009
).
Seasonal and latitudinal acclimatization of cardiac transcriptome responses to thermal stress in porcelain crabs, Petrolisthes cinctipes
.
Mol. Ecol.
18
,
4206
-
4226
.
Sunday
J. M.
,
Calosi
P.
,
Dupont
S.
,
Munday
P. L.
,
Stillman
J. H.
,
Reusch
T. B. H.
(
2014
).
Evolution in an acidifying ocean
.
Trends Ecol. Evol.
29
,
117
-
125
.
Teranishi
K. S.
,
Stillman
J. H.
(
2007
).
A cDNA microarray analysis of the response to heat stress in hepatopancreas tissue of the porcelain crab Petrolisthes cinctipes
.
Comp. Biochem. Physiol.
2D
,
53
-
62
.
Tiedke
J.
,
Cubuk
C.
,
Burmester
T.
(
2013
).
Environmental acidification triggers oxidative stress and enhances globin expression in zebrafish gills
.
Biochem. Biophys. Res. Commun.
441
,
624
-
629
.
Vargas
M.
,
Lagos
M. E.
,
Contreras
D. A.
,
Cáceres
C. W.
(
2010
).
Área de estructuras respiratorias y su efecto en la regulación del equilibrio ácido-base en dos especies de cangrejos porcelánidos intermareales, Petrolisthes laevigatus y Petrolisthes violaceus
.
Rev. Biol. Mar. Oceanogr.
45
,
245
-
253
.
Waldbusser
G. G.
,
Salisbury
J. E.
(
2014
).
Ocean acidification in the coastal zone from an organism's perspective: multiple system parameters, frequency domains, and habitats
.
Ann. Rev. Mar. Sci.
6
,
221
-
247
.
Whiteley
N.
(
2011
).
Physiological and ecological responses of crustaceans to ocean acidification
.
Mar. Ecol. Prog. Ser.
430
,
257
-
271
.
Widdicombe
S.
,
Needham
H. R.
(
2007
).
Impact of CO2-induced seawater acidification on the burrowing activity of Nereis virens and sediment nutrient flux
.
Mar. Ecol. Prog. Ser.
341
,
111
-
122
.
Wootton
J. T.
,
Pfister
C. A.
,
Forester
J. D.
(
2008
).
Dynamic patterns and ecological impacts of declining ocean pH in a high-resolution multi-year dataset
.
Proc. Natl. Acad. Sci. USA
105
,
18848
-
18853
.

Competing interests

The authors declare no competing financial interests.

Supplementary information