The SUGAR Toolbox

The SUGAR Toolbox contains scripts coded in MATLAB for calculating various thermodynamic, kinetic, and geologic properties of substances occurring in the marine environment, particularly gas hydrate and seep systems. Brief descriptions of the toolbox scripts and some notes on the underlying basic theory as well as tables of additional property values can be found in the accompanying documentation.

Classification of coral reef types in Torres Strait (GIS files in zip-archive 319 kB)

Coral reefs represent major accumulations of calcium carbonate (CaCO3). The particularly labyrinthine network of reefs in Torres Strait, north of the Great Barrier Reef (GBR), has been examined in order to estimate their gross CaCO3 productivity. The approach involved a two-step procedure, first characterising and classifying the morphology of reefs based on a classification scheme widely employed on the GBR and then estimating gross CaCO3 productivity rates across the region using a regional census-based approach. This was undertaken by independently verifying published rates of coral reef community gross production for use in Torres Strait, based on site-specific ecological and morphological data. A total of 606 reef platforms were mapped and classified using classification trees. Despite the complexity of the maze of reefs in Torres Strait, there are broad morphological similarities with reefs in the GBR. The spatial distribution and dimensions of reef types across both regions are underpinned by similar geological processes, sea-level history in the Holocene and exposure to the same wind/wave energetic regime, resulting in comparable geomorphic zonation. However, the presence of strong tidal currents flowing through Torres Strait and the relatively shallow and narrow dimensions of the shelf exert a control on local morphology and spatial distribution of the reef platforms. A total amount of 8.7 million tonnes of CaCO3 per year, at an average rate of 3.7 kg CaCO3 m-2 yr-1 (G), were estimated for the studied area. Extrapolated production rates based on detailed and regional census-based approaches for geomorphic zones across Torres Strait were comparable to those reported elsewhere, particularly values for the GBR based on alkalinity-reduction methods. However, differences in mapping methodologies and the impact of reduced calcification due to global trends in coral reef ecological decline and changing oceanic physical conditions warrant further research. The novel method proposed in this study to characterise the geomorphology of reef types based on classification trees provides an objective and repeatable data-driven approach that combined with regional census-based approaches has the potential to be adapted and transferred to different coral reef regions, depicting a more accurate picture of interactions between reef ecology and geomorphology.

Probabilistic models for the prediction of a ship performance in dynamic ice

We introduce two probabilistic, data-driven models that predict a ship's speed and the situations where a ship is probable to get stuck in ice based on the joint effect of ice features such as the thickness and concentration of level ice, ice ridges, rafted ice, moreover ice compression is considered. To develop the models to datasets were utilized. First, the data from the Automatic Identification System about the performance of a selected ship was used. Second, a numerical ice model HELMI, developed in the Finnish Meteorological Institute, provided information about the ice field. The relations between the ice conditions and ship movements were established using Bayesian learning algorithms. The case study presented in this paper considers a single and unassisted trip of an ice-strengthened bulk carrier between two Finnish ports in the presence of challenging ice conditions, which varied in time and space. The obtained results show good prediction power of the models. This means, on average 80% for predicting the ship's speed within specified bins, and above 90% for predicting cases where a ship may get stuck in ice. We expect this new approach to facilitate the safe and effective route selection problem for ice-covered waters where the ship performance is reflected in the objective function.

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.

Last glacial maximum sediment record of the southern Weddell Sea shelf

Abstract: The history of grounded ice-sheet extent on the southern Weddell Sea shelf during the Last Glacial Maximum (LGM) and the timing of post-LGM ice-sheet retreat are poorly constrained. Several glaciological models reconstructed widespread grounding and major thickening of the Antarctic Ice Sheet in the Weddell Sea sector at the LGM. In contrast, recently published onshore data and modelling results concluded only very limited LGM-thickening of glaciers and ice streams feeding into the modern Filchner and Ronne ice shelves. These studies concluded that during the LGM ice shelves rather than grounded ice covered the Filchner and Ronne troughs, two deep palaeo-ice stream troughs eroded into the southern Weddell Sea shelf. Here we review previously published and unpublished marine geophysical and geological data from the southern Weddell Sea shelf. The stratigraphy and geometry of reflectors in acoustic sub-bottom profiles are similar to those from other West Antarctic palaeo-ice stream troughs, where grounded ice had advanced to the shelf break at the LGM. Numerous cores from the southern Weddell Sea shelf recovered sequences with properties typical for subglacially deposited tills or subglacially compacted sediments. These data sets give evidence that grounded ice had advanced across the shelf during the past, thereby grounding in even the deepest parts of the Filchner and Ronne troughs. Radiocarbon dates from glaciomarine sediments overlying the subglacial deposits are limited, but indicate that the ice grounding occurred at the LGM and that ice retreat started before ~15.1 corrected 14C kyrs before present (BP) on the outer shelf and before ~7.7 corrected 14C kyrs BP on the inner shelf, which is broadly synchronous with ice retreat in other Antarctic sectors. The apparent mismatch between the ice-sheet reconstructions from marine and terrestrial data can be attributed to ice streams with very low surface profiles (similar to those of "ice plains") that had advanced through Filchner Trough and Ronne Trough at the LGM. Considering the global sea-level lowstand of ~130 metres below present, a low surface slope of the expanded LGM-ice sheet in the southern Weddell Sea can reconcile grounding-line advance to the shelf break with limited thickening of glaciers and ice streams in the hinterland. This scenario implies that ice-sheet growth in the Weddell Sea sector during the LGM and ice-sheet drawdown throughout the last deglaciation could only have made minor contributions to the major global sea-level fluctuations during these times.

The SUGAR Toolbox - Version 7

The SUGAR Toolbox contains scripts coded in MATLAB for calculating various thermodynamic, kinetic, and geologic properties of substances occurring in the marine environment, particularly gas hydrate and seep systems. Brief descriptions of the toolbox scripts and some notes on the underlying basic theory as well as tables of additional property values can be found in the accompanying documentation.

Stable carbon isotope ratios and accumulation rates of organic carbon and nutrients of DSDP Hole 79-545 (Appendix A)

Future warming is predicted to shift the Earth system into a mode with progressive increase and vigour of extreme climate events possibly stimulating other mechanisms that invigorate global warming. This study provides new data and modelling investigating climatic consequences and biogeochemical feedbacks that happened in a warmer world ~112 Myr ago. Our study focuses on the Cretaceous Oceanic Anoxic Event (OAE) 1b and explores how the Earth system responded to a moderate ~25,000 yr lasting climate perturbation that is modelled to be less than 1 °C in global average temperature. Using a new chronological model for OAE 1b we present high-resolution elemental and bulk carbon isotope records from DSDP Site 545 from Mazagan Plateau off NW Africa and combine this information with a coupled atmosphere-land-ocean model. The simulations suggest that a perturbation at the onset of OAE 1b caused almost instantaneous warming of the atmosphere on the order of 0.3 °C followed by a longer (~45,000 yr) period of ~0.8 °C cooling. The marine records from DSDP Site 545 support that these moderate swings in global climate had immediate consequences for African continental supply of mineral matter and nutrients (phosphorous), subsequent oxygen availability, and organic carbon burial in the eastern subtropical Atlantic, however, without turning the ocean anoxic. The match between modelling results and stratigraphic isotopic data support previous studies [summarized in Jenkyns 2003, doi:10.1098/rsta.2003.1240] in that methane emission from marine hydrates, albeit moderate in dimension, may have been the trigger for OAE 1b, though we can not finally rule out alternative mechanisms. Following the hydrate mechanism a total of 1.15 * 10**18 g methane carbon (delta13C=-60 ?), equivalent to about 10% to the total modern gas hydrate inventory, generated the delta13Ccarb profile recorded in the section. Modelling suggests a combination of moderate-scale methane pulses supplemented by continuous methane emission at elevated levels over ~25,000 yr. The proposed mechanism, though difficult to finally confirm in the geological past, is arguably more likely to occur in a warmer world and apparently perturbs global climate and ocean chemistry almost instantaneously. This study shows that, once set-off, this mechanism can maintain Earth's climate in a perturbed mode over geological time leading to pronounced changes in regional climate.

Corrected ICESat altimetry data, surface mass balance, and firn elevation change on Antarctic ice shelves

Accurate prediction of global sea-level rise requires that we understand the cause of recent, widespread and intensifying glacier acceleration along Antarctic ice-sheet coastal margins. Floating ice shelves buttress the flow of grounded tributary glaciers and their thickness and extent are particularly susceptible to changes in both climate and ocean forcing. Recent ice-shelf collapse led to retreat and acceleration of several glaciers on the Antarctic Peninsula. However, the extent and magnitude of ice-shelf thickness change, its causes and its link to glacier flow rate are so poorly understood that its influence on the future of the ice sheets cannot yet be predicted. Here we use satellite laser altimetry and modelling of the surface firn layer to reveal for the first time the circum-Antarctic pattern of ice-shelf thinning through increased basal melt. We deduce that this increased melt is the primary driver of Antarctic ice-sheet loss, through a reduction in buttressing of the adjacent ice sheet that has led to accelerated glacier flow. The highest thinning rates (~7 m/a) occur where warm water at depth can access thick ice shelves via submarine troughs crossing the continental shelf. Wind forcing could explain the dominant patterns of both basal melting and the surface melting and collapse of Antarctic ice shelves, through ocean upwelling in the Amundsen and Bellingshausen Seas and atmospheric warming on the Antarctic Peninsula. This implies that climate forcing through changing winds influences Antarctic Ice Sheet mass balance, and hence global sea-level, on annual to decadal timescales.