Sampling locations and dates:
Reefs on the southern shore of St. John, United States Virgin Islands (USVI) were targeted for 11 opportunistic sampling events over six years from June 2016 to June 2022. The reefs included (from west to east) Dittlif, Cocoloba, Joel’s Shoal, Europa, Yawzi, Tektite, Booby Rock, and Ram Head, all of which are within the bounds of Virgin Islands National Park, except for Dittlif. The research operations were coastal, predominantly land-based, lasting 1-3 weeks, with seawater sampling and diver-based reef surveys deployed from a power boat. Reef collections included surveys for benthic composition, small volumes of seawater for inorganic and organic nutrients and microbial abundances, larger volumes for microbial biomass and chlorophyll, and CTD casts for temperature and salinity. All collections occurred during daylight hours. Dates of seawater sampling events included June 10-12, 2016, October 28-29, 2016, March 25-28, 2017, July 26-30, 2017, November 27-30, 2017 (occurred aboard R/V Walton Smith), April 11-13, 2018, November 5-9, 2018 (occurred aboard R/V Walton Smith), August 6-10, 2020, January 17-24, 2021, October 20-25, 2021, and June 24-29, 2022.

These sampling points surround two major stressors to St. John, USVI reefs: two category 5 hurricanes, Irma and Maria, which impacted the reefs in September, 2017, and stony coral tissue loss disease, a multi-species disease outbreak that began emerging around St. John, USVI between January-June, 2020. As of August 2020, the disease just began affecting all reefs in the study area, except Europa and Cocoloba. All reefs were impacted by the next sampling in January 2021.

Field collections:
Surveys were conducted via SCUBA. Point intercept benthic reef surveys were conducted to understand the percent cover of organisms and substrates on each reef. Surveys proceeded using four to six 10-meter (m) long transects, with the underlying biological reef organisms or substrate recorded every 10 centimeters (cm). Transects were collected yearly prior to the 2017 hurricane, then were collected at each field sampling event. Following the emergence of stony coral tissue loss disease in 2020, roving diver surveys to quantify that disease were conducted in January and October 2021. Disease surveys lasted thirty minutes within a 100 m² plot and the diver counted all apparently healthy and diseased coral colonies and recorded species level, when possible. At the conclusion of the survey, the area surveyed within the plot was estimated.

At each reef, measurements of temperature (ºC) and salinity (psu) were conducted using a CastAway CTD (SonTek, Xylem, San Diego, CA, USA). Next, small volumes of water were collected from both surface and approximately within 1 m of reef depth (referred to as "benthic" depth) for organic nutrients, inorganic nutrients, and microbial cell abundances. Seawater for total organic carbon (TOC) and total nitrogen (TN) was collected into 40-milliliter (ml) combusted borosilicate glass vials. Seawater for inorganic nutrients (phosphate, ammonium, silicate, nitrite plus nitrate), was collected in acid-clean 30 ml HDPE bottles (Nalgene, ThermoFisher Scientific, Waltham, MA, USA). Bottles for all collections were triple-rinsed with seawater prior to collection. Following collection, organic nutrient samples were fixed with 75 microliters (μl) phosphoric acid and stored at room temperature until analysis. A 1.4 ml aliquot from the inorganic nutrient collection was placed into a cryovial for flow cytometry-based analysis for microbial cell abundances, then fixed with paraformaldehyde (1% final concentration, Electron Microscopy Sciences) in the dark for 20 minutes, then frozen in a liquid nitrogen dry shipper and stored at -80º Celsius (C) until analysis. Bottles for inorganic nutrient analysis were kept frozen at -20ºC until analysis. Inorganic and organic nutrients and microbial cell abundances were collected in biological duplicates for the June and October 2016 sampling events, but not for further events. Duplicates from those events were averaged for comparison to all other timepoints, which were sampled in singlicate.

Seawater for microbial biomass and chlorophyll analysis was collected from surface and benthic depths at each reef. For all benthic seawater collections prior to 2019, a groundwater pump (Mini-Monsoon 12V, Proactive Environmental Products, Bradenton, Florida, USA) was used. Beginning in August 2020, an 8L diver-operated Niskin bottle was used to enable more accurate collection of water directly over the reef habitat. Seawater (4L) for chlorophyll (benthic depth only) and for microbial biomass was collected into acid-clean or bleach-clean, then triple seawater-rinsed Platy® water tank bags (Platypus, Cascade Designs, Seattle, WA, USA) or LDPE Nalgene bottles. Samples were kept in a cooler on ice until filtration within 6 hours of collection. Seawater was filtered for both microbial biomass and chlorophyll via peristalsis using a Masterflex L/S peristaltic pump (Cole-Parmer, Vernon Hills, IL, USA) through acid-clean or 10% bleached silicone tubing (L/S, platinum-cured, #96410-24 size, Cole-Parmer) and a 25-millimeter (mm) filter holder (Swinnex-25, Millipore Corporation) with either a 0.2-micrometer (μm) Supor filter (Pall, Port Washington, New York, USA) or GF/F filter. Seawater was filtered 2 liters at a time for technical duplicates. In some cases, the 0.2 μm filter would clog, so less than the 2 liters were filtered. In those cases, the amount filtered was recorded. Prior to November 2018, the 4L chlorophyll sample was filtered through a single GF/F filter. For consistency with other environmental variables (nutrients and microbial abundances), any technical duplicates of chlorophyll (post-2018) were averaged to yield one variable for environmental analyses. All filters were placed in cryovials and frozen in a liquid nitrogen dry shipper and stored at -80ºC upon return to the Woods Hole Oceanographic Institution until analysis.

Laboratory processing:
Samples for organic nutrients were analyzed on a Shimadzu TOC-VCSH TOC analyzer (Hansell and Carlson, 2001), using a TNM-1 module to generate non-purgeable total organic carbon (TOC) and total nitrogen (TN). The inorganic nutrient samples were analyzed at Oregon State University as described previously (Becker et al., 2020), and used a Technicon AutoAnalyzer II (SEAL Analytical) and Alpkem RFA 300 Rapid Flow Analyzer to generate concentrations of ammonium, phosphate, silicate, nitrite, and nitrate plus nitrite. The microbial abundance samples were analyzed as described in Becker et al. (2020) at the University of Hawaii SOEST Flow Cytometry Facility. The facility used a Beckman-Coulter Altra flow cytometer (Beckman Coulter Life Sciences). The seawater samples were stained with Hoechst 33342 DNA stain and excited with 488 nm and UV wavelengths (Campbell and Vaulot, 1993; Monger and Landry, 1993). Signals were collected and then processed in FlowJo software (Tree Star) to generate populations and abundances (cells per milliliter) of Prochlorococcus, Synechococcus, eukaryotic phytoplankton ("picoeukaryotes"), and non-pigmented bacteria and archaea. Non-pigmented bacteria and archaea are mostly heterotrophic and referred to as "heterotrophic microbes" in this study (Monger and Landry, 1993; Marie et al., 1997).

Chlorophyll was extracted with 90% acetone using standard methods (JGOFS, 1996). Briefly, filters were thawed and placed individually in 5 ml of 90% acetone and capped. If the filter was dark, 10 ml of acetone was used. After a 24-hour extraction in the dark at 4°C, the tubes were vortexed and centrifuged and ~3 ml of solvent was measured on an AquaFluor fluorometer at 664 nanometer (nm) (Turner Designs handheld 800446) fitted with a red-sensitive photomultiplier. Before and after analysis, blanks including air, 90% acetone, and a black standard were run. After each measurement, samples were acidified with two drops of 10% hydrochloric acid, and re-measured to assess phaeopigment concentration. Readings volume-corrected and concentrations generated with a standard curve. Chlorophyll from June and October 2016 was assessed using high-performance liquid chromatography.

DNA extraction of duplicate 0.22 μm filters proceeded as described previously (Becker et al., 2020) using a method that combines physical and chemical lysis with column purification (Santoro et al., 2010). In addition to the seawater filters, DNA from 14 blank filters was extracted as a DNA extraction control. Briefly, the filters were subjected to chemical lysis with a sucrose-EDTA and 10% SDS lysis buffer and physical lysis with a 15 bead-beating step. The DNeasy Blood and Tissue Kit (Qiagen) was then used to purify the lysate. Resulting DNA was diluted 1:100 in UV-sterile PCR-grade water in preparation for PCR. Single, barcoded PCR reactions per sample were run to amplify the small subunit (SSU) ribosomal RNA (rRNA) gene of bacteria and archaea using primers 515F and 806R (Apprill et al., 2015; Parada et al., 2016). In addition to samples, genomic DNA from Microbial Mock Community B (even, low concentration), v5.1L, for 16S rRNA Gene Sequencing, HM-782D, was used as a sequencing control and PCR-grade water was used as a negative PCR control. PCR reactions (50 μl) contained the following: 2 μl DNA template, 0.5 μl of GoTaq DNA Polymerase (Promega), 1 μl each of forward and reverse primers at 10 μM, 1 μl of 10 mM deoxynucleoside triphosphate (dNTP) mix (Promega), 5 μl MgCl2, 10 μl GoTaq 5X colorless flexi buffer (Promega), and 29.5 μl nuclease-free water. Reactions proceeded with: 95°C for 2 minutes; 28 cycles of 95°C for 20 seconds, 55°C for 15 seconds, and 72°C for 5 minutes; and finally 72°C for 10 minutes before holding at 4ºC. PCR products were all purified with the QIAquick 96 PCR purification kit (Qiagen) or MinElute PCR Purification Kit (Qiagen). Concentrations of purified barcoded PCR products were measured with the Qubit 2.0 fluorometer. Each barcoded sample was diluted to 1 nanogram per microliter (ng/μl) and pooled. Samples were sequenced across four runs to maximize read output per sample on an Illumina MiSeq 2 x 250 bp sequencing at the Roy J. Carver Biotechnology Center at the University of Illinois at Urbana-Champaign or the BioMicro center at Massachusetts Institute of Technology. All resulting fastq files were downloaded and used for further analysis. Fastq files were also uploaded to the NCBI sequence read archive under BioProject accession number PRJNA936592.

Known Issues/Problems:
In some cases, data are missing. This is noted as either "not collected", which indicates that the sample does not exist, or "not applicable" meaning the data were not relevant to the sampling design.

usvi_time-series.csv