Body size has large effects on organism physiology, but these effects remain poorly understood in modular animals with complex morphologies. Using two trials of a ∼24 day experiment conducted in 2014 and 2015, we tested the hypothesis that colony size of the coral Pocillopora verrucosa affects the response of calcification, aerobic respiration and gross photosynthesis to temperature (∼26.5 and ∼29.7°C) and PCO2 (∼40 and ∼1000 µatm). Large corals calcified more than small corals, but at a slower size-specific rate; area-normalized calcification declined with size. Whole-colony and area-normalized calcification were unaffected by temperature, PCO2, or the interaction between the two. Whole-colony respiration increased with colony size, but the slopes of these relationships differed between treatments. Area-normalized gross photosynthesis declined with colony size, but whole-colony photosynthesis was unaffected by PCO2, and showed a weak response to temperature. When scaled up to predict the response of large corals, area-normalized metrics of physiological performance measured using small corals provide inaccurate estimates of the physiological performance of large colonies. Together, these results demonstrate the importance of colony size in modulating the response of branching corals to elevated temperature and high PCO2.

Organism size and temperature have profound effects on biological processes (Schmidt-Nielsen, 1984), and these effects play fundamental roles in translating changes in environmental conditions into effects on population dynamics (Nisbet et al., 2000; Brown et al., 2004). The functional relationships between organism size and biological processes, and the extent to which they are modulated by temperature, therefore are of central importance in elucidating the causes of spatio-temporal variation in ecological processes (Nisbet et al., 2000; Brown et al., 2004). While these principles are gaining recognition for their potential to advance understanding of the response of populations to global climate change and ocean acidification (Gaylord et al., 2015; Bruno et al., 2015), the potential remains largely unrealized in colonial modular animals, and corals in particular (Edmunds et al., 2014).

For tropical coral reefs, there is much to be gained by considering organism size and temperature in evaluating their response to changing environmental conditions, because the scleractinian architects of these systems vary greatly in colony size (Barnes, 1973; Pratchett et al., 2015), and the ecosystem is sensitive to climate change and ocean acidification (Hoegh-Guldberg et al., 2007). In part, this sensitivity results from the susceptibility of scleractinians to temperature (Comeau et al., 2016), and their generally negative response to ocean acidification (Chan and Connelly, 2013). The functional relationships between metabolism and temperature are relatively well known for scleractinians (Coles and Jokiel, 1977; Edmunds, 2005; Comeau et al., 2016), particularly in the context of bleaching (Fitt et al., 2000), and similar details are emerging for the effects of ocean acidification (Comeau et al., 2013; Evenhuis et al., 2015). However, there has been little consideration of whether their responses to temperature and seawater pH will be modified by colony size (but see Chan et al., 2016). Perhaps colony size has been overlooked in these studies because the colonial modular design of most scleractinians is thought to confer freedom from allometric constraints (Sebens, 1987; Hughes, 2005). If colony growth occurs while conserving polyp dimensions, it is widely assumed that metabolism will scale isometrically (i.e. proportionally) with colony size (Jackson, 1979; Sebens, 1979, 1987), rather than being constrained through allometric principles (Nisbet et al., 2000; Brown et al., 2004). One consequence of the assumption of isometric scaling is that the field of coral biology has advanced on the premise that studies of small fragments from larger coral colonies (‘nubbins’, after Birkeland, 1976) are representative of a wide variety of corals in the population, regardless of colony size (e.g. Davies, 1984; Edmunds and Davies, 1986; Anthony and Fabricius, 2000; Anthony et al., 2007).

However, as early as 1986 there was evidence that large and small coral colonies were functionally dissimilar (Jokiel and Morrissey, 1986), and throughout the 1980s and early 1990s studies of the interactions of organisms with a fluid environment (Koehl, 1984; Denny et al., 1985; Patterson, 1992a) provided understanding of how colony size could affect physiological performance (Patterson, 1992a,b). Since this work, evidence has emerged that coral polyps can be unequal between colonies differing in size, and within colonies depending on the position of polyps. These effects arise, for example, because biomass varies within colonies (Lough and Barnes, 2000) and among colonies differing in size (Vollmer and Edmunds, 2000), and genotypes of endosymbiotic Symbiodinium can vary among conspecifics differing in age and across the surface of single colonies (Little et al., 2004; Kemp et al., 2008). Light microenvironments differ across the surface of colonies, and especially within branching species (Kaniewska et al., 2011), and the capacity for polyps to capture zooplankton varies among positions within the colony (Sebens et al., 1998). Moreover, for branching corals, corallum morphology interacts with flow to create emergent properties affecting mass transfer as a result of the passage of seawater among the branches (Patterson, 1992b; Reidenbach et al., 2006). Finally, colony size modifies the metabolic rate of non-scleractinian colonial taxa, including bryozoans (Hughes and Hughes, 1986; Peck and Barnes, 2004; White et al., 2011) and ascidians (Nakaya et al., 2003, 2005). Together, the aforementioned effects create a rich landscape in which it is reasonable to expect that coral colony size will affect their response to rising seawater temperature and ocean acidification.

The objective of this study was to test the hypothesis that the response of a branching coral to high PCO2 (i.e. ocean acidification) and temperature is mediated by intraspecific variation in colony size. We focused on the coral Pocillopora verrucosa, which represents a morphological group (thick clustered branches, ‘branching-closed’ after www.coraltraits.org) that is important on many coral reefs, particularly in the Indo-Pacific (Veron, 2000; Pinzón et al., 2013; Schmidt-Roach et al., 2014), and is abundant in Moorea (where the study was conducted) (Edmunds et al., 2016). Colony size was evaluated by planar diameter, and three-dimensional tissue surface area, which typically is proportional to the total tissue biomass of each coral (Edmunds, 2011). A laboratory experiment was designed in which colonies differing in size were exposed during trials in 2014 and 2015 to combinations of temperature and PCO2, with the outcome evaluated through calcification, aerobic respiration and gross photosynthesis.

Overview

Pocillopora verrucosa (Ellis and Solander 1786) was incubated under combinations of two temperature and two PCO2 regimes, using colonies representing the size range of conspecifics found at 10 m depth on the outer reef of Moorea in April 2014 and April 2015. Pocillopora verrucosa was identified morphologically (Veron, 2000), although we cannot be certain that only one species was sampled (Edmunds et al., 2016). Colonies were incubated for 24 days (2014) or 25 days (2015) in flow-through tanks (150 l), creating four treatments in duplicate (8 tanks in 2014, Trial 1) or triplicate (12 tanks in 2015, Trial 2). To prevent the corals from affecting the total alkalinity (AT) of the seawater in the tanks through their metabolism, only one colony of each of four size ranges (described below) was placed in each tank, so that trials started with 32 corals in 2014 and 48 corals in 2015. Sample sizes were determined by results of previous experiments. Colony size was determined by planar diameter when collecting specimens in the field. For analyses, colony size was measured as tissue area and was treated as a continuously distributed variable to evaluate its relationships with physiology in the context of responding to the treatments (see ‘Statistical analyses’, below). All applicable international, national and/or institutional guidelines for the care and use of animals were followed throughout the study.

In our design, we reasoned that if physiological performance of coral modules (i.e. polyps) was independent of colony size, and polyp density was conserved, then area-normalized physiological rates (i.e. whole-colony rates divided by the tissue surface area) would be constant (a slope of zero) across a range of colony sizes, and whole-colony physiological rates would correspondingly increase in direct proportion to surface area. Similar slopes but differences in the intercepts of these relationships between temperature and PCO2 treatments would indicate that treatment effects were independent of colony size. If the slope of the relationship between colony size and physiological performance differed between temperature and PCO2 treatments (i.e. a two- or three-way interaction with colony size), it would indicate that the mean physiological performance of polyps depended on colony size, in a way that varied among treatments. This outcome could arise from one or more of three possibilities, with effects modulated by temperature and/or PCO2: polyps varying in functionality with colony size, polyps differing in functionality within a colony, or polyps varying in density (polyps cm−2) over the colony surface.

Trial 1, 2014

Coral collection

Thirty-two corals were collected on 9 April 2014 from 10 m depth at two sites ∼3 km apart on the north shore of Moorea, and were pooled by site. Mean (±s.e.m. throughout) ambient seawater temperature at the time of collection was 28.7±<0.1°C. Colonies were haphazardly selected in four size classes (1–3, 4–6, 7–8.5 and 8.5–14 cm diameter), removed from the substratum using a hammer and chisel, and returned to the shore submerged and shaded in buckets, where they were placed in a shallow tank supplied with flowing seawater. Colonies were trimmed of rock to allow them to rest upright, identified with a numbered tag, measured for size (as the mean of two planar diameters), and placed in a 1000 l tank for recovery and acclimation to laboratory conditions. The recovery tank was supplied with flowing seawater (∼10 l min−1) at 28.0°C, and illuminated (with SOL White LED, 6500 K, Aqua-Illumination, Ames, IA, USA) at ∼475 µmol photons m−2 s−1 of photosynthetically active radiation (PAR, as measured with a LI-Cor 4-π quantum sensor LI-193, Li-Cor Biosciences, Lincoln, NE, USA). Within the tank, corals were placed upright on a platform rotating at 2 revolutions per day to avoid position effects, and following 3–5 days, were buoyant weighed and allocated to treatment tanks for ∼24 days.

Incubation

Eight 150 l tanks (Aqualogic, San Diego, CA, USA) were randomly assigned to one of four treatments crossing target temperatures of 26.5°C (low temperature; LT) and 29.5°C (high temperature; HT), and targeted PCO2 values of 400 µatm (ambient PCO2; AC) and 1000 µatm (elevated PCO2; HC). The temperatures represented the seasonal range for seawater in Moorea (Edmunds et al., 2010) and the elevated PCO2 represented a pessimistic projection for atmospheric PCO2 by the end of the present century (Representative Concentration Pathway, RCP 8.5) (Riahi et al., 2007; van Vuuren et al., 2011). Seawater in each tank was independently chilled, heated, and circulated with a water pump, and tanks were illuminated (with SOL White LED, 6500 K, Aqua-Illumination) at a target PAR of ∼650 µmol photons m−2 s−1 (measured with a Li-Cor 4-π quantum sensor LI-193). The light intensity simulated the maximum midday light intensity at the collection depth on the outer reef in April (P.J.E., unpublished data). Each tank received fresh seawater at ∼200 ml min−1, which was pumped from ∼8 m depth in Cook's Bay and filtered through sand. One coral from each size class was placed in each tank on 14 April 2014, with corals selected for each tank at random. Thereafter, their position in each tank was randomly changed each day until the experiment ended on 9 May 2014.

CO2 treatments were created by bubbling ambient air (∼400 µatm PCO2), or a mixture of ambient air and pure CO2 (at ∼1000 µatm) into the tanks. The air and CO2 were mixed using a solenoid-controlled valve (Model A352, Qubit Systems, Kingston, ON, Canada) that was regulated dynamically by an infra-red gas analyzer (Model S151, Qubit) to blend the gases in a mixing chamber. Thereafter, the gas mixture was supplied at ∼10–15 l min−1 to the elevated PCO2 tanks. The conditions in the tanks were monitored throughout the experiment, with PAR (Li-Cor 4-π quantum sensor LI-193), temperature (±0.05°C certified, model 15-077, Fisher Scientific, Pittsburgh, PA, USA), salinity (YSI 3100 Conductivity Meter, YSI Inc., Yellow Springs, OH, USA) and pH (Mettler DG115-SC probe fitted to an Orion 3 star meter, Mettler-Toledo, LLC, Columbus, OH, USA) recorded daily using hand-held instruments, and seawater samples (∼50 ml) collected for the determination of seawater chemistry every 3–4 days. Seawater chemistry was evaluated according to standard operating procedures (SOP3b; Dickson et al., 2007), using potentiometric titrations in an open cell, automatic titrator (Model T50, Mettler-Toledo) fitted with a pH probe (Mettler DG115-SC) calibrated on the total scale using Tris/HCl buffers (Dickson et al., 2007). The accuracy of the titrations was evaluated using certified reference materials (CRMs) from A. Dickson (Scripps Institution of Oceanography; batch 122 and 130); measured values departed from certified values by a mean of 0.13±0.47% (N=9). Seawater pH was measured spectrophotometrically using m-Cresol Purple dye (no. 211761, Sigma-Aldrich, St Louis, MO, USA; SOP7, Dickson et al., 2007) and, together with AT determined from the titrations, was used to calculate dissolved inorganic carbon (DIC) parameters using the Seacarb package (https://cran.r-project.org/web/packages/seacarb/index.html) running in R (R Foundation for Statistical Computing). DIC parameters were calculated daily, using the pH recorded with the hand-held meter, salinity, temperature, and AT from the most recent previous analysis of seawater (≤4 days before); preliminary measurements demonstrated that AT was stable among determinations. On days when AT was determined and pH was measured spectrophotometrically, these values were used to calculate DIC parameters.

Response variables

Calcification rate, aerobic dark respiration and gross photosynthesis were used as response variables. Calcification was determined by buoyant weighing corals in seawater (Spencer-Davies, 1989) prior to placing them in the treatments on 16 April 2014, and when the experiment concluded on 8 May 2014. The change in buoyant mass was converted to dry mass using the density of aragonite (2.93 g cm−3) and the empirical density of seawater, and the change in mass was normalized to time and colony (mg day−1). The initial buoyant mass of the corals was converted to dry mass to characterize coral size based on the mass of the skeleton. The live coral tissue area was determined at the end of the experiment by wax dipping (Stimson and Kinzie, 1991), and was used as the measure of colony size, which was treated as a covariate in the statistical design.

Aerobic dark respiration (of the holobiont) and gross photosynthesis (of endosymbiotic Symbiodinium) were determined by measuring O2 flux in two confined chambers designed after Patterson et al. (1991). A limitation of this approach is that the contribution of the Symbiodinium to holobiont aerobic respiration cannot be determined, although classic work suggests this contribution is small (3–9%) (Muscatine et al., 1981), with recent work suggesting these numbers might be ∼31% high (Hawkins et al., 2016). The respiration chambers could not accommodate the largest corals, and measurements were limited to colonies ≤∼8 cm diameter. With two chambers, measuring respiration and gross photosynthesis of 32 corals took multiple days, and this task was completed in the third week of the experiment, with treatment tanks (at a single temperature and PCO2) selected at random for processing of the corals they contained. The chambers had a volume of ∼2130 ml and included a pump (Rule 136 l h−1) that could be operated at varying voltages to control the flow speed within the chamber. Photographs of hydrated brine shrimp eggs were used to evaluate flow speed (Sebens and Johnson, 1991), which was adjusted to ∼8 cm s−1 to create ecologically relevant conditions for a 10 m depth on the outer reef. An acoustic Doppler current profiler (RDI Workhorse, Teledyne RD Instruments, Poway, CA, USA) at 15 m depth on the outer reef (Washburn, 2016) was used to determine in situ flow speeds in a representative year (2013), and for 10 m depth, mean daily flow speed was 8.9±0.1 cm s−1 (N=365), and for April (the month in which the present study was conducted) it was 9.5±0.3 cm s−1 (N=30). Chambers were fitted with a Ruthenium-based optode (Foxy-R, Ocean Optics, Dunedin, FL, USA) connected to a spectrophotometer (USB 2000, Ocean Optics) and light source (USB LS-450) and operated by a PC running the manufacturer’s software (OOISensor). Optodes were calibrated in a zero solution of sodium sulfite and 0.01 mol l−1 sodium tetraborate, and 100% O2 using water-saturated air at the measurement temperature.

Aerobic respiration was measured in darkness, following 5–10 min in the chamber. O2 uptake was recorded continuously at 0.07 Hz for 10–30 min until a stable rate of decline in saturation was achieved. Trials were shortened to maintain O2 saturation ≥80% and, where necessary, aerated water was added for this purpose. Following measurements of respiration, net photosynthesis was recorded at ∼670 µmol photons m−2 s−1 PAR supplied with an LED lamp (SOL White, Aqua-Illumination, 6500 K), which was close to the maximum irradiance found in situ. Each daily set of measurements was accompanied by a control consisting of chambers filled with seawater alone. O2 fluxes in darkness (i.e. respiration) and light (i.e. photosynthesis) were determined by least squares linear regression, adjusted for control values, and converted to changes in O2 concentrations using the solubility of O2 in seawater (N. Ramsing and J. Gundersen, Unisense, Aarhus, Denmark). Aerobic dark respiration was added to net photosynthesis to obtain gross photosynthesis, and normalized to time and colony (in units of µmol O2 h−1). Surface area was determined by wax dipping as described above.

Trial 2 (2015)

The second trial was as similar as possible to Trial 1 to facilitate pooling of the results; any changes were either unavoidable modifications associated with trials in different years or were implemented to improve the experiment. Only the features of the experiment that changed between trials are described.

Coral collection

Forty-eight corals representing the four size classes were collected on 5 April 2015 from two sites at 10 m depth of the fore reef, and were pooled between sites. Mean ambient seawater temperature at the time of collection was 29.0±<0.1°C (for April 2015), which did not differ from the temperature in April 2014 (t=1.006, d.f.=30, P=0.323). Acclimation in the recovery tank took place over 9 days at ∼398 µmol photons m−2 s−1; the corals were buoyant weighed on 14 April 2015, and placed into treatment tanks for 25 days, with the final buoyant mass recorded on 9 May 2015.

Incubation

The experiment was expanded to 12, 150 l tanks that were randomly assigned to four treatments crossing the same target conditions as in Trial 1. High PCO2 treatments were maintained by direct bubbling of pure CO2 through diffusers that were regulated by solenoids (DePaul Submersible) controlled by a pHT-set point measured with an electrode (calibrated with Tris buffers, model PRBPH, Neptune Systems, Morgan Hill, CA, USA) and attached to a microprocessor controller (Apex Aquacontroller, Neptune Systems). Preliminary trials confirmed that the pH in the tanks was stable throughout the day, and did not drift overnight (as determined by analyzing seawater samples collected at 05:30 h). Seawater pH was logged at a frequency of 0.002 Hz and typically varied 0.02–0.04 between consecutive measures. Direct bubbling with CO2 provided more responsive pH regulation than was obtained in 2014, and therefore the flow-through of seawater was increased to 400 ml min−1. Seawater chemistry in the tanks was monitored as described in 2014, except that the pH used to calculate DIC parameters was measured spectrophotometrically using m-Cresol Purple dye. Throughout the experiment, certified reference material (batch 130 from A. Dickson, Scripps Institution of Oceanography) for AT was processed routinely, and measured values departed from certified values by 0.62±0.04% (mean±s.e., N=29).

Response variables

Response variables were measured in a similar way to that in 2014, except that three sizes of chambers were used for the measurement of respiration and photosynthesis (after Patterson et al., 1991), with volumes optimized to the size range of corals. The three chambers had volumes of ∼900 ml, ∼2130 ml (the same chamber used in 2014) and ∼6000 ml, and the pump circulating seawater in each chamber was operated at a unique voltage to obtain ∼8 cm s−1 flow speed in the working section.

Statistical analyses

Physical and chemical conditions were compared between treatment levels and the replicate tanks within each treatment using nested ANOVA in which treatment was the main effect and tanks were nested in each treatment level. These analyses were conducted using Systat 13 software, and the assumptions of ANOVA were tested through graphical analyses of residuals.

Our goals were to estimate the effect of temperature and PCO2 on the corals, using colony size as a covariate, and either colony- or area-normalized physiological rates as dependent variables in separate univariate analyses. The statistical approach involved fitting parametric models and describing the extent to which the quantitative relationships (represented by parameters for the slopes and intercepts) between size and calcification, aerobic respiration and gross photosynthesis differed between the temperature and PCO2 treatments. We used Gaussian linear mixed-effects models where tissue surface area was a continuous predictor, while temperature and PCO2 were categorical predictors (each with two levels, an ANCOVA-like design). Trials (years) and tanks within years were treated as random effects to create 20 groups from eight tanks in 2014, and 12 tanks in 2015. For each response variable, using tissue surface area as the predictor, different slopes and intercepts were estimated in each of the four treatments, while accounting for two sources of variation around the fitted values: (1) variation due to different years and tanks, and (2) variation due to individual corals within each tank.

Models were fitted to data in a Bayesian framework using Markov Chain Monte Carlo (MCMC) methods in the R package MCMCglmm (Hadfield, 2010). Flat (uninformative) inverse-Wishart priors were used. As the group-level variance (i.e. the estimate of variance due to years and tanks) was close to zero, a parameter expansion method was used to improve the mixing properties of the MCMC chain and provide a better description of the posterior distribution (following Hadfield, 2010). In all models, MCMC chains were run for 106 iterations with a burn-in of 3000 and a thinning interval of 10. The mixing properties of all parameters were checked to ensure no systematic trends, and the sensitivity of results to alternative prior specifications was also checked.

Treatment conditions and corals

Trial 1

Treatments were maintained with precision (Table S1), with mean salinities ranging from 33.8 to 34.0, and mean light intensities of 577–600 µmol photons m−2 s−1. Temperature differed between the temperature treatments (F1,6=506.729, P≤0.001) and among tanks (F6,480=2.617, P=0.017), although mean tank temperatures differed ≤0.6°C within each treatment; overall, treatments contrasted 26.5±0.1 with 29.7±≤0.1°C (mean±s.e., N=244). pH and PCO2 differed between PCO2 treatments (F1,60=339.550, P<0.991, and F1,6=17,705.274, P<0.001) and tanks (F6,195=4.709, P<0.001, and F6,204=57,144.154, P=0.001), although mean tank pH differed ≤0.07, and mean PCO2 differed ≤101 µatm; overall, treatment contrasted pH of 7.71±0.01 with 8.03±0.01 (mean±s.e., N=99–104) and 403±6 versus 982±16 µatm PCO2 (mean±s.e., N=100–112). AT did not vary between main effects (F1,6=0.359, P=0.571) and the grand mean (±s.e.) was 2256±5 µmol kg−1 (N=72).

The corals used in Trial 1 ranged in diameter from 2.1 to 13.1 cm, with tissue areas of 7–646 cm2, and initial dry mass of 5–538 g (Table 1). All corals were similar in color at the start of the experiment, and throughout the experiment maintained a natural polyp expansion cycle. None of the colonies died during the incubation, but when the experiment ended there were slight differences in color among the corals. Qualitatively, the colonies formed a gradient from pale to dark for the smallest to largest colonies in all treatments (Fig. 1).

Table 1.

Summary of the characteristics of colonies of Pocillopora verrucosa used in experiments in 2014 (Trial 1) and 2015 (Trial 2)

Summary of the characteristics of colonies of Pocillopora verrucosa used in experiments in 2014 (Trial 1) and 2015 (Trial 2)
Summary of the characteristics of colonies of Pocillopora verrucosa used in experiments in 2014 (Trial 1) and 2015 (Trial 2)
Fig. 1.

Colonies of Pocillopora verrucosa following incubations in 2014 (left) and 2015 (right). Corals were photographed in a single frame for each year; the scale bar for 2014 is shown in C, and that for 2015 is in E. Color cards are portions of Macbeth standards (McCamy et al., 2006), and corals in each quadrant correspond to the four treatments. Within each quadrant, corals are arranged by size class from left to right, with the rows showing the corals assigned to each tank (with two tanks per treatment in 2014, and three tanks per treatment in 2015 except LT-HC). (A) High temperature (HT)-high CO2 (HC), (B) low temperature (LT)-ambient CO2 (AC), (C) LT-HC, (D) HT-AC, (E) LT-HC, (F) LT-AC, (G) HT-HC and (H) HT-AC, where LT ∼26.5°C, HT ∼30.0°C, AC ∼400 µatm PCO2 and HC ∼1000 µatm PCO2 (Table S1).

Fig. 1.

Colonies of Pocillopora verrucosa following incubations in 2014 (left) and 2015 (right). Corals were photographed in a single frame for each year; the scale bar for 2014 is shown in C, and that for 2015 is in E. Color cards are portions of Macbeth standards (McCamy et al., 2006), and corals in each quadrant correspond to the four treatments. Within each quadrant, corals are arranged by size class from left to right, with the rows showing the corals assigned to each tank (with two tanks per treatment in 2014, and three tanks per treatment in 2015 except LT-HC). (A) High temperature (HT)-high CO2 (HC), (B) low temperature (LT)-ambient CO2 (AC), (C) LT-HC, (D) HT-AC, (E) LT-HC, (F) LT-AC, (G) HT-HC and (H) HT-AC, where LT ∼26.5°C, HT ∼30.0°C, AC ∼400 µatm PCO2 and HC ∼1000 µatm PCO2 (Table S1).

Trial 2

Treatments were maintained with precision (Table S1), with mean salinities ranging from 35.4 to 35.6, and mean light intensities of 671–713 µmol photons m−2 s−1. Temperature differed between the temperature treatments (F1,6=1373.247, P<0.001) and tanks (F10,533=2.726, P=0.003), with tanks differing <0.4°C within each treatment; overall treatments contrasted 26.5±<0.1°C and 29.7±<0.1°C (mean±s.e., N=271–274). pH and PCO2 differed between PCO2 treatments (F1,10=783.945, P<0.001, and F1,10=860.239, P<0.001) and among tanks within each treatment (F10,249=6.694, P<0.001, and F10,273=38.339.539, P<0.001), with treatment means of 7.662±0.004 versus 8.010±0.004 pH units and 436±4 versus 1118±11 µatm PCO2 (all means±s.e.). AT did not vary between main effects (F1,108=0.162, P=0.688) or tanks (F9,99=0.850, P=0.572) and the grand mean (±s.e.) was 2305±2 µmol kg−1 (N=120).

The corals used in Trial 2 ranged in diameter from 2.0 to 12.2 cm, with tissue areas of 8–534 cm2, and initial dry mass of 6–417 g (Table 1). All corals appeared similar in color at the start of the experiment, and throughout the experiment they maintained a natural expansion cycle of polyps. Of the 48 corals that started the experiment in 12 tanks, all corals in one tank in the LT-HC treatment died, and one of the smallest corals died in the HT-AC treatment. When the experiment ended, there were slight differences in color among the corals, and qualitatively, the colonies formed a gradient from pale to dark for the smallest to largest colonies in all treatments (Fig. 1).

Relationships between size and physiology as a function of temperature and PCO2

Calcification

All corals grew throughout the experiments, but the mass changes of two colonies (one LT-HC in 2014, and the other HT-HC in 2015) were not available because of logistical constraints. In 2014, the dry mass of colonies increased from 5.7 mg day−1 (a coral 2.3 cm diameter in LT-HC) to 397.0 mg day−1 (a coral 13.7 cm diameter in HT-AC), and in 2015, dry mass increased from 7.3 mg day−1 (a coral 2 cm diameter in LT-AC) to 351.9 mg day−1 (a coral 12.2 cm diameter in HT-HC). Whole-colony calcification increased with tissue surface area (Fig. 2A,B), though this relationship was not affected by temperature, PCO2 or the interaction between the two (Table 2). When normalized to tissue area, calcification (mg cm−2 day−1) was negatively related to colony size (Fig. 3A,B, Table 3), indicating that larger colonies calcified at slower area-normalized rates. Consequently, when scaled up to predict whole-coral calcification rate, area-normalized calcification rate for small corals tended to overestimate whole-colony calcification for larger corals (solid lines versus dashed lines in Fig. 2A,B).

Fig. 2.

Whole-colony responses of P. verrucosa to incubations under factorial combinations of temperature and PCO2, based on trials conducted in 2014 and 2015. Trials were conducted at 26.5 versus 29.7°C and at ambient (∼400 µatm) versus high (∼1000 µatm) PCO2, and values are plotted as a function of size (tissue surface area). (A,B) Calcification, (C,D) aerobic respiration and (E,F) gross photosynthesis. Each symbol represents the results from one coral colony, and the solid lines show the fitted partial regression coefficients from the linear mixed-effects model (Table 2). The dashed lines show the expected whole-colony responses calculated by scaling up the area-normalized responses for the average tissue surface area for corals in our small size class (14.49 cm2) based on parameters in Table 3 [i.e. if the area-normalized regression in each treatment is y=α+βx, where x is colony tissue surface area in cm2, then the expected whole-colony response would be yexpected=α + (β×14.49)x]. Sample sizes in Trial 1: N=8 for all treatment combinations; in Trial 2: N=12 in LT-AC and HT-AC, N=8 in LT-HC and N=11 in HT-HC.

Fig. 2.

Whole-colony responses of P. verrucosa to incubations under factorial combinations of temperature and PCO2, based on trials conducted in 2014 and 2015. Trials were conducted at 26.5 versus 29.7°C and at ambient (∼400 µatm) versus high (∼1000 µatm) PCO2, and values are plotted as a function of size (tissue surface area). (A,B) Calcification, (C,D) aerobic respiration and (E,F) gross photosynthesis. Each symbol represents the results from one coral colony, and the solid lines show the fitted partial regression coefficients from the linear mixed-effects model (Table 2). The dashed lines show the expected whole-colony responses calculated by scaling up the area-normalized responses for the average tissue surface area for corals in our small size class (14.49 cm2) based on parameters in Table 3 [i.e. if the area-normalized regression in each treatment is y=α+βx, where x is colony tissue surface area in cm2, then the expected whole-colony response would be yexpected=α + (β×14.49)x]. Sample sizes in Trial 1: N=8 for all treatment combinations; in Trial 2: N=12 in LT-AC and HT-AC, N=8 in LT-HC and N=11 in HT-HC.

Table 2.

Parameters from Bayesian linear mixed-effects models fit to three response variables measured in P. verrucosa: colony calcification, colony respiration and colony gross photosynthesis

Parameters from Bayesian linear mixed-effects models fit to three response variables measured in P. verrucosa: colony calcification, colony respiration and colony gross photosynthesis
Parameters from Bayesian linear mixed-effects models fit to three response variables measured in P. verrucosa: colony calcification, colony respiration and colony gross photosynthesis
Fig. 3.

Area-normalized responses of P. verrucosa to incubations under factorial combinations of temperature and PCO2, based on trials conducted in 2014 and 2015. Trials were conducted at 26.5 versus 29.7°C and at ambient (∼400 µatm) versus high (∼1000 µatm) PCO2, and values are plotted as a function of size (tissue surface area). (A,B) Calcification, (C,D) respiration and (E,F) gross photosynthesis. Each symbol represents the results from one coral colony, and the lines show the fitted partial regression coefficients from the linear mixed-effects model (Table 3). A flat relationship was expected if colony size did not affect the physiological performance of individual modules (polyps). Sample sizes in Trial 1: N=8 for all treatment combinations; in Trial 2: N=12 in LT-AC and HT-AC, N=8 in LT-HC and N=11 in HT-HC.

Fig. 3.

Area-normalized responses of P. verrucosa to incubations under factorial combinations of temperature and PCO2, based on trials conducted in 2014 and 2015. Trials were conducted at 26.5 versus 29.7°C and at ambient (∼400 µatm) versus high (∼1000 µatm) PCO2, and values are plotted as a function of size (tissue surface area). (A,B) Calcification, (C,D) respiration and (E,F) gross photosynthesis. Each symbol represents the results from one coral colony, and the lines show the fitted partial regression coefficients from the linear mixed-effects model (Table 3). A flat relationship was expected if colony size did not affect the physiological performance of individual modules (polyps). Sample sizes in Trial 1: N=8 for all treatment combinations; in Trial 2: N=12 in LT-AC and HT-AC, N=8 in LT-HC and N=11 in HT-HC.

Table 3.

Parameters from Bayesian linear mixed-effects models fitted to three response variables measured in P. verrucosa and calculated per unit surface area: calcification rate, respiration rate and photosynthesis rate

Parameters from Bayesian linear mixed-effects models fitted to three response variables measured in P. verrucosa and calculated per unit surface area: calcification rate, respiration rate and photosynthesis rate
Parameters from Bayesian linear mixed-effects models fitted to three response variables measured in P. verrucosa and calculated per unit surface area: calcification rate, respiration rate and photosynthesis rate

Aerobic respiration

Whole-colony respiration in 2014 ranged from 1.12 µmol O2 h−1 (LT-HC) to 118.92 µmol O2 h−1 (LT-AC), and in 2015 it ranged from 2.3 µmol O2 h−1 (LT-HC) to 479.7 µmol O2 h−1 (HT-HC). Whole-colony respiration increased with colony size (Fig. 2C,D), and the slopes of these relationships differed between some of the treatments (β4 and β5 in Table 2). Under ambient PCO2 and low temperature, whole-colony respiration increased by 1.0 µmol O2 h−1 (0.9–1.2 µmol O2 h−1, 96% credible interval, CI) for each 1 cm2 increase in tissue area, but was reduced by 0.3 µmol O2 h−1 (0.1–0.5 µmol O2 h−1, 95% CI) for each 1 cm2 increase in tissue area at high temperature (P=0.005) (Fig. 2C,D; β4 in Table 2). At low temperature, the slope of the relationship between whole-colony respiration and tissue surface area was reduced by 0.5 µmol O2 h−1 (0.3–0.7 µmol O2 h−1, 95% CI) for each cm2 increase in tissue area at high PCO2 compared with ambient PCO2 (P<0.001) (Fig. 2C; β5 in Table 2). When standardized to tissue area (Fig. 3C,D), respiration was greater for larger corals, but no differences were detected between temperature and PCO2 treatments (Table 3). When scaled up to predict whole-coral respiration, area-normalized respiration of small corals tended to underestimate whole-colony respiration of larger corals (solid lines versus dashed lines in Fig. 2C,D).

Gross photosynthesis

In 2014, whole-colony gross photosynthesis ranged from 0.7 to 204.1 µmol O2 h−1 (HT-AC), and in 2015 from 9.5 to 762.2 µmol O2 h−1 (HT-HC) (Fig. 2E,F). At ambient PCO2, high temperature increased the slope of the relationship between whole-colony gross photosynthesis and tissue area by 0.3 µmol O2 h−1 (0.0–0.7 µmol O2 h−1, 95% CI) for each cm2 increase in tissue area (P=0.087) compared with the slope at low temperature (Table 2). Area-normalized gross photosynthesis declined with increasing coral size (Fig. 3E,F, Table 3). At low temperature and ambient PCO2, photosynthesis declined by 0.0028 µmol O2 h−1 (0.0001–0.0056 µmol O2 h−1, 95% CI) for each cm2 increase in tissue area, with no detectable differences between treatments (Table 3). When scaled up to predict whole-coral photosynthetic rate, area-normalized estimates of photosynthetic rate for small corals tended to overestimate whole-colony photosynthesis for larger corals (solid lines versus dashed lines in Fig. 2E,F).

Overview

If physiological performance of coral polyps was independent of colony size, then physiological rates per unit surface area should be conserved in corals with a constant polyp density. Whole-colony physiological rates then would increase in proportion to surface area (i.e. dashed lines in Fig. 2), and scaling of physiological traits would be isometric. Our results for calcification, respiration and photosynthesis in P. verrucosa did not conform to this expectation, showing instead that area-normalized physiological rates are dependent on colony size as measured by tissue surface area, and therefore scale allometrically. For example, each unit area of coral surface, on average, grew less, consumed more oxygen and produced less oxygen through photosynthesis in larger versus smaller corals. Therefore, common measures of calcification, respiration and photosynthesis from small nubbins of branching corals are unlikely to accurately scale up to whole-coral physiological rates (compare dashed and solid lines in Fig. 2) that are needed for projections of future change in coral community structure (e.g. Pandolfi et al., 2011; Bramanti et al., 2015). Additionally, some size-dependent effects in P. verrucosa were mediated by temperature and PCO2, which has implications for evaluating the response of populations of this coral to climate change and ocean acidification. Compared with the scaling of whole-colony respiration with size at low temperature and ambient PCO2, whole-colony respiration increased with colony size by a smaller amount when temperature and PCO2 were increased independently, but increased with colony size by a greater amount under the combined effect of elevated temperature and PCO2 (Table 2). Compared with the scaling of whole-colony photosynthesis with size at low temperature, whole-colony photosynthesis increased with size by a larger amount under elevated temperature under both PCO2 conditions. Finally, and in contrast to previous research (Chan and Connolly, 2013; Kroeker et al., 2013; Comeau et al., 2014), our results from the full size range of P. verrucosa colonies found on the outer reefs of Moorea in 2014/2015 revealed no impact of temperature or PCO2 on calcification over 24–25 days. Large corals did, however, deposit less CaCO3 per day than small corals. Together, these results underscore the role of emergent properties of multi-branch coral colonies (versus fragments) in responding to environmental conditions.

Allometry in colonial corals

Our results are consistent with the hypothesis that polyps within a branching colony vary in performance with colony size, and interactions of corallum morphology with seawater flow, temperature and PCO2. This interpretation contradicts the notion that the colonial modular design of scleractinians allows them to escape the allometric constraints of increasing size (Kleiber, 1932; Schmidt-Nielsen, 1984; Hughes and Hughes, 1986) through iterations of independent modules that are conserved in dimensions and function. Classically, coral colonies are thought to increase in size through multiplication of standardized polyps, thereby conferring the biologically rare property of isometry (Jackson, 1979; Sebens, 1979). This principle has had a profound effect on coral biology, by providing a rationale for research relying on the assumption that analyses of small fragments of corals (i.e. nubbins; Birkeland, 1976) are insightful for understanding the biology of a wide range of sizes of coral colonies. Despite the pervasive nature of this assumption, 30 years ago there was evidence that it might not be correct (Jokiel and Morrissey, 1986), and a few years later Patterson (1992b) described allometric scaling in corals as a result of the interactive effects of seawater flow, colony morphology and boundary layers.

A few studies have addressed scaling in corals, and often have revealed allometric scaling through theoretical approaches and mensurative experimentation (Patterson, 1992a,b; Kim and Lasker, 1998; Vollmer and Edmunds, 2000; Edmunds, 2006; Elahi and Edmunds, 2007). Little attention has been paid to the implications of these discoveries, one result of which is that studies investigating the response of corals to ocean acidification (e.g. Chan and Connolly, 2013; Comeau et al., 2014) have relied on experiments with coral nubbins to support inferences applying to large colonies (and higher functional levels), without addressing the assumptions inherent in such extrapolation (Reynaud et al., 2003; Anthony et al., 2008; Edmunds, 2012, and references therein). The mechanisms underpinning an expectation of allometry in colonial corals such as P. verrucosa fall into three non-exclusive groups: effects related to colony size and morphology, effects related to interactions with seawater flow, and effects related to the functional biology of scleractinians.

First, colony size has the potential to be associated with differential allocation of resources (e.g. energy) (Anthony et al., 2002). Further, together with morphology, size can create within-colony heterogeneity that can drive allometry, for example, through variation in tissue biomass (Barnes and Lough, 1992; Edmunds, 2012; P.J.E., unpublished data for P. verrucosa), and the creation of shaded microhabitats that facilitate shade acclimatization (Kaniewska et al., 2011; Wangpraseurt et al., 2014). These effects can be augmented by varying amounts of coenosarc separating polyps, the consequences of within-colony shade acclimatization (Kaniewska et al., 2011; Wangpraseurt et al., 2015) and differential zooplankton capture by peripheral versus inner branches (Sebens et al., 1998), which together are likely to modulate the reliance on autotrophic versus heterotrophic nutrition (Anthony and Fabricius, 2000; Houlbrèque and Ferrier-Pagès, 2009). For branching corals, colony growth creates a complex morphology that modulates the physical and chemical environment within its interstices (Mass et al., 2010; Kaniewska et al., 2011; Chan et al., 2016), making it unlikely that large colonies functionally are the sum of their component branches. Second, consideration of mass transfer between tissue and seawater for organisms varying in size and morphology supports an expectation of allometry (Patterson, 1992a). This expectation arises from the role of size and shape in determining fluid motion around organisms (i.e. the Reynold's number), and the extent to which metabolism is enhanced by water motion (i.e. the Sherwood number) (Patterson, 1992a,b). Finally, fundamental aspects of coral biology favor allometry. For example, as the density and genotypes of the Symbiodinium symbionts of most reef corals can vary among locations within a colony (Rowan et al., 1997; Kemp et al., 2008), and Symbiodinium genotypes can have unique physiological properties (Sampayo et al., 2008; Parkinson and Baums, 2014), large colonies with numerous microhabitats among their branches may be capable of hosting a wider range of Symbiodinium genotypes, so that the physiology of large versus small colonies will not be proportional to their difference in size. Integration among polyps also can favor allometry, and while it contradicts the key assumption of module independence that underpins the expectation of isometry in colonial modular designs (Hughes, 2005), it is facilitated by the common gastrovascular cavity of colonial corals. Evidence of integration associated with chemical transfer among coral polyps is well known (Pearse and Muscatine, 1971; Oren et al., 2001), and has been evoked to explain temperature-associated switching in scaling of growth between isometry and allometry for small corals (Edmunds, 2006). A similar argument has been used to explain transitions between isometric and allometric scaling of metabolism in colonial ascidians (Nakaya et al., 2003).

Colony size and the response to temperature and PCO2

None of the response variables measured for P. verrucosa changed in magnitude in response to the individual or synergistic effects of elevated temperature or high PCO2. A few variables did, however, show a change in slope of response against size, suggesting that treatment effects were size specific. These outcomes generally contrast with the expectation of the responses of tropical reef corals to the treatment levels applied. For temperature, previous studies suggest that coral metabolism should be stimulated by temperature (Coles and Jokiel, 1977; Jokiel and Coles, 1977), at least for the temperatures employed herein and assuming thermal optima of ca. 28–29°C (Pratchett et al., 2015; Comeau et al., 2016). Expectations for the response to high PCO2 are less clear, but it is reasonable to hypothesize depressed calcification (Chan and Connolly, 2013; Comeau et al., 2014) and null effects for respiration and photosynthesis. For instance, PCO2 at a magnitude similar to that employed here (i.e. causing a reduction in aragonite saturation, Ωarag, of ∼1.5) should depress calcification ∼23% based on the results of a meta-analysis (Chan and Connolly, 2013), and previous studies of Pocillopora spp. show declines of ∼55% at 913 µatm PCO2 and 30.1°C (for P. meandrina; Brown and Edmunds, 2016), ∼10% at 986 µatm PCO2 and 27°C (P. verrucosa; Comeau et al., 2014) and 20% at 2072 µatm PCO2 and 27°C (P. damicornis; Comeau et al., 2013). Calcification of P. damicornis was unaffected by PCO2 ≤1000 µatm (at ∼27°C) in one study that contrasted results from Moorea, Hawaii and Okinawa (Comeau et al., 2014). For respiration, despite the thermodynamic prediction that calcification will be energetically more costly at high PCO2 (Erez et al., 2011; Pandolfi et al., 2011), and therefore will cause respiration to increase (McCulloch et al., 2012), most tests of this hypothesis have found no effect (Reynaud et al., 2003; Schneider and Erez, 2006; Comeau et al., 2016), or a decline (Edmunds, 2012; Kaniewska et al., 2012), including for P. damicornis (normalized by biomass, but not area; Comeau et al., 2016). For photosynthesis, the effects of high PCO2 also are equivocal, with reports of increases (Reynaud et al., 2003), decreases (Anthony et al., 2008; Kaniewska et al., 2012) and no effect (Takahashi and Kurihara, 2013; Comeau et al., 2016), which together support a null outcome in a meta-analysis (Kroeker et al., 2013).

While the present experiment does not allow cause and effect to be established for the scaling of relationships, our analyses make a compelling case for allometry. In addition to the scaling principles and examples described above, two recent studies may have relevance in accounting for the present results. First, Jokiel (2011) has posited the ‘two compartment proton flux model’ to reconcile aspects of the morphological design of symbiotic reef corals with the physiological requirements of photosynthesis and calcification. In his model, coral morphology promotes spatial separation within a colony of an inner ‘zone of rapid photosynthesis’ and a peripheral ‘zone of rapid calcification’, that matches the restriction of rapid calcification in branching corals to branch apices (Pearse and Muscatine, 1971; Jokiel, 2011), and facilitates metabolic interchange between the two areas to promote functionality (Jokiel, 2011). Second, Chan et al. (2016) have shown that low flow speeds favor the formation of thick boundary layers in Favites sp. and P. damicornis, in which the metabolic activity of the coral tissue creates pH levels adjacent to the tissue that are elevated relative to that of the ambient seawater. Therefore, under reduced seawater pH as a result of ocean acidification, the interaction of flow and coral morphology has the potential to ameliorate the negative effects of reduced pH on calcification. Together, the hypothesis of Jokiel (2011) and the results of Chan et al. (2016) suggest that morphology in branching corals can attenuate the effects of high seawater PCO2 in a colony size-dependent way, potentially creating organismic results similar to those reported here (Fig. 2A,B).

Our results also reveal size-dependent responses of respiration and photosynthesis to temperature (at ambient PCO2), and of respiration to high PCO2 (at low temperature), and the interactive effects of temperature and PCO2 (HT-HC versus LT-AC; Table 2). None of these effects correspond to a shift in magnitude of the response variable, and thus the response to temperature and PCO2 will depend on the size of colonies considered. Further research will be required to identify the causes and implications of these trends, but given the influence of allometry for P. verrucosa, it is likely that these objectives will be well served by considering the emergent properties of branching morphologies.

Summary and future directions

Our study is the first to address the effects of colony size on the response of a branching coral to elevated temperature and high PCO2, and is among only a handful that have addressed the role of colony size in responding to these conditions individually (Hoeksema, 1991; Shenkar et al., 2005; Brandt, 2009; van Woesik et al., 2012; Chan et al., 2016). There are two important contributions from this study. First, we build from foundational, but overlooked studies (cited above) to underscore the prevalence of allometry in colonial modular cnidarians and, therefore, to challenge the strong tendency to accept isometric scaling for this Bauplan. Our experiments show that area-normalized estimates of physiology from small corals provide misleading estimates of physiological performance of larger colonies. Second, we reveal a null result for the effects of elevated temperature and high PCO2 on the calcification of a range of colony sizes, which we argue reflects properties emerging with the growth of single branches into larger colonies. As investigations of the effects of ocean acidification and climate change on coral reefs expand from coral fragments and small colonies to communities and ecosystems (Edmunds et al., 2016), greater attention will need to be given to emergent properties of coral colonies differing in size.

We thank V. Moriarty and C. Lantz for technical support, A. Ellis, J. Smolenski, A. Estrada and A. Yarid for field assistance, and the staff of the UC Berkeley Richard B. Gump South Pacific Research Station for making our visits to Moorea productive and enjoyable.

Author contributions

P.J.E. and S.C.B. designed the experiments, P.J.E. conducted all empirical research, S.C.B. analyzed primary results, P.J.E. and S.C.B. wrote the manuscript.

Funding

This research was funded by the US National Science Foundation through the Moorea Coral Reef (MCR) LTER (OCE 10-26851 and OCE 12-36905) and grants for research on ocean acidification (OCE 10-41270, 14-15268). This is a product of the MCR-LTER and is contribution number 246 of the Marine Biology Program of California State University, Northridge.

Data availability

Data in support of this manuscript are archived at the LTER Network Catalog https://doi.org/10.6073/PASTA/840E360F30E8028EEECF800548E33A55 and are also available at the MCR LTER site http://mcr.lternet.edu/cgi-bin/showDataset.cgi?docid=knb-lter-mcr.5018.

Anthony
,
K. R. N.
and
Fabricius
,
K. E.
(
2000
).
Shifting roles of heterotrophy and autotrophy in coral energetics under varying turbidity
.
J. Exp. Mar. Biol. Ecol.
252
,
221
-
253
.
Anthony
,
K. R. N.
,
Connolly
,
S. R.
and
Willis
,
B. L.
(
2002
).
Comparative analysis of energy allocation to tissue and skeletal growth in corals
.
Limnol. Oceanogr.
47
,
1417
-
1429
.
Anthony
,
K. R. N.
,
Connolly
,
S. R.
and
Hoegh-Guldberg
,
O.
(
2007
).
Bleaching, energetics, and coral mortality risk: effects of temperature, light, and sediment regime
.
Limnol. Oceanogr.
52
,
716
-
726
.
Anthony
,
K. R. N.
,
Kline
,
D. I.
,
Diaz-Pulido
,
G.
,
Dove
,
S.
and
Hoegh-Guldberg
,
O.
(
2008
).
Ocean acidification causes bleaching and productivity loss in coral reef builders
.
Proc. Natl. Acad. Sci. USA
105
,
17442
-
17446
.
Barnes
,
D. J.
(
1973
).
Growth in colonial Scleractinians
.
Bull. Mar. Sci.
23
,
280
-
298
.
Barnes
,
D. J.
and
Lough
,
J. M.
(
1992
).
Systematic variations in the depth of skeleton occupied by coral tissue in massive colonies of Porites from the Great Barrier Reef
.
J. Exp. Mar. Biol. Ecol.
159
,
113
-
128
.
Birkeland
,
C.
(
1976
).
An experimental method of studying corals during early stages of growth
.
Micronesia
12
,
319
-
322
.
Bramanti
,
L.
,
Iannelli
,
M.
,
Fan
,
T. Y.
and
Edmunds
,
P. J.
(
2015
).
Using demographic models to project the effects of climate change on scleractinian corals: Pocillopora damicornis as a case study
.
Coral Reefs
34
,
505
-
515
.
Brandt
,
M. E.
(
2009
).
The effect of species and colony size on the bleaching response of reef-building corals in the Florida Keys during the 2005 mass bleaching event
.
Coral Reefs
28
,
911
-
924
.
Brown
,
D.
and
Edmunds
,
P. J.
(
2016
).
Differences in the responses of three scleractinians and the hydrocoral Millepora platyphylla to ocean acidification
.
Mar. Biol.
163
,
62
-
72
.
Brown
,
J. H.
,
Gillooly
,
J. F.
,
Allen
,
A. P.
,
Savage
,
V. M.
and
West
,
G. B.
(
2004
).
Toward a metabolic theory of ecology
.
Ecology
85
,
1771
-
1789
.
Bruno
,
J. F.
,
Carr
,
L. A.
and
O'Connor
,
M. I.
(
2015
).
Exploring the role of temperature in the ocean through metabolic scaling
.
Ecology
96
,
3126
-
3140
.
Chan
,
N. C. S.
and
Connolly
,
S. R.
(
2013
).
Sensitivity of coral calcification to ocean acidification: a meta-analysis
.
Glob. Change Biol.
19
,
282
-
290
.
Chan
,
N. C. S.
,
Wangpraseurt
,
D.
,
Kühl
,
M.
and
Connolly
,
S. R.
(
2016
).
Flow and coral morphology control coral surface pH: implications for the effects of ocean acidification
.
Front. Mar. Sci.
3
,
10
.
Coles
,
S. L.
and
Jokiel
,
P. L.
(
1977
).
Effects of temperature on photosynthesis and respiration in hermatypic corals
.
Mar. Biol.
43
,
209
-
216
.
Comeau
,
S.
,
Edmunds
,
P. J.
,
Spindel
,
N. B.
and
Carpenter
,
R. C.
(
2013
).
The responses of eight coral reef calcifiers to increasing partial pressure of CO2 do not exhibit a tipping point
.
Limnol. Oceanogr.
58
,
388
-
398
.
Comeau
,
S.
,
Edmunds
,
P. J.
,
Spindel
,
N. B.
and
Carpenter
,
R. C.
(
2014
).
Fast coral reef calcifiers are more sensitive to ocean acidification in short-term laboratory incubations
.
Limnol. Oceanogr.
59
,
1081
-
1091
.
Comeau
,
S.
,
Carpenter
,
R. C.
,
Lantz
,
C. A.
and
Edmunds
,
P. J.
(
2016
).
Parameterization of the response of calcification to temperature and pCO2 in the coral Acropora pulchra and the alga Lithophyllum kotschyanum
.
Coral Reefs
35
,
929
-
939
.
Davies
,
P. S.
(
1984
).
The role of zooxanthellae in the nutritional energy requirements of Pocillopora eydouxi
.
Coral Reefs
2
,
181
-
186
.
Denny
,
M. W.
,
Daniel
,
T. L.
and
Koehl
,
M. A. R.
(
1985
).
Mechanical limits to size in wave-swept organisms
.
Ecol. Monogr.
55
,
69
-
102
.
Dickson
,
A. G.
,
Sabine
,
C. L.
and
Christian
,
J. R.
(ed.) (
2007
).
Guide to best practices for ocean CO2 measurements
.
PICES Special Publ.
3
,
191
.
Edmunds
,
P. J.
(
2005
).
Effect of elevated temperature on aerobic respiration of coral recruits
.
Mar. Biol.
146
,
655
-
663
.
Edmunds
,
P. J.
(
2006
).
Temperature-mediated transitions between isometry and allometry in a colonial, modular invertebrate
.
Proc. R. Soc. B Biol. Sci.
273
,
2275
-
2281
.
Edmunds
,
P. J.
(
2011
).
Zooplanktivory ameliorates the effects of ocean acidification on the reef coral Porites spp
.
Limnol. Oceanogr.
56
,
2402
-
2410
.
Edmunds
,
P. J.
(
2012
).
Effect of pCO2 on the growth, respiration, and photophysiology of massive Porites spp. in Moorea, French Polynesia
.
Mar. Biol.
159
,
2149
-
2160
.
Edmunds
,
P. J.
and
Davies
,
P. S.
(
1986
).
An energy budget for Porites porites (Scleractinia)
.
Mar. Biol.
92
,
339
-
347
.
Edmunds
,
P. J.
,
Leichter
,
J. J.
and
Adjeroud
,
M.
(
2010
).
Landscape-scale variation in coral recruitment in Moorea, French Polynesia
.
Mar. Ecol. Prog. Ser.
414
,
75
-
89
.
Edmunds
,
P. J.
,
Burgess
,
S. C.
,
Putnam
,
H. M.
,
Baskett
,
M. L.
,
Bramanti
,
L.
,
Fabina
,
N. S.
,
Han
,
X.
,
Lesser
,
M. P.
,
Madin
,
J. S.
,
Wall
,
C. B.
, et al. 
(
2014
).
Evaluating the causal basis of ecological success within the Scleractinia: an integral projection model approach
.
Mar. Biol.
161
,
2719
-
2734
.
Edmunds
,
P. J.
,
Comeau
,
S.
,
Lantz
,
C.
,
Andersson
,
A.
,
Briggs
,
C.
,
Cohen
,
A.
,
Gattuso
,
J.-P.
,
Grady
,
J. M.
,
Gross
,
K.
,
Johnson
,
M.
, et al. 
(
2016
).
Integrating the effects of ocean acidification across functional scales on tropical coral reefs
.
BioScience
66
,
350
-
362
.
Elahi
,
R.
and
Edmunds
,
P. J.
(
2007
).
Determinate growth and the scaling of photosynthetic energy intake in the solitary coral Fungia concinna (Verrill)
.
J. Exp. Mar. Biol. Ecol.
349
,
183
-
193
.
Erez
,
J.
,
Reynaud
,
S.
,
Silverman
,
J.
,
Schneider
,
K.
and
Allemand
,
D.
(
2011
).
Coral calcification under ocean acidification and global change
. In
Coral Reefs: An Ecosystem in Transition
(ed.
Z.
Dubinsky
and
N.
Stambler
), pp.
151
-
176
.
Berlin
:
Springer Science+Business Media
.
Evenhuis
,
C.
,
Lenton
,
A.
,
Cantin
,
N. E.
and
Lough
,
J. M.
(
2015
).
Modelling coral calcification accounting for the impacts of coral bleaching and ocean acidification
.
Biogeosciences
12
,
2607
-
2630
.
Fitt
,
W. K.
,
McFarland
,
F. K.
,
Warner
,
M. E.
and
Chilcoat
,
G. C.
(
2000
).
Seasonal patterns of tissue biomass and densities of symbiotic dinoflagellates in reef corals and relation to coral bleaching
.
Limnol. Oceanogr.
45
,
677
-
685
.
Gaylord
,
B.
,
Kroeker
,
K. J.
,
Sunday
,
J. M.
,
Anderson
,
K. M.
,
Barry
,
J. P.
,
Brown
,
N. E.
,
Connell
,
S. D.
,
Dupont
,
S.
,
Fabricius
,
K. E.
,
Hall-Spencer
,
J. M.
, et al. 
(
2015
).
Ocean acidification through the lens of ecological theory
.
Ecology
96
,
3
-
15
.
Hadfield
,
J. D.
(
2010
).
MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package
.
J. Stat. Software
33
,
1
-
22
.
Hawkins
,
T. D.
,
Hagemeyer
,
J. C. G.
,
Hoadley
,
K. D.
,
Marsh
,
A. G.
and
Warner
,
M. E.
(
2016
).
Partitioning of respiration in an animal-algal symbiosis: implications for different aerobic capacity between Symbiodinium spp
.
Front. Physiol.
7
,
128
.
Hoegh-Guldberg
,
O.
,
Mumby
,
P. J.
,
Hooten
,
A. J.
,
Steneck
,
R. S.
,
Greenfield
,
P.
,
Gomez
,
E.
,
Harvell
,
C. D.
,
Sale
,
P. F.
,
Edwards
,
A. J.
,
Caldeira
,
K.
, et al. 
(
2007
).
Coral reefs under rapid climate change and ocean acidification
.
Science
318
,
1737
-
1742
.
Hoeksema
,
B. W.
(
1991
).
Control of bleaching in mushroom coral populations (Scleractinia: Fungiidae) in the Java Sea: stress tolerance and interference by life history strategy
.
Mar. Ecol. Prog. Ser.
74
,
225
-
237
.
Houlbrèque
,
F.
and
Ferrier-Pagès
,
C.
(
2009
).
Heterotrophy in tropical scleractinian corals
.
Biol. Rev.
84
,
1
-
17
.
Hughes
,
R. N.
(
2005
).
Lessons in modularity: the evolutionary ecology of colonial invertebrates
.
Sci. Mar.
69
,
169
-
179
.
Hughes
,
D. J.
and
Hughes
,
R. N.
(
1986
).
Metabolic implications of modularity: studies on the respiration and growth of Electra pilosa
.
Philos. Trans. R. Soc. B Biol. Sci.
313
,
23
-
29
.
Jackson
,
J. B. C.
(
1979
).
Morphological strategies in sessile animals
. In
Population Biology and Evolution of Clonal Organisms
(ed.
J. B. C.
Jackson
,
L. W.
Buss
and
R. E.
Cook
), pp.
499
-
555
.
New Haven
:
Yale University Press
.
Jokiel
,
P. L.
(
2011
).
The reef coral two compartment proton flux model: a new approach relating tissue-level physiological processes to gross corallum morphology
.
J. Exp. Mar. Biol. Ecol.
409
,
1
-
12
.
Jokiel
,
P. L.
and
Coles
,
S. L.
(
1977
).
Effects of temperature on the mortality and growth of Hawaiian reef corals
.
Mar. Biol.
43
,
201
-
208
.
Jokiel
,
P. L.
and
Morrissey
,
J. I.
(
1986
).
Influence of size on primary production in the reef coral Pocillopora damicornis and the macroalga Acanthophora spicifera
.
Mar. Biol.
91
,
15
-
26
.
Kaniewska
,
P.
,
Magnusson
,
S. H.
,
Anthony
,
K. R. N.
,
Reef
,
R.
,
Kühl
,
M.
and
Hoegh-Guldberg
,
O.
(
2011
).
Importance of macro- versus microstructure in modulating light levels inside coral colonies
.
J. Phycol.
47
,
846
-
860
.
Kaniewska
,
P.
,
Campbell
,
P. R.
,
Kline
,
D. I.
,
Rodriguez-Lanetty
,
M.
,
Miller
,
D. J.
,
Dove
,
S.
and
Hoegh-Guldberg
,
O.
(
2012
).
Major cellular and physiological impacts of ocean acidification on a reef building coral
.
PLoS ONE
7
,
e34659
.
Kemp
,
D. W.
,
Fitt
,
W. K.
and
Schmidt
,
G. W.
(
2008
).
A microsampling method for genotyping coral symbionts
.
Coral Reefs
27
,
289
-
293
.
Kim
,
K.
and
Lasker
,
H. R.
(
1998
).
Allometry of resource capture in colonial cnidarians and constraints on modular growth
.
Funct. Ecol.
12
,
646
-
654
.
Kleiber
,
M.
(
1932
).
Body size and metabolism
.
Hilgardia: J. Agric. Sci.
6
,
315
-
353
.
Koehl
,
M. A. R.
(
1984
).
How do benthic organisms withstand moving water?
Integr. Comp. Biol.
24
,
57
-
70
.
Kroeker
,
K. J.
,
Kordas
,
R. L.
,
Crim
,
R.
,
Hendriks
,
I. E.
,
Ramajo
,
L.
,
Singh
,
G. S.
,
Duarte
,
C. M.
and
Gattuso
,
J.-P.
(
2013
).
Impacts of ocean acidification on marine organisms: quantifying sensitivities and interaction with warming
.
Glob. Change Biol.
19
,
1884
-
1896
.
Little
,
A. F.
,
van Oppen
,
M. J. H.
and
Willis
,
B. L.
(
2004
).
Flexibility in algal endosymbioses shapes growth in reef corals
.
Science
304
,
1492
-
1494
.
Lough
,
J. M.
and
Barnes
,
D. J.
(
2000
).
Environmental controls on growth of the massive coral Porites
.
J. Exp. Mar. Biol. Ecol.
245
,
225
-
243
.
Mass
,
T.
,
Genin
,
A.
,
Shavit
,
U.
,
Grinstein
,
M.
and
Tchernov
,
D.
(
2010
).
Flow enhances photosynthesis in marine benthic autotrophs by increasing the efflux of oxygen from the organism to the water
.
Proc. Natl. Acad. Sci. USA
107
,
2527
-
2531
.
McCamy
,
C. S.
,
Marcus
,
H.
and
Davidson
,
J. G.
(
1976
).
A color-rendition chart
.
J. Appl. Photogr. Eng.
2
,
95
-
99
.
McCulloch
,
M.
,
Falter
,
J.
,
Trotter
,
J.
and
Montagna
,
P.
(
2012
).
Coral resilience to ocean acidification and global warming through pH up-regulation
.
Nat. Clim. Change
2
,
623
-
627
.
Muscatine
,
L.
,
McCloskey
,
L. R.
and
Marian
,
R. E.
(
1981
).
Estimating the daily contribution of carbon from zooxanthellae to coral animal respiration
.
Limnol. Oceanogr.
26
,
601
-
611
.
Nakaya
,
F.
,
Saito
,
Y.
and
Motokawa
,
T.
(
2003
).
Switching of metabolic-rate scaling between allometry and isometry in colonial ascidians
.
Proc. R. Soc. B Biol. Sci.
270
,
1105
-
1113
.
Nakaya
,
F.
,
Saito
,
Y.
and
Motokawa
,
T.
(
2005
).
Experimental allometry: effect of size manipulation on metabolic rate of colonial ascidians
.
Proc. R. Soc. B Biol. Sci.
272
,
1963
-
1969
.
Nisbet
,
R. M.
,
Muller
,
E. B.
,
Lika
,
K.
and
Kooijman
,
S. A. L. M.
(
2000
).
From molecules to ecosystems through dynamic energy budget models
.
J. Anim. Ecol.
69
,
913
-
926
.
Oren
,
U.
,
Benayahu
,
Y.
,
Lubinevsky
,
H.
and
Loya
,
Y.
(
2001
).
Colony integration during regeneration in the stony coral Favia favus
.
Ecology
82
,
802
-
813
.
Pandolfi
,
J. M.
,
Connolly
,
S. R.
,
Marshall
,
D. J.
and
Cohen
,
A. L.
(
2011
).
Projecting coral reef futures under global warming and ocean acidification
.
Science
333
,
418
-
422
.
Parkinson
,
J. E.
and
Baums
,
I. B.
(
2014
).
The extended phenotypes of marine symbioses: ecological and evolutionary consequences of intraspecific genetic diversity in coral-algal associations
.
Front. Microbiol.
5
,
445
.
Patterson
,
M. R.
(
1992a
).
A mass transfer explanation of metabolic scaling relations in some aquatic invertebrates and algae
.
Science
255
,
1421
-
1423
.
Patterson
,
M. R.
(
1992b
).
A chemical engineering view of cnidarian symbioses
.
Integr. Comp. Biol.
32
,
566
-
582
.
Patterson
,
M. R.
,
Sebens
,
K. P.
and
Olson
,
R. R.
(
1991
).
In situ measurements of flow effects on primary production and dark respiration in reef corals
.
Limnol. Oceanogr.
36
,
936
-
948
.
Pearse
,
V. B.
and
Muscatine
,
L.
(
1971
).
Role of symbiotic algae (zooxanthellae) in coral calcification
.
Biol. Bull.
141
,
350
-
363
.
Peck
,
L. S.
and
Barnes
,
D. K. A.
(
2004
).
Metabolic flexibility: the key to long-term evolutionary success in Bryozoa?
Proc. R. Soc B Biol. Sci.
271
Suppl.
,
S18
-
S21
.
Pinzón
,
J. H.
,
Sampayo
,
E.
,
Cox
,
E.
,
Chauka
,
L. J.
,
Chen
,
C. A.
,
Voolstra
,
C. R.
and
Lajeunesse
,
T. C.
(
2013
).
Blind to morphology: genetics identifies several widespread ecologically common species and few endemics among Indo-Pacific cauliflower corals (Pocillopora, Scleractinia)
.
J. Biogeogr.
40
,
1595
-
1608
.
Pratchett
,
M. S.
,
Anderson
,
K. D.
,
Hoogenboom
,
M. O.
,
Widman
,
E.
,
Baird
,
A. H.
,
Pandolfi
,
J. M.
,
Edmunds
,
P. J.
and
Lough
,
J. M.
(
2015
).
Spatial, temporal and taxonomic variation in coral growth–implications for the structure and function of coral reef ecosystems
.
Oceanogr. Mar. Biol. Annu. Rev.
53
,
215
-
296
.
Reidenbach
,
M. A.
,
Koseff
,
J. R.
,
Monismith
,
S. G.
and
Steinbuckc
,
J. V.
(
2006
).
The effects of waves and morphology on mass transfer within branched reef corals
.
Limnol. Oceanogr.
51
,
1134
-
1141
.
Reynaud
,
S.
,
Leclercq
,
N.
,
Romaine-Liuod
,
S.
,
Ferrier-Pagés
,
C.
,
Jaubert
,
J.
and
Gattuso
,
J.-P.
(
2003
).
Interacting effects of CO2 partial pressure and temperature on photosynthesis and calcification in a scleractinian coral
.
Glob. Change Biol.
9
,
1660
-
1668
.
Riahi
,
K.
,
Grübler
,
A.
and
Nakicenovic
,
N.
(
2007
).
Scenarios of long-term socio-economic and environmental development under climate stabilization
.
Technol. Forecast. Soc. Change
74
,
887
-
935
.
Rowan
,
R.
,
Knowlton
,
N.
,
Baker
,
A.
and
Jara
,
J.
(
1997
).
Landscape ecology of algal symbionts creates variation in episodes of coral bleaching
.
Nature
388
,
265
-
269
.
Sampayo
,
E. M.
,
Ridgway
,
T.
,
Bongaerts
,
P.
and
Hoegh-Guldberg
,
O.
(
2008
).
Bleaching susceptibility and mortality of corals are determined by fine-scale differences in symbiont type
.
Proc. Natl. Acad. Sci. USA
105
,
10444
-
10449
.
Schmidt-Nielsen
,
K.
(
1984
).
Scaling: Why is Animal Size So Important?
New York
:
Cambridge University Press
.
Schmidt-Roach
,
S.
,
Miller
,
K. J.
,
Lundgren
,
P.
and
Andreakis
,
N.
(
2014
).
With eyes wide open: a revision of species within and closely related to the Pocillopora damicornis species complex (Scleractinia; Pocilloporidae) using morphology and genetics
.
Zool. J. Linn. Soc.
170
,
1
-
33
.
Schneider
,
K.
and
Erez
,
J.
(
2006
).
The effect of carbonate chemistry on calcification and photosynthesis in the hermatypic coral Acropora eurystoma
.
Limnol. Oceanogr.
51
,
1284
-
1293
.
Sebens
,
K. P.
(
1979
).
The energetics of asexual reproduction and colony formation in benthic marine invertebrates
.
Am. Zool.
19
,
683
-
699
.
Sebens
,
K. P.
(
1987
).
The ecology of indeterminate growth in animals
.
Annu. Rev. Ecol. Syst.
18
,
371
-
407
.
Sebens
,
K. P.
and
Johnson
,
A. S.
(
1991
).
Effects of water movement on prey capture and distribution of reef corals
.
Hydrobiologia
226
,
91
-
101
.
Sebens
,
K. P.
,
Grace
,
S. P.
,
Helmuth
,
B.
,
Maney
,
E. J.
, Jr
and
Miles
,
J. S.
(
1998
).
Water flow and prey capture by three scleractinian corals, Madracis mirabilis, Montastrea cavernosa and Porites porites, in a field enclosure
.
Mar. Biol.
131
,
347
-
360
.
Shenkar
,
N.
,
Fine
,
M.
and
Loya
,
Y.
(
2005
).
Size matters: bleaching dynamics of the coral Oculina patagonica
.
Mar. Ecol. Prog. Ser.
294
,
181
-
188
.
Spencer-Davies
,
P.
(
1989
).
Short-term growth measurements of corals using an accurate buoyant weighing technique
.
Mar. Biol.
101
,
389
-
395
.
Stimson
,
J.
and
Kinzie
,
R. A.
(
1991
).
The temporal pattern and rate of release of zooxanthellae from the reef coral Pocillopora damicornis (Linnaeus) under nitrogen-enrichment and control conditions
.
J. Exp. Mar. Biol. Ecol.
153
,
63
-
74
.
Takahashi
,
A.
and
Kurihara
,
H.
(
2013
).
Ocean acidification does not affect the physiology of the tropical coral Acropora digitifera during a 5-week experiment
.
Coral Reefs
32
,
305
-
314
.
van Vuuren
,
D. P.
,
Edmonds
,
J.
,
Kainuma
,
M.
,
Riahi
,
K.
,
Thomson
,
A.
,
Hibbard
,
K.
,
Hurtt
,
G. C.
,
Kram
,
T.
,
Krey
,
V.
,
Lamarque
,
J.-F.
, et al. 
(
2011
).
The representative concentration pathways: an overview
.
Climatic Change
109
,
5
-
31
.
van Woesik
,
R.
,
Irikawa
,
A.
,
Anzai
,
R.
and
Nakamura
,
T.
(
2012
).
Effects of coral colony morphologies on mass transfer and susceptibility to thermal stress
.
Coral Reefs
31
,
633
-
639
.
Veron
,
J. E. N.
(
2000
).
Corals of the World
. Vol.
1-3
.
Townsville
,
Australia
:
Australian Institute of Marine Science
.
Vollmer
,
S. V.
and
Edmunds
,
P. J.
(
2000
).
Allometric scaling in small colonies of the scleractinian coral Siderastrea siderea (Ellis and Solander)
.
Biol. Bull.
199
,
21
-
28
.
Wangpraseurt
,
D.
,
Polerecky
,
L.
,
Larkum
,
A. W. D.
,
Ralph
,
P. J.
,
Nielsen
,
D. A.
,
Pernice
,
M.
and
Kühl
,
M.
(
2014
).
The in situ light microenvironment of corals
.
Limnol. Oceanogr.
59
,
917
-
926
.
Wangpraseurt
,
D.
,
Pernice
,
M.
,
Guagliardo
,
P.
,
Kilburn
,
M. R.
,
Clode
,
P. L.
,
Polerecky
,
L.
and
Kühl
,
M.
(
2015
).
Light microenvironment and single-cell gradients of carbon fixation in tissues of symbiont-bearing corals
.
ISME J.
10
,
788
-
792
.
Washburn
,
L.
(
2016
).
Moorea Coral Reef LTER: Coral Reef: Ocean Currents and Biogeochemistry: salinity, temperature and current at CTD and ADCP mooring FOR01 from 2004 ongoing. Obtained from knb-lter-mcr.30.31
.
White
,
C. R.
,
Kearney
,
M. R.
,
Matthews
,
P. G. D.
,
Kooijman
,
S. A. L. M.
and
Marshall
,
D. J.
(
2011
).
A manipulative test of competing theories for metabolic scaling
.
Am. Nat.
178
,
746
-
754
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information