The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed. populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (similar to 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

Mesopelagic fish abundance and particle size spectrum analysis from Maria S. Merian cruise MSM26, spring 2013

One of the goals of EU BASIN is to understand variability in production across the Atlantic and the impact of this variability on higher trophic levels. One aspect of these investigations is to examine the biomes defined by Longhurst (2007). These biomes are largely based on productivity measured with remote sensing. During MSM 26, mesopelagic fish and size-spectrum data were collected to test the biome classifications of the north Atlantic. In most marine systems, the size-spectrum is a decay function with more, smaller organisms and fewer larger organisms. The intercept of the size-spectrum has been linked to overall productivity while the slope represents the "rate of decay" of this productivity (Zhou 2006, doi:10.1093/plankt/fbi119). A Laser In-Situ Scattering Transmissometer was used to collect size-spectrum data and net collections were made to capture mesopelagic fish. The relationship among the mesopelagic fish size and abundance distributions will be compared to the estimates of production from the size-spectrum data to evaluate the biomes of the stations occupied during MSM 26.

How optimal life history changes with the community size-spectrum

This paper derives optimal life histories for fishes or other animals in relation to the size spectrum of the ecological community in which they are both predators and prey. Assuming log-linear size-spectra and well known scaling laws for feeding and mortality, we first construct the energetics of the individual. From these we find, using dynamic programming, the optimal allocation of energy between growth and reproduction as well as the trade-off between offspring size and numbers. Optimal strategies were found to be strongly dependent on size spectrum slope. For steep size spectra (numbers declining rapidly with size), determinate growth was optimal and allocation to somatic growth increased rapidly with increasing slope. However, restricting reproduction to a fixed mating season changed optimal allocations to give indeterminate growth approximating a von Bertalanffy trajectory. The optimal offspring size was as small as possible given other restrictions such as newborn starvation mortality. For shallow size spectra, finite optimal maturity size required a decline in fitness for large size or age. All the results are compared with observed size spectra of fish communities to show their consistency and relevance.

The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (~ 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

Particle size spectrum analysis during Maria S. Merian cruise MSM26, spring 2013

A Laser In-Situ Scattering Transmissometer (LISST) was used to collect vertical distribution data of particles from 2.5 to 500 µm in size. The LISST uses a multi-ring detector to measure scattering light of particles from a laser diode. Particles are classified into 32 log-spaced bins and the concentration of each bin is calculated as micro-liters per liter (µl/l). The instrument is rated to a depth of 300 m, and also records temperature and pressure. The sample interval was set to record every second. The LISST was attached to the LOPC frame to conduct casts and allow for particle-size comparisons between the two instruments. The LOPC is rated to a depth of 2000 m, thus a short deployment to a depth of 300 m was first conducted with both instruments. The instruments were then returned to the deck and the LISST removed via a quick release bracket so deep LOPC casts could be continued at a station. Raw LISST size-spectrum data is presented as concentrations for each of the 32 size bins for every second of the cast.

Emerging asymmetric interactions between forage and predator fisheries impose management trade-offs

A size and trait-based marine community model was used to investigate interactions, with potential implications for yields, when a fishery targeting forage fish species (whose main adult diet is zooplankton) co-occurs with a fishery targeting larger-sized predator species. Predicted effects on the size structure of the fish community, growth and recruitment of fishes, and yield from the fisheries were used to identify management trade-offs among the different fisheries. Results showed that moderate fishing on forage fishes imposed only small effects on predator fisheries, whereas predator fisheries could enhance yield from forage fisheries under some circumstances.

Pattern of benthic biomass size spectra from shallow waters in the East China Seas

Benthic biomass size spectra (BSS) and normalized biomass size spectra were constructed, and benthic secondary production was estimated by a size spectrum equation in the shallow waters in the East China Sea, ranging latitudinally from 40A degrees N to 29A degrees N. The BSS patterns were bimodal, two biomass peaks corresponding to meiofauna and macrofauna, respectively, separated by a trough of low biomass at 8-256 mu g individual dry weight which varied in position with median sediment particle size. The BSS also displayed bimodality within meiofauna size ranges, which in most stations was due to the relative proportions of nematodes and other meiofauna taxa. Re-analysis of data from sites in the UK, South Africa, and Antarctic showed a similar bimodality in the adult species body size distribution within the meiofauna size range. Macrofaunal production estimated by the size spectrum equation was very similar to the results of Brey90 empirical equation. However, these production values were much lower than those calculated by Brey01. Different individual dry-to-wet conversion ratios, temperature deviation, and macrofauna taxonomic composition might be responsible for the between-model differences. The macrofaunal P/B ratios calculated by this equation ranged from 0.3 to 3.4 which were in accordance with values from Northern Hemisphere mid-latitudes. Meiofaunal production estimates will need further empirical support.

Zooplankton functional groups on the continental shelf of the yellow sea

Zooplankton plays a vital role in marine ecosystems. Variations in the zooplankton species composition, biomass, and secondary production will change the structure and function of the ecosystem. How to describe this process and make it easier to be modeled in the Yellow Sea ecosystem is the main purpose of this paper. The zooplankton functional groups approach, which is considered a good method of linking the structure of food webs and the energy flow in the ecosystems, is used to describe the main contributors of secondary produciton of the Yellow Sea ecosystem. The zooplankton can be classified into six functional groups: giant crustaceans, large copepods, small copepods, chaetognaths, medusae, and salps. The giant crustaceans, large copepods, and small copepods groups, which are the main food resources for fish, are defined depending on the size spectrum. Medusae and chaetognaths are the two gelatinous carnivorous groups, which compete with fish for food. The salps group, acting as passive filter-feeders, competes with other species feeding on phytoplankton, but their energy could not be efficiently transferred to higher trophic levels. From the viewpoint of biomass, which is the basis of the food web, and feeding activities, the contributions of each functional group to the ecosystem were evaluated; the seasonal variations, geographical distribution patterns, and species composition of each functional group were analyzed. The average zooplankton biomass was 2.1 g dry wt m(-2) in spring, to which the giant crustaceans, large copepods, and small copepods contributed 19, 44, and 26%, respectively. High biomasses of the large copepods and small copepods were distributed at the coastal waters, while the giant crustaceans were mainly located at offshore area. In summer, the mean biomass was 3.1 g dry wt m(-2), which was mostly contributed by the giant crustaceans (73%), and high biomasses of the giant crustaceans, large copepods, and small copepods were all distributed in the central part of the Yellow Sea. During autumn, the mean biomass was 1.8 g dry wt m(-2), which was similarly constituted by the giant crustaceans, large copepods, and small copepods (36, 33, and 23%, respectively), and high biomasses of the giant crustaceans and large copepods occurred in the central part of the Yellow Sea, while the small copepods were mainly located at offshore stations. The giant crustaceans and large copepods dominated the zooplankton biomass (2.9 g dry wt m(-2)) in winter, contributing respectively 57 and 27%, and they, as well as the small copepods, were all mainly located in the central part of the Yellow Sea. The chaetognaths group was mainly located in the northern part of the Yellow Sea during all seasons, but contributed less to the biomass compared with the other groups. The medusae and salps groups were distributed unevenly, with sporadic dynamics, mainly along the coastline and at the northern part of the Yellow Sea. No more than 10 species belonging to the respective functional groups dominated the zooplankton biomass and controlled the dynamics of the zooplankton community. The clear picture of the seasonal and spatial variations of each zooplankton functional group makes the complicated Yellow Sea ecosystem easier to be understood and modeled. (C) 2010 Elsevier Ltd. All rights reserved.

How community ecology links natural mortality, growth, and production of fish populations

Size-spectrum theory is used to show that (i) predation mortality is a decreasing function of individual size and proportional to the consumption rate of predators; (ii) adult natural mortality M is proportional to the von Bertalanffy growth constant K; and (iii) productivity rate P/B is proportional to the asymptotic weight W8 -1/3. The constants of proportionality are specified using individual level parameters related to physiology or prey encounter. The derivations demonstrate how traditional fisheries theory can be connected to community ecology. Implications for the use of models for ecosystem-based fisheries management are discussed.

Factors influencing zooplankton size structure at contrasting temperatures in coastal shallow lakes: Implications for effects of climate change

We assessed the importance of temperature, salinity, and predation for the size structure of zooplankton and provided insight into the future ecological structure and function of shallow lakes in a warmer climate. Artificial plants were introduced in eight comparable coastal shallow brackish lakes located at two contrasting temperatures: cold-temperate and Mediterranean climate region. Zooplankton, fish, and macroinvertebrates were sampled within the plants and at open-water habitats. The fish communities of these brackish lakes were characterized by small-sized individuals, highly associated with submerged plants. Overall, higher densities of small planktivorous fish were recorded in the Mediterranean compared to the cold-temperate region, likely reflecting temperature-related differences as have been observed in freshwater lakes. Our results suggest that fish predation is the major control of zooplankton size structure in brackish lakes, since fish density was related to a decrease in mean body size and density of zooplankton and this was reflected in a unimodal shaped biomass-size spectrum with dominance of small sizes and low size diversity. Salinity might play a more indirect role by shaping zooplankton communities toward more salt-tolerant species. In a global-warming perspective, these results suggest that changes in the trophic structure of shallow lakes in temperate regions might be expected as a result of the warmer temperatures and the potentially associated increases in salinity. The decrease in the density of largebodied zooplankton might reduce the grazing on phytoplankton and thus the chances of maintaining the clear water state in these ecosystems

Testing the influence of small crystals on ice size spectra using Doppler lidar observations

Analysis of the vertical velocity of ice crystals observed with a 1.5micron Doppler lidar from a continuous sample of stratiform ice clouds over 17 months show that the distribution of Doppler velocity varies strongly with temperature, with mean velocities of 0.2m/s at -40C, increasing to 0.6m/s at -10C due to particle growth and broadening of the size spectrum. We examine the likely influence of crystals smaller than 60microns by forward modelling their effect on the area-weighted fall speed, and comparing the results to the lidar observations. The comparison strongly suggests that the concentration of small crystals in most clouds is much lower than measured in-situ by some cloud droplet probes. We argue that the discrepancy is likely due to shattering of large crystals on the probe inlet, and that numerous small particles should not be included in numerical weather and climate model parameterizations.

Estimating drizzle drop size and precipitation rate using two-colour lidar measurements

A method to estimate the size and liquid water content of drizzle drops using lidar measurements at two wavelengths is described. The method exploits the differential absorption of infrared light by liquid water at 905 nm and 1.5 μm, which leads to a different backscatter cross section for water drops larger than ≈50 μm. The ratio of backscatter measured from drizzle samples below cloud base at these two wavelengths (the colour ratio) provides a measure of the median volume drop diameter D0. This is a strong effect: for D0=200 μm, a colour ratio of ≈6 dB is predicted. Once D0 is known, the measured backscatter at 905 nm can be used to calculate the liquid water content (LWC) and other moments of the drizzle drop distribution. The method is applied to observations of drizzle falling from stratocumulus and stratus clouds. High resolution (32 s, 36 m) profiles of D0, LWC and precipitation rate R are derived. The main sources of error in the technique are the need to assume a value for the dispersion parameter μ in the drop size spectrum (leading to at most a 35% error in R) and the influence of aerosol returns on the retrieval (≈10% error in R for the cases considered here). Radar reflectivities are also computed from the lidar data, and compared to independent measurements from a colocated cloud radar, offering independent validation of the derived drop size distributions.

Geochemical and physical sediment properties of sediment cores from the Portuguese shelf

Marine sediments from the Portuguese shelf are influenced by environmental changes in the surrounding continental and marine environment. These are largely controlled by the North Atlantic Oscillation, but additional impacts may arise from episodic tsunamis. In order to investigate these influences, a high resolution multi-proxy study has been carried out on a 5.4 m long gravity core and five box cores from the Tagus prodelta on the western Portuguese margin, incorporating geochemical (Corg/Ntotal ratios, d13Corg, d15N, d18O, Corg and CaCO3 content) and physical sediment properties (magnetic susceptibility, grain-size). Subsurface data of the five box cores indicate no major effect of early postdepositional alteration. Surface data show a higher fraction of terrigenous organic material close to the river mouth and in the southern prodelta. Gravity core GeoB 8903 covers the last 3.2 kyrs with a temporal resolution of at least 0.1 cm/yr. Very high sedimentation rates between 69 and 140 cm core depth indicate a possible disturbance of the record by the AD1755 tsunami, although no evidence for a disturbance is observed in the data. Sea surface temperature and salinity on the prodelta, the local budget of marine NO3- as well as the provenance of organic matter remained virtually constant during the past 3.2 kyrs. A positive correlation between magnetic susceptibility (MS) and North Atlantic Oscillation (NAO) is evident for the past 250 years, coinciding with a negative correlation between mean grain-size and NAO. This is assigned to a constant riverine supply of fine material with high MS, which is diluted by the riverine input of a coarser, low-MS component during NAO negative, high-precipitation phases. End-member modelling of the lithic grain-size spectrum supports this, revealing a third, coarse lithic component. The high abundance of this coarse end-member prior to 2 kyr BP is interpreted as the result of stronger bottom currents, concentrating the coarse sediment fraction by winnowing. As continental climate was more arid prior to 2 kyr BP (Subboreal), the coarse end-member may also consist of dust from local sources. A decrease in grain-size and CaCO3 content after 2 kyr BP is interpreted as a result of decreasing wind strength. The onset of a fining trend and a further decrease in CaCO3 around AD900 occurs simultaneous to climatic variations, reconstructed from eastern North Atlantic records. A strong increase in MS between AD1400 and AD1500 indicates higher lithic terrigenous input, caused by deforestation in the hinterland.

Particle camera and settling chamber measurements at site PC (GeoB11834-3) and SC (GeoB12914-4) off Cape Blanc/Mauritania

Particles sinking out of the euphotic zone are important vehicles of carbon export from the surface ocean. Most of the particles produce heavier aggregates by coagulating with each other before they sink. We implemented an aggregation model into the biogeochemical model of Regional Oceanic Modelling System (ROMS) to simulate the distribution of particles in the water column and their downward transport in the Northwest African upwelling region. Accompanying settling chamber, sediment trap and particle camera measurements provide data for model validation. In situ aggregate settling velocities measured by the settling chamber were around 55 m d**-1. Aggregate sizes recorded by the particle camera hardly exceeded 1 mm. The model is based on a continuous size spectrum of aggregates, characterised by the prognostic aggregate mass and aggregate number concentration. Phytoplankton and detritus make up the aggregation pool, which has an averaged, prognostic and size dependent sinking. Model experiments were performed with dense and porous approximations of aggregates with varying maximum aggregate size and stickiness as well as with the inclusion of a disaggregation term. Similar surface productivity in all experiments has been generated in order to find the best combination of parameters that produce measured deep water fluxes. Although the experiments failed to represent surface particle number spectra, in the deep water some of them gave very similar slope and spectrum range as the particle camera observations. Particle fluxes at the mesotrophic sediment trap site off Cape Blanc (CB) have been successfully reproduced by the porous experiment with disaggregation term when particle remineralisation rate was 0.2 d**-1. The aggregation-disaggregation model improves the prediction capability of the original biogeochemical model significantly by giving much better estimates of fluxes for both upper and lower trap. The results also point to the need for more studies to enhance our knowledge on particle decay and its variation and to the role that stickiness play in the distribution of vertical fluxes.

Bulk rock magnetic measurements and thermomagnetic analyses of sediment core GeoB4901-8 from the southeastern flank of the Niger deep-sea fan

Diagenesis has extensively affected the magnetic mineral inventory of organic-rich late Quaternary sediments in the Niger deep-sea fan. Changes in concentration, grain size, and coercivity document modifications of the primary magnetic mineral assemblages at two horizons. The first front, the modern iron redox boundary, is characterized by a drastic decline in magnetic mineral content, coarsening of the grain size spectrum, and reduction in coercivity. Beneath a second front, the transition from the suboxic to the sulfidic anoxic domain, a further but less pronounced decrease in concentration and bulk grain size occurs. Finer grains and higher coercive magnetic constituents substantially increase in the anoxic environment. Low- and high-temperature experiments were performed on bulk sediments and on extracts which have also been examined by X-ray diffraction. Thermomagnetic analyses proved ferrimagnetic titanomagnetites of terrigenous provenance as the principal primary magnetic mineral components. Their broad range of titanium contents reflects the volcanogenic traits of the Niger River drainage areas. Diagenetic alteration is not only a grain size selective process but also critically depends on titanomagnetite composition. Low-titanium compounds are less resistant to diagenetic dissolution. Intermediate titanium content titanomagnetite thus persists as the predominant magnetic mineral fraction in the sulfidic anoxic sediments. At the Fe redox boundary, precipitation of authigenic, possibly bacterial, magnetite is documented. The presence of hydrogen sulfide in the pore water suggests a formation of secondary magnetic iron sulfides in the anoxic domain. Grain size-specific data argue for a gradual development of a superparamagnetic and single-domain iron sulfide phase in this milieu, most likely greigite.