Gridded bathymetry from multibeam echosounder EM122 data of the cruise M84/2, off Eregli (2011)

We investigated gas bubble emissions at the Don-Kuban paleo-fan in the northeastern Black Sea regarding their geological setting, quantities as well as spatial and temporal variabilities during three ship expeditions between 2007 and 2011. About 600 bubble-induced hydroacoustic anomalies in the water column (flares) originating from the seafloor above the gas hydrate stability zone (GHSZ) at ~700 m water depth were found. At about 890 m water depth a hydrocarbon seep area named "Kerch seep area" was newly discovered within the GHSZ. We propose locally domed sediments ('mounds') discovered during ultra-high resolution bathymetric mapping with an autonomous underwater vehicle (AUV) to result from gas hydrate accumulation at shallow depths. In situ measurements indicated spatially limited temperature elevations in the shallow sediment likely induced by upward fluid flow which may confine the local GHSZ to a few meters below the seafloor. As a result, gas bubbles are suspected to migrate into near-surface sediments and to escape the seafloor through small-scale faults. Hydroacoustic surveys revealed that several flares originated from a seafloor area of about 1 km**2 in size. The highest flare disappeared in about 350 m water depth, suggesting that the released methane remains in the water column. A methane flux estimate, combining data from visual quantifications during dives with a remotely operated vehicle (ROV) with results from ship-based hydroacoustic surveys and gas analysis revealed that between 2 and 87 x 10**6 mol CH4 yr-1 escaped into the water column above the Kerch seep area. Our results show that the finding of the Kerch seep area represents a so far underestimated type of hydrocarbon seep, which has to be considered in methane budget calculations.

Mesozooplankton size data from the fjord branch Kapisigdlit located in the Godthaabsfjord system, West Greenland, 2010

Sampling was conducted from March 24 to August 5 2010, in the fjord branch Kapisigdlit located in the inner part of the Godthåbsfjord system, West Greenland. The vessel "Lille Masik" was used during all cruises except on June 17-18 where sampling was done from RV Dana (National Institute for Aquatic Resources, Denmark). A total of 15 cruises (of 1-2 days duration) 7-10 days apart was carried out along a transect composed of 6 stations (St.), spanning the length of the 26 km long fjord branch. St. 1 was located at the mouth of the fjord branch and St. 6 was located at the end of the fjord branch, in the middle of a shallower inner creek . St. 1-4 was covering deeper parts of the fjord, and St. 5 was located on the slope leading up to the shallow inner creek. Mesozooplankton was sampled by vertical net tows using a Hydrobios Multinet (type Mini) equipped with a flow meter and 50 µm mesh nets or a WP-2 net 50 µm mesh size equipped with a non-filtering cod-end. Sampling was conducted at various times of day at the different stations. The nets were hauled with a speed of 0.2-0.3 m s**-1 from 100, 75 and 50 m depth to the surface at St. 2 + 4, 5 and 6, respectively. The content was immediately preserved in buffered formalin (4% final concentration). All samples were analyzed in the Plankton sorting and identification center in Szczecin (www.nmfri.gdynia.pl). Samples containing high numbers of zooplankton were split into subsamples. All copepods and other zooplankton were identified down to lowest possible taxonomic level (approx. 400 per sample), length measured and counted. Copepods were sorted into development stages (nauplii stage 1 - copepodite stage 6) using morphological features and sizes, and up to 10 individuals of each stage was length measured.

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.

Mesozooplankton abundance data from the fjord branch Kapisigdlit located in the Godthaabsfjord system, West Greenland, 2010

Sampling was conducted from March 24 to August 5 2010, in the fjord branch Kapisigdlit located in the inner part of the Godthåbsfjord system, West Greenland. The vessel "Lille Masik" was used during all cruises except on June 17-18 where sampling was done from RV Dana (National Institute for Aquatic Resources, Denmark). A total of 15 cruises (of 1-2 days duration) 7-10 days apart was carried out along a transect composed of 6 stations (St.), spanning the length of the 26 km long fjord branch. St. 1 was located at the mouth of the fjord branch and St. 6 was located at the end of the fjord branch, in the middle of a shallower inner creek . St. 1-4 was covering deeper parts of the fjord, and St. 5 was located on the slope leading up to the shallow inner creek. Mesozooplankton was sampled by vertical net tows using a Hydrobios Multinet (type Mini) equipped with a flow meter and 50 µm mesh nets or a WP-2 net 50 µm mesh size equipped with a non-filtering cod-end. Sampling was conducted at various times of day at the different stations. The nets were hauled with a speed of 0.2-0.3 m s**-1 from 100, 75 and 50 m depth to the surface at St. 2 + 4, 5 and 6, respectively. The content was immediately preserved in buffered formalin (4% final concentration). All samples were analyzed in the Plankton sorting and identification center in Szczecin (www.nmfri.gdynia.pl). Samples containing high numbers of zooplankton were split into subsamples. All copepods and other zooplankton were identified down to lowest possible taxonomic level (approx. 400 per sample), length measured and counted. Copepods were sorted into development stages (nauplii stage 1 - copepodite stage 6) using morphological features and sizes, and up to 10 individuals of each stage was length measured.

Probability inundation map for a scenario combining storm surge and sea-level rise

Coastal managers require reliable spatial data on the extent and timing of potential coastal inundation, particularly in a changing climate. Most sea level rise (SLR) vulnerability assessments are undertaken using the easily implemented bathtub approach, where areas adjacent to the sea and below a given elevation are mapped using a deterministic line dividing potentially inundated from dry areas. This method only requires elevation data usually in the form of a digital elevation model (DEM). However, inherent errors in the DEM and spatial analysis of the bathtub model propagate into the inundation mapping. The aim of this study was to assess the impacts of spatially variable and spatially correlated elevation errors in high-spatial resolution DEMs for mapping coastal inundation. Elevation errors were best modelled using regression-kriging. This geostatistical model takes the spatial correlation in elevation errors into account, which has a significant impact on analyses that include spatial interactions, such as inundation modelling. The spatial variability of elevation errors was partially explained by land cover and terrain variables. Elevation errors were simulated using sequential Gaussian simulation, a Monte Carlo probabilistic approach. 1,000 error simulations were added to the original DEM and reclassified using a hydrologically correct bathtub method. The probability of inundation to a scenario combining a 1 in 100 year storm event over a 1 m SLR was calculated by counting the proportion of times from the 1,000 simulations that a location was inundated. This probabilistic approach can be used in a risk-aversive decision making process by planning for scenarios with different probabilities of occurrence. For example, results showed that when considering a 1% probability exceedance, the inundated area was approximately 11% larger than mapped using the deterministic bathtub approach. The probabilistic approach provides visually intuitive maps that convey uncertainties inherent to spatial data and analysis.

Crete fault scarp point cloud data

This paper assesses the along strike variation of active bedrock fault scarps using long range terrestrial laser scanning (t-LiDAR) data in order to determine the distribution behaviour of scarp height and the subsequently calculate long term throw-rates. Five faults on Cretewhich display spectacular limestone fault scarps have been studied using high resolution digital elevation model (HRDEM) data. We scanned several hundred square metres of the fault system including the footwall, fault scarp and hanging wall of the investigated fault segment. The vertical displacement and the dip of the scarp were extracted every metre along the strike of the detected fault segment based on the processed HRDEM. The scarp variability was analysed by using statistical and morphological methods. The analysis was done in a geographical information system (GIS) environment. Results show a normal distribution for the scanned fault scarp's vertical displacement. Based on these facts, the mean value of height was chosen to define the authentic vertical displacement. Consequently the scarp can be divided into above, below and within the range of mean (within one standard deviation) and quantify the modifications of vertical displacement. Therefore, the fault segment can be subdivided into areas which are influenced by external modification like erosion and sedimentation processes. Moreover, to describe and measure the variability of vertical displacement along strike the fault, the semi-variance was calculated with the variogram method. This method is used to determine how much influence the external processes have had on the vertical displacement. By combining of morphological and statistical results, the fault can be subdivided into areas with high external influences and areas with authentic fault scarps, which have little or no external influences. This subdivision is necessary for long term throw-rate calculations, because without this differentiation the calculated rates would be misleading and the activity of a fault would be incorrectly assessed with significant implications for seismic hazard assessment since fault slip rate data govern the earthquake recurrence. Furthermore, by using this workflow areas with minimal external influences can be determined, not only for throw-rate calculations, but also for determining samples sites for absolute dating techniques such as cosmogenic nuclide dating. The main outcomes of this study include: i) there is no direct correlation between the fault's mean vertical displacement and dip (R² less than 0.31); ii) without subdividing the scanned scarp into areas with differing amounts of external influences, the along strike variability of vertical displacement is ±35%; iii) when the scanned scarp is subdivided the variation of the vertical displacement of the authentic scarp (exposed by earthquakes only) is in a range of ±6% (the varies depending on the fault from 7 to 12%); iv) the calculation of the long term throw-rate (since 13 ka) for four scarps in Crete using the authentic vertical displacement is 0.35 ± 0.04 mm/yr at Kastelli 1, 0.31 ± 0.01 mm/yr at Kastelli 2, 0.85 ± 0.06 mm/yr at the Asomatos fault (Sellia) and 0.55 ± 0.05 mm/yr at the Lastros fault.

Mesozooplankton biomass data from the Kapisigdlit an inner fjord branch of the Godthaabsfjord system, West Greenland, 2010

Sampling was conducted from March 24 to August 5 2010, in the fjord branch Kapisigdlit located in the inner part of the Godthåbsfjord system, West Greenland. The vessel "Lille Masik" was used during all cruises except on June 17-18 where sampling was done from RV Dana (National Institute for Aquatic Resources, Denmark). A total of 15 cruises (of 1-2 days duration) 7-10 days apart was carried out along a transect composed of 6 stations (St.), spanning the length of the 26 km long fjord branch. St. 1 was located at the mouth of the fjord branch and St. 6 was located at the end of the fjord branch, in the middle of a shallower inner creek . St. 1-4 was covering deeper parts of the fjord, and St. 5 was located on the slope leading up to the shallow inner creek. Mesozooplankton was sampled by vertical net tows using a Hydrobios Multinet (type Mini) equipped with a flow meter and 50 µm mesh nets or a WP-2 net 50 µm mesh size equipped with a non-filtering cod-end. Sampling was conducted at various times of day at the different stations. The nets were hauled with a speed of 0.2-0.3 m s**-1 from 100, 75 and 50 m depth to the surface at St. 2 + 4, 5 and 6, respectively. The content was immediately preserved in buffered formalin (4% final concentration). All samples were analyzed in the Plankton sorting and identification center in Szczecin (www.nmfri.gdynia.pl). Samples containing high numbers of zooplankton were split into subsamples. All copepods and other zooplankton were identified down to lowest possible taxonomic level (approx. 400 per sample), length measured and counted. Copepods were sorted into development stages (nauplii stage 1 - copepodite stage 6) using morphological features and sizes, and up to 10 individuals of each stage was length measured.

Experimental data on photophysiology and growth rates of Pararotalia calcariformata

The eastern Mediterranean is a hotspot of biological invasions. Numerous species of Indo-pacific origin have colonized the Mediterranean in recent times, including tropical symbiont-bearing foraminifera. Among these is the species Pararotalia calcariformata. Unlike other invasive foraminifera, this species has been discovered only two decades ago and is restricted to the eastern Mediterranean coast. Combining ecological, genetic and physiological observations, we attempt to explain the recent invasion of this species in the Mediterranean Sea. Using morphological and genetic data, we confirm the species attribution to P. calcariformata McCulloch 1977 and identify its symbionts as a consortium of diatom species dominated by Minutocellus polymorphus. We document photosynthetic activity of its endosymbionts using Pulse Amplitude Modulated Fluorometry and test the effects of elevated temperatures on growth rates of asexual offspring. The culturing of asexual offspring for 120 days shows a 30-day period of rapid growth followed by a period of slower growth. A subsequent 48-day temperature sensitivity experiment indicates a similar developmental pathway and high growth rate at 28°C, whereas an almost complete inhibition of growth was observed at 20°C and 35°C. This indicates that the offspring of this species may have lower tolerance to cold temperatures than what would be expected for species native to the Mediterranean. We expand this hypothesis by applying a Species Distribution Model (SDM) based on modern occurrences in the Mediterranean using three environmental variables: irradiance, turbidity and yearly minimum temperature. The model reproduces the observed restricted distribution and indicates that the range of the species will drastically expand westwards under future global change scenarios. We conclude that P. calcariformata established a population in the Levant because of the recent warming in the region. In line with observations from other groups of organisms, our results indicate that continued warming of the eastern Mediterranean will facilitate the invasion of more tropical marine taxa into the Mediterranean, disturbing local biodiversity and ecosystem structure.

Pore water (SO4, CH4, alkalinity, HS-, Ca, Sr, Mg, Ba) and solid phase (Ca, Sr, Mg, Ba) data measured in pockmark sediments of the Northern Congo Fan

This study focuses on the vertical distribution of authigenic carbonates (aragonite and high Mg-calcite) in the form of finely disseminated precipitates as well as massive carbonate concretions present in and above gas hydrate bearing sediments of the Northern Congo Fan. Analyses of Ca, Mg, Sr and Ba in pore water, bulk sediments and authigenic carbonates were carried out on gravity cores taken from three pockmark structures (Hydrate Hole, Black Hole and Worm Hole). In addition, a background core was retrieved from an area not influenced by fluid seepage. Pore water Sr/Ca and Mg/Ca ratios are used to reveal the current depths of carbonate formation as well as the mineralogy of the authigenic precipitates. The Sr/Ca and Mg/Ca ratios of bulk sediments and massive carbonate concretions were applied to infer the presence and depth distribution of authigenic aragonite and high Mg-calcite, based on the approach presented by Bayon et al. [Bayon et al. (2007). Sr/Ca and Mg/Ca ratios in Niger Delta sediments: Implications for authigenic carbonate genesis in cold seep environments. Marine Geology 241(1-4), 93-109, doi:10.1016/j.margeo.2007.03.007]. We show that the approach developed by Bayon et al. (2007) for sediments of cold seeps of the Niger Delta is also suitable to identify the mineralogy of authigenic carbonates in pockmark sediments of the Congo Deep-Sea Fan. We expand this approach by combining interstitial with solid phase Sr/Ca and Mg/Ca ratios, which demonstrate that high Mg-calcite is the predominant authigenic carbonate that currently forms at the sulfate/methane reaction zone (SMRZ). This is the first study which investigates both solid phase and pore water signatures typical for either aragonite or high Mg-calcite precipitation for the same sediment cores and thus is able to identify active and fossil carbonate precipitation events. At all investigated pockmark sites fossil horizons of the SMRZ were deduced from high Mg-calcite located above and below the current depths of the SMRZ. Additionally, aragonite enrichments typical for high seepage rates were detected close to the sediment surface at these sites. However, active precipitation of aragonite as indicated by pore water characteristics only occurs at the Black Hole site. Dissolved and solid phase Ba concentrations were used to estimate the time the SMRZ was fixed at the current depths of the diagenetic barite fronts. The combined pore water and solid phase elemental ratios (Mg/Ca, Sr/Ca) and Ba concentrations allow the reconstruction of past changes in methane seepage at the investigated pockmark sites. At the Hydrate Hole and Worm Hole sites the time of high methane seepage was estimated to have ceased at least 600 yr BP. In contrast, a more recent change from a high flux to a more dormant stage must have occurred at the Black Hole site as evidenced by active aragonite precipitation at the sediment surface and a lack of diagenetic Ba enrichments.

Boron contents and boron, carbon, and oxygen isotope data for ODP Sites 160-970 and 160-971

Products of two mud volcanoes from the distal part of the Mediterranean Ridge accretionary complex have been investigated regarding their B, C, and O stable isotope signatures. The mud breccias have been divided into mud matrix, lithified clasts, biogenic deposits, and authigenic cements and crusts related to fluid flow and cementation. Isotope geochemistry is used to evaluate the depth of mobilization of each phase in the subduction zone. B contents and isotope ratios of the mud and mud clasts show a general trend of B enrichment and decreasing d11B values with increasing consolidation (i.e., depth). However, the majority of the clast and matrix samples relate to moderate depths of mobilization within the wedge (1-2 km below seafloor). The carbonate cements of most of these clasts as well as the authigenic crusts, however, provide evidence for a deep fluid influence, probably associated with the décollement at 5-6 km depth. This interpretation is supported by d13C ratios of the crust, which indicate precipitation of C from thermogenic methane, and by the d11B ratios of pore-water samples of mud-breccia drill cores. Clams (Vesicomya sp.) living adjacent to fluid vents have d11B and d18O values corresponding to brines known in the area, which acted as the parent solution for shell precipitation. Such brines are most likely Miocene pore waters trapped at deep levels within the backstop to the accretionary prism, probably prior to desiccation of the Mediterranean in the Messinian (6-5 Ma). Combining all results, deep fluid circulation and expulsion are identified as the main processes triggering mud liquefaction and extrusion, whereas brines contribute only locally. Given the high B contents, mud extrusion has to be considered a major backflux mechanism of B into the hydrosphere.

Ogasawara winter Sr/Ca and U/Ca data of coral core OGA-02-1 and PDO reconstruction (coral-geoduck index)

The Pacific Decadal Oscillation (PDO), the leading mode of sea surface temperature (SST) anomalies in the extratropical North Pacific Ocean, has widespread impacts on precipitation in the Americas and marine fisheries in the North Pacific. However, marine proxy records with a temporal resolution that resolves interannual to interdecadal SST variability in the extratropical North Pacific are extremely rare. Here we demonstrate that the winter Sr/Ca and U/Ca records of an annually-banded reef coral from the Ogasawara Islands in the western subtropical North Pacific are significantly correlated with the instrumental winter PDO index over the last century. The reconstruction of the PDO is further improved by combining the coral data with an existing eastern mid-latitude North Pacific growth ring record of geoduck clams. The spatial correlations of this combined index with global climate fields suggest that SST proxy records from these locations provide potential for PDO reconstructions further back in time.

Stable isotope data and calcification temperature from planktic foraminifera in sediment core BOFS31/1K

Determining the past record of temperature and salinity of ocean surface waters is essential for understanding past changes in climate, such as those which occur across glacial-interglacial transitions. As a useful proxy, the oxygen isotope composition (delta18O) of calcite from planktonic foraminifera has been shown to reflect both surface temperature and seawater delta18O, itself an indicator of global ice volume and salinity (Shackleton, 1974; Rostek et al., 1993, doi:10.1038/364319a0). In addition, magnesium/calcium (Mg/Ca) ratios in foraminiferal calcite show a temperature dependence (Nürnberg, 1995, doi:10.2113/gsjfr.25.4.350; Nürnberg et al., 1996, doi:10.1016/0016-7037(95)00446-7; Lea et al., 1999, doi:10.1016/S0016-7037(99)00197-0) due to the partitioning of Mg during calcification. Here we demonstrate, in a field-based calibration experiment, that the variation of Mg/Ca ratios with temperature is similar for eight species of planktonic foraminifera (when accounting for Mg dissolution effects). Using a multi-species record from the Last Glacial Maximum in the North Atlantic Ocean we found that past temperatures reconstructed from Mg/Ca ratios followed the two other palaeotemperature proxies: faunal abundance (CLIMAP, 1981; Mix et al., 1999, doi:10.1029/1999PA900012) and alkenone saturation (Müller et al., 1998, doi:10.1016/S0016-7037(98)00097-0 ). Moreover, combining Mg/Ca and delta18O data from the same faunal assemblage, we show that reconstructed surface water delta18O from all foraminiferal species record the same glacial-interglacial change-representing changing hydrography and global ice volume. This reinforces the potential of this combined technique in probing past ocean-climate interactions.

E-RUN version 1.1: Observational gridded runoff estimates for Europe, link to data in NetCDF format (69 MB)

River runoff is an essential climate variable as it is directly linked to the terrestrial water balance and controls a wide range of climatological and ecological processes. Despite its scientific and societal importance, there are to date no pan-European observation-based runoff estimates available. Here we employ a recently developed methodology to estimate monthly runoff rates on regular spatial grid in Europe. For this we first assemble an unprecedented collection of river flow observations, combining information from three distinct data bases. Observed monthly runoff rates are first tested for homogeneity and then related to gridded atmospheric variables (E-OBS version 12) using machine learning. The resulting statistical model is then used to estimate monthly runoff rates (December 1950 - December 2015) on a 0.5° x 0.5° grid. The performance of the newly derived runoff estimates is assessed in terms of cross validation. The paper closes with example applications, illustrating the potential of the new runoff estimates for climatological assessments and drought monitoring.