# Sensitivity of sea urchin fertilization to pH varies across a natural pH mosaic

Lydia E Kapsenberg, Daniel E Okamoto, Jessica M. Dutton, Gretchen E.

Hofmann

► To cite this version:

Lydia E Kapsenberg, Daniel E Okamoto, Jessica M. Dutton, Gretchen E. Hofmann. Sensitivity of sea urchin fertilization to pH varies across a natural pH mosaic. Ecology and Evolution, 2017, 10.1002/ece3.2776. hal-01468076

HAL Id: hal-01468076

https://hal.sorbonne-universite.fr/hal-01468076

Submitted on 15 Feb 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Sensitivity of sea urchin fertilization to pH varies across a natural pH mosaic

Lydia Kapsenberg<fnref n="1" />,<fnref n="2" />,† ID | Daniel K. Okamoto<fnref n="1" />,<fnref n="3" />,† | Jessica M. Dutton<fnref n="4" /> |
Gretchen E. Hofmann<fnref n="1" />

**Correspondence**

Lydia Kapsenberg, Sorbonne Universités, Université Pierre et Marie Curie-Paris 6, CNRS-INSU, Laboratoire d'Océanographie de Villefranche, Villefranche-sur-Mer, France.
Email: lydia.kapsenberg@obs-vlfr.fr

**Funding information**

NSF, Grant/Award Number: OCE-1040960; NSF Graduate Research Fellowship; NSF Postdoctoral Research Fellowship, Grant/Award Number: OCE-1521597; NSERC Strategic, Grant/Award Number: STPGP 447247-13; Santa Barbara Coastal LTER, Grant/Award Number: OCE-1232779.

**Abstract**

In the coastal ocean, temporal fluctuations in pH vary dramatically across biogeographic ranges. How such spatial differences in pH variability regimes might shape ocean acidification resistance in marine species remains unknown. We assessed the pH sensitivity of the sea urchin *Strongylocentrotus purpuratus* in the context of ocean pH variability. Using unique male-female pairs, originating from three sites with similar mean pH but different variability and frequency of low pH (pH<sub>T</sub> ≤ 7.8) exposures, fertilization was tested across a range of pH (pH<sub>T</sub> 7.61-8.03) and sperm concentrations. High fertilization success was maintained at low pH via a slight right shift in the fertilization function across sperm concentration. This pH effect differed by site. Urchins from the site with the narrowest pH variability regime exhibited the greatest pH sensitivity. At this site, mechanistic fertilization dynamics models support a decrease in sperm-egg interaction rate with decreasing pH. The site differences in pH sensitivity build upon recent evidence of local pH adaptation in *S. purpuratus* and highlight the need to incorporate environmental variability in the study of global change biology.

**KEYWORDS**

adaptive capacity, fertilization, ocean acidification, pH variability, *Strongylocentrotus*, sea urchin

## 1 | INTRODUCTION

Species exist in spatially and temporally complex climatic environments. The role that such environmental complexity plays in shaping their sensitivity or resistance to anthropogenic climate change is not yet well understood for any biome (Boyd et al., 2016; Coble et al., 2016; Thornton, Ericksen, Herrero, & Challinor, 2014). According to the climate variability hypothesis, exposure to greater environmental variability leads to selection of wider climatic tolerance windows

(Janzen, 1967; Stevens, 1989). This hypothesis is supported for temperature in both terrestrial and marine ectotherms (Sunday, Bates, & Dulvy, 2011). The close relationship between marine species' temperature exposure and thermal tolerance is reflected by their distributions, which maximize the use of the species' thermal niches (Sunday, Bates, & Dulvy, 2012). As such, ocean warming is expected to cause latitudinal shifts in species ranges (Sunday et al., 2012). The climate variability hypothesis provides a framework by which to study the importance of variability exposure in the assessment of species adaptive capacity (sensu Dawson, Jackson, House, Prentice, & Mace, 2011) to other environmental changes.

<footnote n="1">¹Department of Ecology Evolution and Marine Biology, University of California Santa Barbara, Santa Barbara, CA, USA</footnote>

<footnote n="2">²Sorbonne Universités, Université Pierre et Marie Curie-Paris 6, CNRS-INSU, Laboratoire d'Océanographie de Villefranche, Villefranche-sur-Mer, France</footnote>

<footnote n="3">³School of Resource and Environmental Management, Simon Fraser University, Burnaby, BC, Canada</footnote>

<footnote n="4">⁴Wrigley Institute for Environmental Studies, University of Southern California, Los Angeles, CA, USA</footnote>

<footnote>†Authors contributed equally.</footnote>
Compared to warming, predictive effects for other environmental changes, such as ocean acidification (i.e., decrease in seawater pH due to anthropogenic carbon dioxide, CO₂, emissions) are trickier to determine in the absence of a strong latitudinal gradient (Kelly & Hofmann, 2013). Pending climate mitigation, ocean acidification is expected to yield a decrease in global mean seawater pH_T (pH on the total hydrogen ion scale) of 0.13–0.42 (Pörtner et al., 2014). Much like temperature, this mean change is dwarfed by natural variability, which on a global scale ranges between pH_T 7.8 and 8.4 (Rhein et al., 2013). Locally, ocean pH variability regimes arise due to geographic differences in oceanographic and biological features. Variation in these can create natural hot spots of (Hofmann et al., 2014), or refuges from (Kapsenberg & Hofmann, 2016), harmful low pH exposures (e.g., pH_T < 7.7). The potential selection pressure that pH variability envelopes impose is not well studied or understood.

Here, we investigate whether or not pH tolerance is related to local pH variability regimes (defined here as sites with high variance in pH time series observations and frequent exposure to low pH). As the pH variability range increases, the pH range within which organismal physiology must operate widens, creating the environmental regime that would select for fertilization kinetics that are resistant to low pH (Figure 1). Observing this effect in natural populations is extremely valuable as it infers transgenerational plasticity (via maternal provisioning or epigenetic modification Ross, Parker, & Byrne, 2016; Hofmann, 2017) or local adaptation and a potential means for genetic adaptation to future ocean acidification (Hofmann et al., 2014; Kelly & Hofmann, 2013; Kelly, Padilla-Gamiño, & Hofmann, 2013; Pespeni, Chan, Menge, & Palumbi, 2013). Species that can adapt to spatial environmental differences may be better equipped to adapt to temporal changes, such as ocean acidification.

We investigated fertilization success in the purple sea urchin *Strongylocentrotus purpuratus* (Pearse, 2006), across urchin groups spanning radically different pH variability regimes in the eastern boundary California Current Large Marine Ecosystem (CCLME, NE Pacific Ocean). In the CCLME, periodic upwelling brings deep, cold, low-pH (pH_T < 7.7) seawater to the coast (Feely, Sabine, Hernandez-Ayon, Ianson, & Hales, 2008; Hofmann et al., 2011). While extreme low pH events (pH_T < 7.60) are rare, these conditions are likely to become more frequent as upwelling events have increased in strength and duration (Iles et al., 2012) and as ocean acidification progresses (Gruber et al., 2012). Due to spatial variability of upwelling intensity along the CCLME coastline (Feely et al., 2008), the magnitude and frequency of low pH exposures are spatially constrained. These localizes encapsulate the biogeographic range of diverse benthic marine invertebrates.

Despite spanning a large biogeographic range from Alaska to Baja California and colonizing both subtidal and intertidal habitats, *S. purpuratus* lacks a strong genetic structure (Edmands, Moberg, & Burton, 1996; Flowers, Schroeter, & Burton, 2002; Palumbi & Wilson, 1990; Pespeni, Oliver, Manier, & Palumbi, 2010). As such, larvae exhibit little to no differences in temperature sensitivity across sites spanning a ~10°C gradient (Hammond & Hofmann, 2010). Yet, recent studies show that *S. purpuratus* is influenced by

<figure><img src="image_2.png" /><figcaption>**FIGURE 1** Conceptual diagram describing the climate variability hypothesis. In the context of ocean acidification, the environmental parameter is pH. As low pH exposures increase in frequency and the pH variability envelope widens (Site 2), the pH range within which organismal physiology must operate also widens. Through local environmental conditioning or natural selection, this may cause spatial differences in a species' tolerance window across its biogeographic range (dashed lines)</figcaption></figure>

in situ pH exposures where the resistance to ocean acidification is mediated in part by natural selection through differential pH exposures (as shown for larval growth, Kelly et al., 2013; and genomewide allelic frequencies, Pespeni, Chan, et al., 2013). *S. purpuratus* is a broadcast spawner, and thus, fertilization is exposed directly to seawater pH. Given that processes which regulate fertilization are frequently under strong selective pressures (Levitan, 2002, 2004), pH sensitivity in fertilization dynamics has the potential to introduce strong selective pressure for pH tolerance.

Previous studies of urchin fertilization in the context of ocean acidification have revealed a range of pH sensitivities. Responses include reduced fertilization success (Kurihara & Shirayama, 2004; Moulin, Catarino, Claessens, & Dubois, 2011; Reuter, Lotterhos, Crim, Thompson, & Harley, 2011) and sperm motility (Campbell, Levitan, Hosken, & Lewis, 2016; Havenhand, Buttler, Thorndyke, & Williamson, 2008; Schlegel, Havenhand, Gillings, & Williamson, 2012); however, there are some contradictory observations (Byrne, Soars, Selvakumaraswamy, Dworjanyn, & Davis, 2010; Caldwell et al., 2011). In addition, species differences (Frieder, 2014), pair-specific pH sensitivities (Schlegel et al., 2012; Sewell, Millar, Yu, Kapsenberg, & Hofmann, 2014), and sensitivity that scales with tide pool exposure (based on 2 days of observations, Moulin et al., 2011) have been reported. In addition, individuals within populations may be differentially impacted by low pH (Campbell et al., 2016). There is a need to expand fertilization experiments across species' biogeographic ranges and incorporate in situ exposures to the interpretation of results (Havenhand & Schlegel, 2009; Moulin et al., 2011). Using urchins from three sites spanning 1,500 km and different in situ pH exposures, we exposed eggs to a gradient of sperm concentrations under various pH levels to assess fertilization
in the laboratory. As *S. purpuratus* lacks a strong genetic structure across these sites, site differences in fertilization success would be attributable to local environmental effects on adult condition or gamete quality.

2 | MATERIALS AND METHODS

2.1 | Sites and field pH exposures

Three sites with distinct pH variability regimes were chosen (Figure 2a): Fogarty Creek, Oregon (FC, 44°50.200N, 124°03.517W, intertidal), Bodega Marine Reserve, California (BMR, 38°19.110N, 123°04.452W, intertidal), and Goleta Pier in the Santa Barbara Channel, California (SB, 34°24.854N, 119°49.711W, subtidal). Field pH exposures at each site were measured by autonomous, custom-built, Honeywell Durafet®-based pH sensors. These were deployed in the intertidal zone on emergent rocky benches during the upwelling season at FC and BMR, from April to October in 2011–2013. A detailed oceanographic description of pH exposures across the CCLME, including FC and BMR, is forthcoming (Chan et al., in review). For SB, SeaFET pH sensors (Martz, Connery, & Johnson, 2010) were deployed on a subtidal mooring from 2012 to 2015 (Mohawk Reef, 34°23.66N, 119°43.80W), <10 km from the SB urchin collection site (Hofmann & Washburn, 2015).

Mean pH<sub>T</sub> conditions of the time series at each site were similar (pH<sub>T</sub> 7.99–8.01), but variability and frequency of low pH exposures differed (Figure 2a). FC experiences strong, seasonal upwelling events resulting in the most frequent exposures to low pH (18.1% of pH<sub>T</sub> observations were ≤7.80). BMR experiences more intermediate upwelling strength resulting in less frequent low pH exposure and less overall pH variability (5.6% of pH<sub>T</sub> observations were ≤7.80), compared to FC.

SB experiences the weakest upwelling and narrowest pH variability regime with <1% of pH<sub>T</sub> observations ≤7.80. As the study period at FC and BMR targeted the upwelling season, frequency of low pH<sub>T</sub> (≤7.80) exposures is likely to be lower over an annual period than that reported here but still expected to be greater than at SB.

2.2 | Seawater treatments

For fertilization trials, we confined treatment conditions (three per site) to the range of pH observations at the study sites, rationalized by the assumptions that (1) present-day rare pH exposures will become more frequent in the future, and (2) species may not have been exposed to more extreme values than recently observed. Treatments for SB urchins were pH<sub>T</sub> 8.03, 7.87, and 7.76, whereas FC and BMR treatments were pH<sub>T</sub> 8.03, 7.76, and 7.61. Seawater acidification was performed by bubbling mixed dry, CO<sub>2</sub>-free air and pure CO<sub>2</sub> gas at desired pCO<sub>2</sub> levels via venturi injectors in sealed reservoir buckets filled with 0.35 µm filtered, UV-sterilized seawater (FSW) and maintained in temperature-controlled sea tables at ~15°C (Aqua Logic Inc., Digital temperature controller), following modified methods from Fangue et al. (2010).

In order to isolate pH effects, temperature was held constant. The spawning environment of *S. purpuratus* is not homogenous at 15°C and low pH often occurs at low temperatures (Hofmann et al., 2014; Reum et al., 2016). As such, the pH treatments intentionally do not fully represent future habitat exposures. Notably, marine invertebrate fertilization is generally tolerant across a wide range of temperatures (Byrne, 2011).

For each fertilization trial, treatment water was sampled from the reservoirs for salinity (YSI 3100 Conductivity Instrument) and total alkalinity (A<sub>T</sub>, open-cell titration using Mettler-Toledo T50

<figure><img src="image_3.png" /></figure>

**FIGURE 2** Map of study sites in the California Current Large Marine Ecosystem (a) and normal fertilization (NF) functions of *Strongylocentrotus purpuratus* from those sites by pH treatment (b). Sites span >10° latitude from Oregon (FC) to California (BMR, SB) and exhibit similar mean pH values but unique pH variability regimes as described by percent of pH observations ≤7.80 and coefficient of variation (CV) of pH sensor observations (FC and BMR, April-October, 2011–2013; SB, 2012–2015). For the NF functions (b), the gray band with dotted line represents the “global” normal fertilization curve (pooled site and pH treatment) at the 95% confidence interval, while points and thin lines represent mean site–pH treatment combinations. Extreme outliers under the fertilization curve originate largely from two pairs from BMR
titrator) following standard operating procedures (Dickson, Sabine, & Christian, 2007). Temperature and pH were measured immediately at the end of the fertilization trials from control vials that were handled experimentally without sperm or egg additions. pH<sub>T25°C</sub> was determined spectrophotometrically (Dickson et al., 2007) using unpurified m-cresol purple (Sigma-Aldrich®). This method potentially underestimates pH<sub>T</sub> by 0.032 ± 0.006 (Kapsenberg, Kelley, Shaw, Martz, & Hofmann, 2015). Carbonate system parameters were calculated at in situ temperatures from pH<sub>T25°C</sub>, A<sub>T</sub>, and salinity using the program CO2Calc [Version 1.0.1, 2010, U.S. Geological Survey] with dissociation constants from Mehrbach, Culberso, Hawley, and Pytkowic (1973) refit by Dickson and Millero (1987).

## 2.3 | Animal collections and fertilization trials

Adult *S. purpuratus* were collected in 2011 at a subtidal rocky reef at SB in February and from emergent rocky benches and surge channels in the intertidal at BMR and FC in March and April, respectively (Figure 2a). Sampling from tide pools was avoided as this would have introduced undocumented pH exposures that may influence fertilization (Kwiatkowski et al., 2016; Moulin et al., 2011). Collection location (subtidal vs. low intertidal) effects on fertilization have not previously been detected in an echinoderm (Bingham, Bacigalupi, & Johnson, 1997).

Urchins were maintained at the University of California Santa Barbara in flow-through sea tables at ambient temperatures (12–15°C) and fed giant kelp (*Macrocystis pyrifera*) blades ad libidum prior to spawning. In order to follow the natural spawning season at each site and availability of gravid urchins, fertilization trials had to be run sequentially by site: February–March (SB), April (BMR), and May (FC).

All fertilization trials took place in a temperature-controlled aluminum block (~15°C) fitted for 25-ml scintillation vials, for three pH treatments per unique male–female pair. Each male and female were used in only one set of experiments (i.e., males and females were not crossed with different individuals). Fertilization trials were replicated at the pair level to gain inferences on pair-specific pH sensitivity, in sacrifice of technical replicates. Fertilization trials of two pairs were run per day. A total of 12 replicate fertilization trials were conducted on urchins collected from SB and FC. Only 10 replicate fertilization trials were conducted on those collected from BMR, as more urchins failed to spawn with induction.

Individuals were induced to spawn by injecting 0.5 mol/L KCl into the oral surface, at room temperature. Females were inverted on a beaker with ~15°C FSW to collect eggs. Immediately following spawning of each female, eggs were diluted to 0.44 eggs/µl (1,000 eggs per 2.3 ml) in treatment seawater, and allowed to acclimate for 10 min before fertilization at ~15°C. Eggs were checked for quality via shape, size, and color, and sperm was checked for motility, prior to fertilization. For each male, sperm was pipetted dry directly after extrusion from the gonopores, briefly centrifuged (to obtain near-equal concentrations of concentrated sperm across males), and stored on ice while experimental vials were prepared. Just prior to use, concentrated sperm was diluted 1:10 in FSW from which three stock solutions were

made in treatment water (1:100 final dilution). Sperm concentration of each stock solution was determined by preserving a subsample in 2% formaldehyde–seawater that was later scored on a hemocytometer and extrapolated for each serial dilution in the fertilization trial (Sewell et al., 2014).

The pH treatment order for preparation and fertilization was randomized for each fertilization trial. Serial sperm dilutions of 1:10 were set up by adding 2.3 ml of the stock solution to 20.7 ml treatment water in 25-ml scintillation vials for a total of eight sperm concentrations (Levitan, Terhorst, & Fogarty, 2007; Reuter et al., 2011; Sewell et al., 2014). Vials were inverted between dilutions to ensure homogeneity. Eggs (1,000 eggs/2.3 ml) were added to the vials within 4.5 min of sperm activation in treatment water and gently inverted. Sperm-egg contact time was limited to 30 s and controlled by adding 1 ml 0.5 KCl to the vials and inverting to inhibit fertilization (Farley & Levitan, 2001). Embryos were subsequently allowed 2–4 hr for development at ambient temperatures in the sea table. The first 200 embryos encountered on a Sedgewick rafter slide were scored (two-cell to eight-cell stage). The following scoring metrics were used: normal fertilization (smooth raised fertilization membrane and equal cleavage if present), unfertilized (no fertilization membrane), tight fertilization membrane (Sewell et al., 2014; Tyler & Scheer, 1937), and abnormal development (normal fertilization envelope and unequal cleavage). At each sperm concentration, the proportion total eggs fertilized (TF), proportion of fertilized eggs that exhibit abnormal fertilization (AbnF), and proportion of total eggs with normal fertilization (NF) were calculated as follows, where the numerator is the response and denominator is the total count:

$$TF = (\#Tot.Eggs - \#Unfert.)/\#Tot.Eggs \quad (1)$$

$$AbnF = (\#Tight Fert.Membr. + \#Abnorm.Dev.)/ (\#Tot.Eggs - \#Unfert.) \quad (2)$$

$$NF = TF * (1 - AbnF) \quad (3)$$

For AbnF, tight fertilization membrane and abnormal development were grouped, as both phenomena relate to fertilization and tight fertilization membrane can lead to abnormal development (Tyler & Scheer, 1937). Polyspermy is prevalent at high sperm concentrations and was the most likely cause for abnormal early cleavage observed in our data at high sperm concentrations (Franke, Babcock, & Styan, 2002; Levitan, 2004; Levitan & Ferrell, 2006; Levitan et al., 2007; Sewell et al., 2014; Styan, 1998). Cells exhibiting NF but no cellular division were scored as NF, due to the inability to distinguish between delayed cleavage and AbnD. This caveat may result in a slight underestimate of AbnF and overestimate of NF. Likewise, the occasional observation of unfertilized eggs found at high sperm concentrations may be polyspermic due to failure to raise the fertilization envelope, which is known to occur in *Mesocentrotus franciscanus* (Levitan et al., 2007).

## 2.4 | Statistical approach

Two statistical approaches were used to evaluate how site and pH treatments impact fertilization dynamics and provide different forms
of inference about the experimental data. For statistical simplicity, we first used generalized linear mixed effects models (GLMMs) to evaluate how pH treatment and site impact the empirical relationship between sperm concentration and (1) TF and (2) AbnF, in the absence of assumptions of underlying fertilization mechanisms. This approach maintains flexibility for describing statistical relationships without addressing mechanisms. We then employed a mechanistic model of fertilization dynamics (i.e., fertilization kinetics) to generate inference regarding which mechanistic processes (rate of sperm-egg interaction vs. polyspermy block rate) may be affected by pH treatments for each site. The latter approach addresses very specific hypotheses but requires more rigid constraints and assumptions. The model in this case does not allow for more complicated pH effects, such as impacts on postfertilization development, and will only reflect results that are consistent for both TF and AbnF under a single modeling framework. For both approaches, pH treatments were analyzed categorically as pH varied slightly for each fertilization trial (Table 1).

## 2.4.1 | Linear statistical analysis (GLMMs)

For individuals collected from three sites exposed to different in situ pH variability regimes, we estimated the effect of site and pH on fertilization success, separately for TF and AbnF functions in order to isolate processes driving initial fertilization from those influencing

abnormal fertilization and development of fertilized eggs (i.e., polyspermy). To assess the influence of site and pH treatment on TF and AbnF, we estimated models that considered logit-scale probabilities as a function of sperm concentration using a third-order polynomial of $log_{10}(sperm concentration)$. We allowed this relationship to vary by site, pH treatment, and the three-way interaction (fixed effects of site, pH treatment, and sperm concentration) and also allowed the function to vary randomly by pair and by pH treatment within pair. We estimated models using GLMMs with a binomial likelihood estimated with the lme4 R package (Bates, Maechler, Bolker, & Walker, 2015). Because model residuals exhibited overdispersion (larger error variance than accounted for under a binomial distribution), we added an additional Gaussian noise term on the logit-transformed probabilities as in Okamoto (2016). For each model, we tested whether each term was statistically significant using likelihood ratio tests. Main effects were tested in the absence of all interactions containing the focal effect, while two-way interactions were tested in the absence of the three-way interaction. These contrasts were chosen because (1) any nonsignificant interactions may reduce statistical power, (2) one or more significant interactions already suggest a significant effect, and (3) we examine interaction effects in detail using effect size metrics (described below). We restricted data used in analyses to sperm concentrations where either TF or AbnF was not uniformly 0 or 100% ($1 < \text{sperm}/\mu\text{l} < 2 \times 10^5$ for TF and $1 \times 10^2 < \text{sperm}/\mu\text{l} < 2 \times 10^7$ for

<table><caption>**TABLE 1** Experimental conditions for *Strongylocentrotus purpuratus* fertilization trials by site (means ± SD, N = 12 for FC and SB, N = 10 for BMR). Categorical pH treatments were determined as the mean pH treatment across sites: pH~T~ 8.03, 7.87, 7.76, and 7.61. Ω~a~ is aragonite saturation state, A~T~ is total alkalinity</caption><thead><tr><th>Parameter</th><th>pH treatment</th><th>FC</th><th>BMR</th><th>SB</th></tr></thead><tbody><tr><td rowspan="4">Temperature (°C)</td><td>8.03</td><td>14.8 ± 0.2</td><td>14.3 ± 0.3</td><td>15 ± 0.5</td></tr><tr><td>7.87</td><td>-</td><td>-</td><td>15 ± 0.4</td></tr><tr><td>7.76</td><td>14.7 ± 0.2</td><td>14.4 ± 0.3</td><td>15 ± 0.4</td></tr><tr><td>7.61</td><td>14.8 ± 0.3</td><td>14.4 ± 0.3</td><td>-</td></tr><tr><td rowspan="4">pH<sub>T</sub></td><td>8.03</td><td>8.02 ± 0.03</td><td>8.06 ± 0.05</td><td>8.01 ± 0.01</td></tr><tr><td>7.87</td><td>-</td><td>-</td><td>7.87 ± 0.01</td></tr><tr><td>7.76</td><td>7.76 ± 0.02</td><td>7.75 ± 0.01</td><td>7.76 ± 0.02</td></tr><tr><td>7.61</td><td>7.61 ± 0.01</td><td>7.60 ± 0.01</td><td>-</td></tr><tr><td rowspan="4">pCO<sub>2</sub> (μatm)</td><td>8.03</td><td>428 ± 39</td><td>378 ± 46</td><td>435 ± 10</td></tr><tr><td>7.87</td><td>-</td><td>-</td><td>623 ± 11</td></tr><tr><td>7.76</td><td>828 ± 38</td><td>857 ± 18</td><td>832 ± 49</td></tr><tr><td>7.61</td><td>1,213 ± 42</td><td>1,233 ± 38</td><td>-</td></tr><tr><td rowspan="4">Ω<sub>a</sub></td><td>8.03</td><td>2.16 ± 0.15</td><td>2.33 ± 0.22</td><td>2.11 ± 0.04</td></tr><tr><td>7.87</td><td>-</td><td>-</td><td>1.61 ± 0.04</td></tr><tr><td>7.76</td><td>1.28 ± 0.05</td><td>1.22 ± 0.03</td><td>1.27 ± 0.06</td></tr><tr><td>7.61</td><td>0.93 ± 0.03</td><td>0.90 ± 0.03</td><td>-</td></tr><tr><td rowspan="4">A<sub>T</sub> (μmol/kg)</td><td>8.03</td><td>2,246 ± 5</td><td>2,239 ± 2</td><td>2,227 ± 3</td></tr><tr><td>7.87</td><td>-</td><td>-</td><td>2,230 ± 6</td></tr><tr><td>7.76</td><td>2,246 ± 6</td><td>2,240 ± 3</td><td>2,226 ± 6</td></tr><tr><td>7.61</td><td>2,245 ± 5</td><td>2,244 ± 3</td><td>-</td></tr><tr><td rowspan="4">Salinity</td><td>8.03</td><td>33.0 ± 0.2</td><td>33.1 ± 0.1</td><td>33.0 ± 0.0</td></tr><tr><td>7.87</td><td>-</td><td>-</td><td>33.0 ± 0.0</td></tr><tr><td>7.76</td><td>33.0 ± 0.2</td><td>33.0 ± 0.1</td><td>33.0 ± 0.0</td></tr><tr><td>7.61</td><td>33.0 ± 0.2</td><td>33.0 ± 0.1</td><td>-</td></tr></tbody></table>
AbnF) to reduce underdispersion in that region. We estimated variance explained by individual fixed effects and random effects terms as in Johnson's extension (Johnson, 2014) to the method of Nakagawa and Schielzeth (2013) implemented with the MuMln R package (Bartoń, 2016).

We used the resulting TF and AbnF models with full three-way interactions to generate effect size metrics. Specifically, we calculated the sperm concentration required to achieve (1) 50% NF ($S_{NF50}$), (2) optimal NF ($S_{OptNF}$), and (3) 25% AbnF ($S_{AbnF25}$). We generated expected NF function using TF and AbnF (Equation 3). For AbnF, 25% was chosen as it is the approximate observed extent for most treatments so as not to extrapolate beyond observation boundaries. For each site × pH treatment combination within each metric, we used the model fixed effects estimates and covariance matrices to numerically solve for the sperm concentrations required to achieve those metrics. To calculate 95% confidence intervals for estimates and all pairwise contrasts within each site for each metric (i.e., three contrasts for each site and each metric comparing each treatment), we used a fully parametric bootstrap wherein we extracted the mean and variance-covariance matrix of the fixed effects parameters and generated 1,000 random samples from the resulting multivariate normal distribution to generate the confidence intervals and p-values. We inferred homogeneous groups (within site per metric) using p-values with Bonferroni corrections for three comparisons ($\alpha = 0.05/3 = 0.017$). A global NF function was estimated by pooling site and treatments for use as a visual reference to compare fertilization curves in figures.

## 2.4.2 | Mechanistic fertilization dynamics model

We employed the fertilization dynamics model of Okamoto (2016) to evaluate whether any observed changes in the fertilization functions were consistent with (1) shifts in per capita sperm-egg interaction rates (driven by changes in sperm competency or a diversity of other potential processes) or (2) changes in polyspermy block rate (rate at which eggs become invulnerable to a second fertilizer following first fertilization). In brief, the model employs a series of differential equations describing the fertilization process and is the most recent formal extension of the model proposed by Vogel, Czihak, Chang, and Wolf (1982). Other model forms, such as that of Styan (1998) and Millar and Anderson (2003), provide similar results. The model considers the fraction of normal fertilization at time t as:

$$E_N(t) = \frac{\delta}{E_T} \int_0^t \frac{\beta E_T S_0 \gamma [e^{t(\beta E_T - \delta + r)} - 1]}{\exp \left[ t(\beta E_T + r) - \frac{\beta S_0 \gamma (1 - e^{-t(\beta E_T + r)})}{\beta E_T + r} \right] (\beta E_T - \delta + r)} dt \quad (4)$$

where $S_0$ is the initial sperm concentration, $E_T$ is the total egg concentration, $r$ is the sperm viability decay rate, and key parameters are $\beta$ (the rate of sperm-egg interactions per no., per second), $\gamma$ (the dimensionless product of the fraction of sperm that interact with the egg that are acceptable and fertilizable fraction of the egg surface), and $\delta$ (the rate at which eggs fertilized by a single sperm induce a polyspermy block per second). For full model derivation and description, see Okamoto (2016). A similar model with measured egg diameter substituted into

$\beta$ (i.e., $\beta = \text{egg size} \times \beta^*$, where $\beta^*$ is the new estimated parameter), to ensure differences were not due to differences in egg diameter among treatments, yielded no qualitative difference in results (not shown).

We estimated the model parameters in a hierarchical Bayesian framework, where we simultaneously estimated a unique mean value (i.e., $\bar{\beta}_{P,T}$, $\delta_{P,T}$, $\bar{\gamma}_{P,T}$) for each site × pH treatment combination, as well as letting those parameters vary randomly for each pair; we modeled pair level $\beta$ and $\delta$ values as a lognormal [i.e., $\beta_{i,P,T} \sim \text{lognormal}(\bar{\beta}_{P,T}, \sigma_{\beta})$ & $\delta_{i,P,T} \sim \text{lognormal}(\bar{\delta}_{P,T}, \sigma_{\delta})$] and $\gamma$, as truncated normal [$\gamma_i \sim \text{normal}(\bar{\gamma}_{P,T}, \sigma_{\gamma})$, $0.0001 < \gamma_i < 0.15$]. We utilized a beta-binomial likelihood where the response variable is number of normal fertilized eggs given the total number of eggs and a dispersion parameter ($\lambda$) controls for over or under dispersion beyond the constraints of a standard binomial. We utilized vague priors and sampled the model posterior using Hamiltonian Monte Carlo via Stan (Gelman, Lee, & Guo, 2015) using RStan (Stan Development Team, 2016). There is no analytical solution to the integral in Equation 4, so we numerically integrated the function within each sampling iteration using Gauss-Legendre quadrature rules. We used a Bayesian framework for both computational convenience and to integrate over large uncertainty in $\gamma$ and pair-level variation in parameter estimates. There was insufficient information to reliably estimate egg selectivity/sperm viability ($\gamma$), and thus, all estimates are marginalized over the uncertainty in $\gamma$ at all levels. Because we limited sperm-egg contact time to 30 s, the sperm decay rate ($r$) has no influence on parameter estimates (sperm decay over 30 s is negligible, Okamoto, 2016), so we use the value estimated in Okamoto (2016) of $r = 0.0003$ for S. purpuratus from Santa Barbara, California. See Table S1 for the full list of priors and table of model posteriors.

# 3 | RESULTS

## 3.1 | Effect of pH on fertilization

Fertilization response to pH was assessed in urchins from three sites with radically different pH variability regimes: wide with frequent low pH exposure (FC), intermediate (BMR), and narrow with rare low pH exposures (SB). Where pH effects were significant, lower pH increased concentrations of sperm required to achieve a given fertilization rate (i.e., right shifting the curve). The data were analyzed via two approaches. First, the more flexible statistical approach of using GLMMs identified site-specific differences in fertilization metrics wherein urchin fertilization sensitivity to low pH was detected for FC and SB when sperm concentrations were limiting or near optimal (i.e., not so abundant as to produce substantial abnormal fertilization or abnormal development). Under these conditions, SB urchins exhibited a pH sensitivity across a smaller pH range ($pH_T$ 7.76-8.03) compared to FC urchins ($pH_T$ 7.61-8.03). For urchins from BMR, pH effects were only detected at sperm concentrations great enough to produce abnormal fertilization and development (a trend that was observed for SB and statistically significant for FC). Second, the mechanistic models revealed that only experimental results from SB were consistent with a pH-driven shift in per capita rate of interaction among sperm and eggs (rate at which sperm collide with eggs).
## 3.2 | Linear statistical analysis (GLMMs)

Site and pH treatment both exhibited a subtle influence on TF and AbnF. Sperm concentration was by far the greatest determinant of both TF and AbnF in our experiment. In addition, the mean (fixed) effect of sperm concentration alone explained 78% and 20% of variance in TF and AbnF, respectively. The combined addition of site and pH treatment fixed effects yielded an increase in only 1.7% and 2.6% of variance explained for TF and AbnF, respectively. Site and pH treatment accounted for a small, but statistically significant (Table 2), portion of the variation in TF and AbnF functions. The mean effect of pH on TF differed by site, as indicated by a three-way sperm × site × pH treatment interaction (p < .05). In contrast, the AbnF function was impacted by sperm concentration, site × sperm interaction, and an overall pH effect. In addition, among-pair variability in the response accounted for 2.4% and 12% of error variance, with the remaining 18% and 65% of error variance due to within pair variability (i.e., unexplained noise), in TF and AbnF, respectively. Put simply, most of the variation in AbnF was unexplained, while variation in TF was largely explained by sperm concentration alone; however, site and pH still had significant effects despite explaining much smaller percentages of total variation. Figure 2b visually illustrates the aggregate variation about the mean response of fertilization to sperm concentration, while TF and AbnF functions are shown in Figures 3 and 4 along with the estimated global fertilization function (pooled site and pH treatments) for comparison (see Supporting Information for individual fertilization functions).

To quantify the interactive effect of site and pH treatment, pairwise comparisons of fertilization metrics ($S_{NF50}$, $S_{OptNF}$, and $S_{AbnF25}$) were calculated, revealing site-specific pH sensitivities in urchin fertilization (Figure 5). While estimated peak levels of NF were comparable across treatments, both FC and SB urchins exhibited a trend (significant or nonsignificant) in right shifting fertilization curves with decreasing pH (i.e., more sperm required to achieve the same levels of normal, optimal, and abnormal fertilization indicated by $S_{NF50}$, $S_{OptNF}$, and $S_{AbnF25}$, respectively). Changes in sperm concentration with fertilization metrics differed significantly for $S_{NF50}$ and $S_{OptNF}$, at both FC and SB ($\alpha$ = 0.017, following Bonferroni corrections for three comparisons, Table S2). For SB, this meant that slightly higher sperm concentrations were required to reach $S_{NF50}$ and $S_{OptNF}$ at pH$_T$ 7.76, compared to pH$_T$ 8.03. For the FC, the same effect was observed; however, the

lowest pH that differed significantly from control pH$_T$ 8.03 was pH$_T$ 7.61, and not pH$_T$ 7.76 as was the case for SB. Thus, while both SB and FC exhibited sensitivity in $S_{NF50}$ and $S_{OptNF}$ to decreasing pH, the pH range over which this sensitivity was observed was smaller at SB (pH$_T$ 8.03 vs. 7.76) compared to FC (pH$_T$ 8.03 vs. 7.61). In contrast, urchins from BMR exhibited no statistically significant pH sensitivity for either $S_{NF50}$ or $S_{OptNF}$. BMR urchins did exhibit a right shift in $S_{AbnF25}$ at pH$_T$ 7.61 compared to pH$_T$ 8.03 (Figure 5), as did FC urchins but from pH$_T$ 7.76 to pH$_T$ 7.61. Two urchin pairs from BMR exhibited unusually high percentage of AbnD (Figure 4). Excluding these two pairs from analyses neither changed observed patterns for fertilization metrics nor the conclusions (see Supporting Information).

## 3.3 | Mechanistic fertilization dynamics model

To identify the mechanistic underpinnings of the observed site-specific pH sensitivities in fertilization success, we estimated the instantaneous, per capita sperm-egg interaction rate ($\beta$) and polyspermy block rate ($\delta$) for each site × pH treatment (Figure 6). For sperm-egg interaction rate, we only detected a meaningful change across pH treatments in urchins from SB, which was the site with the narrowest pH variability regime. For SB, sperm-egg interaction rate ($\beta$) declined by an estimated 46% from pH$_T$ 8.03 to pH$_T$ 7.76 (upper and lower 95% credible set = 14%-69%, Table S1), with 98% posterior probability of a decline at this treatment. In contrast, we detected no meaningful change in $\beta$ for either FC or BMR (Figure 6, Table S1). No other site-treatment combinations had >95% probability of decline over the control (pH$_T$ 8.03) treatment. For polyspermy block rates, we detected no meaningful effect of pH treatment for any site (Figure 6). This indicates that pH sensitivity of *S. purpuratus* fertilization in this study was related to the pH sensitivity of sperm (and not eggs), which is also shown by the subtle right shift in the fertilization function across sperm concentration (i.e., more sperm required to reach fertilization metrics, Figure 5).

# 4 | DISCUSSION

We investigated the pH sensitivity of fertilization in sea urchins from different coastal ocean pH variability regimes. For *S. purpuratus* collected from three sites, pH sensitivity of fertilization was greatest in urchins from the site exposed to the narrowest pH variability regime

<table><caption>**TABLE 2** Likelihood ratio tests for generalized linear mixed models of total fertilization proportion and abnormal fertilization proportion. Main effects were tested in the absence of interactions involving the focal effect, and two-way interactions were tested in the absence of the three-way interaction (see the Methods section)</caption><thead><tr><th rowspan="2">Source</th><th colspan="3">Total fertilization</th><th colspan="3">Abnormal fertilization</th></tr><tr><th>χ<sup>2</sup></th><th>df</th><th>p</th><th>χ<sup>2</sup></th><th>df</th><th>p</th></tr></thead><tbody><tr><td>Sperm Conc.</td><td>1258.10</td><td>3</td><td>&lt;.001</td><td>490.71</td><td>1</td><td>&lt;.001</td></tr><tr><td>Site</td><td>1.45</td><td>2</td><td>.485</td><td>5.16</td><td>2</td><td>.076</td></tr><tr><td>pH treatment</td><td>5.81</td><td>3</td><td>.121</td><td>8.98</td><td>3</td><td>.030</td></tr><tr><td>Site × Sperm</td><td>27.24</td><td>6</td><td>&lt;.001</td><td>14.41</td><td>2</td><td>&lt;.001</td></tr><tr><td>pH × Sperm</td><td>8.56</td><td>9</td><td>.479</td><td>3.22</td><td>3</td><td>.360</td></tr><tr><td>Site × pH</td><td>17.27</td><td>3</td><td>&lt;.001</td><td>5.24</td><td>3</td><td>.155</td></tr><tr><td>Site × pH × Sperm</td><td>23.05</td><td>9</td><td>&lt;.01</td><td>1.08</td><td>3</td><td>.783</td></tr></tbody></table>
<figure><img src="image_4.png" /></figure>

**FIGURE 3** Proportion total fertilization over a range of sperm concentrations for *Strongylocentrotus purpuratus* from FC, BMR, and SB by pH treatment (colors are the same as in Figure 2b). The dotted line represents the global mean estimate. Solid, thick lines represent the site-pH treatment level estimates with the 95% confidence interval as the colored band. Thin, solid lines represent unique pair estimates (i.e., random effects) for each site-pH treatment. The horizontal line denotes 50% total fertilization. Only three pH treatments were tested per site

<figure><img src="image_5.png" /></figure>

**FIGURE 4** Proportion abnormal fertilization over a range of sperm concentrations for *Strongylocentrotus purpuratus* from FC, BMR, and SB by pH treatment (colors are the same as in Figure 2b). The dotted line represents the global mean estimate. Solid, thick lines represent the site-pH treatment level estimates with the 95% confidence interval as the colored band. Thin, solid lines represent unique pair estimates (i.e., random effects) for each site-pH treatment. The horizontal line denotes 25% abnormal fertilization. Only three pH treatments were tested per site
(SB). Here, the negative pH effect was observed at pH~T~ 7.76 via a right shift in the fertilization function. For urchins from sites with more frequent low pH exposures (FC, BMR), low pH tolerance extended by ≥0.11 units pH~T~. In other words, urchins from sites with frequent low pH exposure required much lower pH treatments to significantly alter fertilization rates. These results are consistent with the climate variability hypothesis and emphasize the need to include spatial and temporal environmental variability in studies of species responses to environmental change (Figure 1).

Across all sites, high mean fertilization success (89%-96%, Figure 5) was achieved at pH~T~ ≥ 7.6 as previously shown for *S. purpuratus* (Frieder, 2014) and other urchin species (Byrne, 2011; Moulin et al., 2011; Reuter et al., 2011). The effect of site and pH treatment on fertilization was small compared to that of sperm concentration or unexplained variation. The following discussion is devoted to the significance of these results in terms of assessing (1) how pH variability might shape adaptive capacity of *S. purpuratus* to ocean acidification, and (2) the ecological implications of the observed pH sensitivity, as ocean acidification might exert only a subtle influence on fertilization dynamics.

## 4.1 | Adaptive capacity

The extended low pH tolerance in fertilization dynamics of urchins from FC and BMR compared to SB could occur as a transgenerational response to in situ low pH exposures (Ross et al., 2016) or stem from natural selection (e.g., postsettlement selection, local adaptation; Kelly et al., 2013; Pespeni, Chan, et al., 2013). Our results are consistent with recent studies showing local adaptation of *S. purpuratus* to pH (Evans, Pespeni, Hofmann, Palumbi, & Sanford, 2017; Kelly et al.,

2013; Pespeni, Chan, et al., 2013; Pespeni, Sanford, et al., 2013). Kelly et al. (2013) found that *S. purpuratus* larvae grew larger at pH~T~ 7.60 when the sires originated from a site exposed to stronger upwelling events and lower pH (Van Damme State Park, ~140 km north of BMR) than when sires originated from a site with narrow pH variability (SB). The potential for local adaptation of *S. purpuratus* in response to pH is further supported by observations of changes in allelic frequencies associated with low-pH adapted alleles of adults across the CCLME (Pespeni, Chan, et al., 2013).

Our results, taken with the context above, indicate that pH tolerance in *S. purpuratus* is likely spatially optimized to current exposures because sensitivities were detected at pH levels that locally occur more rarely (Figure 1). Such spatial fine-tuning suggests that *S. purpuratus* has at least some adaptive capacity (i.e., ability to physiologically adjust via transgenerational plasticity or to evolve) to deal with ocean acidification over temporal scales. In the laboratory, other urchin species have exhibited beneficial transgenerational responses to low pH exposures, when adults were exposed to low pH over a full annual cycle of reproductive conditioning (Dupont, Dorey, Stumpp, Melzner, & Thorndyke, 2013; Suckling et al., 2015).

The cellular and physiological mechanisms that regulate the pH sensitivity of key functional traits, such as fertilization, shape the basis for assessing adaptive capacity but remain underdescribed. Using related mechanistic models, studies on other genera have inferred an impact of pH on polyspermy block rate (Sewell et al., 2014) and sperm-egg interaction rates (Reuter et al., 2011), suggesting potential species-specific sensitivities (Frieder, 2014). For *S. purpuratus*, we did not detect pH effects on block rates. For urchins from SB, however, sperm-egg interaction rates declined with pH, which could be driven by reduced sperm competency. Sperm motility is considered

<figure><img src="image_6.png" /><figcaption>**FIGURE 5** Estimated sperm concentrations required to reach 50% normal fertilization (S~NF50~), optimal normal fertilization (S~OptNF~), and 25% abnormal fertilization (S~AbnF25~) under different pH treatments (y-axis), for *Strongylocentrotus purpuratus* from FC, BMR, and SB. Points within a population-metric combination that do not share a common letter within a metric are significantly different (α = 0.017, following a Bonferroni correction for three comparisons under a parametric bootstrap). Those with no letters indicate no significant differences among treatments (within site, per metric). Percentages are the estimated peak levels of normal fertilization. Error bars are 95% confidence intervals estimated via parametric bootstrap. Colors are the same as in Figure 2b</figcaption></figure>
an important factor contributing to observations of reduced urchin fertilization success in laboratory pH experiments (Bögner, 2016; Havenhand et al., 2008). Sperm motility is activated by a pH-dependent, ATP-hydrolyzing, enzyme complex, axonemal dynein ATPase. Activation of dynein ATPase depends on H⁺ extrusion by Na⁺-H⁺ exchangers that increase intracellular pH (pHᵢ) once the sperm is released into seawater (Bögner, 2016). As such, dynein ATPase activity increases linearly with pHᵢ from 7.4 to 8.0 and is directly related to sperm motility (Trimmer & Vacquier, 1986; but see Caldwell et al., 2011). Seawater pH may influence sperm motility by modulating pHᵢ through direct intracellular acidification via CO₂ diffusion across the cell membrane (i.e., acidosis) or alteration of the effectiveness of trans-membrane proteins that control pHᵢ. Given this mechanism, adult pH exposure could enhance gamete performance via changes in gametic control of pHᵢ. For example, frequent low pH exposure of adults may impact the number or efficiency of transmembrane proteins available for pHᵢ homeostasis in sperm, as discussed by Moulin et al. (2011), potentially through epigenetic changes or selection over time.

Similarly for females, adaptive control on egg pHᵢ could alleviate intracellular acidosis that might underpin the delayed polyspermy block at low pH observed in *M. franciscanus* and *Sterechinus neumayeri* (Reuter et al., 2011; Sewell et al., 2014). We did not observe this effect in *S. purpuratus* (Figure 6). Instead, at the site where pH treatment effects were strongest, the overall sperm-egg interaction rate was impacted with no detectable impact on polyspermy block dynamics (Figure 6). This result may indicate enhanced pH sensitivity of sperm (i.e., enzymes that control sperm motility, viability, binding, and gamete recognition) or factors that contribute to interaction rates (i.e., attractive properties of eggs). Selection effects on males may therefore be the source of varying local pH sensitivities in *S. purpuratus* (Kelly et al., 2013).

The data from BMR, the site of intermediate pH variability but of urchins with greatest pH tolerance, illustrate that pH sensitivity is not

simply a function of frequency of low pH events ≤7.80. Local stress events could influence fertilization dynamics. For example, months prior to the urchin collections at BMR, there had been a rapid mortality event of intertidal *S. purpuratus* (E. Sanford, pers. comm.). This could have influenced surviving individuals differently (e.g., maturation or quality of the gonads) and the BMR urchins that did survive to spawn may have, inadvertently, been more tolerant of environmental stressors than what is representative of the BMR population in general. Urchin densities worldwide are dynamic in space and time with occasional boom and bust cycles (Filbee-Dexter & Scheibling, 2014) and it could be that such short-term perturbations have strong effects in terms of site tolerances at a given location and moment in time. Regardless of the cause, the BMR results contribute valuable information regarding the presence of pH-resistant pairs within a population.

## 4.2 | Ecological implications

From an ecological perspective, pH has a much smaller impact on fertilization dynamics than sperm availability, pair compatibility, or other unidentified processes. First, sperm availability (accounting for 78% of variation in TF) is an obvious determinant of fertilization success. In our experiment, the statistically significant differences in fertilization metrics (e.g., S<sub>OptNF</sub>) occurred within or across one order of magnitude in sperm concentration. However, urchin fertilization rates during broadcast spawning events can vary dramatically in the wild and gradients in sperm concentration can span multiple orders of magnitude (Franke et al., 2002; Levitan, 2002, 2004). Sperm concentration, sperm quality, and rates of sperm-egg encounters are dominant factors determining fertilization success (Levitan, Sewell, & Chia, 1991). In the field, these factors are influenced by urchin density of the spawning population (Gaudette, Wahle, & Himmelman, 2006; Levitan, 2002; Levitan et al., 1991; Wahle & Peckham, 1999), local sex ratios (Levitan, 2004),

<figure><img src="image_7.png" /></figure>

**FIGURE 6** Parameter estimates for the instantaneous per capita sperm-egg interaction rate (per no., per second, a-c) and polyspermy block rate (per second, d-f) by pH treatment, for *Strongylocentrotus purpuratus* from FC (a, d), BMR (b, e), and SB (c, f), using the mechanistic model of Okamoto (2016). Large symbols represent among-pair means with 95% posterior credibility intervals (error bars) and horizontal dashes are the pair-specific mean estimates. Asterisk denotes mean estimates with >95% posterior probability of a decline compared to the ambient (pH<sub>T</sub> = 8.03) treatment. Colors are the same as in Figure 2b
rates of gamete advection (Lauzon-Guay, Scheibling, & Barbeau, 2006; Levitan, 2005), and hydrodynamic mixing (Crimaldi, 2012).

Second, we observed some among-pair and substantial unexplained variation in fertilization responses. Less than 3% of variation in *S. purpuratus* fertilization curves was explained by site and pH treatment, whereas pair effects accounted for up to 12% of variation in AbnF. Pair-specific compatibility varies dramatically in *S. purpuratus* (Levitan & Stapper, 2010; Stapper, Beerli, & Levitan, 2015) as well as in *M. franciscanus* (Levitan, 2012; Levitan & Ferrell, 2006). Even so, the competitive advantage of specific males during a mass spawning event could be pH dependent (Campbell et al., 2016). Unexplained variation could stem from variation in pH dependencies of processes involved to create a fertilization event other than those measured in this study (e.g., proteins involved in gamete fusion and others discussed by Campbell et al., 2016).

Pair-dependent pH sensitivities have been observed in fertilization dynamics of other urchins, *Heliocidaris erythrogramma* (Schlegel et al., 2012) and *S. neumayeri* (Sewell et al., 2014). Our data are consistent with the hypothesis that environmental variability may act on pair-specific sensitivities, thereby increasing the pH tolerance window of local urchin aggregations (Figure 1). Maintaining diverse variation in among-pair compatibility and pH sensitivity (e.g., large, genetically diverse populations) may facilitate the adaptive response to environmental change, despite the fact that among-pair compatibility may be a stronger factor driving fertilization success than pH. However, as fertilization has direct fitness costs, it remains a critical step in population persistence. Small, but widespread shifts in fertilization function due to ocean acidification could still have an important impact, but this process will not likely be a bottleneck for marine invertebrates in the future (Byrne, 2011).

As global ocean change progresses, assessing the adaptive capacity of marine species is of increasing interest to researchers and coastal ocean management groups (Chan et al., 2016). Protecting breeding populations diverse in the pH sensitivities of their functional traits may become an important management approach, especially if such populations are sources to others (Sanford & Kelly, 2011). In addition to urchins (Ross et al., 2016), other marine species show beneficial transgenerational effects from adult pH exposures (oysters, Parker, O’connor, Raftos, Pörtner, & Ross, 2015; Parker et al., 2012; musrels, Fitzer, Cusack, Phoenix, & Kamenos, 2014; fish, Miller, Watson, Donelson, Mccormick, & Munday, 2012; Munday, 2014; Murray, Malvezzi, Gobler, & Baumann, 2014; copepods, Pedersen et al., 2014; Thor & Dupont, 2015; and corals, Putnam & Gates, 2015). Thus, identifying habitats with unique oceanographic features that select for pH-resistant traits among multiple species will be important for conservation efforts (e.g., management of other local stressors) and tracking ecological change. Specifically, it may be that such regions are “nurseries” for genetic diversity. Hot spots of low pH with wide pH variability may, however, at the same time be more sensitive to ocean change compared to sites with narrow pH variability regimes. The interplay of environmental variability and subsequent ecological interactions (Kroeker et al., 2016) presents a new research frontier for understanding the effects of environmental variability and anthropogenic activities on marine ecosystems.

## ACKNOWLEDGMENTS

We thank members of the Ocean Margins Ecosystem Group for Acidification Studies consortium (OMEGAS, a multi-institution collaboration of scientists along the U.S. West Coast) for the oceanographic data collection supporting this study, particularly Dr. Francis Chan who conducted the time series analyses for FC and BMR, and for urchin collections (Menge Lab at Oregon State University and Dr. Eric Sanford at UC Davis Bodega Marine Laboratory). Mohawk Reef pH data are courtesy of the National Science Foundation (NSF) Long Term Ecological Research program (Santa Barbara Coastal LTER, OCE-0620276, OCE-1232779). We thank Christoph Pierre for urchin collections at SB, Evan Hunter, Tullio Rossi, and Ella Bendrick-Chartier for their assistance with the fertilization experiments and sperm concentration determination, and Dr. Kathy Foltz, Dr. Mary Sewell, and two anonymous reviewers for comments on the manuscript. This work was funded by NSF award OCE-1040960 to GEH in support of OMEGAS research consortium; LK was supported by an NSF Graduate Research Fellowship. During the writing of this manuscript, LK was supported by an NSF Postdoctoral Research Fellowship (OCE-1521597), and DKO was supported by an NSERC Strategic grant (STPGP 447247-13) and the Santa Barbara Coastal LTER (OCE-1232779). Data will be made available on the Pangaea database for ocean acidification studies upon publication.

## CONFLICT OF INTEREST

None declared.

## REFERENCES

- Barton, K. (2016). MuMln: Multi-Model Inference. R package version 1.15.6. Retrieved from https://cran.r-project.org/package=MuMln

- Bates, D., Maechler, M., Bolker, B., & Walker, S. (2015). Fitting linear mixed-effects models using lme4. *Journal of Statistical Software*, *67*, 1–48.

- Bingham, B. L., Bacigalupi, M., & Johnson, L. G. (1997). Temperature adaptations of embryos from intertidal and subtidal sand dollars (*Dendraster excentricus*, Eschscholtz). *Oceanographic Literature Review*, *44*, 1313.

- Bögner, D. (2016). Life under climate change scenarios: Sea urchins' cellular mechanisms for reproductive success. *Journal of Marine Science and Engineering*, *4*, 28.

- Boyd, P. W., Cornwall, C. E., Davison, A., Doney, S. C., Fourquez, M., Hurd, C. L., ... McMinn, A. (2016). Biological responses to environmental heterogeneity under future ocean conditions. *Global Change Biology*, *22*, 2633–2650.

- Byrne, M. (2011). Impact of ocean warming and ocean acidification on marine invertebrate life history stages: Vulnerabilities and potential for persistence in a changing ocean. *Oceanography and Marine Biology: An Annual Review*, *49*, 1–42.

- Byrne, M., Soars, N., Selvakumaraswamy, P., Dworjanyn, S. A., & Davis, A. R. (2010). Sea urchin fertilization in a warm, acidified and high pCO₂ ocean across a range of sperm densities. *Marine Environmental Research*, *69*, 234–239.

- Caldwell, G. S., Fitzer, S., Gillespie, C. S., Pickavance, G., Turnbull, E., & Bentley, M. G. (2011). Ocean acidification takes sperm back in time. *Invertebrate Reproduction & Development*, *55*, 217–221.

- Campbell, A. L., Levitan, D. R., Hosken, D. J., & Lewis, C. (2016). Ocean acidification changes the male fitness landscape. *Scientific Reports*, *6*, 31250.
Chan, F., Barth, J. A., Blanchette, C. A., Byrne, R. H., Chavez, F., Cheriton, O., ... Washburn, L. (in review). Evidence for widespread progression of nearshore ocean acidification in the California Current System. *Scientific Reports*.

Chan, F., Boehm, A. B., Barth, J. A., Chornesky, E. A., Dickson, A. G., Feely, R. A., ... Whiteman, E. A. (2016). *The West Coast ocean acidification and hypoxia science panel: Major findings, recommendations, and actions*. Oakland, CA, USA: California Ocean Science Trust.

Coble, A. A., Asch, R. G., Rivero-Calle, S., Heerhartz, S. M., Holding, J. M., Kremer, C. T., ... Strock, K. E. (2016). Climate is variable, but is our science? *Limnology and Oceanography Bulletin*, 25, 71–76.

Crimaldi, J. P. (2012). The role of structured stirring and mixing on gamete dispersal and aggregation in broadcast spawning. *Journal of Experimental Biology*, 215, 1031–1039.

Dawson, T. P., Jackson, S. T., House, J. I., Prentice, I. C., & Mace, G. M. (2011). Beyond predictions: Biodiversity conservation in a changing climate. *Science*, 332, 53–58.

Dickson, A. G., & Millero, F. J. (1987). A comparison of the equilibrium constants for the dissociation of carbonic acid in seawater media. *Deep-Sea Research Part I: Oceanographic Research Papers*, 34, 1733–1743.

Dickson, A. G., Sabine, C. L., & Christian, J. R. (Eds.) (2007). Guide to Best Practices for Ocean CO₂ Measurements. PICES Special Publication 3, 191 pp.

Dupont, S., Dorey, N., Stumpp, M., Melzner, F., & Thorndyke, M. (2013). Long-term and trans-life-cycle effects of exposure to ocean acidification in the green sea urchin *Strongylocentrotus droebachiensis*. *Marine Biology*, 160, 1835–1843.

Edmands, S., Moberg, P. E., & Burton, R. S. (1996). Allozyme and mitochondrial DNA evidence of population subdivision in the purple sea urchin *Strongylocentrotus purpuratus*. *Marine Biology*, 126, 443–450.

Evans, T. G., Pespeni, M. H., Hofmann, G. E., Palumbi, S. R., & Sanford, E. (2017). Transcriptomic responses to seawater acidification among sea urchin populations inhabiting a natural pH mosaic. *Molecular Ecology*, DOI: 10.1111/mec.14038. in press.

Fangue, N. A., O'donnell, M. J., Sewell, M. A., Matson, P. G., Macpherson, A. C., & Hofmann, G. E. (2010). A laboratory-based, experimental system for the study of ocean acidification effects on marine invertebrate larvae. *Limnology and Oceanography–Methods*, 8, 441–452.

Farley, G. S., & Levitan, D. R. (2001). The role of jelly coats in sperm-egg encounters, fertilization success, and selection on egg size in broadcast spawners. *American Naturalist*, 157, 626–636.

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.

Filbee-Dexter, K., & Scheibling, R. E. (2014). Sea urchin barrens as alternative stable states of collapsed kelp ecosystems. *Marine Ecology Progress Series*, 495, 1–25.

Fitzer, S. C., Cusack, M., Phoenix, V. R., & Kamenos, N. A. (2014). Ocean acidification reduces the crystallographic control in juvenile mussel shells. *Journal of Structural Biology*, 188, 39–45.

Flowers, J. M., Schroeter, S. C., & Burton, R. S. (2002). The recruitment sweepstakes has many winners: Genetic evidence from the sea urchin *Strongylocentrotus purpuratus*. *Evolution*, 56, 1445–1453.

Franke, E. S., Babcock, R. C., & Styan, C. A. (2002). Sexual conflict and polyspermy under sperm-limited conditions: In situ evidence from field simulations with the free-spawning marine Echinoid *Evechinus chloroticus*. *The American Naturalist*, 160, 485–496.

Frieder, C. A. (2014). Present-day nearshore pH differentially depresses fertilization in congeneric sea urchins. *The Biological Bulletin*, 226, 1–7.

Gaudette, J., Wahle, R. A., & Himmelman, J. H. (2006). Spawning events in small and large populations of the green sea urchin *Strongylocentrotus droebachiensis* as recorded using fertilization assays. *Limnology and Oceanography*, 51, 1485–1496.

Gelman, A., Lee, D., & Guo, J. (2015). Stan: A probabilistic programming language for bayesian inference and optimization. *Journal of Educational and Behavioral Statistics*, 40, 530–543.

Gruber, N., Hauri, C., Lachkar, Z., Loher, D., Frölicher, T. L., & Plattner, G.-K. (2012). Rapid progression of ocean acidification in the California Current System. *Science*, 337, 220–223.

Hammond, L. M., & Hofmann, G. E. (2010). Thermal tolerance of *Strongylocentrotus purpuratus* early life history stages: Mortality, stress-induced gene expression and biogeographic patterns. *Marine Biology*, 157, 2677–2687.

Havenhand, J. N., Buttler, F. R., Thorndyke, M. C., & Williamson, J. E. (2008). Near-future levels of ocean acidification reduce fertilization success in a sea urchin. *Current Biology*, 18, R651–R652.

Havenhand, J. N., & Schlegel, P. (2009). Near-future levels of ocean acidification do not affect sperm motility and fertilization kinetics in the oyster *Crassostrea gigas*. *Biogeosciences*, 6, 3009–3015.

Hofmann, G. E. (2017). Ecological epigenetics in marine metazoans. *Frontiers in Marine Science*, 4, 4.

Hofmann, G. E., Evans, T. G., Kelly, M. W., Padilla-Gamiño, J. L., Blanchette, C. A., Washburn, L., ... Dutton, J. M. (2014). Exploring local adaptation and the ocean acidification seascape - studies in the California Current Large Marine Ecosystem. *Biogeosciences*, 11, 1053–1064.

Hofmann, G. E., Smith, J. E., Johnson, K. S., Send, U., Levin, L. A., Micheli, F., ... Martz, T. R. (2011). High-frequency dynamics of ocean pH: A multi-ecosystem comparison. *PLoS ONE*, 6, e28983.

Hofmann, G. E., & Washburn, L. (2015). SBC LTER: Ocean: Time-series: Mid-water SeaFET and CO₂ system chemistry at Mohawk Reef (MKO), ongoing since 2012-01-11. Santa Barbara Coastal LTER. doi:10.6073/pasta/826b170f29458104621aa9f0e36c8901.

Iles, A. C., Gouhier, T. C., Menge, B. A., Stewart, J. S., Haupt, A. J., & Lynch, M. C. (2012). Climate-driven trends and ecological implications of event-scale upwelling in the California Current System. *Global Change Biology*, 18, 783–796.

Janzen, D. H. (1967). Why mountain passes are higher in the tropics. *The American Naturalist*, 101, 233–249.

Johnson, P. C. D. (2014). Extension of Nakagawa & Schielzeth's R2GLMM to random slopes models. *Methods in Ecology and Evolution*, 5, 944–946.

Kapsenberg, L., & Hofmann, G. E. (2016). Ocean pH time-series and drivers of variability along the northern Channel Islands, California, USA. *Limnology and Oceanography*, 61, 953–968.

Kapsenberg, L., Kelley, A. L., Shaw, E. C., Martz, T. R., & Hofmann, G. E. (2015). Near-shore Antarctic pH variability has implications for biological adaptation to ocean acidification. *Scientific Reports*, 5, 9638.

Kelly, M. W., & Hofmann, G. E. (2013). Adaptation and the physiology of ocean acidification. *Functional Ecology*, 27, 980–990.

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*. *Global Change Biology*, 19, 2536–2546.

Kroeker, K. J., Sanford, E., Rose, J. M., Blanchette, C. A., Chan, F., Chavez, F. P., ... Washburn, L. (2016). Interacting environmental mosaics drive geographic variation in mussel performance and predation vulnerability. *Ecology Letters*, 19, 771–779.

Kurihara, H., & Shirayama, Y. (2004). Effects of increased atmospheric CO₂ on sea urchin early development. *Marine Ecology Progress Series*, 274, 161–169.

Kwiatkowski, L., Gaylord, B., Hill, T., Hosfelt, J., Kroeker, K. J., Nebuchina, Y., ... Caldeira, K. (2016). Nighttime dissolution in a temperate coastal ocean ecosystem increases under acidification. *Scientific Reports*, 6, 22984.

Lauzon-Guay, J.-S., Scheibling, R. E., & Barbeau, M. A. (2006). Movement patterns in the green sea urchin, *Strongylocentrotus droebachiensis*. *Journal of the Marine Biological Association of the United Kingdom*, 86, 167–174.

Levitan, D. R. (2002). Density-dependent selection on gamete traits in three congeneric sea urchins. *Ecology*, 83, 464–479.

Levitan, D. R. (2004). Density-dependent sexual selection in external fertilizers: Variances in male and female fertilization success along the
continuum from sperm limitation to sexual conflict in the sea urchin *Strongylocentrotus franciscanus*. *The American Naturalist*, **164**, 298-309.

Levitan, D. R. (2005). Sex-specific spawning behavior and its consequences in an external fertilizer. *The American Naturalist*, **165**, 682-694.

Levitan, D. R. (2012). Contemporary evolution of sea urchin gamete-recognition proteins: Experimental evidence of density-dependent gamete performance predicts shifts in allele frequencies over time. *Evolution*, **66**, 1722-1736.

Levitan, D. R., & Ferrell, D. L. (2006). Selection on gamete recognition proteins depends on sex, density, and genotype frequency. *Science*, **312**, 267-269.

Levitan, D. R., Sewell, M. A., & Chia, F. S. (1991). Kinetics of fertilization in the sea urchin *Strongylocentrotus franciscanus*: Interaction of gamete dilution, age, and contact time. *Biological Bulletin*, **181**, 371-378.

Levitan, D. R., & Stapper, A. P. (2010). Simultaneous positive and negative frequency-dependent selection on sperm bindin, a gamete recognition protein in the sea urchin *Strongylocentrotus purpuratus*. *Evolution*, **64**, 785-797.

Levitan, D. R., Terhorst, C. P., & Fogarty, N. D. (2007). The risk of polyspermy in three congenerous sea urchins and its implications for genetic incompatibility and reproductive isolation. *Evolution*, **61**, 2007-2014.

Martz, T. R., Connery, J. G., & Johnson, K. S. (2010). Testing the Honeywell Durafet® for seawater pH applications. *Limnology and Oceanography: Methods*, **8**, 172-184.

Mehrbach, C., Culberso, C. H., Hawley, J. E., & Pytkowic, R. M. (1973). Measurement of apparent dissociation constants of carbonic acid in seawater at atmospheric pressure. *Limnology and Oceanography*, **18**, 897-907.

Millar, R. B., & Anderson, M. J. (2003). The kinetics of monospermic and polyspermic fertilization in free-spawning marine invertebrates. *Journal of Theoretical Biology*, **224**, 79-85.

Miller, G. M., Watson, S.-A., Donelson, J. M., Mccormick, M. I., & Munday, P. L. (2012). Parental environment mediates impacts of increased carbon dioxide on a coral reef fish. *Nature Climate Change*, **2**, 858-861.

Moulin, L., Catarino, A. I., Claessens, T., & Dubois, P. (2011). Effects of seawater acidification on early development of the intertidal sea urchin *Paracentrotus lividus* (Lamarck 1816). *Marine Pollution Bulletin*, **62**, 48-54.

Munday, P. L. (2014). Transgenerational acclimation of fishes to climate change and ocean acidification. *F1000Prime Reports*, **6**, 99.

Murray, C. S., Malvezzi, A., Gobler, C. J., & Baumann, H. (2014). Offspring sensitivity to ocean acidification changes seasonally in a coastal marine fish. *Marine Ecology Progress Series*, **504**, 1-11.

Nakagawa, S., & Schielzeth, H. (2013). A general and simple method for obtaining R2 from generalized linear mixed-effects models. *Methods in Ecology and Evolution*, **4**, 133-142.

Okamoto, D. K. (2016). Competition among eggs shifts to cooperation along a sperm supply gradient in an external fertilizer. *The American Naturalist*, **187**, E129-E142.

Palumbi, S. R., & Wilson, A. C. (1990). Mitochondrial DNA diversity in the sea urchins *Strongylocentrotus purpuratus* and *S. droebachiensis*. *Evolution*, **44**, 403-415.

Parker, L. M., O'connor, W. A., Raftos, D. A., Pörtner, H.-O., & Ross, P. M. (2015). Persistence of positive carryover effects in the oyster, *Saccostrea glomerata*, following transgenerational exposure to ocean acidification. *PLoS ONE*, **10**, e0132276.

Parker, L. M., Ross, P. M., O'connor, W. A., Borysko, L., Raftos, D. A., & Pörtner, H.-O. (2012). Adult exposure influences offspring response to ocean acidification in oysters. *Global Change Biology*, **18**, 82-92.

Pearse, J. S. (2006). Ecological role of purple sea urchins. *Science*, **314**, 940-941.

Pedersen, S. A., Häkedal, O. J., Salaberria, I., Tagliati, A., Gustavson, L. M., Jenssen, B. M., ... Altin, D. (2014). Multigenerational exposure to ocean acidification during food limitation reveals consequences for copepod

scope for growth and vital rates. *Environmental Science & Technology*, **48**, 12275-12284.

Pespeni, M. H., Chan, F., Menge, B. A., & Palumbi, S. R. (2013). Signs of adaptation to local pH conditions across an environmental mosaic in the California Current Ecosystem. *Integrative and Comparative Biology*, **53**, 857-870.

Pespeni, M. H., Oliver, T. A., Manier, M. K., & Palumbi, S. R. (2010). Restriction site tiling analysis: Accurate discovery and quantitative genotyping of genome-wide polymorphisms using nucleotide arrays. *Genome Biology*, **11**, R44.

Pespeni, M. H., Sanford, E., Gaylord, B., Hill, T. M., Hosfelt, J. D., Jaris, H. K., ... Palumbi, S. R. (2013). Evolutionary change during experimental ocean acidification. *Proceedings of the National Academy of Sciences of the United States of America*, **110**, 6937-6942.

Pörtner, H.-O., Karl, D., Boyd, P. W., Cheung, W., Lluch-Cota, S. E., Nojiri, Y., ... Zavialov, P. (2014). Ocean systems. In C. B. Field, V. R. Barros, D. J. Dokken, K. J. Mach, M. D. Mastrandrea, T. E. Bilir, M. Chatterjee, K. L. Ebi, Y. O. Estrada, R. C. Genova, B. Girma, E. S. Kissel, A. N. Levy, S. Maccracken, P. R. Mastrandrea & L. L. White (Eds.), *Climate change 2014: Impacts, adaptation, and vulnerability. Part A: Global and sectoral aspects. Contribution of working group II to the fifth assessment report of the intergovernmental panel on climate change* (pp. 411-484). Cambridge, UK and New York, NY, USA: Cambridge University Press.

Putnam, H. M., & Gates, R. D. (2015). Preconditioning in the reef-building coral *Pocillopora damicornis* and the potential for trans-generational acclimatization in coral larvae under future climate change conditions. *Journal of Experimental Biology*, **218**, 2365-2372.

Reum, J. C. P., Alin, S. R., Harvey, C. J., Bednaršek, N., Evans, W., Feely, R. A., ... Sabine, C. L. (2016). Interpretation and design of ocean acidification experiments in upwelling systems in the context of carbonate chemistry co-variation with temperature and oxygen. *ICES Journal of Marine Science*, **73**, 582-595.

Reuter, K. E., Lotterhos, K. E., Crim, R. N., Thompson, C. A., & Harley, C. D. G. (2011). Elevated pCO₂ increases sperm limitation and risk of polyspermy in the red sea urchin *Strongylocentrotus franciscanus*. *Global Change Biology*, **17**, 163-171.

Rhein, M., Rintoul, S. R., Aoki, S., Campos, E., Chambers, D., Feely, R. A., ... Wang, F. (2013). Observations: Ocean. In D. Qin, G.-K. Plattner, M. Tignor, S. K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex & P. M. Midgley (Eds.), *Climate change 2013: The physical science basis. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change* (pp. 255-315). Cambridge, UK and New York, NY, USA: Cambridge University Press.

Ross, P. M., Parker, L., & Byrne, M. (2016). Transgenerational responses of molluscs and echinoderms to changing ocean conditions. *ICES Journal of Marine Science*, **73**, 537-549.

Sanford, E., & Kelly, M. W. (2011). Local adaptation in marine invertebrates. *Annual Review of Marine Science*, **3**, 509-535.

Schlegel, P., Havenhand, J. N., Gillings, M. R., & Williamson, J. E. (2012). Individual variability in reproductive success determines winners and losers under ocean acidification: A case study with sea urchins. *PLoS ONE*, **7**, e53118.

Sewell, M. A., Millar, R. B., Yu, P. C., Kapsenberg, L., & Hofmann, G. E. (2014). Ocean acidification and fertilization in the Antarctic sea urchin *Sterechinus neumayeri: The importance of polyspermy*. *Environmental Science & Technology*, **48**, 713-722.

Stan Development Team (2016). RStan: the R interface to Stan, Version 2.10.1. Retrieved from http://mc-stan.org/

Stapper, A. P., Beerli, P., & Levitan, D. R. (2015). Assortative mating drives linkage disequilibrium between sperm and egg recognition protein loci in the sea urchin *Strongylocentrotus purpuratus*. *Molecular Biology and Evolution*, **32**, 859-870.

Stevens, G. C. (1989). The latitudinal gradient in geographical range: How so many species coexist in the tropics. *The American Naturalist*, **133**, 240-256.
Styan, C. A. (1998). Polyspermy, egg size, and the fertilization kinetics of free-spawning marine invertebrates. *The American Naturalist*, **152**, 290-297.

Suckling, C. C., Clark, M. S., Richard, J., Morley, S. A., Thorne, M. A. S., Harper, E. M., & Peck, L. S. (2015). Adult acclimation to combined temperature and pH stressors significantly enhances reproductive outcomes compared to short-term exposures. *Journal of Animal Ecology*, **84**, 773-784.

Sunday, J. M., Bates, A. E., & Dulvy, N. K. (2011). Global analysis of thermal tolerance and latitude in ectotherms. *Proceedings of the Royal Society of London B: Biological Sciences*, **278**, 1823-1830.

Sunday, J. M., Bates, A. E., & Dulvy, N. K. (2012). Thermal tolerance and the global redistribution of animals. *Nature Climate Change*, **2**, 686-690.

Thor, P., & Dupont, S. (2015). Transgenerational effects alleviate severe fecundity loss during ocean acidification in a ubiquitous planktonic copepod. *Global Change Biology*, **21**, 2261-2271.

Thornton, P. K., Ericksen, P. J., Herrero, M., & Challinor, A. J. (2014). Climate variability and vulnerability to climate change: A review. *Global Change Biology*, **20**, 3313-3328.

Trimmer, J. S., & Vacquier, V. D. (1986). Activation of sea urchin gametes. *Annual Review of Cell Biology*, **2**, 1-26.

Tyler, A., & Scheer, B. (1937). Inhibition of fertilization in eggs of marine animals by means of acid. *Journal of Experimental Zoology*, **75**, 179-195.

Vogel, H., Czihak, G., Chang, P., & Wolf, W. (1982). Fertilization kinetics of sea urchin eggs. *Mathematical Biosciences*, **58**, 189-216.

Wahle, R. A., & Peckham, S. H. (1999). Density-related reproductive trade-offs in the green sea urchin, *Strongylocentrotus droebachiensis*. *Marine Biology*, **134**, 127-137.

## SUPPORTING INFORMATION

Additional Supporting Information may be found online in the supporting information tab for this article.

**How to cite this article:** Kapsenberg L, Okamoto DK, Dutton JM, Hofmann GE. Sensitivity of sea urchin fertilization to pH varies across a natural pH mosaic. *Ecol Evol*. 2017;00:1-14.
doi:10.1002/ece3.2776.