Population genetic and dispersal modeling data for Bathymodiolus mussels from the Mid-Atlantic Ridge

The zip folder comprises a text file and a gzipped tar archive. 1) The text file contains individual genotype data for 90 SNPs, 9 microsatellites and the mitochondrial ND4 gene that were determined in deep-sea hydrothermal vent mussels from the Mid-Atlantic Ridge (genus Bathymodiolus). Mussel specimens are grouped according to the population (pop)/location from which they have been sampled (first column). The remaining columns contain the respective allele/haplotype codes for the different genetic loci (names in the header line). The data file is in CONVERT format and can be directly transformed into different input files for population genetic statistics. 2) The tar archive contains NetCDF files with larval dispersal probabilities for simulated annual larval releases between 1998 and 2007. For each simulated vent location (Menez Gwen, Lucky Strike, Rainbow, Vent 1-10) two NetCDF files are given, one for an assumed pelagic larval duration of 1 year and the other one for an assumed pelagic larval duration of 6 months (6m).

Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.

Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo 1.0.5-beta)

Sub-ice shelf circulation and freezing/melting rates in ocean general circulation models depend critically on an accurate and consistent representation of cavity geometry. Existing global or pan-Antarctic data sets have turned out to contain various inconsistencies and inaccuracies. The goal of this work is to compile independent regional fields into a global data set. We use the S-2004 global 1-minute bathymetry as the backbone and add an improved version of the BEDMAP topography for an area that roughly coincides with the Antarctic continental shelf. Locations of the merging line have been carefully adjusted in order to get the best out of each data set. High-resolution gridded data for upper and lower ice surface topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier have been carefully merged into the ambient ice and ocean topographies. Multibeam survey data for bathymetry in the former Larsen B cavity and the southeastern Bellingshausen Sea have been obtained from the data centers of Alfred Wegener Institute (AWI), British Antarctic Survey (BAS) and Lamont-Doherty Earth Observatory (LDEO), gridded, and again carefully merged into the existing bathymetry map. The global 1-minute dataset (RTopo-1 Version 1.0.5) has been split into two netCDF files. The first contains digital maps for global bedrock topography, ice bottom topography, and surface elevation. The second contains the auxiliary maps for data sources and the surface type mask. A regional subset that covers all variables for the region south of 50 deg S is also available in netCDF format. Datasets for the locations of grounding and coast lines are provided in ASCII format.

A suite of global, cross-scale topographic variables for environmental and biodiversity modeling, links to files in GeoTIFF format

Topographic variation, the spatial variation in elevation and terrain features, underpins a myriad of patterns and processes in geography and ecology and is key to understanding the variation of life on the planet. The characterization of this variation is scale-dependent, i.e. it varies with the distance over which features are assessed and with the spatial grain (grid cell resolution) of analysis. A fully standardized and global multivariate product of different terrain features has the potential to support many large-scale basic research and analytical applications, however to date, such technique is unavailable. Here we used the digital elevation model products of global 250 m GMTED and near-global 90 m SRTM to derive a suite of topographic variables: elevation, slope, aspect, eastness, northness, roughness, terrain roughness index, topographic position index, vector ruggedness measure, profile and tangential curvature, and 10 geomorphological landform classes. We aggregated each variable to 1, 5, 10, 50 and 100 km spatial grains using several aggregation approaches (median, average, minimum, maximum, standard deviation, percent cover, count, majority, Shannon Index, entropy, uniformity). While a global cross-correlation underlines the high similarity of many variables, a more detailed view in four mountain regions reveals local differences, as well as scale variations in the aggregated variables at different spatial grains. All newly-developed variables are available for download at http://www.earthenv.org and can serve as a basis for standardized hydrological, environmental and biodiversity modeling at a global extent.

Flow conditions, bed form dimensions and shear stress experiments

Large asymmetric bed forms commonly develop in rivers. The turbulence associated with flow separation that develops over their steep lee side is responsible for the form shear stress which can represent a substantial part of total shear stress in rivers. This paper uses the Delft3D modeling system to investigate the effects of bed form geometry and forcing conditions on flow separation length and associated turbulence, and bed form shear stress over angle-of-repose (30 lee side angle) bed forms. The model was validated with lab measurements that showed sufficient agreement to be used for a systematic analysis. The influence of flow velocity, bed roughness, relative height (bed form height/water depth), and aspect ratio (bed form height/length) on the variations of the normalized length of the flow separation zone, the extent of the wake region (where the turbulent kinetic energy (TKE) was more than 70% of the maximum TKE), the average TKE within the wake region and the form shear stress were investigated. Form shear stress was found not to scale with the size of the flow separation zone but to be related to the product of the normalized extent of the wake region (extent of the wake region/extent of water body above the bed form) and the average TKE within the wake region. The results add to understanding of the hydrodynamics of bed forms and may be used for the development of better parameterizations of smallscale processes for application in large-scale studies.

40Ar-39Ar plateau and isochron ages from the Ninetyeast and Broken Ridges (Table 1)

Concordant plateau and isochron ages are obtained from 40Ar-39Ar incremental heating experiments on volcanic rocks recovered by drilling at three Leg 121 sites along the Ninetyeast Ridge and two dredge locations on the southern scarp of the Broken Ridge, eastern Indian Ocean. The new data confirm a northerly increase in the age of volcanism along the Ninetyeast Ridge, from 38 to 82 Ma; this lineament links current hotspot volcanism near the Kerguelen islands with the Rajmahal flood basalt eruptions at M0 time (117 ± 1 Ma). The Broken Ridge was formed over the same hotspot at 88-89 Ma, but later experienced rift-related volcanism in Paleocene time (63 Ma). The geometry and distribution of ages along these prominent volcanic ridges and the Mascarene-Chagos-Laccadive-Maldive ridge system in the western Indian Ocean are most compatible with plate motions over fixed hotspots near Kerguelen and Reunion islands, respectively.

Age determination and stable isotope record of foraminifera of sediment core SO82_5-2

Surface and deepwater paleoclimate records in Irminger Sea core SO82-5 (59°N, 31°W) and Icelandic Sea core PS2644 (68°N, 22°W) exhibit large fluctuations in thermohaline circulation (THC) from 60 to 18 calendar kyr B.P., with a dominant periodicity of 1460 years from 46 to 22 calendar kyr B.P., matching the Dansgaard-Oeschger (D-O) cycles in the Greenland Ice Sheet Project 2 (GISP2) temperature record [Grootes and Stuiver, 1997, doi:10.1029/97JC00880]. During interstadials, summer sea surface temperatures (SSTsu) in the Irminger Sea averaged to 8°C, and sea surface salinities (SSS) averaged to ~36.5, recording a strong Irminger Current and Atlantic THC. During stadials, SSTsu dropped to 2°-4°C, in phase with SSS drops by ~1-2. They reveal major meltwater injections along with the East Greenland Current, which turned off the North Atlantic deepwater convection and hence the heat advection to the north, in harmony with various ocean circulation and ice models. On the basis of the IRD composition, icebergs came from Iceland, east Greenland, and perhaps Svalbard and other northern ice sheets. However, the southward drifting icebergs were initially jammed in the Denmark Strait, reaching the Irminger Sea only with a lag of 155-195 years. We also conclude that the abrupt stadial terminations, the D-O warming events, were tied to iceberg melt via abundant seasonal sea ice and brine water formation in the meltwater-covered northwestern North Atlantic. In the 1/1460-year frequency band, benthic ?18O brine water spikes led the temperature maxima above Greenland and in the Irminger Sea by as little as 95 years. Thus abundant brine formation, which was induced by seasonal freezing of large parts of the northwestern Atlantic, may have finally entrained a current of warm surface water from the subtropics and thereby triggered the sudden reactivation of the THC. In summary, the internal dynamics of the east Greenland ice sheet may have formed the ultimate pacemaker of D-O cycles.

Benthic and planktic foraminiferal stable isotope records combined with Nd isotope ratios of ODP Hole 198-1210B

The Pacific Ocean is the largest water body on Earth, and circulation in the Pacific contributed significantly to climate evolution in the latest Cretaceous, the culmination of a period of long-term cooling. Here, we present new high-resolution late Campanian to Maastrichtian benthic and planktic foraminiferal stable isotope data and a neodymium (Nd) isotope record obtained from sedimentary ferromanganese oxide coatings of Ocean Drilling Program Hole 1210B from the tropical Pacific Ocean (Shatsky Rise). These new records resolve 13 million years in the latest Cretaceous, providing insights into changes in surface and bottom water temperatures and source regions of deep to intermediate waters covering the carbon isotope excursions of the Campanian-Maastrichtian Boundary Event (CMBE) and the Mid-Maastrichtian event (MME). Our new benthic foraminiferal d18O and Nd isotope records together with published Nd isotope data show markedly parallel trends across the studied interval over a broad range of bathyal to abyssal water depths interpreted to reflect changes in the intensity of deep-ocean circulation in the tropical Pacific. In particular, we observe a three-million-year-long period of cooler conditions in the early Maastrichtian (72.5 to 69.5 Ma) when a concomitant change toward less radiogenic seawater Nd isotope signatures probably marks a period of enhanced admixture and northward flow of deep waters with Southern Ocean provenance. We suggest this change to have been triggered by intensified formation and convection of deep waters in the high southern latitudes, a process that weakened during the MME (69.5 to 68.5 Ma). The early Maastrichtian cold interval is closely related to the negative and positive carbon isotope trends of the CMBE and MME. The millions-of-years long duration of these carbon cycle perturbations suggests a tectonic forcing of climatic cooling, possibly related to changes in ocean basin geometry and bathymetry.

Geochemistry of dikes and lavas from the north wall of the Hess Deep Rift

An investigation of ~1-m.y.-old dikes and lavas from the north wall of the Hess Deep Rift (2°15'N, 101°30'W) collected during Alvin expeditions provides a detailed view of the evolution of fast spreading oceanic crust. The study area encompasses 25 km of an east-west flow line, representing ~370,000 years of crustal accretion at the East Pacific Rise. Samples analyzed exhibit depleted incompatible trace element abundances and ratios [(La/Sm)N < 1]. Indices of fractionation (MgO), and incompatible element ratios (La/Sm, Nb/Ti) show no systematic trends along flow line. Rather, over short (<4 m) and long (~25 km) distances, significant variations are observed in major and trace element concentrations and ratios. Modeling of these variations attests to the juxtaposition of dikes of distinct parental magma compositions. These findings, combined with studies of segmentation of the subaxial magma chamber and lateral magma transport in dikes along rift-dominated systems, suggest a more realistic model of the magmatic system underlying the East Pacific Rise relative to the commonly assumed twodimensional model. In this model, melts from a heterogeneous mantle feed distinct portions of a segmented axial magma reservoir. Dikes emanating from these distinct reservoirs transport magma along axis, resulting in interleaved dikes and host lavas with different evolutionary histories. This model suggests the use of axial or flow line lava compositions to infer the evolution of axial magma chambers should be approached with caution because dikes may never erupt lava or may transport magma significant distances along axis and erupt lavas far from their axial magma chamber of origin.