http://lod.bco-dmo.org/id/dataset/878413 eng; USA utf8 dataset Highest level of data collection, from a common set of sensors or instrumentation, usually within the same research project Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. pointOfContact 2022-08-11 ISO 19115-2 Geographic Information - Metadata - Part 2: Extensions for Imagery and Gridded Data ISO 19115-2:2009(E) Results of a study examining how foundation species loss alters multiple ecosystem functions based on community surveys and biogeochemical sampling of tidepools in the Otter Rock Marine Reserve (ORMR) and Marine Garden, Oregon USA from June to August 2019 2022-09-07 publication 2022-09-07 revision Marine Biological Laboratory/Woods Hole Oceanographic Institution Library (MBLWHOI DLA) 2022-09-07 publication https://doi.org/10.26008/1912/bco-dmo.878413.1 Nyssa Silbiger California State University Northridge principalInvestigator Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. publisher Cite this dataset as: Silbiger, N., Fields, J. (2022) Results of a study examining how foundation species loss alters multiple ecosystem functions based on community surveys and biogeochemical sampling of tidepools in the Otter Rock Marine Reserve (ORMR) and Marine Garden, Oregon USA from June to August 2019. Biological and Chemical Oceanography Data Management Office (BCO-DMO). (Version 1) Version Date 2022-09-07 [if applicable, indicate subset used]. doi:10.26008/1912/bco-dmo.878413.1 [access date] Oregon Tide Pools Dataset Description: Methods and Sampling: &lt;p&gt;We selected 32 tide pools within the Oregon Department of Fish and Wildlife’s (ODFW) Otter Rock Marine Reserve (ORMR) and Marine Garden, Oregon USA (44°45'09.1&quot;N 124°03'58.0&quot;W) during the months of June to August 2019. One mussel species (&lt;em&gt;M. califorinanus&lt;/em&gt;) and two surfgrass species (&lt;em&gt;P. scouleri &lt;/em&gt;and &lt;em&gt;P. torreyi&lt;/em&gt;) were present at the site. Of the 32 tide pools, 16 were dominated by the California mussel (&lt;em&gt;M. califorinanus&lt;/em&gt;) and 16 were dominated by surfgrass (&lt;em&gt;Phyllospadix spp&lt;/em&gt;.; 10 pools contained &lt;em&gt;Phyllospadix scouleri &lt;/em&gt;and six pools had both &lt;em&gt;Phyllospadix scouleri &lt;/em&gt;and &lt;em&gt;Phyllospadix torreyi&lt;/em&gt;). Percent cover of foundation species ranged from 45.3 to 98.9% in the mussel pools and 49.5 to 100% in the surfgrass tide pools. Some tide pools had both mussels and surfgrasses; however, the presence of a second foundation species did not exceed 7.2% of the remaining tide pool cover. Tide pools were located in the mid to low intertidal zone ranging from 0.71 meters to 1.77 meters above mean lower-low water. At this tide height, tide pools were isolated for 4-6 hours during summer low tides.&lt;/p&gt; &lt;p&gt;A Before-After-Control-Impact (BACI) design was used to account for changes in ocean chemistry, timing of low tide, and variability within tide pools throughout the experimental period. The BACI design consisted of two 29-day time periods (before and after removal) with control and removal tide pools (Stewart-Oaten et al. 1986) where foundation species were removed from removal tide pools (n = 8 per foundation species) between the before and after removal periods. The before-removal period occurred June to mid-July 2019 and the after-removal period occurred mid-July to August 2019. Tide pools were selected for control or removal groups using a random number generator while accounting for the surface area to percent cover of foundation species. Removal of all foundation species from removal tide pools occurred in mid-July 2019 and tide pools had a one-week recovery period before any surveys or samples were taken. Rhizomes were also removed in surfgrass removal pools due to their ability to alter nutrient cycling (Terrados &amp;amp; Williams 1997). During each time period, we characterized physical parameters (pool volume and tide height), community composition, light, temperature, and biogeochemical fluxes (e.g., dissolved oxygen, pH, nutrients) in each tide pool.&lt;/p&gt; &lt;p&gt;Tide pool physical parameters included tidal height (location within intertidal) and tide pool volume (size of pool). Tide heights for each pool were surveyed with a laser level and stadia rod (DeWalt, Towson, MD, USA). Tide pool volume (V) was determined using a dye method (Pfister 1995) and measured with a SmartSpec3000 spectrophotometer (Bio-Rad Lab, Hercules, CA, USA). Water volume in the tide pools changed by 0 - 20% between the two timepoints due to removal of foundation species. To account for the slight effect of changing biomass on volume, volume was re-measured post-removal and the average between both time periods were used in statistical analyses.&amp;nbsp;&amp;nbsp;&lt;/p&gt; &lt;p&gt;We conducted two rounds of community composition surveys for sessile percent cover and mobile organism counts: three-weeks before removal and one-month post-removal. We temporarily removed seawater from the tide pool and placed a flexible mesh quadrat (Nielsen 2001) with demarcations in 10 x 10 cm squares over the bottom of each pool to survey the entire community (pools ranged from 58 to 754 squares). We measured percent cover for sessile organisms and counted mobile organisms, identifying down to the lowest possible taxonomic unit in the field (usually genus level). We normalized the sum of non-foundation species sessile cover to 100% cover per tide pool, including the opposite foundation species if present (e.g., mussel cover within surfgrass pools was also normalized to 100% cover). Each foundation species cover in their respective tide pools remained the raw percent cover and did not exceed 100%. Both sessile and mobile organisms were grouped into larger functional groups based on their ecological role for data visualizations. All tide pool characterization and community composition methods were completed at least 24 hours before any water sampling event to allow the pools to be flushed at least twice by the ocean before measurements.&lt;/p&gt; &lt;p&gt;Temperature (°C) and light intensity (lumens m-2) were recorded continuously every 15 minutes for 29 days during each time period using HOBO® Pendant loggers bolted facing up on the flattest part of the tide pool on a level platform in the interstitial spaces of the foundation species (Onset® HOBO® Pendent Light Intensity Data Logger MX2202, Bourne, MA, USA). Light intensity (lumens m-2) was converted to photon flux density (PFD; µmol photons m-2 s-1) following Long et al. (2012) field experiment values. We measured the change in maximum temperature between the before and after removal period (i.e., the average of the daily hottest temperature from 7/18/2019 – 8/16/2019 minus the average daily hottest temperature from 6/16/2019 – 7/15/2019) and the percent change in maximum light (i.e., the average maximum light from 7/18/2019 – 8/16/2019 minus the average maximum light from 6/16/2019 – 7/15/2019 divided by the average maximum light from 6/16/2019 – 7/15/2019 multiplied by 100). For the causal model, maximum temperature was extracted from the logger data for the specific dates and times of water collection for comparison with biogeochemistry and ecosystem metabolism measurements. Hourly ocean temperatures over the experimental period were extracted from a nearshore ODFW Marine Reserve mooring sensor at 1 meter depth within ORMR to descriptively compare tide pool temperatures to the local ocean.&lt;/p&gt; &lt;p&gt;To determine biogeochemistry fluxes and ecosystem metabolism (NEC and NEP) before and after removal of foundation species, we collected daytime and nighttime water samples during low tide. We used a block design for water sampling with two daytime and two nighttime sampling events due to the timing of low tide and time restraints of sampling, where n = 16 pools (n=8 per foundation species) were measured on separate day and night sampling events. Each sampling event had an equal number of pools per foundation species and treatment group (removal or control: n = 8 pools). In situ temperature, DO, salinity, pH, and discrete samples for dissolved inorganic nutrients (NH4+, NO2− + NO3−, PO43 −) were collected hourly over a four-hour period in each pool and the adjacent ocean following methods by Silbiger &amp;amp; Sorte (2018). Temperature, DO, and salinity were measured with a calibrated multi-parameter pro meter directly in each pool (YSI Pro 2030, Lot #18B100763, Yellowsprings, OH, USA). For pH, nutrients, and TA, we collected 400 mL discrete water samples into a sealed Erlenmeyer flask using a vacuum hand pump (Mityvac, St. Louis, MO, USA) from the deepest part of the pool. Discrete samples (~250mL) for total alkalinity (TA) were taken four times over the low tide period. To compare tide pool conditions to the open ocean, ocean measurements were taken from the surface adjacent to the site.&lt;/p&gt; &lt;p&gt;pHT was measured within one hour of water collection from the sealed Erlenmeyer flask using an Orion Star Multiparameter Meter with a ROSS Ultra glass electrode (Thermo Scientific, USA, accuracy = ±0.2 mV, resolution = ±0.1, drift &amp;lt; 0.005 pH units per day) and a traceable digital thermometer (5-077-8, accuracy = 0.05 °C, resolution = 0.001 °C; Control Company, Friendswood, TX, USA) following Dickson SOP 6 (Dickson et al. 2007). The glass electrode measured millivolts (mV) and was calibrated within 48 h of each sampling event using a multipoint calibration to a tris standard solution from the Dickson Lab at Scripps Institution of Oceanography following Dickson SOP 6a (Dickson et al. 2007). TA seawater samples were placed in 250 mL Nalgene bottles with 100 µl of 50% saturated HgCl2 to preserve the water within five hours of collection. Seawater samples for nutrient analysis were filtered through GF/F filters (0.7µm) with a syringe into designated 50 ml centrifuge tubes and frozen within five hours of collection. All sampling and storage containers were soaked in 10% HCl for 24 hours, rinsed with MilliQ water, and rinsed three times with sample water before sampling events.&lt;/p&gt; &lt;p&gt;In situ pH was calculated using the seacarb package in R (Gattuso et al. 2018) by correcting for the in situ temperature in each tide pool from a multi-parameter pro meter. Total alkalinity seawater samples were processed using open-cell potentiometric titrations on a Mettler-Toledo T5 auto-titrator (Columbus, OH, USA) following Dickson SOP 3b (Dickson et al. 2007). A certified reference material (CRM) from the Dickson Lab at the Scripps Institution of Oceanography was used at the beginning of each sample group run. The accuracy of the CRM never exceeded ± 0.79% (precision = 5 µmol kg-1) from the standard value and TA samples were corrected for deviations from the standard value. Dissolved inorganic nutrients (NH4+, NO2− + NO3−, PO43 −) were analyzed at Moss Landing Marine Laboratory using a Lachat Quickchem 8000 Flow Injection Analyzer (± 1.21% NH4+, ± 0.26% NO2− + NO3−, ± 3.57% PO43 − instrument precision; Hach, Loveland, CO, USA). After processing nutrients, one mussel control pool (Pool ID 30) was removed from all resource flux analyses due to abnormally high values of ammonium compared to other pools (4290 ± 188.1 µmol g-1), likely due to contamination, on one sampling day leaving n = 7 control pools.&lt;/p&gt; &lt;p&gt;We used the total alkalinity anomaly technique (Chisholm &amp;amp; Gattuso 1991) to calculate net ecosystem calcification ( mmol CaCO3 m−2 hr−1) where total alkalinity values were divided by 1000 to convert from μmol kg−1 to mmol kg−1. ΔTA/2 is the salinity-normalized and nutrient-corrected TA (mmol kg−1) between each time point (n = 2−3 values per pool) divided by 2, where one mole of CaCO3 is formed per 2 moles of TA times ρ is the density of seawater (1023 kg m−3) times V is the volume of water in the pool at each time point (m3) divided by SA is the bottom surface area of the tide pool (m2) and t is the time between sampling points (h). and net ecosystem production rates (mmol C m-2 hr-1) were calculated from differences in dissolved inorganic carbon (DIC), calculated from pHT and TA (error propagation = 7.09 ± 0.06 mmol kg -1) with the seacarb package in R (Gattuso et al. 2018), using Gattuso et al. 1999’s equation: NEP = ΔDIC is the difference in salinity-normalized DIC (mmol kg−1) between each time point (around n = 3 values per pool) times density of seawater times the volume of the tide pool divided by the surface area of the tide pool and the time between sampling. NEC is subtracted to account for changes in DIC by the precipitation or dissolution of CaCO3, and FCO2 (mmol m−2 h−1) is the air−sea flux of CO2, which was subtracted to account for the flux in CO2 from the air−sea exchange. FCO2 is k is the gas transfer velocity (m h−1) calculated from wind speed (Ho et al. 2006) using the closest weather station, around 10 miles south of Otter Rock (NOAA Station NWPO3: 44° 36’ 46.8” N, 124° 04’ 01.2” W) multiplied by s (the solubility of CO2 in seawater calculated from in situ temperature and salinity (Weiss 1974) times CO2 (μatm) in water is calculated from pHT and TA values minus CO2 in the air was 410 μatm based on concurrent measurements at the Mauna Loa Observatory (Tans &amp;amp; Keeling 2019).&lt;/p&gt; &lt;p&gt;&lt;strong&gt;Known Issues: &lt;/strong&gt;&lt;/p&gt; &lt;ul&gt; &lt;li&gt;Tide pool 30 was removed because of high nutrient values in the after period. All time period 5 data were removed from analysis due to issues with sampling with the incoming tide.&lt;/li&gt; &lt;li&gt;Missing 8 samples from nutrient analysis out of 612 samples: N_17_5_BC was lost before taking nutrients, D_24_1_BC, N_24_1_BC, N2_Ocean_5_BC, D_1_2_AI, D_16_2_AI, D_4_3_AI, N_20_3_AI were missing from lab analyses. We removed time points that nutrients were missing. N_26_1_BC has very low phosphate that yields it to be an outlier in later analysis could be a typo or the tide pool was just exposed from the ocean during sampling so time point was removed.&lt;/li&gt; &lt;li&gt;Tide pool 18 was removed from mobile community composition analysis.&lt;/li&gt; &lt;/ul&gt; Funding provided by NSF Division of Ocean Sciences (NSF OCE) Award Number: OCE-2044837 Award URL: https://www.nsf.gov/awardsearch/show-award?AWD_ID=2044837 completed Nyssa Silbiger California State University Northridge 818-677-2901 HI USA silbiger@hawaii.edu pointOfContact asNeeded Dataset Version: 1 Unknown Pool_ID Id_code Time_Point Foundation_spp Before_After Removal_Control Day_Night Group TADeltaTime DICDeltaTime DeltaTime DeltaNN DeltaNH4 DeltaPO DeltapH DeltapCO2 DeltaDO SurfaceArea Vol SAtoV TideHeight SamplingVolume AdjSamplingVolume DeltaTA_N_Norm NEC Temp pCO2 Salinity Wind pH DO NN NH4 PO K dpCO2 a F_CO2 NEP AdjMusselCover AdjSurfgrassCover MusselCover SurfgrassCover allCCA macroalgae macrophytes macroCCA consumers allconsumers prodphyllodom allproddom FisherBrand™ Traceable™; Model 5-077-8 Lachat Quickchem 8000 Flow Injection Analyzer (Hach) Orion Star Multiparameter Meter Onset. HOBO. Pendant light intensity data logger MX2202, Bourne ROSS Ultra glass electrode Mettler-Toledo T5 auto-titrator Multi-parameter pro meter (YSI Pro 2030, Lot #18B100763, Yellowsprings, OH, USA theme None, User defined Site_ID station id sample time dominant_species comments time_of_day treatment total alkalinity (TA) dissolved inorganic Carbon time difference nitrate plus nitrite Ammonium phosphate pH Partial pressure of CO2 dissolved Oxygen surface_area volume ratio tide volume of water filtered calcification water temperature at measurement depth salinity wind_speed No BCO-DMO term carbon dioxide flux production percent coverage featureType BCO-DMO Standard Parameters digital thermometer Flow Injection Analyzer Multi Parameter Bench Meter Onset HOBO Pendant Temperature/Light Data Logger pH Sensor Titrator YSI Professional Plus Multi-Parameter Probe instrument BCO-DMO Standard Instruments otherRestrictions otherRestrictions Access Constraints: none. Use Constraints: Please follow guidelines at: http://www.bco-dmo.org/terms-use Distribution liability: Under no circumstances shall BCO-DMO be liable for any direct, incidental, special, consequential, indirect, or punitive damages that result from the use of, or the inability to use, the materials in this data submission. If you are dissatisfied with any materials in this data submission your sole and exclusive remedy is to discontinue use. CAREER: Predicting ecosystem metabolism of rocky intertidal communities in warming and acidifying oceans. https://osprey.bco-dmo.org/project/876576 CAREER: Predicting ecosystem metabolism of rocky intertidal communities in warming and acidifying oceans. <p><strong>NSF Award Abstract:</strong><br /> This award is funded in whole or in part under the American Rescue Plan Act of 2021 (Public Law 117-2).</p> <p>The devastating impacts of global warming and ocean acidification (OA) on rocky intertidal ecosystems are expected to increase as the oceans continue to warm and acidify. Further, loss of critical foundation species as a result of anthropogenic stressors lead to changes in local conditions that alter ecosystem functioning. While data exist on the physiological response of individual organisms to OA and warming in rocky systems, far less research has been conducted on community and ecosystem-scale metabolic responses. The proposed work provides a critical step in understanding how altered environmental conditions affect ecosystem functioning in rocky intertidal systems through the combination of controlled laboratory studies, field experiments, and a synthetic modeling approach integrating experimental results with pre-existing time-series data. The proposed work is focused on the overarching question: How does shifting environmental variability and loss of foundation species interact to affect ecosystem functioning in rocky intertidal communities? The PI is integrating research with ecological and quantitative educational opportunities including classroom and hands-on training in lab/field methods in marine ecology, and data science and coding bootcamps for undergraduate students. Importantly, these opportunities financially and educationally support traditionally underrepresented students at one of the largest minority-serving institutions in the country. This project provides support for intensive mentoring and training for 35 undergraduate students (25 in data science and 10 in field/lab science, all paid), 2-3 masters students, and hands-on marine ecology opportunities in the classroom (~125 students). In addition to formal education, the PI is collaborating with an artist-in-residence to communicate science to the broader public through interactive and immersive art installations in Los Angeles.</p> <p>Because rocky intertidal systems provide important ecosystem services including food production, coastal protection, and tourism, it is critical to understand how warming (both air and ocean), acidification, and altered community states affect reef-scale ecosystem metabolism. While information exists on responses of a variety of individual intertidal taxa to temperature stress, ocean acidification, and habitat loss, there is notably less on the response of ecosystem metabolism (e.g. NEP and NEC) to these stressors. This proposed work is focused on a series of conceptual knowledge gaps and tests mechanisms through which different warming regimes, lowered pH, and community disturbance lead to altered community metabolism and ultimately affect ecosystem function. Specifically, this study: 1) Describes community thermal performance curves of multiple ecosystem functions under differing pH conditions in experimental mesocosms, 2) characterizes drivers of ecosystem functioning (e.g. Net Ecosystem Production [NEP] and Net Ecosystem Calcification [NEC]) in situ using natural changes in environmental variability before and after a disturbance, and 3) integrates 16 years of publicly available community composition and environmental time-series data with lab and field data to hindcast ecosystem metabolic rates and predict how ecosystem metabolism may change in the future.</p> <p>This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.</p> TIDES largerWorkCitation project eng; USA oceans -124.066111 -124.066111 44.752528 44.752528 2019-06-01 2019-08-31 California, USA 0 BCO-DMO catalogue of parameters from Results of a study examining how foundation species loss alters multiple ecosystem functions based on community surveys and biogeochemical sampling of tidepools in the Otter Rock Marine Reserve (ORMR) and Marine Garden, Oregon USA from June to August 2019 Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. pointOfContact http://lod.bco-dmo.org/id/dataset-parameter/878432.rdf Name: Pool_ID Units: unitless Description: <p>tide pool id</p> http://lod.bco-dmo.org/id/dataset-parameter/878433.rdf Name: Id_code Units: unitless Description: <p>day or night, tide pool id, sampling period, and if it&#039;s in the before (BC) or after (AI) foundation species removal period</p> http://lod.bco-dmo.org/id/dataset-parameter/878434.rdf Name: Time_Point Units: unitless Description: <p>Time point of sampling (1-4)</p> http://lod.bco-dmo.org/id/dataset-parameter/878435.rdf Name: Foundation_spp Units: unitless Description: <p>mussel (Mytilus) or surfgrass (Phyllospadix) tide pool</p> http://lod.bco-dmo.org/id/dataset-parameter/878436.rdf Name: Before_After Units: unitless Description: <p>Before or after foundation species removal</p> http://lod.bco-dmo.org/id/dataset-parameter/878437.rdf Name: Removal_Control Units: unitless Description: <p>Removal pools that had foundation species removed in the after period or control</p> http://lod.bco-dmo.org/id/dataset-parameter/878438.rdf Name: Day_Night Units: unitless Description: <p>Day or night sampling</p> http://lod.bco-dmo.org/id/dataset-parameter/878439.rdf Name: Group Units: unitless Description: <p>Combination or Before/After and Removal/Control</p> http://lod.bco-dmo.org/id/dataset-parameter/878440.rdf Name: TADeltaTime Units: mmol kg_1 Description: <p>Change in total alkalinity over sampling time points (2-1, 3-2, 4-3)</p> http://lod.bco-dmo.org/id/dataset-parameter/878441.rdf Name: DICDeltaTime Units: mmol kg_1 Description: <p>Change in dissolved inorganic carbon over sampling time points (2-1, 3-2, 4-3)</p> http://lod.bco-dmo.org/id/dataset-parameter/878442.rdf Name: DeltaTime Units: hrs Description: <p>change in time between sampling periods</p> http://lod.bco-dmo.org/id/dataset-parameter/878443.rdf Name: DeltaNN Units: µmol/L Description: <p>Change in NO2_ + NO3_ over sampling time points (2-1, 3-2, 4-3)</p> http://lod.bco-dmo.org/id/dataset-parameter/878444.rdf Name: DeltaNH4 Units: µmol/L Description: <p>Change in NH4+ over sampling time points (2-1, 3-2, 4-3)</p> http://lod.bco-dmo.org/id/dataset-parameter/878445.rdf Name: DeltaPO Units: µmol/L Description: <p>Change in PO43 _ over sampling time points (2-1, 3-2, 4-3)</p> http://lod.bco-dmo.org/id/dataset-parameter/878446.rdf Name: DeltapH Units: unitless Description: <p>Change in pH over sampling time points (2-1, 3-2, 4-3)</p> http://lod.bco-dmo.org/id/dataset-parameter/878447.rdf Name: DeltapCO2 Units: µatm Description: <p>Change in pCO2 over sampling time points (2-1, 3-2, 4-3)</p> http://lod.bco-dmo.org/id/dataset-parameter/878448.rdf Name: DeltaDO Units: mg/L Description: <p>Change in dissolved oxygen over sampling time points (2-1, 3-2, 4-3)</p> http://lod.bco-dmo.org/id/dataset-parameter/878449.rdf Name: SurfaceArea Units: m2 Description: <p>surface area of tide pool</p> http://lod.bco-dmo.org/id/dataset-parameter/878450.rdf Name: Vol Units: m3 Description: <p>volume of tide pool via dye method</p> http://lod.bco-dmo.org/id/dataset-parameter/878451.rdf Name: SAtoV Units: m Description: <p>ratio of surface area to volume</p> http://lod.bco-dmo.org/id/dataset-parameter/878452.rdf Name: TideHeight Units: m Description: <p>tide height of tide pools</p> http://lod.bco-dmo.org/id/dataset-parameter/878453.rdf Name: SamplingVolume Units: mL Description: <p>sampling volume from water sampling</p> http://lod.bco-dmo.org/id/dataset-parameter/878454.rdf Name: AdjSamplingVolume Units: mL Description: <p>sampling volume from water sampling included collection error</p> http://lod.bco-dmo.org/id/dataset-parameter/878455.rdf Name: DeltaTA_N_Norm Units: mmol kg_1 Description: <p>Change in total alkalinity over sampling time points (2-1, 3-2, 4-3) normalized with nutrients</p> http://lod.bco-dmo.org/id/dataset-parameter/878456.rdf Name: NEC Units: mmol CaCO3 m_2 hr_1 Description: <p>net ecosystem calcification rates</p> http://lod.bco-dmo.org/id/dataset-parameter/878457.rdf Name: Temp Units: degrees celsius (°C) Description: <p>mean temperature over time points</p> http://lod.bco-dmo.org/id/dataset-parameter/878458.rdf Name: pCO2 Units: µatm Description: <p>mean of pCO2 between time points</p> http://lod.bco-dmo.org/id/dataset-parameter/878459.rdf Name: Salinity Units: ppt Description: <p>mean salinity between time points</p> http://lod.bco-dmo.org/id/dataset-parameter/878460.rdf Name: Wind Units: m h_1 Description: <p>mean wind between time points</p> http://lod.bco-dmo.org/id/dataset-parameter/878461.rdf Name: pH Units: unitless Description: <p>mean pH between time points</p> http://lod.bco-dmo.org/id/dataset-parameter/878462.rdf Name: DO Units: mg/L Description: <p>mean dissolved oxygen between time points</p> http://lod.bco-dmo.org/id/dataset-parameter/878463.rdf Name: NN Units: µmol/L Description: <p>mean NO2_ + NO3_ between time points</p> http://lod.bco-dmo.org/id/dataset-parameter/878464.rdf Name: NH4 Units: µmol/L Description: <p>mean NH4+ between time points</p> http://lod.bco-dmo.org/id/dataset-parameter/878465.rdf Name: PO Units: µmol/L Description: <p>mean PO43- between time points</p> http://lod.bco-dmo.org/id/dataset-parameter/878466.rdf Name: K Units: m h_1 Description: <p>gas transfer velocity</p> http://lod.bco-dmo.org/id/dataset-parameter/878467.rdf Name: dpCO2 Units: µatm Description: <p>difference btween air and seawater pCO2 pCO2_water - pCO2_atm</p> http://lod.bco-dmo.org/id/dataset-parameter/878468.rdf Name: a Units: mmol L^-1 atm^-1 or mmol m^-3 uatm^-1 Description: <p>CO2 solubility constant according to Weiss (1974) solubility</p> http://lod.bco-dmo.org/id/dataset-parameter/878469.rdf Name: F_CO2 Units: mmol m_2 hr_1 Description: <p>CO2 flux</p> http://lod.bco-dmo.org/id/dataset-parameter/878470.rdf Name: NEP Units: mmol C m-2 hr-1 Description: <p>Net ecosystem production rates</p> http://lod.bco-dmo.org/id/dataset-parameter/878471.rdf Name: AdjMusselCover Units: percent Description: <p>percent cover of mussels adjusted to be proportional to 100%</p> http://lod.bco-dmo.org/id/dataset-parameter/878472.rdf Name: AdjSurfgrassCover Units: percent Description: <p>percent cover of surfgrass adjusted to be proportional to 100%</p> http://lod.bco-dmo.org/id/dataset-parameter/878473.rdf Name: MusselCover Units: percent Description: <p>raw percent cover of mussels</p> http://lod.bco-dmo.org/id/dataset-parameter/878474.rdf Name: SurfgrassCover Units: percent Description: <p>raw percent cover of surfgrass</p> http://lod.bco-dmo.org/id/dataset-parameter/878475.rdf Name: allCCA Units: percent Description: <p>percent cover crustose coralline and coralline algae</p> http://lod.bco-dmo.org/id/dataset-parameter/878476.rdf Name: macroalgae Units: percent Description: <p>percent cover of algae</p> http://lod.bco-dmo.org/id/dataset-parameter/878477.rdf Name: macrophytes Units: percent Description: <p>percent cover algae + surfgrass</p> http://lod.bco-dmo.org/id/dataset-parameter/878478.rdf Name: macroCCA Units: percent Description: <p>percent cover all macroalgae + all CCA</p> http://lod.bco-dmo.org/id/dataset-parameter/878479.rdf Name: consumers Units: count Description: <p>percent cover consumers (e.g. barnacles)</p> http://lod.bco-dmo.org/id/dataset-parameter/878480.rdf Name: allconsumers Units: percent Description: <p>percent cover consumers (e.g. barnacles) + mussel cover</p> http://lod.bco-dmo.org/id/dataset-parameter/878481.rdf Name: prodphyllodom Units: percent Description: <p>percent cover macroalgae - all consumers</p> http://lod.bco-dmo.org/id/dataset-parameter/878482.rdf Name: allproddom Units: percent Description: <p>percent cover producers - all consumers</p> GB/NERC/BODC > British Oceanographic Data Centre, Natural Environment Research Council, United Kingdom Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. pointOfContact 131901 https://darchive.mblwhoilibrary.org/bitstream/1912/29320/1/dataset-878413_oregon-tide-pools__v1.tsv download https://doi.org/10.26008/1912/bco-dmo.878413.1 download onLine dataset &lt;p&gt;We selected 32 tide pools within the Oregon Department of Fish and Wildlife’s (ODFW) Otter Rock Marine Reserve (ORMR) and Marine Garden, Oregon USA (44°45'09.1&quot;N 124°03'58.0&quot;W) during the months of June to August 2019. One mussel species (&lt;em&gt;M. califorinanus&lt;/em&gt;) and two surfgrass species (&lt;em&gt;P. scouleri &lt;/em&gt;and &lt;em&gt;P. torreyi&lt;/em&gt;) were present at the site. Of the 32 tide pools, 16 were dominated by the California mussel (&lt;em&gt;M. califorinanus&lt;/em&gt;) and 16 were dominated by surfgrass (&lt;em&gt;Phyllospadix spp&lt;/em&gt;.; 10 pools contained &lt;em&gt;Phyllospadix scouleri &lt;/em&gt;and six pools had both &lt;em&gt;Phyllospadix scouleri &lt;/em&gt;and &lt;em&gt;Phyllospadix torreyi&lt;/em&gt;). Percent cover of foundation species ranged from 45.3 to 98.9% in the mussel pools and 49.5 to 100% in the surfgrass tide pools. Some tide pools had both mussels and surfgrasses; however, the presence of a second foundation species did not exceed 7.2% of the remaining tide pool cover. Tide pools were located in the mid to low intertidal zone ranging from 0.71 meters to 1.77 meters above mean lower-low water. At this tide height, tide pools were isolated for 4-6 hours during summer low tides.&lt;/p&gt; &lt;p&gt;A Before-After-Control-Impact (BACI) design was used to account for changes in ocean chemistry, timing of low tide, and variability within tide pools throughout the experimental period. The BACI design consisted of two 29-day time periods (before and after removal) with control and removal tide pools (Stewart-Oaten et al. 1986) where foundation species were removed from removal tide pools (n = 8 per foundation species) between the before and after removal periods. The before-removal period occurred June to mid-July 2019 and the after-removal period occurred mid-July to August 2019. Tide pools were selected for control or removal groups using a random number generator while accounting for the surface area to percent cover of foundation species. Removal of all foundation species from removal tide pools occurred in mid-July 2019 and tide pools had a one-week recovery period before any surveys or samples were taken. Rhizomes were also removed in surfgrass removal pools due to their ability to alter nutrient cycling (Terrados &amp;amp; Williams 1997). During each time period, we characterized physical parameters (pool volume and tide height), community composition, light, temperature, and biogeochemical fluxes (e.g., dissolved oxygen, pH, nutrients) in each tide pool.&lt;/p&gt; &lt;p&gt;Tide pool physical parameters included tidal height (location within intertidal) and tide pool volume (size of pool). Tide heights for each pool were surveyed with a laser level and stadia rod (DeWalt, Towson, MD, USA). Tide pool volume (V) was determined using a dye method (Pfister 1995) and measured with a SmartSpec3000 spectrophotometer (Bio-Rad Lab, Hercules, CA, USA). Water volume in the tide pools changed by 0 - 20% between the two timepoints due to removal of foundation species. To account for the slight effect of changing biomass on volume, volume was re-measured post-removal and the average between both time periods were used in statistical analyses.&amp;nbsp;&amp;nbsp;&lt;/p&gt; &lt;p&gt;We conducted two rounds of community composition surveys for sessile percent cover and mobile organism counts: three-weeks before removal and one-month post-removal. We temporarily removed seawater from the tide pool and placed a flexible mesh quadrat (Nielsen 2001) with demarcations in 10 x 10 cm squares over the bottom of each pool to survey the entire community (pools ranged from 58 to 754 squares). We measured percent cover for sessile organisms and counted mobile organisms, identifying down to the lowest possible taxonomic unit in the field (usually genus level). We normalized the sum of non-foundation species sessile cover to 100% cover per tide pool, including the opposite foundation species if present (e.g., mussel cover within surfgrass pools was also normalized to 100% cover). Each foundation species cover in their respective tide pools remained the raw percent cover and did not exceed 100%. Both sessile and mobile organisms were grouped into larger functional groups based on their ecological role for data visualizations. All tide pool characterization and community composition methods were completed at least 24 hours before any water sampling event to allow the pools to be flushed at least twice by the ocean before measurements.&lt;/p&gt; &lt;p&gt;Temperature (°C) and light intensity (lumens m-2) were recorded continuously every 15 minutes for 29 days during each time period using HOBO® Pendant loggers bolted facing up on the flattest part of the tide pool on a level platform in the interstitial spaces of the foundation species (Onset® HOBO® Pendent Light Intensity Data Logger MX2202, Bourne, MA, USA). Light intensity (lumens m-2) was converted to photon flux density (PFD; µmol photons m-2 s-1) following Long et al. (2012) field experiment values. We measured the change in maximum temperature between the before and after removal period (i.e., the average of the daily hottest temperature from 7/18/2019 – 8/16/2019 minus the average daily hottest temperature from 6/16/2019 – 7/15/2019) and the percent change in maximum light (i.e., the average maximum light from 7/18/2019 – 8/16/2019 minus the average maximum light from 6/16/2019 – 7/15/2019 divided by the average maximum light from 6/16/2019 – 7/15/2019 multiplied by 100). For the causal model, maximum temperature was extracted from the logger data for the specific dates and times of water collection for comparison with biogeochemistry and ecosystem metabolism measurements. Hourly ocean temperatures over the experimental period were extracted from a nearshore ODFW Marine Reserve mooring sensor at 1 meter depth within ORMR to descriptively compare tide pool temperatures to the local ocean.&lt;/p&gt; &lt;p&gt;To determine biogeochemistry fluxes and ecosystem metabolism (NEC and NEP) before and after removal of foundation species, we collected daytime and nighttime water samples during low tide. We used a block design for water sampling with two daytime and two nighttime sampling events due to the timing of low tide and time restraints of sampling, where n = 16 pools (n=8 per foundation species) were measured on separate day and night sampling events. Each sampling event had an equal number of pools per foundation species and treatment group (removal or control: n = 8 pools). In situ temperature, DO, salinity, pH, and discrete samples for dissolved inorganic nutrients (NH4+, NO2− + NO3−, PO43 −) were collected hourly over a four-hour period in each pool and the adjacent ocean following methods by Silbiger &amp;amp; Sorte (2018). Temperature, DO, and salinity were measured with a calibrated multi-parameter pro meter directly in each pool (YSI Pro 2030, Lot #18B100763, Yellowsprings, OH, USA). For pH, nutrients, and TA, we collected 400 mL discrete water samples into a sealed Erlenmeyer flask using a vacuum hand pump (Mityvac, St. Louis, MO, USA) from the deepest part of the pool. Discrete samples (~250mL) for total alkalinity (TA) were taken four times over the low tide period. To compare tide pool conditions to the open ocean, ocean measurements were taken from the surface adjacent to the site.&lt;/p&gt; &lt;p&gt;pHT was measured within one hour of water collection from the sealed Erlenmeyer flask using an Orion Star Multiparameter Meter with a ROSS Ultra glass electrode (Thermo Scientific, USA, accuracy = ±0.2 mV, resolution = ±0.1, drift &amp;lt; 0.005 pH units per day) and a traceable digital thermometer (5-077-8, accuracy = 0.05 °C, resolution = 0.001 °C; Control Company, Friendswood, TX, USA) following Dickson SOP 6 (Dickson et al. 2007). The glass electrode measured millivolts (mV) and was calibrated within 48 h of each sampling event using a multipoint calibration to a tris standard solution from the Dickson Lab at Scripps Institution of Oceanography following Dickson SOP 6a (Dickson et al. 2007). TA seawater samples were placed in 250 mL Nalgene bottles with 100 µl of 50% saturated HgCl2 to preserve the water within five hours of collection. Seawater samples for nutrient analysis were filtered through GF/F filters (0.7µm) with a syringe into designated 50 ml centrifuge tubes and frozen within five hours of collection. All sampling and storage containers were soaked in 10% HCl for 24 hours, rinsed with MilliQ water, and rinsed three times with sample water before sampling events.&lt;/p&gt; &lt;p&gt;In situ pH was calculated using the seacarb package in R (Gattuso et al. 2018) by correcting for the in situ temperature in each tide pool from a multi-parameter pro meter. Total alkalinity seawater samples were processed using open-cell potentiometric titrations on a Mettler-Toledo T5 auto-titrator (Columbus, OH, USA) following Dickson SOP 3b (Dickson et al. 2007). A certified reference material (CRM) from the Dickson Lab at the Scripps Institution of Oceanography was used at the beginning of each sample group run. The accuracy of the CRM never exceeded ± 0.79% (precision = 5 µmol kg-1) from the standard value and TA samples were corrected for deviations from the standard value. Dissolved inorganic nutrients (NH4+, NO2− + NO3−, PO43 −) were analyzed at Moss Landing Marine Laboratory using a Lachat Quickchem 8000 Flow Injection Analyzer (± 1.21% NH4+, ± 0.26% NO2− + NO3−, ± 3.57% PO43 − instrument precision; Hach, Loveland, CO, USA). After processing nutrients, one mussel control pool (Pool ID 30) was removed from all resource flux analyses due to abnormally high values of ammonium compared to other pools (4290 ± 188.1 µmol g-1), likely due to contamination, on one sampling day leaving n = 7 control pools.&lt;/p&gt; &lt;p&gt;We used the total alkalinity anomaly technique (Chisholm &amp;amp; Gattuso 1991) to calculate net ecosystem calcification ( mmol CaCO3 m−2 hr−1) where total alkalinity values were divided by 1000 to convert from μmol kg−1 to mmol kg−1. ΔTA/2 is the salinity-normalized and nutrient-corrected TA (mmol kg−1) between each time point (n = 2−3 values per pool) divided by 2, where one mole of CaCO3 is formed per 2 moles of TA times ρ is the density of seawater (1023 kg m−3) times V is the volume of water in the pool at each time point (m3) divided by SA is the bottom surface area of the tide pool (m2) and t is the time between sampling points (h). and net ecosystem production rates (mmol C m-2 hr-1) were calculated from differences in dissolved inorganic carbon (DIC), calculated from pHT and TA (error propagation = 7.09 ± 0.06 mmol kg -1) with the seacarb package in R (Gattuso et al. 2018), using Gattuso et al. 1999’s equation: NEP = ΔDIC is the difference in salinity-normalized DIC (mmol kg−1) between each time point (around n = 3 values per pool) times density of seawater times the volume of the tide pool divided by the surface area of the tide pool and the time between sampling. NEC is subtracted to account for changes in DIC by the precipitation or dissolution of CaCO3, and FCO2 (mmol m−2 h−1) is the air−sea flux of CO2, which was subtracted to account for the flux in CO2 from the air−sea exchange. FCO2 is k is the gas transfer velocity (m h−1) calculated from wind speed (Ho et al. 2006) using the closest weather station, around 10 miles south of Otter Rock (NOAA Station NWPO3: 44° 36’ 46.8” N, 124° 04’ 01.2” W) multiplied by s (the solubility of CO2 in seawater calculated from in situ temperature and salinity (Weiss 1974) times CO2 (μatm) in water is calculated from pHT and TA values minus CO2 in the air was 410 μatm based on concurrent measurements at the Mauna Loa Observatory (Tans &amp;amp; Keeling 2019).&lt;/p&gt; &lt;p&gt;&lt;strong&gt;Known Issues: &lt;/strong&gt;&lt;/p&gt; &lt;ul&gt; &lt;li&gt;Tide pool 30 was removed because of high nutrient values in the after period. All time period 5 data were removed from analysis due to issues with sampling with the incoming tide.&lt;/li&gt; &lt;li&gt;Missing 8 samples from nutrient analysis out of 612 samples: N_17_5_BC was lost before taking nutrients, D_24_1_BC, N_24_1_BC, N2_Ocean_5_BC, D_1_2_AI, D_16_2_AI, D_4_3_AI, N_20_3_AI were missing from lab analyses. We removed time points that nutrients were missing. N_26_1_BC has very low phosphate that yields it to be an outlier in later analysis could be a typo or the tide pool was just exposed from the ocean during sampling so time point was removed.&lt;/li&gt; &lt;li&gt;Tide pool 18 was removed from mobile community composition analysis.&lt;/li&gt; &lt;/ul&gt; Specified by the Principal Investigator(s) &lt;p&gt;&lt;strong&gt;Data Processing:&lt;/strong&gt;&lt;br /&gt; We conducted separate statistical analyses for surfgrass and mussel tide pools for all response variables due to the difference in effects that each foundation species has on the ecosystem. Change in foundation species loss between the before and after removal periods varied continuously across treatment groups from -27.5% to 100% change in cover in the surfgrass tide pools and -10.4 to 98.5% change in cover in the mussel pools (Fig. S1, Fields &amp;amp; Silbinger 2022), with negative values indicating there was an increase in foundation species cover between time points. Therefore, we investigated foundation species loss as a continuous gradient rather than using categorical ANOVA-style design typical for BACI studies (Smith 2002).&lt;/p&gt; &lt;p&gt;Multivariate plots and analyses were used to assess how community composition and biogeochemistry parameters changed between time periods. We used a principal coordinate analysis (PCoA) in the vegan package (Oksanen et al. 2007) to visualize how the gradient of percent foundation species loss altered the change (after – before removal periods) in the cover of sessile organisms and counts in mobile organisms by functional group. Mobile community data were square root transformed to reduce the effect of rare values. One surfgrass control tide pool (Pool ID 18) was removed from the mobile community analysis because of outlier values of limpets and littorine snails. A PERMANOVA was used to test if foundation species loss altered functional community composition, with tide height and tide pool size included as covariates. For the biogeochemistry data, principal component analyses (PCA) were used to reduce the dimensionality and visualize the multivariate data. Data was centered and standardized before visualizing the PCA. Specifically, shifts between the before and after-removal period in biogeochemistry (DO, nutrients, and pH) and temperature mean, maximum, and variance of n = 8 (four day and four night) time points were visualized to observe differences for local tide pool conditions and the adjacent ocean from water sampling events.&lt;/p&gt; &lt;p&gt;We used general linear models to test the effect of foundation species loss on change (after – before removal periods) in species richness, percent maximum light, and maximum temperature in tide pools, with tide height and tide pool volume as covariates. Model assumptions for general linear models were verified by investigating residual plots of the model for homogeneity of variance and normality and data that did not meet assumptions were transformed. Multi-collinearity was tested with Pearson correlations between each predictor. All data was processed and analyzed using R statistical program (v 4.0.2; R Core Team 2020).&lt;/p&gt; &lt;p&gt;We conducted a piecewise structural equation model (SEM) using PiecewiseSEM 2.1.0 (Lefcheck 2016) understand how mussel and surfgrass percent loss and tide pool physical characteristics affected community changes, physical environmental conditions, biogeochemistry, and ecosystem metabolism of tide pool communities. Piecewise SEMs combine multiple linear models into a single causal framework (Shipley 2009, Lefcheck 2016). Because the model is pieced together with local estimation, it is more flexible than traditional SEMs and allows for inclusion of non-normal distributions, random effects, nested models, and smaller sample size within the causal model network (Lefcheck &amp;amp; Duffy 2015). Specifically, we tested direct and indirect pathways of surfgrass and mussel percent loss on ecosystem metabolism (NEC and NEP) mediated by changes in micro/macroalgae cover, temperature, nutrients (dissolved nitrogen to phosphate ratio), pH, and physical parameters of the tide pool. We conducted a multigroup analysis to test if each path varied by day and night, as expected from previous tide pool studies (Kwiatkowski et al. 2016, Bracken et al. 2018, Silbiger &amp;amp; Sorte 2018, Wolfe et al. 2020). The fit of each model was determined with a Shipley’s test of d-separation, which verifies there are no missing inferences from the model using a Fisher’s C statistic, where p&amp;gt;0.05 indicates that the model is a good fit for the data (Shipley 2009). However, since our study’s ratio of total number of samples (16) to the number of variables (9) is 1.78 versus the recommended 5 (Grace et al. 2015, Lefcheck 2016), the smaller sample size could influence the goodness of fit of the model. Path coefficients were standardized allowing for comparisons of the magnitude of effect among groups.&lt;/p&gt; &lt;p&gt;For all model components, we averaged values over the low tide period (n=3 - 4 values per tide pool for pH, nutrients, and temperature; n= 2 - 3 values for NEC and NEP), and then calculated the difference between the before and after removal time period (n = 16 pools for the surfgrass model and n =15 pools for the mussel model). Positive values indicate an increase, zero indicates no change, and negative indicates a decrease in that parameter after the removal period. For the community composition components of the model, we used the sessile community metrics that changed the most between the before and after removal period, including, surfgrass, mussel, and non-surfgrass, non-calcifying producers (all fleshy micro and macroalgae). We excluded coralline algae from the analysis due their different functional roles within the tide pool (e.g., calcification and production and to reduce the number of pathways in our analysis). Tide pool size, represented as volume, and tide pool height was used as a covariate in all models because physical parameters can affect temperature and biogeochemistry (Legrand et al. 2018, Wolfe et al. 2020).&lt;/p&gt; &lt;p&gt;Due to the collinearity between nutrient species, the dissolved nitrogen to phosphate ratio ([nitrite + nitrate + ammonium]: phosphate) was used as the measurement for nutrients in the SEM models. Because temperature and light were highly correlated and change in light was zero in most pools at night, only maximum temperature was used with the SEM framework (Fig. S3, Fields &amp;amp; Silbinger 2022). We used existing knowledge of intertidal systems to create hypothesized paths within each of the models (Grace 2008, Grace et al. 2012):&lt;/p&gt; &lt;p&gt;Micro/macroalgal cover ~ Foundation species loss + Volume + Tide height [1]&lt;/p&gt; &lt;p&gt;Model 1 represents the hypothesis that loss of surfgrass and mussel cover would increase producer cover by increasing free space for fast colonizers, with the volume of the pool and tide height as a covariates (Dethier 1981, 1984). &amp;nbsp;&lt;/p&gt; &lt;p&gt;Maximum Temperature ~ Foundation species loss + Volume + Tide height [2]&lt;/p&gt; &lt;p&gt;Model 2 tests the understanding that temperature will increase with surfgrass and mussel loss (Stephens &amp;amp; Bertness 1991, Shelton 2010), with size of the pool and tide height as covariates.&lt;/p&gt; &lt;p&gt;N:P ~ Foundation species loss + Volume + Tide height [3]&lt;/p&gt; &lt;p&gt;Model 3 represents the hypothesis that the N:P ratio will increase with surfgrass loss because the loss of dominant producers will lead to less nitrogen uptake (Terrados &amp;amp; Williams 1997, Ramirez-Garcia et al. 2002). We hypothesized that the N:P ratio will decrease with mussel loss due to the mussels’ ability to recycle nitrogen (Nielsen 2003, Bracken &amp;amp; Nielsen 2004, Pfister 2007, Pfister &amp;amp; Altabet 2019). N:P ratio may be affected by the size of the tide pool due to diffusive properties of tide pools (Hurd 2000) and tide height due to emersion time.&lt;/p&gt; &lt;p&gt;NEP ~ Maximum temperature + Micro/Macroalgal Cover + N:P + Tide Height;&lt;/p&gt; &lt;p&gt;NEP~ N:P + Micro/Macroalgal Cover + Tide Height [4]&lt;/p&gt; &lt;p&gt;Model 4 represents the different mechanisms that influence NEP. NEP is increased by producer cover, temperature, and nutrients (N:P) and affected by tide height (Pfister 2007, Kwiatkowski et al. 2016, Takeshita et al. 2016, Bracken et al. 2018, Duarte &amp;amp; Krause-Jensen 2018, Silbiger &amp;amp; Sorte 2018, Pfister &amp;amp; Altabet 2019, Wolfe et al. 2020). Surfgrass and mussel loss were both correlated with most of the components of the NEP model and were therefore excluded from the NEP model. Since volume was already accounted for in the calculation of NEP, it was removed from the model equation. Maximum temperature and N:P were correlated within mussel pools, so N:P ratio was used because there is more background knowledge on mussels’ effect on production via altering nutrients (Pfister 2007, Pfister &amp;amp; Altabet 2019).&lt;/p&gt; &lt;p&gt;pH ~ NEP + Foundation species loss + Volume + Tide height [5]&lt;/p&gt; &lt;p&gt;Model 5 represents our understanding that there is a positive relationship between NEP and pH (Silbiger &amp;amp; Sorte 2018). Production increases pH, while respiration decreases pH due to the uptake and release of CO2 respectively. Foundation species loss, size of tide pool, and tide height may also affect pH. Foundation species loss was also included within the pH model as a missing link determined by tests of directed separation.&lt;/p&gt; &lt;p&gt;NEC ~ pH + Maximum temperature + Tide Height [6]&lt;/p&gt; &lt;p&gt;Model 6 represents our understanding that NEC increases with pH and temperature and may be affected by the tide height (Kwiatkowski et al. 2016, Silbiger &amp;amp; Sorte 2018, Wolfe et al. 2020). Since volume was already accounted for in the calculation of NEC, it was removed from the model equation. Foundation species loss was correlated with both pH and maximum temperature and therefore was excluded from the model.&lt;br /&gt; R statistical program (v 4.0.2; R Core Team 2020)&lt;/p&gt; &lt;p&gt;&amp;nbsp;&lt;/p&gt; &lt;p&gt;&lt;strong&gt;BCO-DMO Processing Description:&lt;/strong&gt;&lt;br /&gt; - Adjusted field/parameter names to comply with BCO-DMO naming conventions&lt;br /&gt; - Added a conventional header with dataset name, PI names, version date&lt;br /&gt; &amp;nbsp;&lt;/p&gt; Specified by the Principal Investigator(s) asNeeded 7.x-1.1 Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. pointOfContact FisherBrand™ Traceable™; Model 5-077-8 FisherBrand™ Traceable™; Model 5-077-8 PI Supplied Instrument Name: FisherBrand™ Traceable™; Model 5-077-8 Instrument Name: digital thermometer Instrument Short Name: Instrument Description: An instrument that measures temperature digitally. Lachat Quickchem 8000 Flow Injection Analyzer (Hach) Lachat Quickchem 8000 Flow Injection Analyzer (Hach) PI Supplied Instrument Name: Lachat Quickchem 8000 Flow Injection Analyzer (Hach) Instrument Name: Flow Injection Analyzer Instrument Short Name:FIA Instrument Description: An instrument that performs flow injection analysis. Flow injection analysis (FIA) is an approach to chemical analysis that is accomplished by injecting a plug of sample into a flowing carrier stream. FIA is an automated method in which a sample is injected into a continuous flow of a carrier solution that mixes with other continuously flowing solutions before reaching a detector. Precision is dramatically increased when FIA is used instead of manual injections and as a result very specific FIA systems have been developed for a wide array of analytical techniques. Community Standard Description: http://vocab.nerc.ac.uk/collection/L05/current/LAB36/ Orion Star Multiparameter Meter Orion Star Multiparameter Meter PI Supplied Instrument Name: Orion Star Multiparameter Meter PI Supplied Instrument Description:Orion Star Multiparameter Meter with a ROSS Ultra glass electrode (Thermo Scientific; accuracy = +/- 0.2 mV, resolution = +/-0.1 mV, drift < 0.005 pH units d−1)   Instrument Name: Multi Parameter Bench Meter Instrument Short Name:Multi Parameter Bench Meter Instrument Description: An analytical instrument that can measure multiple parameters, such as pH, EC, TDS, DO and Temperature with one device. Onset. HOBO. Pendant light intensity data logger MX2202, Bourne Onset. HOBO. Pendant light intensity data logger MX2202, Bourne PI Supplied Instrument Name: Onset. HOBO. Pendant light intensity data logger MX2202, Bourne Instrument Name: Onset HOBO Pendant Temperature/Light Data Logger Instrument Short Name:Onset HOBO Pendant Temp/Light Instrument Description: The Onset HOBO (model numbers UA-002-64 or UA-001-64) is an in-situ instrument for wet or underwater applications. It supports light intensity, soil temperature, temperature, and water temperature. A two-channel logger with 10-bit resolution can record up to approximately 28,000 combined temperature and light measurements with 64K bytes memory. It has a polypropylene housing case. Uses an optical USB to transmit data. A solar radiation shield is used for measurement in sunlight. Temperature measurement range: -20 deg C to 70 deg C (temperature). Light measurement range: 0 to 320,000 lux. Temperature accuracy: +/- 0.53 deg C from 0 deg C to 50 deg C. Light accuracy: Designed for measurement of relative light levels. Water depth rating: 30 m. ROSS Ultra glass electrode ROSS Ultra glass electrode PI Supplied Instrument Name: ROSS Ultra glass electrode PI Supplied Instrument Description:(Thermo Scientific; accuracy = +/- 0.2 mV, resolution = +/-0.1 mV, drift < 0.005 pH units d−1)   Instrument Name: pH Sensor Instrument Short Name:pH Sensor Instrument Description: An instrument that measures the hydrogen ion activity in solutions. The overall concentration of hydrogen ions is inversely related to its pH.  The pH scale ranges from 0 to 14 and indicates whether acidic (more H+) or basic (less H+).  Mettler-Toledo T5 auto-titrator Mettler-Toledo T5 auto-titrator PI Supplied Instrument Name: Mettler-Toledo T5 auto-titrator Instrument Name: Titrator Instrument Short Name:Titrator Instrument Description: Titrators are instruments that incrementally add quantified aliquots of a reagent to a sample until the end-point of a chemical reaction is reached. Community Standard Description: http://vocab.nerc.ac.uk/collection/L05/current/LAB12/ Multi-parameter pro meter (YSI Pro 2030, Lot #18B100763, Yellowsprings, OH, USA Multi-parameter pro meter (YSI Pro 2030, Lot #18B100763, Yellowsprings, OH, USA PI Supplied Instrument Name: Multi-parameter pro meter (YSI Pro 2030, Lot #18B100763, Yellowsprings, OH, USA Instrument Name: YSI Professional Plus Multi-Parameter Probe Instrument Short Name:YSI ProPlus Instrument Description: The YSI Professional Plus handheld multiparameter meter provides for the measurement of a variety of combinations for dissolved oxygen, conductivity, specific conductance, salinity, resistivity, total dissolved solids (TDS), pH, ORP, pH/ORP combination, ammonium (ammonia), nitrate, chloride and temperature. More information from the manufacturer.