Introduction

The planktonr package provides access to plankton and biogeochemical data collected by the Integrated Marine Observing System (IMOS) around Australia. This vignette introduces the main data sources and demonstrates how to visualise sampling effort and taxonomic coverage.

Data Sources

National Reference Stations (NRS)

The NRS network consists of nine coastal monitoring stations around Australia that collect plankton, biogeochemical, and physical oceanographic data. Most stations sample monthly using:

  • Niskin bottles - For phytoplankton (0.5-5 L samples from 0-50 m depths)
  • Vertical net tows - For zooplankton (100 µm mesh, 0-50 m integrated samples)

Stations include: Darwin (DAR), Yongala (YON), North Stradbroke Island (NSI), Port Hacking (PHB), Maria Island (MAI), Rottnest Island (ROT), Esperance (ESP), Ningaloo (NIN), and Kangaroo Island (KAI).

Some stations have limited sampling: * Darwin samples quarterly to capture monsoonal variability * Kangaroo Island samples 5-6 times per year due to weather constraints * Esperance and Ningaloo were discontinued after initial sampling periods

Continuous Plankton Recorder (CPR)

The CPR is a high-speed plankton sampler towed behind ships of opportunity at ~10 m depth. It samples plankton continuously, with each sample representing approximately 3 nautical miles (~5.6 km).

The CPR primarily samples in: * The Great Barrier Reef region * The East Australian Current region * Bass Strait and the South-east region

CPR sampling is quarterly in main areas of interest and more ad-hoc in other regions. Unlike NRS samples which are collected during daylight hours, CPR samples are collected day and night as ships operate continuously.

Southern Ocean Time Series (SOTS)

The SOTS mooring south of Tasmania provides high-frequency physical and biogeochemical measurements, with intermittent biological sampling similar to NRS protocols.

Examine the frequency of sampling

The sampling at the National Reference Stations is predominantly monthly. Darwin samples on a quarterly basis and completes multiple trips across the tidal cycle. Kangaroo Island sampling is restricted to 5-6 times per year by sea conditions.

The Continuous Plankton Recorder samples quarterly through our main areas of interest (the Great Barrier Reef and the East Australian Current) and on a more ad-hoc basis in other regions, depending on ship availability.

datCPR <- pr_get_CPRTrips()
datNRS <- pr_get_NRSTrips()

ggCPR <- planktonr::pr_plot_Gantt(datCPR)
ggNRS <- planktonr::pr_plot_Gantt(datNRS)

wrap_plots(ggCPR, ggNRS, ncol = 1) &
  theme(text = ggplot2::element_text(size = 12, face = "bold"))

Plot the spatial distribution of NRS and CPR samples

The National Reference Stations are situated throughout the coastal region of Australia, with most Marine Bioregions having at least one station. The Western Australian stations at Ningaloo and Esperance were discontinued after an initial sampling period due to logistical constraints.

The CPR tracks are also more sparse in Western Australia as shipping options are limited. CPR sampling is concentrated along major shipping routes, particularly the east coast and Bass Strait.

df <- pr_get_ProgressMapData(Survey = c("NRS", "CPR"))

pr_plot_ProgressMap(df)

What functional groups do we commonly see?

We can learn a lot about the community structure of plankton by examining functional groups. These pie charts also show how sampling methods impact what we find in the samples.

Functional groups represent taxa with similar ecological roles or morphological characteristics. For example: * Phytoplankton: Centric diatoms, pennate diatoms, dinoflagellates, cyanobacteria * Zooplankton: Copepods, appendicularians, chaetognaths, molluscs

Phytoplankton

Centric diatoms dominate in the NRS samples, whilst pennate diatoms become more important in the CPR samples. This likely reflects differences in sampling depth and method:

  • NRS: Niskin bottle samples from discrete depths (0-50 m), collected in daylight
  • CPR: Continuous sampling at ~10 m depth, collected day and night
p1 <- pr_plot_PieFG(pr_get_FuncGroups(Survey = "NRS", Type = "Phytoplankton"))
p2 <- pr_plot_PieFG(pr_get_FuncGroups(Survey = "CPR", Type = "Phytoplankton"))

(wrap_plots(p1, p2, guides = "collect", ncol = 2) &
    theme(text = element_text(size = 16), legend.position = "bottom", 
          plot.title = element_text(face = "bold")))

Zooplankton

Copepods dominate both the CPR and the NRS net samples, which is typical for marine zooplankton communities. However, the relative abundance of other groups differs:

  • Chaetognaths (arrow worms) and cladocerans (water fleas) are sampled more effectively with the CPR
  • Molluscs (primarily pteropods) are caught more frequently in the NRS net samples

These differences reflect the different mesh sizes and sampling methods: * NRS: 100 µm mesh net, vertical tows (0-50 m) * CPR: Silk mesh (~270 µm), horizontal tow at ~10 m

p3 <- pr_plot_PieFG(pr_get_FuncGroups(Survey = "NRS", Type = "Zooplankton"))
p4 <- pr_plot_PieFG(pr_get_FuncGroups(Survey = "CPR", Type = "Zooplankton"))

(wrap_plots(p3, p4, guides = "collect", ncol = 2) &
    theme(text = element_text(size = 16), legend.position = "bottom", 
          plot.title = element_text(face = "bold")))

How many taxa are we identifying?

Over the course of the surveys since 2009, taxonomic resolution has improved through analyst training and exposure to diverse assemblages. We are now identifying:

  • >600 phytoplankton taxa from the CPR
  • >1000 zooplankton taxa from the CPR, including life stage and sex for many taxa

Zooplankton diversity is generally higher than phytoplankton diversity because: * Zooplankton are larger and morphological differences between species are more definitive * Zooplankton identification includes life stages (nauplii, copepodites, adults) and sex * Zooplankton sampling captures a wider size range of organisms

Phytoplankton

These accumulation curves show the cumulative number of taxa identified over time as sampling effort increases. The curves are not yet asymptotic, suggesting more taxa will be discovered with continued sampling.

PSpNRSAccum <- pr_get_TaxaAccum(Survey = "NRS", Type = "Phytoplankton")
PSpCPRAccum <- pr_get_TaxaAccum(Survey = "CPR", Type = "Phytoplankton")

p1 <- pr_plot_TaxaAccum(PSpNRSAccum) 
p2 <- pr_plot_TaxaAccum(PSpCPRAccum)


wrap_plots(p1, p2, ncol = 2) &
  theme(text = element_text(size = 16, face = "bold"),
        aspect.ratio = 1)

Zooplankton

The higher taxonomic diversity in zooplankton is evident in these accumulation curves. The CPR has now identified over 1000 zooplankton taxa, while the NRS continues to add new taxa regularly.

ZSpNRSAccum <- pr_get_TaxaAccum(Survey = "NRS", Type = "Zooplankton")
ZSpCPRAccum <- pr_get_TaxaAccum(Survey = "CPR", Type = "Zooplankton")

p3 <- pr_plot_TaxaAccum(ZSpNRSAccum)
p4 <- pr_plot_TaxaAccum(ZSpCPRAccum)

wrap_plots(p3, p4, ncol = 2) &
  theme(text = element_text(size = 16, face = "bold"),
        aspect.ratio = 1)