Three species distribution models for the marine diatom Fragilariopsis kerguelensis of the Southern Ocean, links to NetCDF files

Fragilariopsis kerguelensis, a dominant diatom species throughout the Antarctic Circumpolar Current, is coined to be one of the main drivers of the biological silicate pump. Here, we study the distribution of this important species and expected consequences of climate change upon it, using correlative species distribution modeling and publicly available presence-only data. As experience with SDM is scarce for marine phytoplankton, this also serves as a pilot study for this organism group. We used the maximum entropy method to calculate distribution models for the diatom F. kerguelensis based on yearly and monthly environmental data (sea surface temperature, salinity, nitrate and silicate concentrations). Observation data were harvested from GBIF and the Global Diatom Database, and for further analyses also from the Hustedt Diatom Collection (BRM). The models were projected on current yearly and seasonal environmental data to study current distribution and its seasonality. Furthermore, we projected the seasonal model on future environmental data obtained from climate models for the year 2100. Projected on current yearly averaged environmental data, all models showed similar distribution patterns for F. kerguelensis. The monthly model showed seasonality, for example, a shift of the southern distribution boundary toward the north in the winter. Projections on future scenarios resulted in a moderately to negligibly shrinking distribution area and a change in seasonality. We found a substantial bias in the publicly available observation datasets, which could be reduced by additional observation records we obtained from the Hustedt Diatom Collection. Present-day distribution patterns inferred from the models coincided well with background knowledge and previous reports about F. kerguelensis distribution, showing that maximum entropy-based distribution models are suitable to map distribution patterns for oceanic planktonic organisms. Our scenario projections indicate moderate effects of climate change upon the biogeography of F. kerguelensis.

Magnetic observation of ODP Site 152-919 in the Irminger Basin, North Atlantic

Component natural remanent magnetizations derived from u-channel and 1-qcm discrete samples from ODP Site 919 (Irminger Basin) indicate the existence of four intervals of negative inclinations in the upper Brunhes Chronozone. According to the age model based on planktic oxygen isotope data, these "excursional" intervals occur in sediments deposited during the following time intervals: 32-34 ka, 39-41 ka, 180-188 ka and 205-225 ka. These time intervals correspond to polarity excursions detected elsewhere, known as Mono Lake, Laschamp, Iceland Basin and Pringle Falls. The isotope-based age model is supported by the normalized remanence (paleointensity) record that can be correlated to other calibrated paleointensity records for the 0-500 ka interval, such as that from ODP Site 983. For the intervals associated with the Mono Lake and Laschamp excursions, virtual geomagnetic poles (VGPs) reach equatorial latitudes and mid-southerly latitudes, respectively. For intervals associated with the Iceland Basin and Pringle Falls excursions, repeated excursions of VGPs to high southerly latitudes indicate rapid directional swings rather than a single short-lived polarity reversal. The directional instability associated with polarity excursions is not often recorded, probably due to smoothing of the sedimentary record by the process of detrital remanence (DRM) acquisition.

Temperature observation of CORK installation during ODP Leg 168

Four boreholes, drilled a few tens of meters into igneous basement on the eastern flank of the Juan de Fuca Ridge during ODP Leg 168, were sealed and instrumented for long-term monitoring to observe the hydrologic state of young sediment-sealed oceanic crust. The thermal regime is dominated by the effects of rapid fluid circulation in uppermost igneous basement driven by very small non-hydrostatic pressure gradients. Upper basement temperatures are uniform laterally between pairs of holes over distances of hundreds of meters to kilometers. In the case of two holes drilled into a sediment-buried basement ridge and adjacent valley, basement temperatures differ by less than 2 K despite the 2.2 km lateral separation of the sites and the 2.5:1 contrast in sediment cover thickness. Under conductive conditions, upper basement temperatures would differ by roughly 50 K. By comparison with modeling results, the observed degree of isothermality suggests a fluid flux of at least 10**-6 m/s (30 m/yr), and an effective permeability in the range of 10**-10 to 10**-9 m**2 in the uppermost igneous crust. The pressure difference available to drive this rapid flux between the ridge and valley, estimated by comparing the observed pressures via the isothermal upper basement hydrostat that is inferred to connect the two sites, is small (~2 kPa) and also suggests high permeability. Relative to the hydrostats defined by the local conductive sediment geotherms, substantial super-hydrostatic pressure (+18 kPa) is present within the buried basement ridge, and sub-hydrostatic pressure is present in the adjacent valley (-26 kPa). Such pressure differentials are the direct consequence of the advection-dominated thermal regime and small pressure losses in high-permeability basement, and are available to drive fluid seepage through sediment sections vertically up above and horizontally away from buried ridges, and down above valleys. No constraints are provided by any of the observations on the depth in the crust to which thermally or chemically significant flow might extend, although just as in the overlying sediments, the pattern of deep flow may be affected by the near-isothermal and near-hydrostatic conditions present in the permeable uppermost crustal section.

(Table 1) Snow and ice thickness and layering at observation points of the North Pole 34 drifting station

Variability of total alkalinity in sea ice of the high-latitudinal Arctic from November 2005 to May 2006 is considered. For the bulk of one- and two-year sea ice, alkalinity dependence on salinity is described as TA = k x Sal, where k is salinity/alkalinity ratio in under-ice water. The given relationship is valid within a wide range of salinity from 0.1 psu in desalinated fraction of two-year ice to 36 psu in snow on the young ice surface. Geochemically significant deviations from the relationship noted were observed exclusively in snow and the upper layer of one-year ice. In the upper layer of one-year ice, deficiency of alkalinity is observed ( delta TA ~= -0.07 mEq/kg, or -15%). In snow on the surface of the one-year ice, alkalinity excess is formed under desalination ( delta TA is as high as 1.3 mEq/kg, or 380%). Deviations registered are caused by possibility of carbonate precipitation in form of CaCO3 x 6H2O under seawater freezing. It is shown that ice formation and the following melting might cause losses of atmospheric CO2 of up to 3 x 10**12 gC/year.

Global distributions of epipelagic macrozooplankton abundance and biomass - Gridded data product (NetCDF) - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

Macrozooplankton are an important link between higher and lower trophic levels in the oceans. They serve as the primary food for fish, reptiles, birds and mammals in some regions, and play a role in the export of carbon from the surface to the intermediate and deep ocean. Little, however, is known of their global distribution and biomass. Here we compiled a dataset of macrozooplankton abundance and biomass observations for the global ocean from a collection of four datasets. We harmonise the data to common units, calculate additional carbon biomass where possible, and bin the dataset in a global 1 x 1 degree grid. This dataset is part of a wider effort to provide a global picture of carbon biomass data for key plankton functional types, in particular to support the development of marine ecosystem models. Over 387 700 abundance data and 1330 carbon biomass data have been collected from pre-existing datasets. A further 34 938 abundance data were converted to carbon biomass data using species-specific length frequencies or using species-specific abundance to carbon biomass data. Depth-integrated values are used to calculate known epipelagic macrozooplankton biomass concentrations and global biomass. Global macrozooplankton biomass has a mean of 8.4 µg C l-1, median of 0.15 µg C l-1 and a standard deviation of 63.46 µg C l-1. The global annual average estimate of epipelagic macrozooplankton, based on the median value, is 0.02 Pg C. Biomass is highest in the tropics, decreasing in the sub-tropics and increasing slightly towards the poles. There are, however, limitations on the dataset; abundance observations have good coverage except in the South Pacific mid latitudes, but biomass observation coverage is only good at high latitudes. Biomass is restricted to data that is originally given in carbon or to data that can be converted from abundance to carbon. Carbon conversions from abundance are restricted in the most part by the lack of information on the size of the organism and/or the absence of taxonomic information. Distribution patterns of global macrozooplankton biomass and statistical information about biomass concentrations may be used to validate biogeochemical models and Plankton Functional Type models.

Global Ocean Colour for Carbon Cycle Research (full product set)

In 2005, the International Ocean Colour Coordinating Group (IOCCG) convened a working group to examine the state of the art in ocean colour data merging, which showed that the research techniques had matured sufficiently for creating long multi-sensor datasets (IOCCG, 2007). As a result, ESA initiated and funded the DUE GlobColour project (http://www.globcolour.info/) to develop a satellite based ocean colour data set to support global carbon-cycle research. It aims to satisfy the scientific requirement for a long (10+ year) time-series of consistently calibrated global ocean colour information with the best possible spatial coverage. This has been achieved by merging data from the three most capable sensors: SeaWiFS on GeoEye's Orbview-2 mission, MODIS on NASA's Aqua mission and MERIS on ESA's ENVISAT mission. In setting up the GlobColour project, three user organisations were invited to help. Their roles are to specify the detailed user requirements, act as a channel to the broader end user community and to provide feedback and assessment of the results. The International Ocean Carbon Coordination Project (IOCCP) based at UNESCO in Paris provides direct access to the carbon cycle modelling community's requirements and to the modellers themselves who will use the final products. The UK Met Office's National Centre for Ocean Forecasting (NCOF) in Exeter, UK, provides an understanding of the requirements of oceanography users, and the IOCCG bring their understanding of the global user needs and valuable advice on best practice within the ocean colour science community. The three year project kicked-off in November 2005 under the leadership of ACRI-ST (France). The first year was a feasibility demonstration phase that was successfully concluded at a user consultation workshop organised by the Laboratoire d'Océanographie de Villefranche, France, in December 2006. Error statistics and inter-sensor biases were quantified by comparison with insitu measurements from moored optical buoys and ship based campaigns, and used as an input to the merging. The second year was dedicated to the production of the time series. In total, more than 25 Tb of input (level 2) data have been ingested and 14 Tb of intermediate and output products created, with 4 Tb of data distributed to the user community. Quality control (QC) is provided through the Diagnostic Data Sets (DDS), which are extracted sub-areas covering locations of in-situ data collection or interesting oceanographic phenomena. This Full Product Set (FPS) covers global daily merged ocean colour products in the time period 1997-2006 and is also freely available for use by the worldwide science community at http://www.globcolour.info/data_access_full_prod_set.html. The GlobColour service distributes global daily, 8-day and monthly data sets at 4.6 km resolution for, chlorophyll-a concentration, normalised water-leaving radiances (412, 443, 490, 510, 531, 555 and 620 nm, 670, 681 and 709 nm), diffuse attenuation coefficient, coloured dissolved and detrital organic materials, total suspended matter or particulate backscattering coefficient, turbidity index, cloud fraction and quality indicators. Error statistics from the initial sensor characterisation are used as an input to the merging methods and propagate through the merging process to provide error estimates for the output merged products. These error estimates are a key component of GlobColour as they are invaluable to the users; particularly the modellers who need them in order to assimilate the ocean colour data into ocean simulations. An intensive phase of validation has been undertaken to assess the quality of the data set. In addition, inter-comparisons between the different merged datasets will help in further refining the techniques used. Both the final products and the quality assessment were presented at a second user consultation in Oslo on 20-22 November 2007 organised by the Norwegian Institute for Water Research (NIVA); presentations are available on the GlobColour WWW site. On request of the ESA Technical Officer for the GlobColour project, the FPS data set was mirrored in the PANGAEA data library.