Dome C age scale EDC1 for ice core EDC99

A tentative age scale (EDC1) for the last 45 kyr is established for the new 788 m EPICA Dome C ice core using a simple ice flow model. The age of volcanic eruptions, the end of the Younger Dryas event, and the estimated depth and age of elevated 10Be, about 41 kyr ago were used to calibrate the model parameters. The uncertainty of EDC1 is estimated to ±10 yr for 0 to 700 yr BP, up to ±200 yr back to 10 kyr BP, and up to ±2 kyr back to 41 kyr BP. The age of the air in the bubbles is calculated with a firn densification model. In the Holocene the air is about 2000 yr younger than the ice and about 5500 yr during the last glacial maximum.

Framework for risk assessment for maritime transportation systems

Maritime accidents involving ships carrying passengers may pose a high risk with respect to human casualties. For effective risk mitigation, an insight into the process of risk escalation is needed. This requires a proactive approach when it comes to risk modelling for maritime transportation systems. Most of the existing models are based on historical data on maritime accidents, and thus they can be considered reactive instead of proactive. This paper introduces a systematic, transferable and proactive framework estimating the risk for maritime transportation systems, meeting the requirements stemming from the adopted formal definition of risk. The framework focuses on ship-ship collisions in the open sea, with a RoRo/Passenger ship (RoPax) being considered as the struck ship. First, it covers an identification of the events that follow a collision between two ships in the open sea, and, second, it evaluates the probabilities of these events, concluding by determining the severity of a collision. The risk framework is developed with the use of Bayesian Belief Networks and utilizes a set of analytical methods for the estimation of the risk model parameters. The model can be run with the use of GeNIe software package. Finally, a case study is presented, in which the risk framework developed here is applied to a maritime transportation system operating in the Gulf of Finland (GoF). The results obtained are compared to the historical data and available models, in which a RoPax was involved in a collision, and good agreement with the available records is found.

(Table A1) Chemical composition of abyssal peridotites from ODP Hole 153-920

The compatibility of vanadium (V) during mantle melting is a function of oxygen fugacity (fO2): at high fO2's, V becomes more incompatible. The prospects and limitations of using the V content of peridotites as a proxy for paleo-fO2 at the time of melt extraction were investigated here by assessing the uncertainties in V measurements and the sensitivity of V as a function of degree of melt extracted and fO2. V-MgO and V-Al2O3 systematics were found to be sensitive to fO2 variations, but consideration of the uncertainties in measurements and model parameters indicates that V is sensitive only to relative fO2 differences greater than ~2 log units. Post-Archean oceanic mantle peridotites, as represented by abyssal peridotites and obducted massif peridotites, have V-MgO and -Al2O3 systematics that can be modeled by 1.5 GPa melting between FMQ - 3 and FMQ - 1. This is consistent with fO2's of the mantle source for mid-ocean ridge basalts (MORBs) as determined by the Fe3+ activity of peridotitic minerals and basaltic glasses. Some arc-related peridotites have slightly lower V for a given degree of melting than oceanic mantle peridotites, and can be modeled by 1.5 GPa melting at fO2's as high as FMQ. However, the majority of arc-related peridotites have V-MgO systematics overlapping that of oceanic mantle peridotites, suggesting that although some arc mantle may melt under slightly oxidizing conditions, most arc mantle does not. The fact that thermobarometrically determined fO2's in arc peridotites and lavas can be significantly higher than that inferred from V systematics, suggests that V retains a record of the fO2 during partial melting, whereas the activity of Fe3+ in arc peridotitic minerals and lavas reflect subsequent metasomatic overprints and magmatic differentiation/emplacement processes, respectively. Peridotites associated with middle to late Archean cratonic mantle are characterized by highly variable V-MgO systematics. Tanzanian cratonic peridotites have V systematics indistinguishable from post-Archean oceanic mantle and can be modeled by 3 GPa partial melting at ~FMQ - 3. In contrast, many South African and Siberian cratonic peridotites have much lower V contents for a given degree of melting, suggesting at first glance that partial melting occurred at high fO2's. More likely, however, their unusually low V contents for a given degree of melting may be artifacts of excess orthopyroxene, a feature that pervades many South African and Siberian peridotites but not the Tanzanian peridotites. This is indicated by the fact that the V contents of South African and Siberian peridotites are correlated with increases in SiO2 content, generating data arrays that cannot be modeled by partial melting but can instead be generated by the addition of orthopyroxene through processes unrelated to primary melt depletion. Correction for orthopyroxene addition suggests that the South African and Siberian peridotites have V-MgO systematics similar to those of Tanzanian peridotites. Thus, if the Tanzanian peridotites represent the original partial melting residues, and if the South African and Siberian peridotites have been modified by orthopyroxene addition, then there is no indication that Archean cratonic mantle formed under fO2's significantly greater than that of modern oceanic mantle. Instead, the fO2's inferred from the V systematics in these three cratonic peridotite suites are within range of modern oceanic mantle. This also suggests that the transition from a highly reducing mantle in equilibrium with a metallic core to the present oxidized state must have occurred by late Archean times.

Bathymetric map of Heron Reef, Australia, derived from airborne hyperspectral data at 1 m resolution

A simple method for efficient inversion of arbitrary radiative transfer models for image analysis is presented. The method operates by representing the shape of the function that maps model parameters to spectral reflectance by an adaptive look-up tree (ALUT) that evenly distributes the discretization error of tabulated reflectances in spectral space. A post-processing step organizes the data into a binary space partitioning tree that facilitates an efficient inversion search algorithm. In an example shallow water remote sensing application, the method performs faster than an implementation of previously published methodology and has the same accuracy in bathymetric retrievals. The method has no user configuration parameters requiring expert knowledge and minimizes the number of forward model runs required, making it highly suitable for routine operational implementation of image analysis methods. For the research community, straightforward and robust inversion allows research to focus on improving the radiative transfer models themselves without the added complication of devising an inversion strategy.

Strong shift from HCO3- to CO2 uptake in Emiliania huxleyi with acidification: new approach unravels acclimation versus short-term pH effects

Effects of ocean acidification on Emiliania huxleyi strain RCC 1216 (calcifying, diploid life-cycle stage) and RCC 1217 (non-calcifying, haploid life-cycle stage) were investigated by measuring growth, elemental composition, and production rates under different pCO2 levels (380 and 950 µatm). In these differently acclimated cells, the photosynthetic carbon source was assessed by a (14)C disequilibrium assay, conducted over a range of ecologically relevant pH values (7.9-8.7). In agreement with previous studies, we observed decreased calcification and stimulated biomass production in diploid cells under high pCO2, but no CO2-dependent changes in biomass production for haploid cells. In both life-cycle stages, the relative contributions of CO2 and HCO3 (-) uptake depended strongly on the assay pH. At pH values =< 8.1, cells preferentially used CO2 (>= 90 % CO2), whereas at pH values >= 8.3, cells progressively increased the fraction of HCO3 (-) uptake (~45 % CO2 at pH 8.7 in diploid cells; ~55 % CO2 at pH 8.5 in haploid cells). In contrast to the short-term effect of the assay pH, the pCO2 acclimation history had no significant effect on the carbon uptake behavior. A numerical sensitivity study confirmed that the pH-modification in the (14)C disequilibrium method yields reliable results, provided that model parameters (e.g., pH, temperature) are kept within typical measurement uncertainties. Our results demonstrate a high plasticity of E. huxleyi to rapidly adjust carbon acquisition to the external carbon supply and/or pH, and provide an explanation for the paradoxical observation of high CO2 sensitivity despite the apparently high HCO3 (-) usage seen in previous studies.

SeaWiFS (version R2010.0) Monthly Global 9 km Microplankton Fraction, links to NetCDF files

Phytoplankton cell size is important to biogeochemical and food web processes. The goal of this study is to estimate phytoplankton cell size distribution from satellite imagery of spectral remote sensing reflectance (Rrs(lambda)). Previous studies have indicated phytoplankton size classes have distinctive absorption spectra despite the physiological and taxonomic variability within an assemblage. For this study, the chlorophyll specific absorption spectra for phytoplankton size class extremes, pico- and microphytoplankton, are weighted by the percent microplankton (Sfm) and are the basis of phytoplankton size retrieval from SeaWiFS imagery. Satellite retrievals of Sfm are done through implementation of a forward optical model look-up table (LUT) that incorporates the range of absorption and scattering variability due to phytoplankton size, chlorophyll concentration ([Chl]) and dissolved and detrital matter (acdm(443)) in the global ocean from which Rrs(lambda) is calculated by the radiative transfer software, Hydrolight. The Hydrolight modeled Rrs(lambda) options for a given combination of [Chl] and acdm(443) within the LUT vary only due to Sfm. For a given pixel, the LUT search space was limited by satellite imagery of [Chl] and acdm(443). Within the narrowed search space, SeaWiFS Rrs(lambda) was matched with the closest LUT Rrs(lambda) option and the associated Sfm was assigned. Thresholds at which changes in Rrs(lambda) due to Sfm could be discerned were established in terms of [Chl] and acdm(443). In situ high-precision liquid chromatography-derived estimates of cell size are used in conjunction with matched daily satellite estimates of Sfm for validation and agree well. A single month is displayed as an example of the Sfm retrieval.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.