knitr::opts_chunk$set(echo = F, warning = FALSE, fig.height=4, fig.width=7.5, fig.align="center", fig.pos='H')

require(tibble)
require(knitr)
require(DT)
# devtools::install_github("sbashevkin/deltareportr", "v0.2.0")
require(deltareportr)

create_dt <- function(x, caption=NULL){
  DT::datatable(x,
                rownames = FALSE,
                caption = caption,
                extensions = 'Buttons',
                options = list(dom = 'Blfrtip',
                               buttons = c('csv'),
                               lengthMenu = list(c(10,25,50,-1),
                                                 c(10,25,50,"All")),
                               columnDefs = list(list(className = 'dt-right', targets="_all"))))
}

End_year=2018

Data <- DeltaDater()
Data$Water_quality<-dplyr::filter(Data$Water_quality, Year<=2018)
WQ <- DeltaWQer(Data$Water_quality, End_year=End_year)
DF <- DeltaDayFlower(End_year=End_year)
PH <- DeltaPhyter(Data$Phytoplankton, End_year=End_year)
FISH <- DeltaSmelter(End_year=End_year)
META <- DeltaMetadater(Data)
ZOOP <- readRDS("ZOOP.Rds") # Created using most recent version of deltareportr, but with data filtered to 2018 or earlier. 
BIV <- DeltaBivalver(Data$Bivalves, End_year=End_year)

if(params$Map=="Yes"){
  DeltaMapper(Save=TRUE,
            Save_location="../vignettes/figures/map.png")
}


Accessibility

For any assistance accessing the map, graphs, images, or other material in this report, please contact us at accessibility@deltacouncil.ca.gov.

Erratum

A prior version of this report contained an incorrect zooplankton graph and dataset. This version of the report has now been corrected but the original report with the incorrect zooplankton data is here for reference.

Introduction

This report is intended for managers and scientists interested in tracking environmental conditions considered important for Delta Smelt (Hypomesus transpacificus). The variables shown in this report represent a selected subset of the variables that have been incorporated into conceptual models of Delta Smelt life history and analyses of the Delta Smelt population and its habitat. The variables selected have also been consistently measured as part of the long term monitoring efforts of the Interagency Ecological Program and other entities. Data for selected variables are summarized from 2002 to present, to reflect conditions since the beginning of the Pelagic Organism Decline. The report does not include data from short-term special studies. This report also does not provide detailed assessments of trends or the outcomes of specific management actions. Such assessments require specialized analysis and synthesis by expert teams and are beyond the scope of this report.

Delta Smelt is a pelagic species that prefers cooler temperatures (< 22 °C) in low salinity (0 - 6), turbid habitats. Some additional habitat concerns are food availability, invasive species, contaminants and harmful algal blooms.

Delta Smelt is generally described as a semi-anadromous species, moving from brackish low salinity (0.5-6 PSU) habitats to tidal fresh water habitats to spawn; however, the species expresses a variety of life histories. These spawning movements generally occur in the winter, after the “first-flush” of turbid freshwater following the first major precipitation event. A variable proportion of the population also exhibit resident life histories, residing in freshwater or brackish water for their entire life cycle.

Since the POD (Pelagic Organism Decline, 2002-2003), Delta Smelt have been most commonly found in the North Delta region (Cache Slough/Liberty Island and Sacramento River Deep Water Shipping Channel) during the winter spawning season. Larval and juvenile fish of the migratory life history move into low salinity habitat (Suisun Bay and Marsh, Lower Sacramento River) during the spring and early summer and remain there until their winter spawning movements into fresh water. The Lower San Joaquin River and Southern Delta are highly modified habitats characterized by warm summer temperatures and low turbidities with low abundance of Delta Smelt. Proximity to the water export facilities, which may result in some probability of entrainment under some environmental conditions and operations, may also be a factor.


(ref:DeltaSmeltphoto) Juvenile Delta Smelt Dale Kolke / DWR.

include_graphics("../vignettes/figures/DBK_Delta_Smelt-2463edit.jpg")

Methods

Most variables are divided among 7 of the 8 regions from the Enhanced Delta Smelt Monitoring (EDSM) 2018-19 phase I strata, with the Western Delta region excluded because it is almost never occupied by Delta Smelt. Missing data are denoted by vertical dashed lines. Data are plotted from 2002 (or the start of data collection) until present. The graphs presented in this report are for the season when a selected variable is most critical in the life history of Delta Smelt, as determined from literature. The graphs also highlight in red the most recent year, to facilitate visual comparisons of recent conditions with previous years. This report will be produced annually.

(ref:Deltamap) The Sacramento San Joaquin Delta divided into the EDSM 2018-19 phase I strata (Suisun Bay, Suisun Marsh, Lower Sacramento River, Sacramento Deep Water Shipping Channel, Cache Slough/Liberty Island, Lower San Joaquin River, and the Southern Delta).

include_graphics("../vignettes/figures/map.png")

Sampling programs

knitr::kable(tibble::tibble(Acronym = c("EMP", "FMWT", "SKT", "STN", "20mm", "EDSM", "Dayflow"), `Dataset name` = c("Environmental Monitoring Program", "Fall Midwater Trawl", "Spring Kodiak Trawl", "Summer Townet", "20 mm Survey", "Enhanced Delta Smelt Monitoring", "Dayflow"), Agency = c("California Department of Fish and Wildlife; Department of Water Resources", "California Department of Fish and Wildlife", "California Department of Fish and Wildlife", "California Department of Fish and Wildlife", "California Department of Fish and Wildlife", "United States Fish and Wildlife Service", "Department of Water Resources")), caption="Acronyms and responsible agencies for each dataset used in this report.")

Definitions of seasons

knitr::kable(tibble(Season = c("Winter", "Spring", "Summer", "Fall"), Months = c("December (of prior year), January, February", "March, April, May", "June, July, August", "September, October, November")), caption = "Definitions of seasons as used in this report.")

Sampling effort

(ref:Metadatafigure) Average yearly sampling effort for each variable, region, season, and survey. Cyanobacteria data are included in the phytoplankton category. Delta Outflow and X2 are derived from the Dayflow model while Delta Smelt abundance is represented by the official indices produced by FMWT, SKT, STN, 20mm, and EDSM. Raw data are available in table 3.

META$Plot

Abiotic drivers

Delta outflow

(ref:Outflow) Monthly modeled Delta outflow from the Dayflow model. Fall months are highlighted in orange. Raw data are available in table 4.

DF$Plots$Out

Delta Smelt abiotic habitat (as defined by salinity, turbidity, and water temperature) is determined largely by freshwater flow because Delta Smelt occupy the fresh, and low salinity portion of the SFE during the summer and fall. High Delta outflow in the summer-fall increases the amount of low salinity habitat available in Suisun Bay and Suisun Marsh. This is thought to be favorable for Delta Smelt. High outflow may also improve other abiotic environmental conditions and reduce harmful algal blooms.

X2

(ref:X2) Monthly modeled X2 from the Dayflow model. Fall months are highlighted in orange. Raw data are available in table 4.

DF$Plots$X2

X2 is a measure of the location of low salinity habitat. X2 is defined as the horizontal distance in kilometers from the Golden Gate up the axis of the estuary to where tidally averaged near-bottom salinity is 2. X2 is used as a proxy for the amount of low salinity habitat in the Delta; lower X2 in the fall generally means more low salinity habitat is available.

Summer temperature

(ref:Temperature) Mean summer water temperature from discrete measurements collected by the STN, EMP, and EDSM surveys, shaded by suitability for Delta Smelt (bottom blue section under 20 °C = good, middle yellow section between 20 and 22 °C = marginal, top red section over 22 °C = bad). Gray shading represents the standard deviation across all measurements from the summer of each year. Vertical dotted lines represent years with missing data. Raw data are available in table 5.

WQ$Temperature$Plot

Delta Smelt are sensitive to high water temperatures. Results from laboratory and observational studies suggest that Delta Smelt growth is optimal at water temperatures ≤ 20 °C, growth likely declines at 20-22 °C with poor growth at 22-24 °C, and the onset of physiological stress around 24 °C. Chronic exposure to water temperatures above 26 °C are likely fatal. Delta Smelt are rarely captured above 22 °C.

Fall Secchi depth

(ref:Secchi) Mean fall Secchi depth from the FMWT, EMP, and EDSM surveys. Means were calculated by pooling discrete measurements of the three surveys for the fall of each year. Gray shading represents the standard deviation across all measurements from the fall of each year. Raw data are available in table 6.

WQ$Secchi$Plot

Secchi depth is a measure of turbidity. Lower Secchi depth indicates higher turbidity, which is preferred by Delta Smelt. Delta Smelt distributions are clustered around areas with low Secchi depth (high turbidity) in the fall. Delta Smelt evolved in the historically turbid San Francisco Estuary and the larvae rely on this turbidity to effectively feed and hide from predators.

Fall salinity

(ref:Salinity) Mean fall salinity from the FMWT and EMP surveys. Means were calculated by pooling discrete measurements of the two surveys for the fall of each year. Gray shading represents the standard deviation across all measurements from the fall of each year. Raw data are available in table 7.

WQ$Salinity$Plot

During the fall, Delta Smelt prefer fresh to brackish water are most abundant in salinities less than about 6 PSU, decline in abundance at higher salinities, and become rare at salinities above about 14 PSU. Juveniles rely on freshwater and low salinity habitat in the fall to feed, grow, and mature into adults.

Biotic drivers

Summer chlorophyll concentrations

(ref:Chlorophyll) Mean summer chlorophyll concentrations from EMP. Means were calculated by pooling discrete measurements for the summer of each year. No data are available from the Sac Deep Water Shipping Channel or Cache/Slough/Liberty Island regions. Gray shading represents the standard deviation across all measurements from the summer of each year. Raw data are available in table 8.

WQ$Chlorophyll$Plot

Chlorophyll is a measure of productivity at the base of the food web. Higher chlorophyll is an indicator of more food availability for zooplankton, which are important prey for many fish, including Delta Smelt.

Summer phytoplankton composition

(ref:Phytoplankton) Summer phytoplankton community composition from EMP. Cyanobacteria are presented in figure 11. Means were calculated by pooling discrete measurements for the summer of each year. No data are available from the Sac Deep Water Shipping Channel or Cache/Slough/Liberty Island regions. Data are only shown from 2008 - 2018 to reflect the time frame of consistent sample processing methodology. Vertical dotted lines represent years with missing data. Raw data are available in table 9.

PH$Plots$Phytoplankton

Phytoplankton are the base of the aquatic food web. They provide food for zooplankton, which are important prey for Delta Smelt. Diatoms and cryptophytes are currently considered the best quality zooplankton food, although other taxa may be important as well.

Summer cyanobacteria

(ref:Cyanobacteria) Summer cyanobacteria abundance from EMP. Means were calculated by pooling discrete measurements for the summer of each year. No data are available from the Sac Deep Water Shipping Channel or Cache/Slough/Liberty Island regions. Data are only shown from 2008 - 2018 to reflect the time frame of consistent sample processing methodology. Vertical dotted lines represent years with missing data. Raw data are available in table 9.

PH$Plots$Cyanobacteria

Cyanobacteria are an unfavorable type of phytoplankton Cyanobacteria are considered poor quality zooplankton food and can produce toxins, such as microcystins.

Summer Microcystis

(ref:Microcystis) Summer Microcystis bloom intensity from EMP and STN. Means were calculated by pooling discrete measurements from the two surveys for the summer of each year. Microcystis bloom presence and intensity are measured on a qualitative scale with 5 categories: absent, low (widely scattered colonies), medium (adjacent colonies), high (contiguous colonies), and very high (concentration of contiguous colonies forming mats/scum). Vertical dotted lines represent years with missing data. Raw data are available in table 10.

WQ$Microcystis$Plot

Microcystis is a toxin-producing cyanobacteria. The toxins, microcystins, are harmful to human and animal health. Microcystis toxins can directly harm zooplankton populations which may affect food resources for Delta Smelt. Blooms occur yearly in the summer and fall.

Fall zooplankton

(ref:Zooplankton) Fall zooplankton biomass and composition from EMP. Means were calculated by pooling discrete measurements for the fall of each year. No data are available from the Sac Deep Water Shipping Channel or Cache/Slough/Liberty Island regions. Raw data are available in table 11.

ZOOP$Plot

Delta Smelt rely on zooplankton for food throughout their entire lives. Calanoid copepods and mysids are particularly important in the diet of Delta Smelt.

Fall invasive bivalve abundance

(ref:Bivalves) Fall abundance of two invasive bivalve species (the overbite clam Potamocorbula amurensis and the freshwater clam Corbicula fluminea) from EMP. Means were calculated by pooling discrete measurements for the fall of each year. No data are available from the Sac Deep Water Shipping Channel or Cache/Slough/Liberty Island regions. Raw data are available in table 12.

BIV$Plot

Invasive bivalves (clams) consume phytoplankton and zooplankton, reducing the amount of food available for fishes. The overbite clam Potamocorbula amurensis invaded the estuary in 1987, prefers brackish waters, and feeds on both zooplankton and phytoplankton. The freshwater clam Corbicula fluminea invaded sometime before 1945, prefers freshwater, and feeds primarily on phytoplankton but can also feed on organic material in the substrate when phytoplankton are rare.

Delta Smelt

IEP Delta Smelt index values

(ref:DeltasmeltIEP) Official Delta Smelt indices from four Interagency Ecological Program (IEP) surveys: Spring Kodiak Trawl (SKT), the Summer Townet Survey (STN), the 20mm survey, and the Fall Midwater Trawl (FMWT). The 20mm index could not be calculated in 2018 due to low catch. For all surveys, indices are calculated from selected set of stations that have been consistently sampled over the full time series of the survey. Raw data are available in table 13.

FISH$IEP$Plot

Delta Smelt abundance is estimated by 4 IEP surveys that target different life stages. The Spring Kodiak Trawl (SKT) estimates spawning adult Delta Smelt abundance from January through May, the Summer Townet Survey (STN) estimates juvenile Delta Smelt abundance in June, the 20mm survey estimates larval and juvenile Delta Smelt abundance from March through August, and the Fall Midwater Trawl (FMWT) estimates juvenile and adult Delta Smelt abundance from September through December.

EDSM Delta Smelt abundance

(ref:DeltasmeltEDSM) Monthly Delta Smelt abundance estimates (with 95% confidence intervals) from the Enhanced Delta Smelt Monitoring survey (EDSM). No Delta Smelt were collected from the Southern Delta. Vertical dotted lines represent months with missing data. Note that the y-axis is on the log scale and the x axis starts in mid 2017. The y axis starts at the lowest estimated detection value; all values below the limit are 0s. Raw data are available in table 14.

FISH$EDSM$Plot

EDSM is a new survey that calculates Delta Smelt abundance estimates throughout the year.

Recent publications

Data appendix

Detailed methods

Data collected from 2002 - present within our regional strata were included in this report. We used regional strata from the Enhanced Delta Smelt Monitoring (EDSM) 2018-19 phase I strata. The Western Delta region was excluded because it is almost never occupied by Delta Smelt. Temperature, Secchi depth, salinity, chlorophyll, phytoplankton, cyanobacteria, zooplankton, and bivalve data were first filtered to only include data within the season of interest, then annual means within each region were calculated. Monthly means were calculated for outflow and X2. The relative frequencies of Microcystis index values (1-5) were calculated for each year and region. Official Delta Smelt indices were used from the Fall Midwater Trawl (FMWT), Spring Kodiak Trawl (SKT), Summer Townet (STN), and 20-mm (20mm) surveys. Official abundance estimates from EDSM were averaged monthly within each region and confidence intervals were calculated based on the lognormal distribution as described in Polansky et al. (2019).

Data sources

Sources bolded below were used in the FLOAT MAST report. This section lists data sources available for all seasons. For data sources included in specific focal seasons, as well as sampling effort in each region, see figure 3.

Outflow

  1. Dayflow

X2

  1. Dayflow

Water Temperature

  1. EMP
  2. FMWT
  3. STN
  4. EDSM

Secchi depth

  1. EMP
  2. FMWT
  3. STN
  4. EDSM

Salinity

  1. EMP
  2. FMWT
  3. STN

Chlorophyll-a

  1. EMP

Phytoplankton

  1. EMP

Cyanobacteria

  1. EMP

Microcystis

  1. EMP

Zooplankton

  1. EMP

Bivalve biomass

  1. EMP

Delta Smelt Abundance

  1. FMWT
  2. SKT
  3. STN
  4. 20mm
  5. EDSM

Raw data

Sampling effort

cat("<table>",paste0("<caption>", "(#tab:Metadatadata)", "Data from <a href='#fig:Metadatafigure'> figure 3.</a>", "</caption>"),"</table>", sep ="\n")

create_dt(META$Data)

Outflow and X2

cat("<table>",paste0("<caption>", "(#tab:Dayflowdata)", "Data from <a href='#fig:Outflow'> figure 4</a> and <a href='#fig:X2'> figure 5.</a> Outflow is in ft<sup>3</sup>/s and X2 is in km.", "</caption>"),"</table>", sep ="\n")

create_dt(DF$Data)

Summer temperature

cat("<table>",paste0("<caption>", "(#tab:Temperaturedata)", "Data from <a href='#fig:Temperature'> figure 6.</a> Standard deviation and temperature are in &deg;C. Temperatures under 20 &deg;C are good, between 20 and 22 &deg;C are marginal, and over 22 &deg;C are bad for Delta Smelt.", "</caption>"),"</table>", sep ="\n")

create_dt(WQ$Temperature$Data)

Fall Secchi depth

cat("<table>",paste0("<caption>", "(#tab:Secchidata)", "Data from <a href='#fig:Secchi'> figure 7.</a> Standard deviation and Secchi depth are in cm.", "</caption>"),"</table>", sep ="\n")

create_dt(WQ$Secchi$Data)

Fall salinity

cat("<table>",paste0("<caption>", "(#tab:Salinitydata)", "Data from <a href='#fig:Salinity'> figure 8.</a>", "</caption>"),"</table>", sep ="\n")

create_dt(WQ$Salinity$Data)

Summer chlorophyll concentration

cat("<table>",paste0("<caption>", "(#tab:Chlorophylldata)", "Data from <a href='#fig:Chlorophyll'> figure 9.</a> Standard deviation and chlorophyll are in &micro;g/L.", "</caption>"),"</table>", sep ="\n")

create_dt(WQ$Chl$Data)

Summer phytoplankton and cyanobacteria

cat("<table>",paste0("<caption>", "(#tab:Phytodata)", "Data from <a href='#fig:Phytoplankton'> figure 10</a> and <a href='#fig:Cyanobacteria'> figure 11.</a> Count per unit effort represents the number of cells, colonies, or filaments / ml.", "</caption>"),"</table>", sep ="\n")

create_dt(PH$Data)

Summer Microcystis

cat("<table>",paste0("<caption>", "(#tab:Microcystisdata)", "Data from <a href='#fig:Microcystis'> figure 12.</a>", "</caption>"),"</table>", sep ="\n")

create_dt(WQ$Microcystis$Data)

Fall zooplankton

cat("<table>",paste0("<caption>", "(#tab:Zooplanktondata)", "Data from <a href='#fig:Zooplankton'> figure 13.</a> Biomass per unit effort is in &micro;g/m<sup>3</sup>.", "</caption>"),"</table>", sep ="\n")

create_dt(ZOOP$Data)

Fall invasive bivalve abundance

cat("<table>",paste0("<caption>", "(#tab:Bivalvedata)", "Data from <a href='#fig:Bivalves'> figure 14.</a> Count per unit effort represents the number of clams/m<sup>2</sup>.", "</caption>"),"</table>", sep ="\n")

create_dt(BIV$Data)

IEP Delta Smelt Index Values

cat("<table>",paste0("<caption>", "(#tab:DeltasmeltIEPdata)", "Data from <a href='#fig:DeltasmeltIEP'> figure 15.</a>", "</caption>"),"</table>", sep ="\n")

create_dt(FISH$IEP$Data)

EDSM Delta Smelt Index Values

cat("<table>",paste0("<caption>", "(#tab:DeltasmeltEDSMdata)", "Data from <a href='#fig:DeltasmeltEDSM'> figure 16.</a> CI indicates confidence intervals.", "</caption>"),"</table>", sep ="\n")

create_dt(FISH$EDSM$Data)

Report version

The current Git commit details are:

if (all(c("git2r", "here") %in% installed.packages()) & git2r::in_repository(path = ".")) git2r::repository(here::here())  


sbashevkin/deltareportr documentation built on Oct. 25, 2022, 2:03 p.m.