Euphausiid respiration model revamped, link to model results

Euphausiids constitute major biomass component in shelf ecosystems and play a fundamental role in the rapid vertical transport of carbon from the ocean surface to the deeper layers during their daily vertical migration (DVM). DVM depth and migration patterns depend on oceanographic conditions with respect to temperature, light and oxygen availability at depth, factors that are highly dependent on season in most marine regions. Changes in the abiotic conditions also shape Euphausiid metabolism including aerobic and anaerobic energy production. Here we introduce a global krill respiration model which includes the effect of latitude (LAT), the day of the year of interest (DoY), and the number of daylight hours on the day of interest (DLh), in addition to the basal variables that determine ectothermal oxygen consumption (temperature, body mass and depth) in the ANN model (Artificial Neural Networks). The newly implemented parameters link space and time in terms of season and photoperiod to krill respiration. The ANN model showed a better fit (r**2=0.780) when DLh and LAT were included, indicating a decrease in respiration with increasing LAT and decreasing DLh. We therefore propose DLh as a potential variable to consider when building physiological models for both hemispheres. We also tested for seasonality the standard respiration rate of the most common species that were investigated until now in a large range of DLh and DoY with Multiple Linear Regression (MLR) or General Additive model (GAM). GAM successfully integrated DLh (r**2= 0.563) and DoY (r**2= 0.572) effects on respiration rates of the Antarctic krill, Euphausia superba, yielding the minimum metabolic activity in mid-June and the maximum at the end of December. Neither the MLR nor the GAM approach worked for the North Pacific krill Euphausia pacifica, and MLR for the North Atlantic krill Meganyctiphanes norvegica remained inconclusive because of insufficient seasonal data coverage. We strongly encourage comparative respiration measurements of worldwide Euphausiid key species at different seasons to improve accuracy in ecosystem modelling.

Global habitat suitability of deep-water Antipatharia

Antipatharia are a diverse group of corals with many species found in deep water. Many Antipatharia are habitat for associates, have extreme longevity and some species can occur beyond 8500 m depth. As they are major constituents of 'coral gardens', which are Vulnerable Marine Ecosystems (VMEs), knowledge of their distribution and environmental requirements is an important pre-requisite for informed conservation planning particularly where the expense and difficulty of deep-sea sampling prohibits comprehensive surveys. This study uses a global database of Antipatharia distribution data to perform habitat suitability modelling using the Maxent methodology to estimate the global extent of black coral habitat suitability. The model of habitat suitability is driven by temperature but there is notable influence from other variables of topography, surface productivity and oxygen levels. This model can be used to predict areas of suitable habitat, which can be useful for conservation planning. The global distribution of Antipatharia habitat suitability shows a marked contrast with the distribution of specimen observations, indicating that many potentially suitable areas have not been sampled, and that sampling effort has been disproportionate to shallow, accessible areas inside marine protected areas (MPAs). Although 25% of Antipatharia observations are located in MPAs, only 7-8% of predicted suitable habitat is protected, which is short of the Convention on Biological Diversity target to protect 10% of ocean habitats by 2020.

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

This study is a first effort to compile the largest possible body of data available from different plankton databases as well as from individual published or unpublished datasets regarding diatom distribution in the world ocean. The data obtained originate from time series studies as well as spatial studies. This effort is supported by the Marine Ecosystem Data (MAREDAT) project, which aims at building consistent data sets for the main PFTs (Plankton Functional Types) in order to help validate biogeochemical ocean models by using converted C biomass from abundance data. Diatom abundance data were obtained from various research programs with the associated geolocation and date of collection, as well as with a taxonomic information ranging from group down to species. Minimum, maximum and average cell size information were mined from the literature for each taxonomic entry, and all abundance data were subsequently converted to biovolume and C biomass using the same methodology.

Modelled global ocean temperature from 6 ka to present

We compare the ocean temperature evolution of the Holocene as simulated by climate models and reconstructed from marine temperature proxies. This site provides informations about the Holocene temperature trends as simulated by the models. We use transient simulations from a coupled atmosphere-ocean general circulation model, as well as an ensemble of time slice simulations from the Paleoclimate Modelling Intercomparison Project. The general pattern of sea surface temperature (SST) in the models shows a high latitude cooling and a low latitude warming. The proxy dataset comprises a global compilation of marine alkenone- and Mg/Ca-derived SST estimates. Independently of the choice of the climate model, we observe significant mismatches between modelled and estimated SST amplitudes in the trends for the last 6000 years. Alkenone-based SST records show a similar pattern as the simulated annual mean SSTs, but the simulated SST trends underestimate the alkenone-based SST trends by a factor of two to five. For Mg/Ca, no significant relationship between model simulations and proxy reconstructions can be detected. We tested if such discrepancies can be caused by too simplistic interpretations of the proxy data. We tested different seasons and depths in the model to compare the proxy data trends, and can reconcile only part of the mismatches on a regional scale. We therefore considered the additional environmental factor changes in the planktonic organisms' habitat depth and a time-shift in the recording season to diagnose whether invoking those environmental factors can help reconciling the proxy records and the model simulations. We find that invoking shifts in the living season and habitat depth can remove some of the model-data discrepancies in SST trends. Regardless whether such adjustments in the environmental parameters during the Holocene are realistic, they indicate that when modeled temperature trends are set up to allow drastic shifts in the ecological behavior of planktonic organisms, they do not capture the full range of reconstructed SST trends. Our findings indicate that climate model and reconstructed temperature trends are to a large degree only qualitatively comparable, thus providing a challenge for the interpretation of proxy data as well as the models' sensitivity to orbital forcing.

Part of the global DOC versus AOU (dissolved organic carbon/apparent oxygen utilization) data compilation, Pacific Ocean

Total organic carbon (TOC) was analyzed on four transects along 140°W in 1992 using a high temperature combustion/discrete injection (HTC/DI) analyzer. For two of the transects, the analyses were conducted on-board ship. Mixed-layer concentrations of organic carbon varied from about 80 µM C at either end of the transect (12°N and 12°S) to about 60 µM C at the equator. Total organic carbon concentrations decreased rapidly below the mixed-layer to about 38-40 µM C at 1000 m across the transect. Little variation was observed below this depth; deep water concentrations below 2000 m were virtually monotonic at about 36 µM C. Repeat measurements made on subsequent cruises consistently found the same concentrations at 1000 m or deeper, but substantial variations were observed in the mixed-layer and the upper water column above 400 m depth. Linear mixing models of total organic carbon versus sigmaT exhibited zones of organic carbon formation and consumption. TOC was found to be inversely correlated with apparent oxygen utilization (AOU) in the region between the mixed-layer and the oxygen minimum. In the mixed-layer, TOC concentrations varied seasonally. Part of the variations in TOC at the equator was driven by changes in the upwelling rate in response to variations in physical forcing related to an El Niño and to the passage of tropical instability waves. TOC export fluxes, calculated from simple box models, averaged 8±4 mmol C/m**2/day at the equator and also varied seasonally. These export fluxes account for 50-75% of the total carbon deficit and are consistent with other estimates and model predictions.

Stable carbon and oxygen isotope ratios and age model of ODP Hole 182-1127B

The aim of this study was to evaluate the potential of constructing an oxygen and carbon isotope stratigraphy for the late Pleistocene succession from Hole 1127B drilled on the Great Australian Bight. Stable isotope analyses were performed on bulk- and fine-fraction (<38 µm) sediment samples. The oxygen isotope variations are generally smaller in magnitude than expected from global pelagic records. This is most likely due to the neriticly dominated sediment composition. Correlation of the oxygen isotope data with carbonate mineralogy and downhole logging data shows simultaneous variations and trends, which are particularly evident in the mid-Pleistocene sediments. Correlation of the oxygen isotope data with the classic SPECMAP curve is used to evaluate the stratigraphic potential of the Site 1127 sediments. This study indicates that an isotope stratigraphy based on planktonic and benthic foraminifers is needed to fully evaluate the response of cool-water carbonates deposited in a margin setting to global ice-volume fluctuations and, hence, the associated sea level variations.

Static GPS measurements and elevation model in the hinterland of Neumayer

As part of the CryoSat Cal/Val activities and the pre-site survey for an ice core drilling contributing to the International Partnerships in Ice Core Sciences (IPICS), ground based kinematic GPS measurements were conducted in early 2007 in the vicinity of the German overwintering station Neumayer (8.25° W and 70.65° S). The investigated area comprises the regions of the ice tongues Halvfarryggen and Søråsen, which rise from the Ekströmisen to a maximum of about 760 m surface elevation, and have an areal extent of about 100 km x 50 km each. Available digital elevation models (DEMs) from radar altimetry and the Antarctic Digital Database show elevation differences of up to hundreds of meters in this region, which necessitated an accurate survey of the conditions on-site. An improved DEM of the Ekströmisen surroundings is derived by a combination of highly accurate ground based GPS measurements, satellite derived laser altimetry data (ICESat), airborne radar altimetry (ARA), and radio echo sounding (RES). The DEM presented here achieves a vertical accuracy of about 1.3 m and can be used for improved ice dynamic modeling and mass balance studies.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Static GPS measurements and elevation model in the vicinity of the EDML deep-drilling site

Interpretation of ice-core records requires accurate knowledge of the past and present surface topography and stress-strain fields. The European Project for Ice Coring in Antarctica (EPICA) drilling site (0.0684° E and 75.0025° S, 2891.7 m) in Dronning Maud Land, Antarctica, is located in the immediate vicinity of a transient and splitting ice divide. A digital elevation model is determined from the combination of kinematic GPS measurements with the GLAS12 data sets from the ICESat satellite. Based on a network of stakes, surveyed with static GPS, the velocity field around the EDML drilling site is calculated. The annual mean velocity magnitude of 12 survey points amounts to 0.74 m/a. Flow directions mainly vary according to their distance from the ice divide. Surface strain rates are determined from a pentagon-shaped stake network with one center point, close to the drilling site. The strain field is characterised by along flow compression, lateral dilatation, and vertical layer thinning.

Age model for the time between 3.45 and 2.45 Ma of ODP Site 189-1172

Calcareous nannoplankton assemblages and benthic d18O isotopes of Pliocene deep-sea sediments of ODP site 1172 (East of Tasmania) have been studied to improve our knowledge of the Southern Ocean paleoceanography. Our study site is located just north of the Subtropical Front (STF), an ideal setting to monitor migrations of the STF during our study period, between 3.45 and 2.45 Ma. The assemblage identified at ODP site 1172 has been interpreted as characteristic for the transitional zone water mass, located south of the STF, based on: (i) the low abundances (< 1%) of subtropical taxa, (ii) relatively high percentages of Coccolithus pelagicus, a subpolar type species, (iii) abundances from 2-10% of Calcidiscus leptoporus, a species that frequently inhabits the zone south of the STF and (iv) the high abundances of small Noelaerhabdaceae which at present dominates the zone south of the STF. Across our interval the calcareous nannoplankton manifests glacial-interglacial variability. We have identified cold events, characterized by high abundances of C. pelagicus which coincide with glacial periods, except during G7. After 3.1 Ma cold events are more frequent, in concordance with global cooling trends. Around 2.75 Ma, the interglacial stage G7 is characterized by anomalous low temperatures which most likely are linked to definite closure of the Central American Seaway (CAS), an event that is believed to have had global consequences. A gradual increase of very small Reticulofenestra across our section marks a significant trend in the small Noelaerhabdaceae species group and has been linked to a general enhanced mixing of the water column in agreement with previous studies. It is suggested that a rapid decline of small Gephyrocapsa after isotopic stage G7 might be related to the cooling observed in our study site after the closure of the CAS.