Latitudinal phytoplankton distribution and the neutral theory of biodiversity

Aim Recent studies have suggested that global diatom distributions are not limited by dispersal, in the case of both extant species and fossil species, but rather that environmental filtering explains their spatial patterns. Hubbell's neutral theory of biodiversity provides a framework in which to test these alternatives. Our aim is to test whether the structure of marine phytoplankton (diatoms, dinoflagellates and coccolithophores) assemblages across the Atlantic agrees with neutral theory predictions. We asked: (1) whether intersite variance in phytoplankton diversity is explained predominantly by dispersal limitation or by environmental conditions; and (2) whether species abundance distributions are consistent with those expected by the neutral model. Location Meridional transect of the Atlantic (50 degrees N50 degrees S). Methods We estimated the relative contributions of environmental factors and geographic distance to phytoplankton composition using similarity matrices, Mantel tests and variation partitioning of the species composition based upon canonical ordination methods. We compared the species abundance distribution of phytoplankton with the neutral model using Etienne's maximum-likelihood inference method. Results Phytoplankton communities are slightly more determined by niche segregation (24%), than by dispersal limitation and ecological drift (17%). In 60% of communities, the assumption of neutrality in species' abundance distributions could not be rejected. In tropical zones, where oceanic gyres enclose large stable water masses, most communities showed low species immigration rates; in contrast, we infer that communities in temperate areas, out of oligotrophic gyres, have higher rates of species immigration. Conclusions Phytoplankton community structure is consistent with partial niche assembly and partial dispersal and drift assembly (neutral processes). The role of dispersal limitation is almost as important as habitat filtering, a fact that has been largely overlooked in previous studies. Furthermore, the polewards increase in immigration rates of species that we have discovered is probably caused by water mixing conditions and productivity.

Long-term responses of North Atlantic calcifying plankton to climate change

The global increase in atmospheric carbon dioxide concentration is potentially threatening marine biodiversity in two ways. First, carbon dioxide and other greenhouse gases accumulating in the atmosphere are causing global warming1. Second, carbon dioxide is altering sea water chemistry, making the ocean more acidic2. Although temperature has a cardinal influence on all biological processes from the molecular to the ecosystem level3, acidification might impair the process of calcification or exacerbate dissolution of calcifying organisms4. Here, we show however that North Atlantic calcifying plankton primarily responded to climate-induced changes in temperatures during the period 1960–2009, overriding the signal from the effects of ocean acidification. We provide evidence that foraminifers, coccolithophores, both pteropod and nonpteropod molluscs and echinoderms exhibited an abrupt shift circa 1996 at a time of a substantial increase in temperature5 and that some taxa exhibited a poleward movement in agreement with expected biogeographical changes under sea temperature warming6,7. Although acidification may become a serious threat to marine calcifying organisms, our results suggest that over the study period the primary driver of North Atlantic calcifying planktonwas oceanic temperature.

Identifying four phytoplankton functional types from space: An ecological approach

Deriving maps of phytoplankton taxa based on remote sensing data using bio-optical properties of phytoplankton alone is challenging. A more holistic approach was developed using artificial neural networks, incorporating ecological and geographical knowledge together with ocean color, bio-optical characteristics, and remotely sensed physical parameters. Results show that the combined remote sensing approach could discriminate four major phytoplankton functional types (diatoms, dinoflagellates, coccolithophores, and silicoflagellates) with an accuracy of more than 70%. Models indicate that the most important information for phytoplankton functional type discrimination is spatio-temporal information and sea surface temperature. This approach can supply data for large-scale maps of predicted phytoplankton functional types, and an example is shown.

Development and validation of a molecular technique for the analysis of archived formalin-preserved phytoplankton samples permits retrospective assessment of Emiliania huxleyi communities

The coccolithophores, particularly the species Emiliania huxleyi (Lohmann) Hay & Mohler, account for the bulk of global calcium carbonate production and as such play a fundamental role in global CO2 cycling and the carbonate chemistry of the oceans. To evaluate the response of this functional group to the effects of climate change, we undertook a feasibility study to determine whether a retrospective approach could be used on archived coccolithophore datasets. We demonstrate for the first time a technique for the extraction of E. huxleyi nucleic acids from archived formalin-fixed samples of the long-term Continuous Plankton Recorder. Molecular analysis of a nine year old formalin-fixed sample reveals the presence of a diverse population of E. huxleyi genotypes within a developing coccolithophore bloom. In addition, E. huxleyi sequences were amplified from a number of formalin-fixed samples, the earliest of which was collected in August 1972. This molecular assay promises the possibility of studying global variations in the distribution and genetic make-up of E. huxleyi communities over extensive periods of time. (c) 2008 Elsevier B.V. All rights reserved.

Coccolithophore bloom detection in the north east Atlantic using SeaWiFS: Algorithm description, application and sensitivity analysis

Coccolithophores are the largest source of calcium carbonate in the oceans and are considered to play an important role in oceanic carbon cycles. Current methods to detect the presence of coccolithophore blooms from Earth observation data often produce high numbers of false positives in shelf seas and coastal zones due to the spectral similarity between coccolithophores and other suspended particulates. Current methods are therefore unable to characterise the bloom events in shelf seas and coastal zones, despite the importance of these phytoplankton in the global carbon cycle. A novel approach to detect the presence of coccolithophore blooms from Earth observation data is presented. The method builds upon previous optical work and uses a statistical framework to combine spectral, spatial and temporal information to produce maps of coccolithophore bloom extent. Validation and verification results for an area of the north east Atlantic are presented using an in situ database (N = 432) and all available SeaWiFS data for 2003 and 2004. Verification results show that the approach produces a temporal seasonal signal consistent with biological studies of these phytoplankton. Validation using the in situ coccolithophore cell count database shows a high correct recognition rate of 80% and a low false-positive rate of 0.14 (in comparison to 63% and 0.34 respectively for the established, purely spectral approach). To guide its broader use, a full sensitivity analysis for the algorithm parameters is presented.

Coccolithophore surface distributions in the North Atlantic and their modulation of the air-sea flux of CO2 from 10 years of satellite Earth observation data

Coccolithophores are the primary oceanic phytoplankton responsible for the production of calcium carbonate (CaCO3). These climatically important plankton play a key role in the oceanic carbon cycle as a major contributor of carbon to the open ocean 5 carbonate pump (�50%) and their formation can affect the atmosphere-to-ocean (airsea) uptake of carbon dioxide (CO2) through increasing the seawater partial pressure of CO2 (pCO2). Here we document variations in the areal extent of surface blooms of the globally important coccolithophore, Emiliania huxleyi, in the North Atlantic over a 10-year period (1998–2007), using Earth observation data from the Sea-viewing Wide 10 Field of view Sensor (SeaWiFS).We calculate the annual mean surface areal coverage of E. huxleyi in the North Atlantic to be 474 000±119 000km2 yr−1, which results in a net CaCO3 production of 0.62±0.15 Tg CaCO3 carbon per year. However, this surface coverage and net production can fluctuate by −54/+81% about these mean values and are strongly correlated with the El Ni˜no/Southern Oscillation (ENSO) climate os15 cillation index (r =0.75, p<0.02). Our analysis evaluates the spatial extent over which the E. huxleyi blooms in the North Atlantic can increase the pCO2 and thus decrease the localised sink of atmospheric CO2. In regions where the blooms are prevalent, the average reduction in the monthly CO2 sink can reach 12 %. The maximum reduction of the monthly CO2 sink in the time series is 32 %. This work suggests that the high 20 variability, frequency and distribution of these calcifying plankton and their impact on pCO2 should be considered within modelling studies of the North Atlantic if we are to fully understand the variability of its air-to-sea CO2 flux.

Exploiting satellite earth observation to quantify current global oceanic DMS flux and its future climate sensitivity

We used coincident Envisat RA2 and AATSR temperature and wind speed data from 2008/2009 to calculate the global net sea-air flux of dimethyl sulfide (DMS), which we estimate to be 19.6 Tg S a21. Our monthly flux calculations are compared to open ocean eddy correlation measurements of DMS flux from 10 recent cruises, with a root mean square difference of 3.1 lmol m22 day21. In a sensitivity analysis, we varied temperature, salinity, surface wind speed, and aqueous DMS concentration, using fixed global changes as well as CMIP5 model output. The range of DMS flux in future climate scenarios is discussed. The CMIP5 model predicts a reduction in surface wind speed and we estimate that this will decrease the global annual sea-air flux of DMS by 22% over 25 years. Concurrent changes in temperature, salinity, and DMS concentration increase the global flux by much smaller amounts. The net effect of all CMIP5 modelled 25 year predictions was a 19% reduction in global DMS flux. 25 year DMS concentration changes had significant regional effects, some positive (Southern Ocean, North Atlantic, Northwest Pacific) and some negative (isolated regions along the Equator and in the Indian Ocean). Using satellite-detected coverage of coccolithophore blooms, our estimate of their contribution to North Atlantic DMS emissions suggests that the coccolithophores contribute only a small percentage of the North Atlantic annual flux estimate, but may be more important in the summertime and in the northeast Atlantic.

Nanoplankton and picoplankton in the Western English Channel: abundance and seasonality from 2007–2013

The nano- and picoplankton community at Station L4 in the Western English Channel was studied between 2007 and 2013 by flow cytometry to quantify abundance and investigate seasonal cycles within these communities. Nanoplankton included both photosynthetic and heterotrophic eukaryotic single-celled organisms while the picoplankton included picoeukaryote phytoplankton, Synechococcus sp. cyanobacteria and heterotrophic bacteria. A Box–Jenkins Transfer Function climatology analysis of surface data revealed that Synechococcus sp., cryptophytes, and heterotrophic flagellates had bimodal annual cycles. Nanoeukaryotes and both high and low nucleic acid-containing bacteria (HNA and LNA, respectively) groups exhibited unimodal annual cycles. Phaeocystis sp., whilst having clearly defined abundance maxima in spring was not detectable the rest of the year. Coccolithophores exhibited a weak seasonal cycle, with abundance peaks in spring and autumn. Picoeukaryotes did not exhibit a discernable seasonal cycle at the surface. Timings of maximum group abundance varied through the year. Phaeocystis sp. and heterotrophic flagellates peaked in April/May. Nanoeukaryotes and HNA bacteria peaked in June/July and had relatively high abundance throughout the summer. Synechococcus sp., cryptophytes and LNA bacteria all peaked from mid to late September. The transfer function model techniques used represent a useful means of identifying repeating annual cycles in time series data with the added ability to detect trends and harmonic terms at different time scales from months to decades.

Feeding selectivity of bivalve larvae on natural plankton assemblages in the Western English Channel

Meroplankton, including bivalve larvae, are an important and yet understudied component of coastal marine food webs. Understanding the baseline of meroplankton ecology is imperative to establish and predict their sensitivity to local and global marine stressors. Over an annual cycle (October 2009–September 2010), bivalve larvae were collected from the Western Channel Observatory time series station L4 (50°15.00′N, 4°13.02′W). The morphologically similar larvae were identified by analysis of the 18S nuclear small subunit ribosomal RNA gene, and a series of incubation experiments were conducted to determine larval ingestion rates on natural plankton assemblages. Complementary gut content analysis was performed using a PCR-based method for detecting prey DNA both from field-collected larvae and those from the feeding experiments. Molecular identification of bivalve larvae showed the community composition to change over the course of the sampling period with domination by Phaxas in winter and higher diversity in autumn. The larvae selected for nanoeukaryotes (2–20 µm) including coccolithophores (<20 µm) which together comprised >75 % of the bivalve larvae diet. Additionally, a small percentage of carbon ingested originated from heterotrophic ciliates (<30 µm). The molecular analysis of bivalve larvae gut content provided increased resolution of identification of prey consumed and demonstrated that the composition of prey consumed established through bottle incubations conferred with that established from in situ larvae. Despite changes in bivalve larvae community structure, clearance rates of each prey type did not change significantly over the course of the experiment, suggesting different bivalve larvae species may consume similar prey.

Metais em águas e sedimentos da plataforma dos Açores

A presente tese tem por objetivo principal contribuir para o conhecimento da geoquímica sedimentar da zona oceânica da crista da Terceira e montanhas submarinas a sul (região entre 29-39ºN e 27-32ºW), integrando também a caraterização dos metais e nutrientes na coluna de água e propondo concentrações para servirem de referência nesta região do Atlântico Central. Para o efeito foram realizadas amostragens na coluna de água em sete locais e de sedimento em cinco locais, durante a campanha oceanográfica designada por EMEPC/AÇORES/G3/2007 a bordo do navio SV Kommandor Jack, no âmbito do projeto da Estrutura de Missão para a Extensão da Plataforma Continental (EMEPC). Os perfis de CTD da coluna de água na região estudada revelam a presença de massas de água distintas: a Western North Atlantic Central Water (WNACW), a Eastern North Atlantic Central Water tropical (ENACWt), a Eastern North Atlantic Central Water polar (ENACWp), a Mediterranean Overflow Water (MOW), a Deep Mediterranean Water (DMW) e a North Eastern Atlantic Deep Water (NEADW). Observou-se nos perfis de temperatura e salinidade, referentes aos primeiros 200 m da coluna de água, um gradiente meridional negativo entre as estações localizadas na crista da Terceira e as estações localizadas mais a sul. Observou-se nas águas superficiais valores de oxigénio dissolvido de 93% e de pH de 8,1, assim como que as concentrações dos nutrientes NOx, PO4 e SiO2 variam de acordo com a atividade biológica, tendo-se registado concentrações medianas mais baixas, respetivamente de 6,5, 0,23 e 1,3 mol L-1, que aumentam com a profundidade devido à ausência de produção primária (respetivamente 31, 1,4 e 22 mol L-1). As concentrações de NH4 e de SO4 não variam significativamente nas massas de água, sendo os valores medianos mínimos e máximos de 0,69 a 0,79 mol L-1 para o NH4 e de 30 a 32 mol L-1 para o SO4. São propostas concentrações de referência para as massas de água, para os elementos cobre, cádmio, chumbo e arsénio. Os perfis de sedimento analisados permitem distinguir os sedimentos na crista da Terceira (core A) dos restantes (cores B a E). A grande variabilidade textural encontrada no core A, que contrasta com os outros cores analisados, deve-se a importantes contribuições terrígenas, originadas pela erosão sub-aérea e pela atividade vulcânica das ilhas próximas. iv resumo (continuação) A análise mineralógica, efetuada à fração areia e à fração fina (< 63 μm), confirma que os sedimentos do core A derivam de rochas vulcânicas formadas maioritariamente por piroxenas, olivinas, anfíbolas, biotite, alterites e ainda calcite, plagióclase e magnetite, tendo-se identificado ao microscópio a glauconite e o vidro vulcânico. De acordo com a composição química destes minerais o core A apresenta valores mais elevados de Al, Fe, K, P, Mg, Si, Na, Zn, V, Cr e Mn relativamente aos cores B a E. Os cores B a E apresentam grandes quantidades de calcite (>80%) formada maioritariamente por foraminíferos e nanoplâncton calcário (cocolitóforos). A fração areia confirma a composição maioritariamente carbonatada com grande abundância de material biogénico formado por oozes de foraminíferos (planctónicos e bentónicos) com raras espículas de espongiários e restos de conchas. Os cores B a E apresentam valores muito mais elevados que o core A para os elementos Ca e Sr. Os resultados para o Al, Fe, K, P, Si, Na, As, Cu, Ni, Zn, V, Cr, Li, Pb, Cd e Co presentes nos locais B, C, D e E sugerem que estes cores são comparáveis aos sedimentos de fundo carbonatados. Propõe-se concentrações de referência para a região do Atlântico compreendida entre 29-39ºN e 27-32ºW considerando a primeira camada colhida em cada core. Para o core A as concentrações são normalizadas a 5% de Al e CaCO3, enquanto que para os cores B a E são normalizadas a 2% de Al e CaCO3. Assim as concentrações de referência para o core A são: As – 18 mg kg-1, Cr – 91 mg kg-1, Cu – 127 mg kg-1, Ni – 84 mg kg-1, Pb – 41 mg kg-1, Hg – 41 ng g-1 e Zn – 482 mg kg-1. Para os cores B a E as concentrações de referência são: As – 3 mg kg-1, Cr – 10 mg kg-1, Cu – 36 mg kg-1 Ni – 12 mg kg -1, Hg – 3 ng g-1 e Zn – 20 mg kg-1. Para os restantes metais as concentrações de referência para o core A são: Al – 9%, Si – 25%, Fe – 6%, Ca – 13%, K – 2%, Mg – 2%, Na – 3%, P – 0,4%, Sr – 900 mg kg-1, Li – 10 mg kg-1, Mn – 1200 mg kg-1, Ba – 700 mg kg-1 e V – 140 mg kg-1. Para os cores B a E as concentrações de referência são: Al – 0,9%, Si – 2%, Fe – 0,2%, Ca – 95%, K – 0,3%, Mg – 0,4%, Na – 0,3%, P – 0,04%, Sr – 2600 mg kg-1, Li – 5 mg kg-1, Mn – 240 mg kg-1, Ba – 345 mg kg-1, Co – 2 mg kg-1 e V – 6 mg kg-1. Os resultados da presente tese constituem um contributo para a caraterização geoquímica da região e podem servir de referência à monitorização futura do mar dos Açores e montes submarinos a sul.

(Table 1) Calcium isotope composition, Mg/Ca, Sr/Ca and U/Ca ratios of calcareous cysts of Thoracosphaera heimii

We investigated the influences of temperature, salinity and pH on the calcium isotope as well as trace and minor element (uranium, strontium, magnesium) to Ca ratios on calcium carbonate cysts of the calcareous dinoflagellate species Thoracosphaera heimii grown in laboratory cultures. The natural habitat of this species is the photic zone (preferentially at the chlorophyll maximum depth) of temperate to tropical oceans, and it is abundant in deep-sea sediments over the entire Cenozoic. In our experiments, temperatures ranged from 12 to 30 °C, salinity from 36.5 to 38.8 and pH from 7.9 to 8.4. The delta44/40Ca of T. heimii cysts resembles that of other marine calcifiers, including coccolithophores, foraminifers and corals. However, its temperature sensitivity is considerably smaller and statistically insignificant, and T. heimii might serve as a recorder of changes in seawater delta44/40Ca over geologic time. The Sr/Ca ratios of T. heimii cysts show a pronounced temperature sensitivity (0.016 mmol/mol °C**-1) and have the potential to serve as a palaeo-sea surface temperature proxy. No clear temperature- and pH-dependences were observed for Mg/Ca. U/Ca seems to be influenced by temperature and pH, but the correlations change sign at 23 °C and pH 8.2, respectively.

Quantitative coccolithophore analyses in sediment core MD01-2446

Quantitative coccolithophore analyses were performed in core MD01-2446, located in the mid-latitude North Atlantic, to reconstruct climatically induced sea-surface water conditions throughout Marine Isotope Stages (MIS) 14-9. The data are compared to new and available paleoenvironmental proxies from the same site as well as other nearby North Atlantic records that support the coccolithophore signature at glacial-interglacial to millennial climate scale. Total coccolithophore absolute abundance increases during interglacials but abruptly drops during the colder glacial phases and deglaciations. Coccolithophore warm-water taxa (wwt) indicate that MIS11c and MIS9e experienced warmer and more stable conditions throughout the whole photic zone compared to MIS13. MIS11 was a long-lasting warmer and stable interglacial characterized by a climate optimum during MIS11c when a more prominent influence of the subtropical front at the site is inferred. The wwt pattern also suggests distinct interstadial and stadial events lasting about 4-10 kyr. The glacial increases of Gephyrocapsa margereli-G. muellerae 3-4 µm along with higher values of Corg, additionally supported by the total alkenone abundance at Site U1313, indicate more productive surface waters, likely reflecting the migration of the polar front into the mid-latitude North Atlantic. Distinctive peaks of G. margereli-muellerae (> 4 µm), C. pelagicus pelagicus, Neogloboquadrina pachyderma left coiling, and reworked nannofossils, combined with minima in total nannofossil accumulation rate, are tracers of Heinrich-type events during MIS12 and MIS10. Additional Heinrich-type events are suggested during MIS12 and MIS14 based on biotic proxies, and we discuss possible iceberg sources at these times. Our results improve the understanding of mid-Brunhes paleoclimate and the impact on phytoplankton diversity in the mid-latitude North Atlantic region.

Northwards and eastwards velocities extracted from a historical run (for the year 2000) of the UVic Earth System Climate Model of Weaver et al. (2001)

This study provides a theoretical assessment of the potential bias due to differential lateral transport on multi-proxy studies based on a range of marine microfossils. Microfossils preserved in marine sediments are at the centre of numerous proxies for paleoenvironmental reconstructions. The precision of proxies is based on the assumption that they accurately represent the overlying watercolumn properties and faunas. Here we assess the possibility of a syn-depositional bias in sediment assemblages caused by horizontal drift in the water column, due to differential settling velocities of sedimenting particles based on their shape, size and density, and due to differences in current velocities. Specifically we calculate the post-mortem lateral transport undergone by planktic foraminifera and a range of other biological proxy carriers (diatoms, radiolaria and fecal pellets transporting coccolithophores) in several regions with high current velocities. We find that lateral transport of different planktic foraminiferal species is minimal due to high settling velocities. No significant shape- or size-dependent sorting occurs before reaching the sediment, making planktic foraminiferal ideal proxy carriers. In contrast, diatoms, radiolaria and fecal pellets can be transported up to 500km in some areas. For example in the Agulhas current, transport can lead to differences of up to 2°C in temperature reconstructions between different proxies in response to settling velocities. Therefore, sediment samples are likely to contain different proportions of local and imported particles, decreasing the precision of proxies based on these groups and the accuracy of the temperature reconstruction.

Quaternary calcareous nannofossils from Periantarctic basins

Currently, coccolithophores have a widespread oceanic distribution and are reported from most latitudes, but not those higher than 65°S. Fifteen piston cores were sampled with the aim of investigating the distribution and abundance variation of Quaternary calcareous nannofossils of the Antarctic region, south of the Antarctic Divergence (>65°S), particularly from Maud Rise, Bausan Bank, and from Weddell, Ross and Bellingshausen Seas. A calcareous nannofossil cold-taxa association is present in most of the cores examined and their discontinuous occurrence is thought to indicate key environmental relationships. The presence of calcareous nannofossils is correlated with interglacial intervals with warmer SSTs and may indicate high productivity and an open-ocean environment. Our results confirm that, during short periods of the late Quaternary, coccolithophorids occurred at southern high latitudes, in the western Antarctic basins, while in the eastern Antarctic basins they are nearly absent, suggesting more variable SSTs near West Antarctic Ice Sheet.

Current velocities and settling experimental results

This study provides a theoretical assessment of the potential bias due to differential lateral transport on multi-proxy studies based on a range of marine microfossils. Microfossils preserved in marine sediments are at the centre of numerous proxies for paleoenvironmental reconstructions. The precision of proxies is based on the assumption that they accurately represent the overlying watercolumn properties and faunas. Here we assess the possibility of a syn-depositional bias in sediment assemblages caused by horizontal drift in the water column, due to differential settling velocities of sedimenting particles based on their shape, size and density, and due to differences in current velocities. Specifically we calculate the post-mortem lateral transport undergone by planktic foraminifera and a range of other biological proxy carriers (diatoms, radiolaria and fecal pellets transporting coccolithophores) in several regions with high current velocities. We find that lateral transport of different planktic foraminiferal species is minimal due to high settling velocities. No significant shape- or size-dependent sorting occurs before reaching the sediment, making planktic foraminiferal ideal proxy carriers. In contrast, diatoms, radiolaria and fecal pellets can be transported up to 500 km in some areas. For example in the Agulhas current, transport can lead to differences of up to 2°C in temperature reconstructions between different proxies in response to settling velocities. Therefore, sediment samples are likely to contain different proportions of local and imported particles, decreasing the precision of proxies based on these groups and the accuracy of the temperature reconstruction.