An Island Mass Effect Resolved Near Mo’orea, French Polynesia

We sought to resolve the extent, variability, and magnitude of productivity enrichment around a high tropical island consistent with the phenomenon of an Island Mass Effect (IME). Key biogeochemical constituents and physical oceanographic parameters were measured offshore over the upper 500 m from July 27 to August 7, 2014 around the Society Island of Mo’orea in French Polynesia in association with the nearshore measurements of the Mo’orea Coral Reef Long Term Ecological Research program. High-resolution synoptic sampling in a rectangular grid around the island revealed vertical and horizontal patterns in hydrographic conditions, inorganic nutrients, rates of productivity, and concentrations of organic matter that are characteristic of oligotrophic gyre ecosystems. Within the upper euphotic zone (0–75 m), levels of net primary productivity (NPP), chlorophyll a (Chl), heterotrophic bacterioplankton productivity (BP), and particulate organic carbon (POC) exhibited concurrent enhancement at stations located within 5–15 km of shore, relative to stations farther offshore. These observations of enhanced productivity near an island are consistent with an IME. Particulate organic matter nitrogen isotopes (POM-δ15N) were significantly lower near the island than at stations farther offshore, further emphasizing spatial differences in water column biogeochemistry consistent with an IME. Vertical profiles suggest thermocline shoaling and mixing associated with the pycnocline impinging on the island’s submerged flanks and coral reef slope may have contributed to the decreasing depth and increasing intensity of chlorophyll-a concentration in the DCM at nearshore stations relative to farther offshore. Shipboard measurements of an anticyclonic flow within the upper 75 m of the water column in the vicinity of Mo’orea suggest that retention of inorganic nutrients and organic matter near Mo’orea may also have contributed to the patterns in NPP, Chl, BP, POC, and POM-δ15N, providing a potential mechanistic understanding of the processes driving an IME.

We sought to resolve the extent, variability, and magnitude of productivity enrichment around a high tropical island consistent with the phenomenon of an Island Mass Effect (IME). Key biogeochemical constituents and physical oceanographic parameters were measured offshore over the upper 500 m from July 27 to August 7, 2014 around the Society Island of Mo'orea in French Polynesia in association with the nearshore measurements of the Mo'orea Coral Reef Long Term Ecological Research program. High-resolution synoptic sampling in a rectangular grid around the island revealed vertical and horizontal patterns in hydrographic conditions, inorganic nutrients, rates of productivity, and concentrations of organic matter that are characteristic of oligotrophic gyre ecosystems. Within the upper euphotic zone (0-75 m), levels of net primary productivity (NPP), chlorophyll a (Chl), heterotrophic bacterioplankton productivity (BP), and particulate organic carbon (POC) exhibited concurrent enhancement at stations located within 5-15 km of shore, relative to stations farther offshore. These observations of enhanced productivity near an island are consistent with an IME. Particulate organic matter nitrogen isotopes (POM-δ 15 N) were significantly lower near the island than at stations farther offshore, further emphasizing spatial differences in water column biogeochemistry consistent with an IME. Vertical profiles suggest thermocline shoaling and mixing associated with the pycnocline impinging on the island's submerged flanks and coral reef slope may have contributed to the decreasing depth and increasing intensity of chlorophyll-a concentration in the DCM at nearshore stations relative to farther offshore. Shipboard measurements of an anticyclonic flow within the upper 75 m of the water column in the vicinity of Mo'orea suggest that retention of inorganic nutrients and organic matter near Mo'orea may also have contributed to the patterns in NPP, Chl, BP, POC, and POM-δ 15 N, providing a potential mechanistic understanding of the processes driving an IME.

INTRODUCTION
High productivity and biomass observed on coral reefs surrounding islands such as Mo'orea, French Polynesia, are thought to be supported by the transport and retention of key, waterborne biogeochemical constituents including organic and inorganic nutrients and biological particles (Odum and Odum, 1955;Johannes et al., 1972;Leichter et al., 2013). Mo'orea, located in the westward flowing South Equatorial Current, is surrounded by nutrient-depleted waters that form the northern portion of the South Pacific Subtropical Gyre (Rougerie and Rancher, 1994). Surface waters in the South Pacific Subtropical Gyre are among the least biologically productive waters in the global oceans (Longhurst et al., 1995). However, the shallow waters near Mo'orea are distinct from offshore waters (Nelson et al., 2011;Leichter et al., 2013), indicating a potential for islands to alter the chemical and biological composition of water in their vicinity through a phenomenon known as the Island Mass Effect (IME). The IME was first described over 60 years ago as a mechanistic explanation of enhanced phytoplankton biomass observed near island-reef ecosystems (Doty and Oguri, 1956). IMEs have since been documented for numerous islandreef ecosystems, and the mechanisms that promote enhanced productivity have been shown to be ecologically important in coral reef development and persistence (e.g., Signorini et al., 1999;Palacios, 2002;Andrade et al., 2014;Gove et al., 2016 and citations therein).
Mechanisms contributing to enhanced production and the modification of water column biogeochemical characteristics near islands in otherwise stable oceanic environments include both inputs from the island terrestrial environments, and a suite of oceanographic processes contributing to local mixing, transport, and retention of water masses near islands (Gilmartin and Revelante, 1974;Heywood et al., 1996;Palacios, 2002). Transport of terrigenous material through riverine outflow and both surface and groundwater runoff can increase concentrations of nutrients and organic material nearshore (Anderson et al., 2002;Garrison et al., 2003;Knee et al., 2016). Local human impacts such as urban development and agricultural land-use changes, as well as waste water deposition may also enhance nutrient and organic matter inputs to nearshore waters (e.g., Anderson et al., 2002;Gove et al., 2016). In addition to direct inputs from land, multiple oceanographic mechanisms can increase nutrient inputs to waters surrounding oceanic islands. Coral reef organisms are known to biogeochemically transform ocean water flowing over reefs (e.g., Wyatt et al., 2010Wyatt et al., , 2012Nelson et al., 2011). Current wakes formed through the interaction of oceanic currents with island bathymetry, may form eddies and induce upwelling of cool, deep, nutrient rich water (e.g., Takahashi et al., 1981;Wolanski and Hamner, 1988). Island wakes can also influence particle accumulation through partial trapping of water near islands (Rissik et al., 1997). Internal waves generated along and impinging on the submerged slopes of islands result in lateral and vertical transport of deep water masses as well as near-bottom mixing (Wolanski and Pickard, 1983;Merrifield et al., 2001), influencing nutrient and particle concentrations and the thermal environments of island reef ecosystems (Leichter et al., 2003(Leichter et al., , 2012Wyatt et al., 2020). Because multiple biophysical drivers work to influence the prevalence and magnitude of IMEs, detailed observations of these processes can enhance understanding of the environmental contexts in which coral reef ecosystems develop and persist around oceanic islands.
Here we consider Mo'orea, French Polynesia as a model system, and provide one of the most comprehensive synoptic oceanographic surveys to date of the waters surrounding an oceanic island. We focus on in situ measurements of key biogeochemical parameters that support the presence of an IME surrounding the island. We find evidence suggesting enhanced productivity and potentially distinct nutrient sources in the vicinity of the island relative to farther offshore. Previous evidence of the prevalence of internal waves around Mo'orea (Leichter et al., 2012) and a net anticyclonic flow along the 15 m isobath surrounding the island (Leichter et al., 2013) indicates the potential for island-altered physical flow patterns that may magnify IMEs through localized upwelling, and the retention and accumulation of laterally advected nutrients and organic material near Mo'orea. Our spatial survey results were consistent with patterns extracted from satellite time series of chlorophyll a, suggesting that IMEs like that observed near Mo'orea may have long-term impacts on biogeochemical processes near island coral reef ecosystems in general.

Site Description
Mo'orea is a high volcanic island and part of an island chain surrounded by barrier reefs within 0.5-1 km of shore. Fringing reefs around the island (∼10 m deep) are separated from the barrier reef by shallow lagoons (<3 m). The lagoons are connected to the open ocean by reef-passes that are present around the island, and spaced approximately 5-10 km apart. The forereef slopes steeply from the reef crest offshore to >500 m within 1-2 km of the reef. Local ocean circulation near Mo'orea is influenced by the South Equatorial Current, which generally flows westward and constitutes the northern portion of the anticyclonic South Pacific Subtropical Gyre (Rougerie and Rancher, 1994). Evidence of net anticyclonic flow at a local scale along the 15 m isobath at Mo'orea has also been observed using current meter records from fixed locations around the island that were vertically averaged and low-pass filtered on 2-week or longer timescales (Leichter et al., 2013).

Sampling Scheme
High resolution hydrographic and biogeochemical data were sampled during a 3-week research cruise from July 27 to August 7, 2014 aboard the R/V Kilo Moana in the waters surrounding the island of Mo'orea, French Polynesia, in the South Pacific Subtropical Gyre (17 • 36 S, 149 • 43 W). The cruise included a sampling grid comprising 40 stations around the island (Figure 1). Stations near Mo'orea, as well as those near Tahiti, are referred to as "nearshore" stations and ranged from roughly 5 to 15 km offshore of the islands. Stations farther offshore ranged from roughly 15 to 30 km from Mo'orea and are referred to as "offshore stations." At each station, Conductivity, Temperature, Depth (CTD) profiles and bottle samples were collected to determine the hydrographic and biogeochemical properties of the water column. Vertical profiles of water samples to measure biogeochemical variables including inorganic nutrients (nitrogen and phosphorous), net primary production (NPP), chlorophyll a (Chl), heterotrophic bacterial productivity (BP), bacterial abundance (BA), dissolved organic carbon (DOC), particulate organic carbon (POC) concentrations, and isotopic composition of POM were collected via 12 L Niskin bottles. Sample bottles were attached to a rosette and CTD profiler (Sea-Bird Electronics SBE 911) which had sensors for oxygen (SBE 43), PAR (Biospherical QSP-2300), Chl fluorescence (Wetlabs ECO FLNTU), an optical nitrate sensor (Satlantic ISUS), and optical transmittance (Wetlabs C-Star 25 cm transmissometer, 650 nm light source). Water samples were collected at ten depths from the surface to 500 m (nominally 10, 25, 50, 75, 100, 150, 200, 250, 300, and 500 m), with the depth of the nominal 100 m bottle adjusted and targeted to the deep chlorophyll fluorescence maximum (DCM). In addition, oceanographic currents were measured using two shipboard broadband ADCP as described below.

CTD Measurements
CTD profile data were processed using standard Sea-Bird software (Sea-Bird Electronics, 2014) and averaged into 1-m depth bins to form water property profiles that extended from the surface to about 500 m at all stations. Mixed layer depth (MLD) was calculated as the depth where density was 0.01 kg m −3 greater than the local near-surface value.

Sample Processing
Inorganic nutrients -Unfiltered seawater samples for inorganic nutrient analysis were collected in 20 mL HDPE vials, frozen at −20 • C, and analyzed using flow injection analysis on a QuickChem 8000 (Lachat Instruments, Zellweger Analytics, Inc.) by the University of California, Santa Barbara Marine Science Institute Analytical Laboratory (detection limits: NO 2 − + NO 3 − , 0.2 µmol L −1 ; PO 4 3− , 0.1 µmol L −1 ). Net primary productivity (NPP) -Primary production was measured using the 14 C method according to USJGOFS protocols (Knap et al., 1996) with the following modifications (detection limit: 0.01 mg m −3 d −1 ). Seawater was collected from six depths between 10 and 150 m (nominally, 10, 25, 50, 75, 100, and 150) into 500 mL acid-washed polycarbonate bottles. One light and one dark bottle (each containing 250 mL seawater) from each depth were incubated on deck for 24-25 h in a 1 m 3 clear plexiglass incubator under continuous flow with seawater at ambient upper MLD temperature and at light levels corresponding to 75%, 50%, 16%, 3%, 0.6%, and 0.02% of full sunlight based on a clear water optical model (Morel et al., 2007) and previous local Chl estimates from the MCR LTER dataset 10 1 . 100 µL of 14 C at a specific activity of 80 uCi mL −1 was added to each bottle. At the end of the incubation 100 µL was pipetted from each bottle to determine the total 14 C present and then the entire bottle was filtered onto 25 mm Whatman GF/F filters and placed into scintillation vials. These were acidified with 0.25 mL of 0.5 M HCl and allowed to sit for 4-6 h to remove inorganic 14 C. 10 mL of scintillation cocktail (Ultima Gold LLT, PerkinElmer) was then added to each vial and the radioactivity was counted on a Packard TriCarb Liquid Scintillation Counter. An average alkalinity value of 23611 mg C m −3 (standard deviation 0.006%) determined over four seasons from 2005 to 2012 at an offshore station 5 km north of Mo'orea by previous LTER research 2 was used to calculate the carbon fixed. Gradients in alkalinity across our grid were not expected and slight variations in alkalinity would have insignificant impacts on NPP comparisons among stations.
Chlorophyll -Chlorophyll a (Chl) concentrations were determined by filtering 1035 mL of seawater from each depth and station through a glass fiber filter (GF/F) and stored at −80 • C (detection limit: 0.01 µg L −1 ). Samples were analyzed within 3 weeks of collection by extraction in 90% acetone in the dark at −4 • C for 36 h and measured with a fluorometer (Turner 10AU) according to Smith et al. (1981). While GF/F filters may not retain all of the smallest picoplankton, JGOFS protocols were followed allowing comparisons across stations to be internally consistent.
Heterotrophic bacterioplankton productivity -Samples for bacterioplankton production (BP) were measured via 3H-leucine incorporation (14-19 nmol L −1 H-leucine; specific activity 54.1 Ci mmol −1 ; PerkinElmer, Boston, MA, United States), using a modified microcentrifuge method (Halewood et al., 2012;Wear et al., 2015). All samples were processed at sea with incubations conducted using water drawn directly into foam coolers (up to 4) and samples incubated in the dark within 2 • C of in situ temperatures for 2-3 h. Radioactivity was analyzed using a liquid scintillation counter (Packard TriCarb) and corrected with external standard and quench curve (detection limit: roughly 0.2 pmol L −1 hr −1 ). Leucine incorporation was converted to bacterial carbon production by applying a carbon conversion factor of 1.5 kg C mol Leu −1 (Simon and Azam, 1989).
Bacterioplankton abundance -Samples for bacterioplankton abundance were estimated with flow cytometry (BD Bioscience LSR II with SYBR Green I) according to sample preparation and instrument settings described in Nelson et al. (2011). Flow cytometry analysis enumerates total prokaryotic abundance. Thus, we were unable to differentiate between bacterial and archaeal domains and refer to the combined cell densities as bacterioplankton abundance (Glöckner et al., 1999).
Dissolved organic carbon -DOC samples were filtered through combusted glass fiber filters (GF/F) placed in polycarbonate cartridges and attached directly to the Niskin bottle via silicone tubing. Samples were collected into 60 mL high-density polyethylene bottles and frozen at −20 • C until analysis on shore at UC Santa Barbara. Samples were analyzed via high temperature combustion method on a modified Shimadzu TOC-V or Shimadzu TOC-L using the standardization and referencing approaches described in Carlson et al. (2010); detection limit roughly 1 µmol C L −1 .
Particulate organic carbon -Positive pressure displacement was used to concentrate POC from 4 L of seawater onto precombusted GF/F (25 mm) placed in high density polyethylene filter holders. Sample bottles were pressurized at low pressure (approximately 4-7 PSI of nitrogen gas). Immediately following filtration, filters were placed into pre-combusted glass scintillation vials and stored at −20 • C until prepped for analysis at UC Santa Barbara. To prepare samples for analysis, glass vials containing filters were removed from −20 • C storage and placed in an acid fumer (deep covered glass dish containing a 25 mL beaker of concentrated hydrochloric acid in center) and allowed to sit overnight. The following day, vials were loosely capped and placed in a 65 • C drying oven for 48 h. When dry, the vials were removed from the oven, capped tightly, and submitted to the Marine Science Institute analytics laboratory for analysis. POC was quantified by combustion elemental analysis (CE440 elemental analyzer, Exeter Analytical Inc.) with a detection limit of 2-7 µg C.
Particulate organic matter stable isotope compositions -Samples for the determination of stable isotope composition of POM were collected from Niskin bottles at depths of 10 m (near the surface), 50 m, and in the DCM layer (typically 90-110 m as described above). For each sample, 8L of seawater were filtered as described above through pre-combusted GF/F (25 mm) filters. Isotope ratios were measured on the preserved filters using a continuous flow Elemental Analyzer-Isotope Ratio Mass Spectrometer (EA-IRMS) system consisting of a Flash 2000 EA coupled to a Delta V Advantage IRMS via a Conflo IV (Thermo Fisher Scientific, Inc., Bremen, Germany). Dried filters were acidified to remove inorganic carbon by dropwise addition of 1 mol L −1 HCl before being re-dried at 60 • C for 24 h. Redried filters were then combusted to N 2 and CO 2 in tin capsules (Säntis Analytical AG, Teufen, Switzerland) using the EA, and the gases were purified by chromatography and injected into the IRMS. N isotope ratios (δ 15 N) are reported in parts per thousand ( ) relative to air and C isotope ratios (δ 13 C) relative to standard Pee Dee Belemnite (V-PDB). Isotope values were normalized to an international scale using a working standard, L-Alanine (SI Science Co., Ltd., AZ101, Lot: SS13; δ 13 C = −19.6 , δ 15 N = 13.7 ), regularly spaced (every 2-3 samples) and analyzed across a range of C and N masses encompassing the range of C and N in the filters for which the precision was 0.05 for δ 13 C and 0.1 for δ 15 N. As described elsewhere, the working standard was periodically checked in-house against primary international standards (e.g., USGS) using dual-inlet IRMS (Wyatt et al., 2019). Precision (1 SD) based on repeated analysis of two different working standards treated as unknowns was always better than 0.05 for δ 13 C and 0.1 for δ 15 N. Current measurements -Horizontal currents were measured using two shipboard broadband ADCPs operating simultaneously at 300 kHz (Workhorse) and 38 kHz (Ocean Surveyor, Teledyne RD Instruments). Velocity profiles from the 300 kHz ADCP were collected in 2-m bins with a ping interval of 0.8 s and 150 pings were ensemble-averaged over 120 s. Ship motion was removed from the ADCP velocity data using an integrated inertial navigation and GPS system. Velocity errors for both ADCPs were about 0.01 m s −1 . Shipboard ADCP processing was conducted according to the recommendation of Firing and Hummon (2010). The shallowest bin that consistently provided usable data was centered at 11 m. Low concentrations of scattering particles limited the range of the 300 kHz ADCP to about 75 m depth. Velocity profiles from the 38 kHz ADCP were collected in 12-m bins with a ping interval of 4 s and 75 pings were ensemble-averaged over 300 s. The shallowest bin that consistently yielded data was centered at 37 m and good data extended down to about 1100 m.
To visualize current patterns and compare with sampling on the grid stations, the ADCP data were interpolated onto a uniform grid with 4 km spacing using optimal interpolation (e.g., Thompson and Emery, 2014). In the optimal interpolation, a spatial autocorrelation function was calculated from the data and fitted to a Gaussian curve with a standard deviation of 15 km.

Satellite Chl Observations
Subsurface Chl concentrations from the NASA Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua sensor (level-3 binned operational Chl data product, R2018.0, approximately 4.6 km resolution) were extracted for the Mo'orea-Tahiti region covered by the cruise (Figure 1). Data from nine 40 to 50 km long transects distributed around Mo'orea (6) and Tahiti (3), starting at least one-half pixel away from the reef crest and running perpendicular to it were extracted. Data from the transects were interpolated to a common ∼ 3.25 km grid (equal to half the length of a pixel diagonal, resulting in each transect data point representing an approximately 10.6 km 2 surface area) so the data could be aggregated from the 9 transects as a function of the distance from the reef crest (Gove et al., 2016). Data for two different time scales were analyzed: (1) A monthly composite dataset matching the time of the cruise (July and August 2014), and (2) A long-term average that covers most of the current MODIS-Aqua mission (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017).

Data and Statistical Analyses
Statistical analysis of the difference in biogeochemical rates and concentrations between nearshore and offshore stations (Figure 1) was evaluated using a one-way analysis of variance (ANOVA; Scheffe, 1959). Ranged Major Axis (RMA) model II simple linear regression was conducted to assess relationships between parameters according to Legendre and Legendre (2012). Statistical analyses were conducted on log 10 -transformed values.

Vertical Variability
We compared vertical variability in hydrographic conditions and biogeochemical variables by comparing direct measurements of temperature, salinity, rates of productivity, and concentrations of inorganic nutrients and organic matter across the upper 300 m of the water column in vertical profile plots. In this way, we were able to evaluate vertical patterns within distinct depth layers including the upper euphotic zone (0-75 m), the lower euphotic zone which contained the deep chlorophyll maximum (75-150 m), and the upper mesopelagic zone (150-300 m).

Hydrographic Conditions
Sampling near Mo'orea during July 2014 (austral winter) revealed characteristics of stratification in the upper 300 m of the water column typical of Pacific subtropical gyre ecosystems Church, 2014, 2017;Wilson et al., 2015). Relatively warm surface temperatures (mean, 26.8 ± standard deviation 0.2 • C, range: 25.9-27.2 • C) and low salinity (mean: 35.7 ± 0.1, range: 35.6-36.1) resulted in a low-density layer of water that extended from the surface to approximately 75 m and that separated the upper euphotic zone from the remainder of the water column (Supplementary Figure S1). In the lower euphotic zone, between 75 and 150 m (approximately 0.1% light level), mean temperatures decreased to 25.1 ± 1 • C and mean salinity increased to 36.0 ± 0.2, leading to the depth of maximum vertical density stratification occurring between approximately 75 and 150 m across stations. These observations of stratification are similar to those previously observed near Mo'orea during the austral summer (Leichter et al., 2012). The MLD, defined as above, was variable across stations, ranging from 5 to 99 m with a mean and standard deviation of 54 ± 24 m.  Figure S2). This observation of separation between the upper portion of the nutricline and the MLD is consistent with a previous study near Mo'orea during the austral summer (Leichter et al., 2012) and indicates that vertical mixing at the time of sampling was not deep enough to entrain nutrients from the nutricline into the upper euphotic zone.

Net Primary Productivity and Chlorophyll a
Rates of NPP were highest and most variable in the upper euphotic zone, ranging from 0 to 1.5 µmol C L −1 d −1 with a mean value of 0.3 ± 0.3 µmol C L −1 d −1 (Figure 2A). The depth of maximum NPP was consistently shallower than the depth of the upper portion of the nutricline and was approximately 65-90 m shallower than the DCM (compare Figures 2A,B). This observation is consistent with oligotrophic gyre ecosystems in the North Pacific and North Atlantic, where the DCM is generally deeper in the water column than maximum rates of NPP due to photoadaptation of phytoplankton cells over depth (Steinberg et al., 2001;Karl and Church, 2014;Cullen, 2015). Rates of NPP decreased by an order of magnitude in the lower euphotic zone and were below limits of detection below approximately 120 m for all stations (Figure 2A).
In contrast to NPP, mean concentrations of Chl were approximately 33% lower in the upper euphotic zone (0.12 ± 0.06 µg Chl L −1 ), compared with the lower euphotic zone (0.18 ± 0.09 µg Chl L −1 ; Figure 2B). The DCM was consistently present in the lower euphotic zone between the 23.5 and 24.5 kg m −3 potential density isopycnals, and was bound by sampling depths of approximately 75-150 m ( Figure 2B). The depth of maximum concentrations of Chl varied by station, ranging from 74 to 140 m, with a mean depth of 107 ± 16 m. Six stations exhibited mean or above mean concentrations of Chl throughout the upper 75 m of the water column ( Figure 2B). All of these stations exhibited MLDs that were shallower (≤76 m, Figure 2B) than the depth of maximum Chl concentration, and in most cases the upper portion of the DCM, suggesting that vertical mixing at the time of sampling was not sufficiently deep to entrain high concentrations of Chl from the DCM into the upper euphotic zone. We also note the depth of the DCM in Figure 2 was often shallower at stations near the island than stations farther offshore. This appears to be consistent with shoaling of the density stratification near the island, and with nutrient mixing associated with interactions between the pycnocline and the submerged flanks of the island.

Bacterioplankton Dynamics
Rates of BP and concentrations of BA decreased over the upper 300 m of the water column and ranged from 0 to 0.1 µmol C L −1 d −1 , and 0.6-8.5 × 10 5 cells mL −1 , respectively. To evaluate the contribution of bacterioplankton biomass to carbon dynamics, we converted BA to bacterioplankton carbon (BC) using a carbon conversion factor of 10 fg C cell −1 to reflect small cells in oligotrophic regions (Caron et al., 1995;Fukuda et al., 1998;Gundersen et al., 2002). Concentrations of BC, and rates of BP, were elevated in the upper and lower euphotic zones and substantially decreased in the upper mesopelagic (Figures 2C,D).
Rates of BP were substantially reduced and demonstrated low variability in the upper mesopelagic ( Figure 2D). In contrast, concentrations of BC were variable and exhibited elevated concentrations that extended throughout the upper mesopelagic at stations offshore, as well as near Mo'orea ( Figure 2C). These elevated concentrations of BC in the upper mesopelagic may be indicative of previous vertical export of organic material and subsequent utilization of this material by extant bacterioplankton communities (e.g., Hansell and Carlson, 2001;Letscher et al., 2015).

Organic Carbon
Concentrations of POC and DOC were highest in the upper and lower euphotic zones and decreased in the mesopelagic (Figures 2E,F). Concentrations of POC in the upper euphotic zone were low (<3 µmol C L −1 ), and similar in amount to those observed in the oligotrophic waters of the North  Pacific (Karl and Church, 2014;Wilson et al., 2015) and North Atlantic Subtropical Gyres (Lomas et al., 2013). Concentrations of DOC in the upper euphotic zone reflected the accumulation of DOC within surface waters of subtropical gyres (Hansell et al., 2009) and ranged from 66 to 77 µmol C L −1 , consistent with previously reported ranges from these waters (Nelson et al., 2011;Carlson and Hansell, 2015).
Despite lower concentrations of POC and DOC in the upper mesopelagic, several stations exhibited elevated concentrations of POC ( Figure 2E) and DOC ( Figure 2F) similar to those observed for concentrations of BC ( Figure 2C). In contrast to BC, comparison of offshore and nearshore stations in the upper mesopelagic revealed significantly elevated POC (p value = 0.02) and DOC (p value = 0.03) nearshore (one-way ANOVA; see red box indicating nearshore stations in Figures 2E,F). Additionally, statistical analysis revealed a positive correlation between logtransformed mean POC and DOC concentrations in the upper mesopelagic (integrated and depth normalized; p value < 0.01, R 2 adjusted 0.18, RMA), indicating coincident enhancement of particulate and dissolved organic carbon within this depth layer. These coincident elevated concentrations extended through isopycnal surfaces and may have occurred in the absence of deep mixing. Though we lack measurements of the amount of POC exported from the upper and lower euphotic zone, as well as sinking rates into the upper mesopelagic, solubilization of POC to DOC (Smith et al., 1992) may have partially contributed to contemporaneous enhancement of organic carbon in the upper mesopelagic. Submesoscale dynamics arising from eddy interactions may also help to explain some of the elevated pulses of DOC, POC observed in the mesopelagic at depths well below the MLD (Mahadevan, 2016;Erickson and Thompson, 2018).
There were two clear patterns in the POM stable isotope data. First, the carbon and nitrogen isotope compositions were distinct in the DCM relative to the shallower layers (Figure 3). In the DCM mean δ 13 C was lower (−23.90 ± 0.09 %) and mean δ 15 N higher (7.19 ± 0.08 %) than in the surface (−21.50 ± 0.06 % and 5.87 ± 0.10 %, respectively) and mid layers (−21.60 ± 0.06 % and 5.85 ± 0.10 %, respectively). Second, for the surface and mid layers there was evidence of significantly lower POM-δ 15 N values at the stations near the island relative to stations farther offshore, e.g., 5.58 ± 0.09 % at near stations compared to 6.23 ± 0.25 % at far stations, respectively (p < 0.01, one-way ANOVA; Figure 3). A variety of factors including temperature and irradiance can lead to differential C fractionation in marine phytoplankton (Hinga et al., 1994). The differences in δ 13 C with depth could reflect light availability with samples in the DCM likely containing phytoplankton growing in the lower euphotic zone, whereas surface and mid layers may have had phytoplankton growing at higher irradiance. The higher δ 15 N in the DCM relative to shallower depths could represent the role of deeper water nitrate supporting primary production in the DCM. Bottle samples showed that nitrate concentrations increased markedly below the pycnocline and DCM and deep-water nitrate in oceanic environments including the central South Pacific Gyre would be expected to have elevated δ 15 N relative to shallow water nutrients (Sigman et al., 1997;Altabet, 2001). Lower POM-δ 15 N in the surface layers could also reflect higher assimilation of nitrogen from nitrate (fractionation ∼5 %) relative to ammonium (fractionation ≤ 20 %), or a greater contribution from nitrogen fixation (with expected δ 15 N ∼ 0 %) (Sigman et al., 2009). The significant differences in δ 15 N values between near and far stations at the surface and mid depths support the observations of spatial differences in biogeochemical properties suggestive of an IME. We can speculate regarding the causes of the direction of the differences. If there were a strong signal of terrestrial nutrients including agricultural or waste water runoff we would expect elevated δ 15 N values in nearshore samples, however, that is not the direction of the observed differences. Another possibility is that the lower δ 15 N values observed in stations closer to the island reflect organic material with a nutrient source from nitrogen fixation on the reef, and/or trophic alteration of the POM by reef organisms. Coral reef communities also contain exceptional abundance of suspension feeding organisms which are known to consume large quantities of zooplankton from overlying waters (Hamner et al., 2007;Wyatt et al., 2013). A result of extensive consumption of zooplankton by the reef communities might be that the remaining POM pool could represent material with an overall lower trophic position and average δ 15 N than the zooplankton. Again, while these potential mechanistic explanations would require additional research to confirm or refute, the spatial differences in isotope compositions between nearshore and offshore stations is clearly suggestive of an overall influence of island proximity reflecting an IME.
Estimates of apparent oxygen utilization (AOU) account for the effects of temperature and salinity on dissolved oxygen concentrations, and reflect remineralization of organic carbon and the removal of dissolved oxygen by respiration (Doval and Hansell, 2000;Carlson et al., 2010). To estimate the contribution of DOC remineralization to AOU, we converted AOU to carbon equivalents (AOU-C eq ) using a molar ratio of − C: O 2 of 0.72 (Anderson, 1995). Comparing direct measurements of logtransformed concentrations of DOC in the upper mesopelagic to direct estimates of log-transformed AOU, we found DOC remineralization accounted for approximately 18% of AOU within the upper mesopelagic. This observation is at the lower end of previous reports of the contribution of DOC oxidation to AOU in the upper mesopelagic of subtropical gyre ecosystems (e.g., Carlson et al., 2010;Pan et al., 2014). The remaining AOU signal is likely due to remineralization of sinking particle flux over an unknown period of time.
Upper Euphotic Zone (0-75 m) Variability and Evidence of an IME To assess the presence of an IME, we compared vertical profiles of rates of productivity, as well as concentrations and , and dissolved organic carbon (DOC; F). Dark blue points indicated means ± SD for offshore stations, while red represents means ± SD for stations near Mo'orea. Asterisks denote significant differences (p values < 0.05) in log(10)-transformed measured biogeochemical rates and concentrations as a function of station location (offshore versus nearshore). Statistical analysis was evaluated using a one-way analysis of variance (ANOVA).
compositions of organic matter within the upper euphotic zone (0-75 m; Figures 3, 4). Stations were grouped by their relative proximity to island shores (Figure 1) into offshore (# 1-24) and nearshore (# 25-41) stations. Statistical differences between the two groups were analyzed using log-transformed values of each measured parameter.
In the upper euphotic zone, grouping stations by proximity to Mo'orea ("nearshore") revealed coincident horizontal variability in NPP, Chl, BP, POC, and POM-δ 15 N (0-75 m; Figures 3B, 4A,B,D,E and Supplementary Table S1). Collectively, nearshore stations exhibited significantly elevated log-transformed mean concentrations of Chl (p value < 0.01) and POC (p value < 0.01), as well as lower values of POMδ 15 N (p value < 0.01) throughout the upper euphotic zone (0-75 m; Figures 4B,E), compared with offshore stations (Supplementary Table S1). In addition, nearshore stations exhibited elevated log-transformed mean rates of NPP (p value < 0.01) and BP (p value ≤ 0.01) over the upper 50 m of the water column, compared with offshore stations (Figures 4A,D and Supplementary Table S1). These observations in the upper euphotic zone indicate that key biogeochemical constituents were enhanced at nearshore stations, compared with offshore Black vectors show optimally interpolated currents from a shipboard ADCP during 28 July-9 August 2014 as described in the text. Red current sticks near Mo'orea show vertically averaged currents from three bottom-mounted ADCPs on the 15 m isobath during 28 July-9 August 2014. Velocity scales at lower left correspond to ADCP vectors (black arrows) and moored currents (red sticks). Gray line shows ship track during ADCP data collection. stations, during our sampling period and may be indicative of an IME that resulted in distinct biogeochemical characteristics near the island (e.g., Takahashi et al., 1981;Xie et al., 2001;Caldeira et al., 2005;Gove et al., 2016 and citations therein).
An IME may have contributed to elevated rates of productivity, as well as concentrations of POC and Chl, through the introduction of limiting nutrients to stations nearshore, thus stimulating productivity and production of organic matter. Enhanced rates of BP in the upper euphotic zone may indicate a rapid response of bacterioplankton to this freshly produced organic matter. However, in contrast to previous observations of upwelling induced IMEs near deep-sea islands (Furuya et al., 1986;Heywood et al., 1990;Hasegawa et al., 2004), our observations of MLDs suggest that vertical mixing was not sufficiently deep to entrain nutrients from the upper portion of the nutricline into the upper euphotic zone during our sampling period (see MLD in white line in Supplementary Figure S2).
Observations by satellite remote sensing, shipboard ADCPs, and nearshore moorings were all consistent with a general anticyclonic flow pattern around Mo'orea. Satellite altimetry data 3 averaged over the sampling period revealed generally positive sea surface height anomalies (SSHAs) at large scales around Mo'orea. Geostrophic flow vectors derived from the SSHAs indicated the predominance of anticyclonic flow in the region (Figure 5A). A few areas of cyclonic flow were also observed at these scales 200-300 km southeast and southwest of Mo'orea. Shipboard ADCP measurements in the upper ocean (averaged over depths ≤75 m) indicated a general pattern of anticyclonic flow at smaller scales around Mo'orea, consistent with the SSHA observations ( Figure 5B). It should be noted that the ADCP data were collected over 2 weeks and changes in the current field around Mo'orea likely occurred over the measurement period. Despite this, the orientations of most current vectors measured by the ADCPs suggested an anticyclonic flow pattern. Nearshore moored current observations averaged over the 2-week sampling period (vertically averaged in 15-m water depths) also indicated anticyclonic flow around Mo'orea ( Figure 5B). Anticyclonic flow at mesoscales results in elevated sea surface height and downward displacement of the pycnocline, often depressing the nutricline and reducing nutrient supply to the upper euphotic zone (McGillicuddy et al., 1999). Consequently, the observed IME in the upper euphotic zone around Mo'orea is most likely not driven by nutrient supply from below the nutricline. We therefore suggest that lateral transport of inorganic nutrients originating from the islands, rather than vertical transport via mixing, stimulated productivity at these stations. Observations of offshore flow along western Tahiti, as well as anticyclonic flow around Mo'orea ( Figure 5B), during the sampling period may have led to the introduction of nutrients near Mo'orea via mechanisms such as submarine groundwater discharge, local human impacts, atmospheric deposition, or previous mixing induced by upwelling near the island (e.g., Furuya et al., 1986;Hasegawa et al., 2004). Furthermore, entrainment of inorganic 3 http://marine.copernicus.eu/ nutrients and organic matter in the anticyclonic flow around Mo'orea could lead to retention of this material near the island.

Satellite Chl Observations
Enhanced productivity near the island is also evidenced from the MODIS Chl imagery. The satellite data show higher Chl concentrations in the close vicinity of both Mo'orea and Tahiti compared to offshore waters over the time period July-August 2014 ( Figure 6A). On average, Chl values within ∼10-12 km of the reef crest are higher than in more offshore waters (Figure 6B). This is true for the two time-scales that were tested in this study: July to August 2014 and 2003-2017 ( Figure 6B). The enhanced Chl pattern and its magnitude as observed around Mo'orea and Tahiti agree well with those presented in Gove et al. (2016) for other islands in the tropical Pacific. The fact that increased Chl close to the islands is observed at very different time scales suggests that this is a persistent pattern in the Mo'orea and Tahiti island system.

Summary
Overall, our observations of physical and biogeochemical variables in the upper 300 m of the water column near Mo'orea provide evidence of enhancement of key biogeochemical variables (i.e., NPP, Chl, BP, and POC) in the upper euphotic zone within 5 to 15 km of shore, compared with farther offshore. This, along with distinct physical dynamics and the presence of POM with distinct isotopic characteristics at stations near Mo'orea, provides evidence of an IME which remote sensing data also appear to support. While the work presented here was primarily observational in nature, it does lend support to prior considerations of the mechanisms supporting IMEs. Mechanisms contributing to the IME surrounding Mo'orea likely include inputs to the marine environment from the relatively high (∼1,200 m highest elevation) volcanic island, potentially including both surface and groundwater run off (Knee et al., 2016), as well as from alteration of oceanic water by the island's reefs and oceanic sources associated with internal waves (Wolanski and Delesalle, 1995;Leichter et al., 2012) and islandinduced mixing and potential retention processes (Leichter et al., 2013). The presence of an IME likely augments the spatial variability in key biogeochemical variables in an otherwise stable oceanic environment and ultimately enhances the role of vertical and lateral transport along with retention of inorganic nutrients and organic matter in supporting the high biomass and productivity on oceanic island coral reefs such as Mo'orea.

AUTHOR CONTRIBUTIONS
All authors designed the research, conducted the research, and contributed to writing the manuscript.