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.

Wet tundra polygon centers and distance between polygon centroids on Samoylov Island, Lena Delta, Siberia, Russia, with links to shape files

Subgrid processes occur in various ecosystems and landscapes but, because of their small scale, they are not represented or poorly parameterized in climate models. These local heterogeneities are often important or even fundamental for energy and carbon balances. This is especially true for northern peatlands and in particular for the polygonal tundra, where methane emissions are strongly influenced by spatial soil heterogeneities. We present a stochastic model for the surface topography of polygonal tundra using Poisson-Voronoi diagrams and we compare the results with available recent field studies. We analyze seasonal dynamics of water table variations and the landscape response under different scenarios of precipitation income. We upscale methane fluxes by using a simple idealized model for methane emission. Hydraulic interconnectivities and large-scale drainage may also be investigated through percolation properties and thresholds in the Voronoi graph. The model captures the main statistical characteristics of the landscape topography, such as polygon area and surface properties as well as the water balance. This approach enables us to statistically relate large-scale properties of the system to the main small-scale processes within the single polygons.

Digital surface model, hillshade and orthophoto of the world's largest fossil oyster reef, links to GeoTIFFs

The world's largest fossil oyster reef, formed by the giant oyster Crassostrea gryphoides and located in Stetten (north of Vienna, Austria) is studied by Harzhauser et al., 2015, 2016; Djuricic et al., 2016. Digital documentation of the unique geological site is provided by terrestrial laser scanning (TLS) at the millimeter scale. Obtaining meaningful results is not merely a matter of data acquisition with a suitable device; it requires proper planning, data management, and postprocessing. Terrestrial laser scanning technology has a high potential for providing precise 3D mapping that serves as the basis for automatic object detection in different scenarios; however, it faces challenges in the presence of large amounts of data and the irregular geometry of an oyster reef. We provide a detailed description of the techniques and strategy used for data collection and processing in Djuricic et al., 2016. The use of laser scanning provided the ability to measure surface points of 46,840 (estimated) shells. They are up to 60-cm-long oyster specimens, and their surfaces are modeled with a high accuracy of 1 mm. In addition to laser scanning measurements, more than 300 photographs were captured, and an orthophoto mosaic was generated with a ground sampling distance (GSD) of 0.5 mm. This high-resolution 3D information and the photographic texture serve as the basis for ongoing and future geological and paleontological analyses. Moreover, they provide unprecedented documentation for conservation issues at a unique natural heritage site.

A global seasonal surface ocean climatology of phytoplankton types based on CHEMTAX analysis of HPLC pigments

Much advancement has been made in recent years in field data assimilation, remote sensing and ecosystem modeling, yet our global view of phytoplankton biogeography beyond chlorophyll biomass is still a cursory taxonomic picture with vast areas of the open ocean requiring field validations. High performance liquid chromatography (HPLC) pigment data combined with inverse methods offer an advantage over many other phytoplankton quantification measures by way of providing an immediate perspective of the whole phytoplankton community in a sample as a function of chlorophyll biomass. Historically, such chemotaxonomic analysis has been conducted mainly at local spatial and temporal scales in the ocean. Here, we apply a widely tested inverse approach, CHEMTAX, to a global climatology of pigment observations from HPLC. This study marks the first systematic and objective global application of CHEMTAX, yielding a seasonal climatology comprised of ~1500 1°x1° global grid points of the major phytoplankton pigment types in the ocean characterizing cyanobacteria, haptophytes, chlorophytes, cryptophytes, dinoflagellates, and diatoms, with results validated against prior regional studies where possible. Key findings from this new global view of specific phytoplankton abundances from pigments are a) the large global proportion of marine haptophytes (comprising 32 ± 5% of total chlorophyll), whose biogeochemical functional roles are relatively unknown, and b) the contrasting spatial scales of complexity in global community structure that can be explained in part by regional oceanographic conditions. These publicly accessible results will guide future parameterizations of marine ecosystem models exploring the link between phytoplankton community structure and marine biogeochemical cycles.

(Table 1) Stratigraphic position, abundance, shape, color, and refractive indices of glass shards of selected samples from DSDP Leg 60 Holes vitric muds

Precisely determined refractive indices of glass shards from 32 ash-rich, volcaniclastic sediments, mostly turbidites interbedded with nonvolcanic sediments in the Mariana Trough, range from 1.480 to 1.585 (corresponding to SiO2 ca. 75 to 49%), with most in the range 1.500 to 1.540 (SiO2 ca. 70-62%) and a second, smaller mode between ca. 1.560 and 1.585 (57 to 49% SiO2). Shards are almost exclusively colorless from 1.480 to ca. 1.530, light brown with minor colorless and green tones between 1.530 and 1.560, and dominantly brown at higher refractive indices. Tubular pumice shards are more common at higher silica percentages and non- to poorly-vesicular cuniform shards at low SiO2 values, but there is no clear correlation between shape and composition of shards. About half of the samples have bimodal shard populations with silica differences ranging up to 20 percent; unimodal layers have a range of up to about 7 percent SiO2. Of 21 samples in which one type of shard dominates, seven have the main mode in the rhyolitic composition (>69% SiO2), eight in the intermediate range (56 to 69% SiO2), and five in mafic composition (SiO2 <53%). These unusually abundant mafic shards occur mainly in site survey piston cores, SP-IA and 4E, and in Holes 454, 456, 458, and 459B. These are the sites closest to the present arc. Hole 453, containing by far the most vitric tuff turbidites, shows a gradual increase in silica content of ash layers upward to the hole from Cores 36 to 19 (about 4.6 to 3.0 Ma). A drastic decrease in ash-rich beds in the younger (Pleistocene) part of this hole was noted by the shipboard party (see site chapter, Site 453) and was interpreted by them as indicating increasing distance from the arc volcanoes as the trough opened. The increase in silica in ashes from the early to the late Pliocene at Site 453 could be interpreted in the same way and might indicate that the trough started to open in early Pliocene time.

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.

Heterocyst glycolipids in surface sediments and water samples of the Amazon shelf

Marine endosymbiotic heterocystous cyanobacteria make unique heterocyst glycolipids (HGs) containing pentose (C5) moieties. Functionally similar HGs with hexose (C6) moieties found in free-living cyanobacteria occur in the sedimentary record, but C5 HGs have not been documented in the natural environment. Here we developed a high performance liquid chromatography multiple reaction monitoring (MRM) mass spectrometry (HPLC-MS2) method specific for trace analysis of long chain C5HGs and applied it to cultures of Rhizosolenia clevei Ostenfeld and its symbiont Richelia intracellularis which were found to contain C5 HGs and no C6 HGs. The method was then applied to suspended particulate matter (SPM) and surface sediment from the Amazon plume region known to harbor marine diatoms carrying heterocystous cyanobacteria as endosymbionts. C5 HGs were detected in both marine SPM and surface sediments, but not in SPM or surface sediment from freshwater settings in the Amazon basin. Rather, the latter contained C6 HGs, established biomarkers for free-living heterocystous cyanobacteria. Our results indicate that the C5 HGs may be potential biomarkers for marine endosymbiotic heterocystous cyanobacteria.

Temperature, salinity and ikaite crystal characteristics in snow, slush and sea ice during R/V Lance cruise in April-May 2011

We identified ikaite crystals (CaCO3 · 6H2O) and examined their shape and size distribution in first-year Arctic pack ice, overlying snow and slush layers during the spring melt onset north of Svalbard. Additional measurements of total alkalinity (TA) were made for melted snow and sea-ice samples. Ikaite crystals were mainly found in the bottom of the snowpack, in slush and the surface layers of the sea ice where the temperature was generally lower and salinity higher than in the ice below. Image analysis showed that ikaite crystals were characterized by a roughly elliptical shape and a maximum caliper diameter of 201.0±115.9 µm (n = 918). Since the ice-melting season had already started, ikaite crystals may already have begun to dissolve, which might explain the lack of a relationship between ikaite crystal size and sea-ice parameters (temperature, salinity, and thickness of snow and ice). Comparisons of salinity and TA profiles for melted ice samples suggest that the precipitation/dissolution of ikaite crystals occurred at the top of the sea ice and the bottom of the snowpack during ice formation/melting processes.

(Table 1) Bulk and grain density, porosity, and p- and s-wave velocity characteristics of ODP Hole 163-990A sediments

Seismic velocities in rocks are influenced by the properties of the solid, the pore fluid, and the pore space. Cracks dramatically affect seismic velocities in rocks; their influence on the effective elastic moduli of rocks depends on their shape and concentration. Thin cracks (or fractures) substantially lower the moduli of a rock relative to the effect of spherical voids (or vesicles), and lower moduli are reflected by lower P- and S-wave velocities. The objective of this research is to determine the types and concentrations of cracks and their influence on the seismic properties of subaerially erupted basalts drilled from Hole 990A on the Southeast Greenland margin during Ocean Drilling Program Leg 163. Ellipsoidal cracks are used to model the voids in the rocks. The elastic moduli of the solid (grains) are also free parameters in the inverse modeling procedure. The apparent grain moduli reflect a weighted average of the moduli of the constituent minerals (e.g., plagioclase, augite, and clay minerals). The results indicate that (1) there is a strong relationship between P-wave velocity and porosity, suggesting a similarity of pore shape distributions, (2) the distribution of crack types within the massive, central region of aa flows from Hole 990A is independent of total porosity, (3) thin cracks are the first to be effectively sealed by alteration products, and (4) grain densities (an alteration index) and apparent grain moduli of the basalt samples are directly related.

Static GPS measurements in the vicinity of EDML and an elevation model for Dronning Maud Land

The knowledge of ice sheet surface topography and the location of the ice divides are essential for ice dynamic modeling. An improved digital elevation model (DEM) of Dronning Maud Land (DML), Antarctica, is presented in this paper. It is based on ground-based kinematic GPS profiles, airborne radar altimetry, and data of the airborne radio-echo sounding system, as well as spaceborne laser altimetry from NASA's Ice, Cloud and land Elevation Satellite (ICESat). The accuracy of ICESat ice sheet altimetry data in the area of investigation is discussed. The location of the ice divides is derived from aspect calculation of the topography and is verified with several velocity data derived from repeated static GPS measurements.

(Table 1) Sea-surface temperatures, intermediate-water temperatures, and latitudinal gradient at three time slices at the end of the Cretaceous

Climatic and oceanographic variations during the last 2 m.y. of the Maastrichtian inferred from high-resolution (10 k.y.) stable isotope analysis of the mid-latitude South Atlantic Deep Sea Drilling Project Site 525 reveal a major warm pulse followed by rapid cooling prior to the Cretaceous-Tertiary boundary. Between 66.85 and 65.52 Ma, cool but fluctuating temperatures average 9.9 and 15.4°C in intermediate and surface waters, respectively. This interval is followed by an abrupt short-term warming between 65.45 and 65.11 Ma, which increased temperatures by 2-3°C in intermediate waters, and decreased the vertical thermal gradient to an average of 2.7°C. This warm pulse may be linked to increased atmospheric pCO2, increased poleward heat transport, and the switch of an intermediate water source from high to low-middle latitudes. During the last 100 k.y. of the Maastrichtian, intermediate and surface temperatures decreased by an average of 2.1 and 1.4°C, respectively, compared to the maximum temperature between 65.32 and 65.24 Ma.

Ice thickness distribution and glacier bed of Hallstätter Gletscher

Hallstätter Glacier is the northernmost glacier of Austria. Appendant to the northern Limestone Alps, the glacier is located at 47°28'50'' N, 13°36'50'' E in the Dachstein-region. At the same time with its advance linked to the Little Ice Age (LIA), research on changes in size and mass of Hallstätter glacier was started in 1842 by Friedrich Simony. He observed and documented the glacier retreat related to its last maximum extension in 1856. In addition, Hallstätter Glacier is a subject to scientific research to date. In this thesis methods and results of ongoing mass balance measurements are presented and compared to long term volume changes and meteorological observations. The current mass balance monitoring programm using the direct glaciological method was started 2006. In this context, 2009 the ice thickness was measured with ground penetrating radar. The result are used with digital elevation models reconstucted from historical maps and recent digital elevation models to calculate changes in shape and volume of Hallstätter Glacier. Based on current meteorological measurements near the glacier and longtime homogenized climate data provided by HISTALP, time series of precipitation and temperature beginning at the LIA are produced. These monthly precipitation and monthly mean temperature data are used to compare results of a simple degree day model with the volume change calculated from the difference of the digital elevation models. The two years of direct mass balance measurements are used to calibrate the degree day model. A number of possible future scenarios are produced to indicate prospective changes. Within the 150-year-period between 1856 and 2007 the Hallstätter Glacier lost 1940 meters of its length and 2.23 km**2 in area. 37% of the initial volume of 1856 remained. This retreat came along with a change in climate. The application of a running avarage of 30 years shows an increase in precipitation of 18.5% and a warming of 1.3°C near the glacier between 1866 and 1993. The mass loss was continued in the hydrological years 2006/2007 and 2007/2008 showing mean specific mass balance of -376 mm and -700 mm, respectively. Applying a temperature correction for the different minimum elevations of the glacier, the degree day approach based on the two measured mass balances can reproduce sign and order of magnitude of the volume change of Hallstätter Glacier since 1856. Nevertheless, the relative deviation is significant. Future scenarios show, that 30% of the entire glacier volume remains after subtracting the elevation changes between the digital elevation models of 2002 and 2007 ten times from the surface of 2007. The past and present mass changes of Hallstätter Glacier are showing a retreating glacier as a consequence of rising temperatures. Due to high precepitation, increased with previous warming, the Hallstätter Glacier can and will exist in lower elevation compared to inner alpine glaciers.

Siliceous microplankton in surface sediments and the water column of the Red Sea

We studied the siliceous microplankton assemblages (mainly diatoms) from plankton tows (mesh size 20 µm) and surface sediment samples collected along a N-S transect in the northern Red Sea (28-21°N). In addition, we analyzed differences/similarities between plankton and sediment assemblages within a brine-filled basin (the southern basin) of the Shaban Deep and compared these assemblages with those from outside the brine. Plankton samples revealed the overwhelming dominance of diatoms over other siliceous groups. Diatoms accounted for ca. 97% of all biosiliceous particles at 120-20 m (vs. 2.9% silicoflagellates and 0.4% radiolarians), and ca. 94% at 200-120 m (vs. 4.5% silicoflagellates and 1.6% radiolarians). In general, a marine, warm-water (tropical/subtropical) diatom assemblage characterizes the plankton samples. Representatives of the Nitzschia bicapitata group are by far the most abundant contributors at both depth intervals (average=43%), ranging from ca. 30% in the North to ca. 60% in the South. Biogenic opal content in non-brine surface sediments is very low, (below 0.2 wt.% SiO2); and concentration of siliceous microorganisms is also low and of the order of 5*10**3-10**4 microorganisms/g dry sediment. Diatoms are the main contributors to the opal signal in the 20-40 µm fraction, while they share dominance with radiolarians in the >40 µm fraction. Total diatom concentrations average 1.2*10**4 valves/g in the 20-40 µm fraction and 4*10**3 valves/g in the >40 µm fraction. Robust taxa of warm water affinity (Alveus marinus, Azpeitia neocrenulata, Azpeitia nodulifera and Roperia tesselata) characterize the surface sediments. In contrast, biogenic opal content in brine surface sediment samples is much higher than in the non-brine samples, ranging from 2.8 to 3.8 wt.% SiO2, and concentration of siliceous microorganisms is 3-4 orders of magnitude higher. In addition here, diatoms dominate the opal signal. The taxa found in these samples are a mixture of non-brine and plankton samples, and fragile forms (e.g., N. bicapitata group, Neodelphineis indica) are well preserved in these sediments. Thus, brine sediments in this region seem to offer a great potential for palaeoenvironmental studies.