The distribution of sulfur in oceanic crust with respect to oxidative diagenesis (Table 1)

A detailed geochemical-petrological examination of layer 2 basalts recovered during Leg 37 of the DSDP has revealed that the original distribution, form and abundance of igneous sulfide have been profoundly altered during low-grade oxidative diagenesis. The net result appears to have been a rather pervasive remobilization of igneous sulfide to form secondary pyrite accompanied by a bulk loss of sulfur equivalent to about 50-60% of the original igneous value, assuming initial saturation. It is suggested that during infiltration of seawater into the massive crystalline rock, igneous sulfide has experienced pervasive oxidation, under conditions of limited oxidation potential, to form a series of unstable, soluble sulfur species, primarily in the form of SO3[2-] and S2O3[2-]. Spontaneous decomposition of these intermediate compounds through disproportionation has resulted in partial reconstitution of the sulfur as secondary pyrite and the generation of SO4[2-] ion, which, due to its kinetic stability, has been lost from the basalt system and ultimately transferred to the ocean. This model not only satisfies the geochemical and petrological observations but also provides a suitable explanation for the highly variable delta34S values which characterize secondary sulfides in deep ocean floor basalts.

(Table 1) Porewater analysis with respect to gypsum from different ODP holes

The stability of gypsum in marine sediments has been investigated through the calculation of its saturation index at the sediment in situ temperature and pressure, using the entire ODP/IODP porewater composition database (14416 samples recovered from sediments collected during 95 ODP and IODP Legs). Saturation is reached in sediment porewaters of 26 boreholes drilled at 23 different sites, during 12 ODP/IODP Legs. As ocean bottom seawater is largely undersaturated with respect to gypsum, the porewater Ca content or its SO4 concentration, or both, must increase in order to reach equilibrium. At several sites equilibrium is reached either through the presence of evaporitic gypsum layers found in the sedimentary sequence, and/or through a salinity increase due to the presence of evaporitic brines with high concentrations of Ca and SO4. Saturation can also be reached in porewaters of seawater-like salinity (~ 35 per mil), provided sulfate reduction is limited. In this case, saturation is due to the alteration of volcanogenic material which releases large amounts of Ca to the porewaters, where the Ca concentration can reach 55 times its seawater value as for example at ODP Leg 134 site 833. At a few sites, saturation is reached in hydrothermal environments, or as a consequence of the alteration of the basaltic basement. In addition to the well known influence of brines on the formation of gypsum, these results indicate that the alteration of sediments rich in volcanogenic material is a major process leading to gypsum saturation in marine sediment porewaters. Therefore, the presence of gypsum in ancient and recent marine sediments should not be systematically interpreted as due to hypersaline waters, especially if volcanogenic material is present.

30 years of synoptic observations from Neumayer Station with links to datasets

The analysis of research data plays a key role in data-driven areas of science. Varieties of mixed research data sets exist and scientists aim to derive or validate hypotheses to find undiscovered knowledge. Many analysis techniques identify relations of an entire dataset only. This may level the characteristic behavior of different subgroups in the data. Like automatic subspace clustering, we aim at identifying interesting subgroups and attribute sets. We present a visual-interactive system that supports scientists to explore interesting relations between aggregated bins of multivariate attributes in mixed data sets. The abstraction of data to bins enables the application of statistical dependency tests as the measure of interestingness. An overview matrix view shows all attributes, ranked with respect to the interestingness of bins. Complementary, a node-link view reveals multivariate bin relations by positioning dependent bins close to each other. The system supports information drill-down based on both expert knowledge and algorithmic support. Finally, visual-interactive subset clustering assigns multivariate bin relations to groups. A list-based cluster result representation enables the scientist to communicate multivariate findings at a glance. We demonstrate the applicability of the system with two case studies from the earth observation domain and the prostate cancer research domain. In both cases, the system enabled us to identify the most interesting multivariate bin relations, to validate already published results, and, moreover, to discover unexpected relations.

A fractional vegetation cover remote sensing product on pan-arctic scale with link to geotiff image

The paper presents first results of a pan-boreal scale land cover harmonization and classification. A methodology is presented that combines global and regional vegetation datasets to extract percentage cover information for different vegetation physiognomy and barren for the pan-arctic region within the ESA Data User Element Permafrost. Based on the legend description of each land cover product the datasets are harmonized into four LCCS (Land Cover Classification System) classifiers which are linked to the MODIS Vegetation Continuous Field (VCF) product. Harmonized land cover and Vegetation Continuous Fields products are combined to derive a best estimate of percentage cover information for trees, shrubs, herbaceous and barren areas for Russia. Future work will concentrate on the expansion of the developed methodology to the pan-arctic scale. Since the vegetation builds an isolation layer, which protects the permafrost from heat and cold temperatures, a degradation of this layer due to fire strongly influences the frozen conditions in the soil. Fire is an important disturbance factor which affects vast processes and dynamics in ecosystems (e.g. biomass, biodiversity, hydrology, etc.). Especially in North Eurasia the fire occupancy has dramatically increased in the last 50 years and has doubled in the 1990s with respect to the last five decades. A comparison of global and regional fire products has shown discrepancies between the amounts of burn scars detected by different algorithms and satellite data.

Circumarctic wetland dataset based on ENVISAT ASAR Global Monitoring Mode with links to geotiff image and NetCDF files (2005-01 to 2011-12)

A circumpolar representative and consistent wetland map is required for a range of applications ranging from upscaling of carbon fluxes and pools to climate modelling and wildlife habitat assessment. Currently available data sets lack sufficient accuracy and/or thematic detail in many regions of the Arctic. Synthetic aperture radar (SAR) data from satellites have already been shown to be suitable for wetland mapping. Envisat Advanced SAR (ASAR) provides global medium-resolution data which are examined with particular focus on spatial wetness patterns in this study. It was found that winter minimum backscatter values as well as their differences to summer minimum values reflect vegetation physiognomy units of certain wetness regimes. Low winter backscatter values are mostly found in areas vegetated by plant communities typically for wet regions in the tundra biome, due to low roughness and low volume scattering caused by the predominant vegetation. Summer to winter difference backscatter values, which in contrast to the winter values depend almost solely on soil moisture content, show expected higher values for wet regions. While the approach using difference values would seem more reasonable in order to delineate wetness patterns considering its direct link to soil moisture, it was found that a classification of winter minimum backscatter values is more applicable in tundra regions due to its better separability into wetness classes. Previous approaches for wetland detection have investigated the impact of liquid water in the soil on backscatter conditions. In this study the absence of liquid water is utilized. Owing to a lack of comparable regional to circumpolar data with respect to thematic detail, a potential wetland map cannot directly be validated; however, one might claim the validity of such a product by comparison with vegetation maps, which hold some information on the wetness status of certain classes. It was shown that the Envisat ASAR-derived classes are related to wetland classes of conventional vegetation maps, indicating its applicability; 30% of the land area north of the treeline was identified as wetland while conventional maps recorded 1-7%.

Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.

Seawater carbonate chemistry, calcification and respiration duirng experiments with a pteropod Creseis acicula, 2010

Among marine calcifiers, shelled pteropods are expected to be particularly sensitive to ocean acidification, generated by the uptake of anthropogenic CO2 by the ocean, and the associated decrease of the seawater saturation state with respect to aragonite (omega aragonite). The few available studies have mostly focused on polar species although pteropods are also important components of temperate and tropical ecosystems. It is also unknown which parameter of the carbonate system controls calcification. Specimens of the temperate Mediterranean species Creseis acicula were maintained under seven different conditions of the carbonate chemistry, obtained by manipulating pH and total alkalinity, with the goal to disentangle the effects of pH and omega aragonite. Respiration, excretion as well as rates of net and gross calcification were not directly affected by a decrease in pH but decreased significantly with a decrease of omega aragonite. The decrease of gross calcification rates is consistent with that reported for polar species. Although the organisms were apparently able to maintain gross calcification rates under slightly undersaturated aragonite conditions, the clear net dissolution signal observed below saturation suggests that they are not able to build a shell in seawater corrosive to aragonite. The decrease in respiration and excretion, and the low O:N molar ratio, could be due to the short time that the organisms were allowed to acclimatize to their new environment.

A geomorphological seabed classification for the Weddell Sea, Antarctica, with links to ArcGIS project files

Sea floor morphology plays an important role in many scientific disciplines such as ecology, hydrology and sedimentology since geomorphic features can act as physical controls for e.g. species distribution, oceanographically flow-path estimations or sedimentation processes. In this study, we provide a terrain analysis of the Weddell Sea based on the 500 m × 500 m resolution bathymetry data provided by the mapping project IBCSO. Seventeen seabed classes are recognized at the sea floor based on a fine and broad scale Benthic Positioning Index calculation highlighting the diversity of the glacially carved shelf. Beside the morphology, slope, aspect, terrain rugosity and hillshade were calculated. Applying zonal statistics to the geomorphic features identified unambiguously the shelf edge of the Weddell Sea with a width of 45-70 km and a mean depth of about 1200 m ranging from 270 m to 4300 m. A complex morphology of troughs, flat ridges, pinnacles, steep slopes, seamounts, outcrops, and narrow ridges, structures with approx. 5-7 km width, build an approx. 40-70 km long swath along the shelf edge. The study shows where scarps and depressions control the connection between shelf and abyssal and where high and low declination within the scarps e.g. occur. For evaluation purpose, 428 grain size samples were added to the seabed class map. The mean values of mud, sand and gravel of those samples falling into a single seabed class was calculated, respectively, and assigned to a sediment texture class according to a common sediment classification scheme.

At surface behaviour, dive depth and dive duration and associated hydrographic data of adult male southern elephant seals from King George Island from expedition JUB2010 with links to datasets

Understanding the distribution and foraging ecology of major consumers within pelagic systems, specifically in relation to physical parameters, can be important for the management of bentho-pelagic systems undergoing rapid change associated with global climate change and other anthropogenic disturbances such as fishing (i.e., the Antarctic Peninsula and Scotia Sea). We tracked 11 adult male southern elephant seals (Mirounga leonina), during their five-month post-moult foraging migrations from King George Island (Isla 25 de Mayo), northern Antarctic Peninsula, using tags capable of recording and transmitting behavioural data and in situ temperature and salinity data. Seals foraged mostly within the Weddell-Scotia Confluence, while a few foraged along the western Antarctic Peninsula shelf of the Bellingshausen Sea. Mixed model outputs suggest that the at-sea behaviour of seals was associated with a number of environmental parameters, especially seafloor depth, sea-ice concentrations and the temperature structure of the water column. Seals increased dive bottom times and travelled at slower speeds in shallower areas and areas with increased sea-ice concentrations. Changes in dive depth and durations, as well as relative amount of time spent during the bottom phases of dives, were observed in relation to differences in overall temperature gradient, likely as a response to vertical changes in prey distribution associated with temperature stratification in the water column. Our results illustrate the likely complex influences of bathymetry, hydrography and sea ice on the behaviour of male southern elephant seals in a changing environment and highlight the need for region-specific approaches to studying environmental influences on behaviour.

Towards convection-resolving, global atmospheric simulations with the Model for Prediction Across Scales (MPAS) v3.1: an extreme scaling experiment, links to sourcecode of model and grides

The Model for Prediction Across Scales (MPAS) is a novel set of Earth system simulation components and consists of an atmospheric model, an ocean model and a land-ice model. Its distinct features are the use of unstructured Voronoi meshes and C-grid discretisation to address shortcomings of global models on regular grids and the use of limited area models nested in a forcing data set, with respect to parallel scalability, numerical accuracy and physical consistency. This concept allows one to include the feedback of regional land use information on weather and climate at local and global scales in a consistent way, which is impossible to achieve with traditional limited area modelling approaches. Here, we present an in-depth evaluation of MPAS with regards to technical aspects of performing model runs and scalability for three medium-size meshes on four different high-performance computing (HPC) sites with different architectures and compilers. We uncover model limitations and identify new aspects for the model optimisation that are introduced by the use of unstructured Voronoi meshes. We further demonstrate the model performance of MPAS in terms of its capability to reproduce the dynamics of the West African monsoon (WAM) and its associated precipitation in a pilot study. Constrained by available computational resources, we compare 11-month runs for two meshes with observations and a reference simulation from the Weather Research and Forecasting (WRF) model. We show that MPAS can reproduce the atmospheric dynamics on global and local scales in this experiment, but identify a precipitation excess for the West African region. Finally, we conduct extreme scaling tests on a global 3?km mesh with more than 65 million horizontal grid cells on up to half a million cores. We discuss necessary modifications of the model code to improve its parallel performance in general and specific to the HPC environment. We confirm good scaling (70?% parallel efficiency or better) of the MPAS model and provide numbers on the computational requirements for experiments with the 3?km mesh. In doing so, we show that global, convection-resolving atmospheric simulations with MPAS are within reach of current and next generations of high-end computing facilities.