phytoplankton :: Ocean Carbon & Biogeochemistry

Studying marine ecosystems and biogeochemical cycles in the face of environmental change

Archive for phytoplankton

A new proxy for ocean iron bioavailability

Posted by mmaheigan

· Monday, July 26th, 2021

In many oceanic regions, iron exerts strong control on phytoplankton growth, ecosystem structure and carbon cycling. Yet, iron bioavailability and uptake rates by phytoplankton in the ocean are poorly constrained.

Recently, Shaked et al. (2020) (see GEOTRACES highlight), established a new approach for quantifying the availability of dissolved Fe (dFe) in natural seawater based on its uptake kinetics by Fe-limited cultured phytoplankton. In a follow up study published in GBC, this approach was extended to in situ phytoplankton, establishing a standardized proxy for dFe bioavailability in low-Fe ocean regions.

As explained in the short video lecture above, Yeala Shaked, Ben Twining, and their colleagues have analyzed large datasets collected during 10 research cruises (including 3 GEOTRACES section and process cruises) in multiple ocean regions. Dissolved Fe bioavailability was estimated through single cell Fe uptake rates, calculated by combining measured Fe contents of individual phytoplankton cells collected with concurrently-measured dFe concentrations, as well as modeled growth rates (Figure). Then the authors applied this proxy for: a) comparing dFe bioavailability among organisms and regions; b) calculating dFe uptake rates and residence times in low-Fe oceanic regions; and c) constraining Fe uptake parameters of earth system models to better predict ocean productivity in response to climate-change.

The data suggest that dFe species are highly available in low-Fe settings, likely due to photochemical reactions in sunlit waters.

Figure 1: The new bioavailability proxy (an uptake rate constant-kin-app) was calculated for ~1000 single cells from multiple ocean regions. For each cell, the iron quota was measured with synchrotron x-ray fluorescence (left panel), a growth rate was estimated from the PISCES model for the corresponding phytoplankton group (right panel), and the dissolved Fe concentration was measured concurrently (middle panel).

Authors:
Y. Shaked (Hebrew University and Interuniversity Institute for Marine Sciences)
B.S. Twining (Bigelow Lab)
A. Tagliabue (University of Liverpool)
M.T. Maldonado (University of British Columbia)
K.N. Buck (University of South Florida)
T. Mellett (University of South Florida)

References:
Shaked, Y., Twining, B. S., Tagliabue, A., & Maldonado, M. T. (2021). Probing the bioavailability of dissolved iron to marine eukaryotic phytoplankton using in situ single cell iron quotas. Global Biogeochemical Cycles, e2021GB006979. https://doi.org/10.1029/2021GB006979

Shaked, Y., Buck, K. N., Mellett, T., & Maldonado, M. T. (2020). Insights into the bioavailability of oceanic dissolved Fe from phytoplankton uptake kinetics. The ISME Journal, 1–12. https://doi.org/10.1038/s41396-020-0597-3

Joint highlight with GEOTRACES – read here.

Using BGC-Argo to obtain depth-resolved net primary production

Posted by mmaheigan

· Friday, July 23rd, 2021

Net primary production (NPP)—the organic carbon produced by the phytoplankton minus the organic carbon respired by phytoplankton themselves—serves as a major energy source of the marine ecosystem. Traditional methods for measuring NPP rely on ship-based discrete sampling and bottle incubations (e.g., ¹⁴C incubation), which introduce potential artifacts and limit the spatial and temporal data coverage of the global ocean. The global distribution of NPP has been estimated using satellite observations, but the satellite remote sensing approach cannot provide direct information at depth.

Figure 1. Panel A. Trajectories of 5 BGC-Argo and 1 SOS-Argo with the initial float deployment locations denoted by filled symbols. The dash-line at 47° N divided the research area into the northern (temperate) and southern (subtropical) regions. Stars indicate ship stations where ¹⁴C NPP values were measured during NAAMES cruises and compared with NPP from nearby Argo floats. Panels B and C. Monthly climatologies of net primary production (NPP, mmol m^-3 d^-1) profiles in the northern and southern regions of the research area, derived from BGC-Argo measurements using the PPM model. The shadings indicate one standard deviation. The red dotted line indicates mixed layer depth (MLD, m), and the yellow dashed line shows euphotic depth (Z_1%, m).

To fill this niche, a recent study in Journal of Geophysical Research: Biogeosciences, applied bio-optical measurements from Argo profiling floats to study the year-round depth-resolved NPP of the western North Atlantic Ocean (39° N to 54° N). The authors calculated NPP with two bio-optical models (Carbon-based Productivity Model, CbPM; and Photoacclimation Productivity Model, PPM). A comparison with NPP profiles from ¹⁴C incubation measurements showed advantages and limitations of both models. CbPM reproduced the magnitude of NPP in most cases, but had artifacts in the summer (a large NPP peak in the subsurface) due to the subsurface chlorophyll maximum caused by photoacclimation. PPM avoided the artifacts in the summer from photoacclimation, but the magnitude of PPM-derived NPP was smaller than the ¹⁴C result. Latitudinally varying NPP were observed, including higher winter NPP/lower summer NPP in the south, timing differences in NPP seasonal phenology, and different NPP depth distribution patterns in the summer months. With a 6-month record of concurrent oxygen and bio-optical measurements from two Argo floats, the authors also demonstrated the ability of Argo profiling floats to obtain estimates of the net community production (NCP) to NPP ratio (f-ratio), ranging from 0.3 in July to -1.0 in December 2016.

This work highlights the utility of float bio-optical profiles in comparison to traditional measurements and indicates that environmental conditions (e.g. light availability, nutrient supply) are major factors controlling the seasonality and spatial (horizontal and vertical) distributions of NPP in the western North Atlantic Ocean.

Authors:
Bo Yang (University of Virginia, UM CIMAS/NOAA AOML)
James Fox (Oregon State University)
Michael J. Behrenfeld (Oregon State University)
Emmanuel S. Boss (University of Maine)
Nils Haëntjens (University of Maine)
Kimberly H. Halsey (Oregon State University)
Steven R. Emerson (University of Washington)
Scott C. Doney (University of Virginia)

Wildfire impacts on coastal ocean phytoplankton

Posted by mmaheigan

· Wednesday, February 24th, 2021

Wildfire frequency, size, and destructiveness has increased over the last two decades, particularly in coastal regions such as Australia, Brazil, and the western United States. While the impact of fire on land, plants, and people is well documented, very few studies have been able to evaluate the impact of fires on ocean ecosystems. A serendipitously planned research cruise one week after the Thomas Fire broke out in California in December 2017 allowed the authors of this study and their colleagues to sample the adjacent Santa Barbara Channel during this devastating extreme fire event.

In a recent paper published in Journal of Geophysical Research: Oceans, the authors describe the phytoplankton community in the Santa Barbara Channel during the Thomas Fire. Phytoplankton community composition was described using a combination of images of phytoplankton from the Imaging FlowCytobot (McLane Labs) and phytoplankton pigments. Dinoflagellates were the dominant phytoplankton group in the surface ocean during the Thomas Fire, according to both methods (Figure 1).

Figure 1. (A) The fraction of total particle volume imaged by the Imaging FlowCytobot (IFCB) comprised of phytoplankton (green) and detritus (brown). Example IFCB images of ash (counted as part of detritus) particles are outlined in brown. (B) The phytoplankton fraction is then further divided by taxonomy, showing the abundance of nano-sized phytoplankton and especially dinoflagellates during the week of sampling. Example IFCB images of Gonyaulax (outlined in dark green), Prorocentrum (outlined in light green), and Umbilicosphaera (outlined in purple) cells are also shown.

While this study was not able to demonstrate a causal relationship between the Thomas Fire and the presence of dinoflagellates, this result is quite different from previous winters in the Santa Barbara Channel, when picophytoplankton and diatoms typically dominate the winter community. The incidence of dinoflagellates in the Santa Barbara Channel in December 2017 was correlated with the warmer-than-average water temperature during this study, which matched observations from other areas along the Central California coast that winter.

At the time this study was conducted, the Thomas Fire was the largest wildfire in California history. Since then, California fires have increased in danger, destruction, and human mortality; the Mendocino Fire complex (summer 2018) and five separate wildfires in summer 2020 exceeded the impacts of the Thomas Fire. With wildfire severity and frequency increasing not only in California but in coastal regions worldwide, this study gives an important first look at the impact of wildfire smoke and ash on oceanic primary productivity and community composition.

Authors:
Sasha Kramer (University of California Santa Barbara)
Kelsey Bisson (Oregon State University)
Alexis Fischer (University of California Santa Cruz)

Species loss alters ecosystem function in plankton communities

Posted by mmaheigan

· Monday, February 8th, 2021

Climate change impacts on the ocean such as warming, altered nutrient supply, and acidification will lead to significant rearrangement of phytoplankton communities, with the potential for some phytoplankton species to become extinct, especially at the regional level. This leads to the question: What are phytoplankton species’ redundancy levels from ecological and biogeochemical standpoints—i.e. will other species be able to fill the functional ecological and/or biogeochemical roles of the extinct species? Authors of a paper published recently in Global Change Biology explored these ideas using a global three-dimensional computer model with diverse planktonic communities, in which single phytoplankton types were partially or fully eliminated. Complex trophic interactions such as decreased abundance of a predator’s predator led to unexpected “ripples” through the community structure and in particular, reductions in carbon transfer to higher trophic levels. The impacts of changes in resource utilization extended to regions beyond where the phytoplankton type went extinct. Redundancy appeared lowest for types on the edges of trait space (e.g., smallest) or those with unique competitive strategies. These are responses that laboratory or field studies may not adequately capture. These results suggest that species losses could compound many of the already anticipated outcomes of changing climate in terms of productivity, trophic transfer, and restructuring of planktonic communities. The authors also suggest that a combination of modeling, field, and laboratory studies will be the best path forward for studying functional redundancy in phytoplankton.

Figure caption: Examples of the modelled ecological and biogeochemical responses to the extinction of different phytoplankton species.Figure caption: Examples of the modelled ecological and biogeochemical responses to the extinction of different phytoplankton species.

Authors:
Stephanie Dutkiewicz (Massachusetts Institute of Technology)
Philip W. Boyd (Institute for Marine and Antarctic Studies, University of Tasmania)
Ulf Riebesell (GEOMAR Helmholtz Centre for Ocean Research Kiel)

How environmental drivers regulated the long-term evolution of the biological pump

Posted by mmaheigan

· Friday, January 22nd, 2021

The marine biological pump (BP) plays a crucial role in regulating earth’s atmospheric oxygen and carbon dioxide levels by transferring carbon fixed by primary producers into the ocean interior and marine sediments, thereby controlling the habitability of our planet. The rise of multicellular life and eukaryotic algae in the ocean about 700 million years ago would likely have influenced the physical characteristics of oceanic aggregates (e.g., sinking rate), yet the magnitude of the impact this biological innovation had on the efficiency of BP is unknown.

Figure. 1. The impact of biological innovations (left) and environmental factors (atmospheric oxygen level and seawater temperature; right) on the efficiency of marine biological pump (BP). Temperatures are ocean surface temperatures (SST), and atmospheric pO2 is shown relative to the present atmospheric level (PAL). The BP efficiency is calculated as the fraction of carbon exported from the surface ocean that is delivered to the sediment-water interface. The results indicate that evolution of larger sized algae and zooplanktons has little influence on the long-term evolution of biological pump (left panel). The change in the atmospheric oxygen level and seawater surface temperature as environmental factors, on the other hand, have a stronger leverage on the efficiency of biological pump (right panel).

The authors of a recent paper in Nature Geoscience constructed a particle-based stochastic model to explore the change in the efficiency of the BP in response to biological and physical changes in the ocean over geologic time. The model calculates the age of organic particles in each aggregate based on their sinking rates, and considers the impact of primary producer cell size, aggregation, temperature, dust flux, biomineralization, ballasting by mineral phases, oxygen, and the fractal geometry (porosity) of aggregates. The model results demonstrate that while the rise of larger-sized eukaryotes led to an increase in the average sinking rate of oceanic aggregates, its impact on BP efficiency was minor. The evolution of zooplankton (with daily vertical migration in the water column) had a larger impact on the carbon transfer into the ocean interior. But results suggest that environmental factors most strongly affected the marine carbon pump efficiency. Specifically, increased ocean temperatures and greater atmospheric oxygen abundance led to a significant decrease in the efficiency of the BP. Cumulatively, these results suggest that while major biological innovations influenced the efficiency of BP, the long-term evolution of the marine carbon pump was primarily controlled by environmental drivers such as climate cooling and warming. By enhancing the rate of heterotrophic microbial degradation, our results suggest that the anthropogenically-driven global warming can result in a less efficient BP with reduced power of marine ecosystem in sequestering carbon from the atmosphere.

Authors:
Mojtaba Fakhraee (Yale University, Georgia Tech, and NASA Astrobiology Institute)
Noah J. Planavsky (Yale University, and NASA Astrobiology Institute)
Christopher T. Reinhard (Georgia Tech, and NASA Astrobiology Institute)

Water clarity impacts temperature and biogeochemistry in Chesapeake Bay

Posted by mmaheigan

· Thursday, December 3rd, 2020

Estuarine water clarity is determined by suspended materials in the water, including colored dissolved organic matter, phytoplankton, sediment, and detritus. These constituents directly affect temperature because when water is opaque, sunlight heats only the shallowest layers near the surface, but when water is clear, sunlight can penetrate deeper, warming the waters below the surface. Despite the importance of accurately predicting temperature variability, many numerical modeling studies do not adequately parameterize this fundamental relationship between water clarity and temperature.

In a recent study published in Estuaries and Coasts, the authors quantified the impact of a more realistic representation of water clarity in a hydrodynamic-biogeochemical model of the Chesapeake Bay by comparing two simulations: (1) water clarity is constant in space and time for the calculation of solar heating vs. (2) water clarity varies with modeled concentrations of light-attenuating materials. In the variable water clarity simulation (2), the water is more opaque, particularly in the northern region of the Bay. During the spring and summer months, the lower water clarity in the northern Bay is associated with warmer surface temperatures and colder bottom temperatures. Warmer surface temperatures encourage phytoplankton growth and nutrient uptake near the head of the Bay, thus fewer nutrients are transported downstream. These conditions are exacerbated during high-river flow years, when differences in temperature, nutrients, phytoplankton, and zooplankton extend further seaward.

Figure 1: Top row: Difference in the light attenuation coefficient for shortwave heating, kh[m-1] (variable minus constant light attenuation simulation). June, July, and August average for (A) 2001, (B) average of 2001-2005, and (C) 2003; difference in bottom temperatures [oC] (variable minus constant). Bottom row: Difference in June, July, and August average bottom temperature for (D) 2001, (E) average of 2001-2005, and (F) 2003. Data for 2001 are representative of low river discharge, and 2003 are representative high river discharge years.

This work demonstrates that a constant light attenuation scheme for heating calculations in coupled hydrodynamic-biogeochemical models underestimates temperature variability, both temporally and spatially. This is an important finding for researchers who use models to predict future temperature variability and associated impacts on biogeochemistry and species habitability.

Authors:
Grace E. Kim (NASA, Goddard Space Flight Center)
Pierre St-Laurent (VIMS, William & Mary)
Marjorie A.M. Friedrichs (VIMS, William & Mary)
Antonio Mannino (NASA, Goddard Space Flight Center)

Tiny phytoplankton seen from space

Posted by mmaheigan

· Thursday, November 19th, 2020

Picophytoplankton, the smallest phytoplankton on Earth, are dominant in over half of the global surface ocean, growing in low-nutrient “ocean deserts” where diatoms and other large phytoplankton have difficult to thrive. Despite their small size, picophytoplankton collectively account for well over 50% of primary production in oligotrophic waters, thus playing a major role in sustaining marine food webs.

In a recent paper published in Optics Express, the authors use satellite-detected ocean color (namely remote-sensing reflectance, R_rs(λ)) and sea surface temperature to estimate the abundance of the three picophytoplankton groups—the cyanobacteria Prochlorococcus and Synechococcus, and autotrophic picoeukaryotes. The authors analysed R_rs(λ) spectra using principal component analysis, and principal component scores and SST were used in the predictive models. Then, they trained and independently evaluated the models with in-situ data from the Atlantic Ocean (Atlantic Meridional Transect cruises). This approach allows for the satellite detection of the succession of species across ocean oligotrophic ecosystem boundaries, where these cells are most abundant (Figure 1).

Figure 1. Cell abundances of the three major picophytoplankton groups (the cyanobacteria Prochlorococcus and Synechococcus, and a collective group of autotrophic picoeukaryotes) in surface waters of the Atlantic Ocean. Abundances are shown for the dominant group in terms of total biovolume (converted from cell abundance).

Since these organisms can be used as proxies for marine ecosystem boundaries, this method can be used in studies of climate and ecosystem change, as it allows a synoptic observation of changes in picophytoplankton distributions over time and space. For exploring spectral features in hyperspectral R_rs(λ) data, the implementation of this model using data from future hyperspectral satellite instruments such as NASA PACE’s Ocean Color Instrument (OCI) will extend our knowledge about the distribution of these ecologically relevant phytoplankton taxa. These observations are crucial for broad comprehension of the effects of climate change in the expansion or shifts in ocean ecosystems.

Authors:
Priscila K. Lange (NASA Goddard Space Flight Center / Universities Space Research Association / Blue Marble Space Institute of Science)
Jeremy Werdell (NASA Goddard Space Flight Center)
Zachary K. Erickson (NASA Goddard Space Flight Center)
Giorgio Dall’Olmo (Plymouth Marine Laboratory)
Robert J. W. Brewin (University of Exeter)
Mikhail V. Zubkov (Scottish Association for Marine Science)
Glen A. Tarran (Plymouth Marine Laboratory)
Heather A. Bouman (University of Oxford)
Wayne H. Slade (Sequoia Scientific, Inc)
Susanne E. Craig (NASA Goddard Space Flight Center / Universities Space Research Association)
Nicole J. Poulton (Bigelow Laboratory for Ocean Sciences)
Astrid Bracher (Alfred-Wegener-Institute Helmholtz Center for Polar and Marine Research / University of Bremen)
Michael W. Lomas (Bigelow Laboratory for Ocean Sciences)
Ivona Cetinić (NASA Goddard Space Flight Center / Universities Space Research Association)

Marine heatwave implications for future phytoplankton blooms

Posted by mmaheigan

· Thursday, October 15th, 2020

Ocean temperature extreme events such as marine heatwaves are expected to intensify in coming decades due to anthropogenic warming. Although the effects of marine heatwaves on large plants and animals are becoming well documented, little is known about how these warming events will impact microbes that regulate key biogeochemical processes such as ocean carbon uptake and export, which represent important feedbacks on the global carbon cycle and climate.

Figure caption: Relationship between phytoplankton bloom response to marine heatwaves and background nitrate concentration in the 23 study regions. X-axis denotes the annual-mean sea-surface nitrate concentration based on the model simulation (1992-2014; OFAM3, blue) and the in situ climatology (WOA13, orange). Y-axis denotes the mean standardised anomalies (see Equation 1 of the paper) of simulated sea-surface phytoplankton nitrogen biomass (1992-2014; OFAM3, blue) and observed sea-surface chlorophyll a concentration (2002-2018; MODIS, orange) during the co-occurrence of phytoplankton blooms and marine heatwaves.

In a recent study published in Global Change Biology, authors combined model simulations and satellite observations in tropical and temperate oceanographic regions over recent decades to characterize marine heatwave impacts on phytoplankton blooms. The results reveal regionally‐coherent anomalies depicted by shallower surface mixed layers and lower surface nitrate concentrations during marine heatwaves, which counteract known light and nutrient limitation effects on phytoplankton growth, respectively (Figure 1). Consequently, phytoplankton bloom responses are mixed, but derive from the background nutrient conditions of a study region such that blooms are weaker (stronger) during marine heatwaves in nutrient-poor (nutrient-rich) waters.

Given the projected expansion of nutrient-poor waters in the 21st century ocean, the coming decades are likely to see an increased occurrence of weaker blooms during marine heatwaves, with implications for higher trophic levels and biogeochemical cycling of key elements.

Authors:
Hakase Hayashida (University of Tasmania)
Richard Matear (CSIRO)
Pete Strutton (University of Tasmania)

Profiling floats reveal fate of Southern Ocean phytoplankton stocks

Posted by mmaheigan

· Tuesday, September 1st, 2020

More observations are needed to constrain the relative roles of physical (advection), biogeochemical (downward export), and ecological (grazing and biological losses) processes in driving the fate of phytoplankton blooms in Southern Ocean waters. In a recent paper published in Nature Communications, authors used seven Biogeochemical Argo (BGC-Argo) floats that vertically profiled the upper ocean every ten days as they drifted for three years across the remote Sea Ice Zone of the Southern Ocean. Using the floats’ biogeochemical sensors (chlorophyll, nitrate, and backscattering) and regional ratios of nitrate consumption:chlorophyll synthesis, the authors developed a new approach to remotely estimate the fate of the phytoplankton stocks, enabling calculations of herbivory and of downward carbon export. The study revealed that the major fate of phytoplankton biomass in this region is grazing, which consumes ~90% of stocks. The remaining 10% is exported to depth. This pattern was consistent throughout the entire sea ice zone where the floats drifted, from 60°-69° South.

Figure Caption: Southern Ocean Chlorophyll a climatology and floats’ trajectories (top panel). Total losses of Chlorophyll a (including grazing and phytodetritus export, left panel). Phytodetritus export (right panel).

This study region comprises two of the three major krill growth and development areas—the eastern Weddell and King Haakon VII Seas and Prydz Bay and the Kerguelen Plateau—so the observed grazing was probably due to Antarctic krill, underscoring their pivotal importance in this ecosystem. Building upon the greater understanding of ocean ecosystems via satellite ocean colour development in the 1990s, BGC-Argo floats and this new approach will allow remote monitoring of the different fates of phytoplankton stocks and insights into the status of the ecosystem.

Authors:
Sebastien Moreau (Norwegian Polar Institute, Tromsø, Norway)
Philip Boyd (Institute for Marine and Antarctic Studies, Hobart, Australia)
Peter Strutton (Institute for Marine and Antarctic Studies, Hobart, Australia)

A close-up view of biomass controls in Southern Ocean eddies

Posted by mmaheigan

· Thursday, August 20th, 2020

Southern Ocean biological productivity is instrumental in regulating the global carbon cycle. Previous correlative studies associated widespread mesoscale activity with anomalous chlorophyll levels. However, eddies simultaneously modify both the physical and biogeochemical environments via several competing pathways, making it difficult to discern which mechanisms are responsible for the observed biological anomalies within them. Two recently published papers track Southern Ocean eddies in a global, eddy-resolving, 3-D ocean simulation. By closely examining eddy-induced perturbations to phytoplankton populations, the authors are able to explicitly link eddies to co-located biological anomalies through an underlying mechanistic framework.

Figure caption: Simulated Southern Ocean eddies modify phytoplankton division rates in different directions of depending on the polarity of the eddy and background seasonal conditions. During summer anticyclones (top right panel) deliver extra iron from depth via eddy-induced Ekman pumping and fuel faster phytoplankton division rates. During winter (bottom right panel) the extra iron supply is eclipsed by deeper mixed layer depths and elevated light limitation resulting in slower division rates. The opposite occurs in cyclones.

In the first paper, the authors observe that eddies primarily affect phytoplankton division rates by modifying the supply of iron via eddy-induced Ekman pumping. This results in elevated iron and faster phytoplankton division rates in anticyclones throughout most of the year. However, during deep mixing winter periods, exacerbated light stress driven by anomalously deep mixing in anticyclones can dominate elevated iron and drive division rates down. The opposite response occurs in cyclones.

The second paper tracks how eddy-modified division rates combine with eddy-modified loss rates and physical transport to produce anomalous biomass accumulation. The biomass anomaly is highly variable, but can exhibit an intense seasonal cycle, in which cyclones and anticyclones consistently modify biomass in different directions. This cycle is most apparent in the South Pacific sector of the Antarctic Circumpolar Current, a deep mixing region where the largest biomass anomalies are driven by biological mechanisms rather than lateral transport mechanisms such as eddy stirring or propagation.

It is important to remember that the correlation between chlorophyll and eddy activity observable from space can result from a variety of physical and biological mechanisms. Understanding the nuances of how these mechanisms change regionally and seasonally is integral in both scaling up local observations and parameterizing coarser, non-eddy resolving general circulation models with embedded biogeochemistry.

Authors:
Tyler Rohr (Australian Antarctic Partnership Program, previously at MIT/WHOI)
Cheryl Harrison (University of Texas Rio Grande Valley)
Matthew Long (National Center for Atmospheric Research)
Peter Gaube (University of Washington)
Scott Doney (University of Virginia)

Funding for the Ocean Carbon & Biogeochemistry Project Office is provided by the National Science Foundation (NSF) and the National Aeronautics and Space Administration (NASA). The OCB Project Office is housed at the Woods Hole Oceanographic Institution.

Archive for phytoplankton

Filter by Keyword