DOI: 10.5281/zenodo.11186903
title: Data and Analysis for Kaplanis, Denny, and Raimondi 2024, "Vertical distribution of rocky intertidal organisms shifts with sea-level variability on the Northeast Pacific Coast".

description: This repository contains all the data and R scripts used to produce all analyses and figures for Kaplanis, Denny, and Raimondi 2024, as well as all intermediate outputs and final figures. To access this content, download and unzip the intertidalvertdist folder (for intertidal vertical distribution). The R Project is titled "intertidalvertdist". All pertinent information needed to access data, replicate the analyses, and produce figures is contained within the README file, but a brief desciption is below.

Directory Architecture:

Data:Contains all data. Within this folder are two subdirectories - Raw Data, and Processed Data. Raw Data are unmanipulated, straight from the data source. Processed Data are outputs from scripted data wrangling and transformations.   

Within each of these folders are two more subdirectories: Tide Gauge Data, and MARINe Data. These are the two data sources used in this manuscript - monthly sea-level data from The National Oceanic and Atmospheric Administration Center for Operational Oceanographic Products and Services (NOAA CO-OPS) tide gauge stations, and long-term rocky intertidal biological monitoring data from Multi-Agency Rocky Intertidal Network (MARINe) survey sites.

Scripts:All R scripts are contained within the Scripts folder. The scripts have the prefix IVD (for intertidal vertical distribution), then a name that indicates the major function of the code. The scripts either downloads data, manipulates data, conducts analyses, and/or produces a figure.

Outputs:Any figures and tables from preliminary analyses, but that are not used in the final manuscript, are saved in Outputs.

Figures:All final figures and tables are contained in the Figures folder. All figures are produced by scripts, except Figs. 1 and 2, which are schematics produced manually in a graphics editor. This folder contains two other folders: Supplemenatary Figures, and Partial Regression Plots. Partial Regression plots are the same as the final Figures 8-12, except they are grouped by taxa rather than by explanatory variable.

Data Processing Workflow - Overview: Tide Gauge Data (Data/Raw Data/Tide Gauge Data/individual stations) were downloaded using the NOAA Co-Ops API URL Builder (https://tidesandcurrents.noaa.gov/api-helper/url-generator.html), merged, then analyzed. Three MARINe data sets from the Coastal Biodiversity Survey (CBS) were accessed via data requests (https://marine.ucsc.edu/explore-the-data/contact/data-request-form.html). The first MARINe dataset (Data/Raw Data/MARINe Data/CBS_Percent Cover Data, both First Sample and Full Sample) was used to determine the top ten most abundant taxa (hereafter termed “dominant taxa”) across CBS survey sites during the monitoring period of 2001-01-01 to 2021-09-30. The second MARINe dataset (Data/Raw Data/MARINe Data/CBS_Elevation Data) was used to describe the upper limits of vertical distribution of dominant taxa through time. The third MARINe dataset (Data/Raw Data/MARINe Data/CBS_Presence Data) was used to visualize latitudinal distribution of taxa.

Location information for Tide Gauge Stations and CBS Survey Sites were assembled into a table (Data/Raw Data/CBS_Tide Gauge_Data.csv)

Tide Gauge Data were processed first, then MARINe Data. To replicate this workflow follow the steps described in the README file, in order.