Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

Additional CaCO3 measurements from IODP Site 306-U1313

We present further %CaCO3 data from Site U1313 across the Pliocene-Pleistocene intensification of Northern Hemisphere glaciation. This data was measured on the U1313 secondary splice. We also present tie points between the primary and secondary splice for this interval based on graphical tuning of L* (sediment lightness).

Quantitative component analysis of the coarse fraction of surface sediments from the Persian Gulf

In the Persian Gulf and the Gulf of Oman marl forms the primary sediment cover, particularly on the Iranian side. A detailed quantitative description of the sediment components > 63 µ has been attempted in order to establish the regional distribution of the most important constituents as well as the criteria governing marl sedimentation in general. During the course of the analysis, the sand fraction from about 160 bottom-surface samples was split into 5 phi° fractions and 500 to 800 grains were counted in each individual fraction. The grains were cataloged in up to 40 grain type catagories. The gravel fraction was counted separately and the values calculated as weight percent. Basic for understanding the mode of formation of the marl sediment is the "rule" of independent availability of component groups. It states that the sedimentation of different component groups takes place independently, and that variation in the quantity of one component is independent of the presence or absence of other components. This means, for example, that different grain size spectrums are not necessarily developed through transport sorting. In the Persian Gulf they are more likely the result of differences in the amount of clay-rich fine sediment brought in to the restricted mouth areas of the Iranian rivers. These local increases in clayey sediment dilute the autochthonous, for the most part carbonate, coarse fraction. This also explains the frequent facies changes from carbonate to clayey marl. The main constituent groups of the coarse fraction are faecal pellets and lumps, the non carbonate mineral components, the Pleistocene relict sediment, the benthonic biogene components and the plankton. Faecal pellets and lumps are formed through grain size transformation of fine sediment. Higher percentages of these components can be correlated to large amounts of fine sediment and organic C. No discernable change takes place in carbonate minerals as a result of digestion and faecal pellet formation. The non-carbonate sand components originate from several unrelated sources and can be distinguished by their different grain size spectrum; as well as by other characteristics. The Iranian rivers supply the greatest amounts (well sorted fine sand). Their quantitative variations can be used to trace fine sediment transport directions. Similar mineral maxima in the sediment of the Gulf of Oman mark the path of the Persian Gulf outflow water. Far out from the coast, the basin bottoms in places contain abundant relict minerals (poorly sorted medium sand) and localized areas of reworked salt dome material (medium sand to gravel). Wind transport produces only a minimal "background value" of mineral components (very fine sand). Biogenic and non-biogenic relict sediments can be placed in separate component groups with the help of several petrographic criteria. Part of the relict sediment (well sorted fine sand) is allochthonous and was derived from the terrigenous sediment of river mouths. The main part (coarse, poorly sorted sediment), however, was derived from the late Pleistocene and forms a quasi-autochthonous cover over wide areas which receive little recent sedimentation. Bioturbation results in a mixing of the relict sediment with the overlying younger sediment. Resulting vertical sediment displacement of more than 2.5 m has been observed. This vertical mixing of relict sediment is also partially responsible for the present day grain size anomalies (coarse sediment in deep water) found in the Persian Gulf. The mainly aragonitic components forming the relict sediment show a finely subdivided facies pattern reflecting the paleogeography of carbonate tidal flats dating from the post Pleistocene transgression. Standstill periods are reflected at 110 -125m (shelf break), 64-61 m and 53-41 m (e.g. coare grained quartz and oolite concentrations), and at 25-30m. Comparing these depths to similar occurrences on other shelf regions (e. g. Timor Sea) leads to the conclusion that at this time minimal tectonic activity was taking place in the Persian Gulf. The Pleistocene climate, as evidenced by the absence of Iranian river sediment, was probably drier than the present day Persian Gulf climate. Foremost among the benthonic biogene components are the foraminifera and mollusks. When a ratio is set up between the two, it can be seen that each group is very sensitive to bottom type, i.e., the production of benthonic mollusca increases when a stable (hard) bottom is present whereas the foraminifera favour a soft bottom. In this way, regardless of the grain size, areas with high and low rates of recent sedimentation can be sharply defined. The almost complete absence of mollusks in water deeper than 200 to 300 m gives a rough sedimentologic water depth indicator. The sum of the benthonic foraminifera and mollusca was used as a relative constant reference value for the investigation of many other sediment components. The ratio between arenaceous foraminifera and those with carbonate shells shows a direct relationship to the amount of coarse grained material in the sediment as the frequence of arenaceous foraminifera depends heavily on the availability of sand grains. The nearness of "open" coasts (Iranian river mouths) is directly reflected in the high percentage of plant remains, and indirectly by the increased numbers of ostracods and vertebrates. Plant fragments do not reach their ultimate point of deposition in a free swimming state, but are transported along with the remainder of the terrigenous fine sediment. The echinoderms (mainly echinoids in the West Basin and ophiuroids in the Central Basin) attain their maximum development at the greatest depth reached by the action of the largest waves. This depth varies, depending on the exposure of the slope to the waves, between 12 to 14 and 30 to 35 m. Corals and bryozoans have proved to be good indicators of stable unchanging bottom conditions. Although bryozoans and alcyonarian spiculae are independent of water depth, scleractinians thrive only above 25 to 30 m. The beginning of recent reef growth (restricted by low winter temperatures) was seen only in one single area - on a shoal under 16 m of water. The coarse plankton fraction was studied primarily through the use of a plankton-benthos ratio. The increase in planktonic foraminifera with increasing water depth is here heavily masked by the "Adjacent sea effect" of the Persian Gulf: for the most part the foraminifera have drifted in from the Gulf of Oman. In contrast, the planktonic mollusks are able to colonize the entire Persian Gulf water body. Their amount in the plankton-benthos ratio always increases with water depth and thereby gives a reliable picture of local water depth variations. This holds true to a depth of around 400 m (corresponding to 80-90 % plankton). This water depth effect can be removed by graphical analysis, allowing the percentage of planktonic mollusks per total sample to be used as a reference base for relative sedimentation rate (sedimentation index). These values vary between 1 and > 1000 and thereby agree well with all the other lines of evidence. The "pteropod ooze" facies is then markedly dependent on the sedimentation rate and can theoretically develop at any depth greater than 65 m (proven at 80 m). It should certainly no longer be thought of as "deep sea" sediment. Based on the component distribution diagrams, grain size and carbonate content, the sediments of the Persian Gulf and the Gulf of Oman can be grouped into 5 provisional facies divisions (Chapt.19). Particularly noteworthy among these are first, the fine grained clayey marl facies occupying the 9 narrow outflow areas of rivers, and second, the coarse grained, high-carbonate marl facies rich in relict sediment which covers wide sediment-poor areas of the basin bottoms. Sediment transport is for the most part restricted to grain sizes < 150 µ and in shallow water is largely coast-parallel due to wave action at times supplemented by tidal currents. Below the wave base gravity transport prevails. The only current capable of moving sediment is the Persian Gulf outflow water in the Gulf of Oman.

Tritium concentrations of hydrography in the North Atlantic around 40°N

We report, numerically and in graphical form, meaured tritium concentrations from five hydrographic stations in the North Atlantic. Fairly homogeneous concentrations are observed in a surface layer typically 400 m deep. In the thermocline, concentrations decrease steadily down to the sigma-theta = 27.3 density horizon, and are more variable further down. The tritium in the lower part of the thermocline originates from the Subarctic Intermediate Water and the Mediterranean Water. There is a relative tritium maximum associated with the Mediterranean Water on the easternmost station of the section. In the deep water (sigma-theta > 27.8), concentrations east of the Midatlantic Ridge are close to the limit of detection down to 2500m, and undetectable further down, while west of the ridge tritium is found throughout the water column. The deep water tritium is associated with the deep-water advective cores of Arctic origin. The present tritium data can serve as northern boundary values in attempts to use tritium in studies of the North Atlantic main thermocline dynamics. The present data together with data from the literature point to a general division of the North Atlantic main thermocline into two layers separated by an isopycnal surface near sigma-theta = 27.3.

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.

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.

Characteristics of the Vernagtferner mass balance for the period from 1964 to 2014

The main characteristics of the Vernagtferner mass balance are sumarized in the table below. The mass balance years from 1964/65 to 2003/2004 are listed. The table includes the total area of the glacier (basis for the calculations), the equilibrium line altitude (ELA), percentage of the accumulation area in relation to the total area (AAR) and the specific net mass balance in mm w.e. (water equivalent) per year. It becomes clear that, after a rather minor growth period in the mid 1970's, the glacier continually lost mass since the beginning of the 1980's. Besides that, a clear increase of mass balance years with extreme mass losses could be observed in the last decade. The "glacier-friendly" summer with a well-balanced mass balance in 1999 could only interrupt the series of years with extreme mass losses, but this means no change in the trend. The minor mass loss in 1999 was caused by a winter snow cover above average, which prevented the glacier from becoming snow free over large areas and thus resulted in a lower ice melt. Although real summer conditions in 2000 were mainly restricted to August and produced a snow free area only slightly larger than in 1999, there have been further ice losses. This trend of negative mass balance continued also in the years 2001 and 2002. Nevertheless, the losses are moderate because a smaller part of the glacier became ice free until autumn (appr. 50 %). The summer 2003 caused a loss of ice in a dimension never seen since the beginning of the scientific investigations. This resulted from a combination of different factors: after only a moderate winter snowcover the glacier became snow free very early. For the first time the ablation area spanned over the entire glacier (blue fields in the mass balance tables!). Only one short snowfall event interrupted the ablation period, which lasted twice as long as in the years of large losses in the 1990's. The extreme mass loss in 2003 will also influence the mass balance in the following year 2004. The graphical representation of the elevation distribution of the specific mass balance together with the absolute mass balance can be found individually for each year by choosing one of the mass balance values from the table. These diagrams also include the area-height-distribution of the glacier and the ablation area. A tabular version of the numeric values in dependence of the elevation, provided separately for the accumulation area, the ablation area and the total glacier, can be found in colums "Persistent Identifier". The tables include the results for three different parts of the glacier and for the total glacier.

Large-Area Photo-Mosaicking tool (LAPM tool)

This paper presents a new tool for large-area photo-mosaicking (LAPM tool). This tool was developed specifically for the purpose of underwater mosaicking, and it is aimed at providing end-user scientists with an easy and robust way to construct large photo-mosaics from any set of images. It is notably capable of constructing mosaics with an unlimited number of images on any modern computer (minimum 1.30 GHz, 2 GB RAM). The mosaicking process can rely on both feature matching and navigation data. This is complemented by an intuitive graphical user interface, which gives the user the ability to select feature matches between any pair of overlapping images. Finally, mosaic files are given geographic attributes that permit direct import into ArcGIS. So far, the LAPM tool has been successfully used to construct geo-referenced photo-mosaics with photo and video material from several scientific cruises. The largest photo-mosaic contained more than 5000 images for a total area of about 105,000 m**2. This is the first article to present and to provide a finished and functional program to construct large geo-referenced photo-mosaics of the seafloor using feature detection and matching techniques. It also presents concrete examples of photo-mosaics produced with the LAPM tool.

Grain size investigations of sediments off the coast of Istria

Composition, grain-size distribution, and areal extent of Recent sediments from the Northern Adriatic Sea along the Istrian coast have been studied. Thirty one stations in four sections vertical to the coast were investigated; for comparison 58 samples from five small bays were also analyzed. Biogenic carbonate sediments are deposited on the shallow North Adriatic shelf off the Istrian coast. Only at a greater distance from the coast are these carbonate sediments being mixed with siliceous material brought in by the Alpine rivers Po, Adige, and Brenta. Graphical analysis of grain-size distribution curves shows a sediment composition of normally three, and only in the most seaward area, of four major constituents. Constituent 1 represents the washed-in terrestrial material of clay size (Terra Rossa) from the Istrian coastal area. Constituent 2 consists of fine to medium sand. Constituent 3 contains the heterogeneous biogenic material. Crushing by organisms and by sediment eaters reduces the coarse biogenic material into small pieces generating constituent 2. Between these two constituents there is a dynamic equilibrium. Depending upon where the equilibrium is, between the extremes of production and crushing, the resulting constituent 2 is finer or coarser. Constituent 4 is composed of the fine sandy material from the Alpine rivers. In the most seaward area constituents 2 and 4 are mixed. The total carbonate content of the samples depends on the distance from the coast. In the near coastal area in high energy environments, the carbonate content is about 80 %. At a distance of 2 to 3 km from the coast there is a carbonate minimum because of the higher rate of sedimentation of clay-sized terrestrial, noncarbonate material at extremely low energy environments. In an area between 5 and 20 km off the coast, the carbonate content is about 75 %. More than 20 km from the shore, the carbonate content diminishes rapidly to values of about 30 % through mixing with siliceous material from the Alpine rivers. The carbonate content of the individual fractions increases with increasing grain-size to a maximum of about 90 % within the coarse sand fractions. Beyond 20 km from the coast the samples show a carbonate minimum of about 13 % within the sand-size classes from 1.5 to 0.7 zeta¬? through mixing with siliceous material from the alpine rivers. By means of grain-size distribution and carbonate content, four sediment zones parallel to the coast were separated. Genetically they are closely connected with the zonation of the benthic fauna. Two cores show a characteristic vertical distribution of the sediment. The surface zone is inversely graded, that means the coarse fractions are at the top and the fine fractions are at the bottom. This is the effect of crushing of the biogenic material produced at the surface by predatory organisms and by sediment eaters. lt is proposed that at a depth of about 30 cm a chemical solution process begins which leads to diminution of the original sediment from a fine to medium sand to a silt. The carbonate content decreases from about 75 % at the surface to 65 % at a depth of 100 cm. The increase of the noncarbonate components by 10 % corresponds to a decrease in the initial amount of sediment (CaC03=75 %) by roughly 30 % through solution. With increasing depth the carbonate content of the individual fractions becomes more and more uniform. At the surface the variation is from 30 % to 90 %, at the bottom it varies only between 50 % and 75 %. Comparable investigations of small-bay sediments showed a c1ear dependence of sediment/faunal zonation from the energy of the environment. The investigations show that the composition and three-dimensional distribution of the Istrian coastal sediments can not be predicted only from one or a few measurable factors. Sedimentation and syngenetic changes must be considered as a complex interaction between external factors and the actions of producing and destroying organisms that are in dynamic equilibrium. The results obtained from investigations of these recent sediments may be of value for interpreting fossil sediments only with strong limitations.

Granulometry and Bulk Geochemistry of Nearshore Sediments, Herschel Island, Yukon Territory, Canada

This dataset contains the results of granulometric and bulk geochemical analyses of Van Veen surface samples obtained by the Alfred Wegener Institute (AWI) in the course of the 2012 and 2013 summer field seasons. The sampling was performed along transects in depths generally <13 m, to a distance of about <5 km off Herschel Island. In 2012, 75 samples in Pauline Cove and in the vicinity of Simpson Point were obtained. Sample collection was expanded in 2013, on transects established the previous year, with additional locations in Tetris Bay and Workboat Passage. Samples consisted of approximately 100 g of the top 3-6 cm of sediment, and were frozen in the field and freeze dried at the AWI before undergoing analytical procedures. Sample locations were recorded with the onboard global positioning system (GPS) unit. Grain size distributions in our study were obtained using laser diffractometry at the AWI (Beckman Coulter LS200) on the <1 mm fraction of samples oxidized with 30% H2O2 until effervescence ceased to remove organics. Some samples were also sieved using a sieve stack with 1 phi intervals. GRADISTAT (Blott and Pye, 2001) was used to calculate graphical grain size statistics (Folk and Ward, 1957). Grain diameters were logarithmically transformed to phi values, calculated as phi=-log2d, where d is the grain diameter in millimeters (Blott and Pye, 2001; Krumbein, 1934). Freeze dried samples were ground and ground using an Elemetar Vario EL III carbon-nitrogen-sulphur analyzer at the AWI to measure total carbon (TC) and total nitrogen (TN). Tungsten oxide was added to the samples as a catalyst to the pyrolysis. Following this analysis, total organic carbon (TOC) was determined using an Elementar VarioMax. Stable carbon isotope ratios of 13C/12C of 118 samples were determined on a DELTAplusXL mass spectrometer (ThermoFisher Scientific, Bremen) at the German Research Centre for Geosciences (GFZ) in Potsdam, Germany . An additional analysis on 69 samples was carried out at the University of Hamburg with an isotope ratio mass spectrometer (Delta V, Thermo Scientific, Germany) coupled to an elemental analyzer (Flash 2000, Thermo Scientific, Germany). Prior to analysis, soil samples were treated with phosphoric acid (43%) to release inorganic carbon. Values are expressed relative to Vienna Peedee belemnite (VPDB) using external standards (USGS40, -26.4 per mil VPDB and IVA soil 33802153, -27.5 per mil VPDB).

Kiel fjord pCO2 datasets between 2012 (July) and 2015 (January) measured using a HydroC® pCO2 sensor

The HydroC® CO2 sensor was deployed from a pontoon at the waterfront of the GEOMAR west shore building into Kiel Fjord, Western Baltic Sea (Kiel, Germany; 54°19'48.78"N, 010° 8'59.44"E). Since the pontoon is floating the deployment depth of the sensor was constant at 1m. Data of three deployment intervals are published here: 1) July 2012 - December 2012 2) April 2013 - June 2013 3) November 2013 - January 2015 Data are processed and corrected, for documentation and graphical overview see further details.