Global Lithological Map Database v1.0 (gridded to 0.5° spatial resolution)

Lithology describes the geochemical, mineralogical, and physical properties of rocks. It plays a key role in many processes at the Earth surface, especially the fluxes of matter to soils, ecosystems, rivers, and oceans. Understanding these processes at the global scale requires a high resolution description of lithology. A new high resolution global lithological map (GLiM) was assembled from existing regional geological maps translated into lithological information with the help of regional literature. The GLiM represents the rock types of the Earth surface using 1,235,400 polygons. The lithological classification consists of three levels. The first level contains 16 lithological classes comparable to previously applied definitions in global lithological maps. The additional two levels contain 12 and 14 subclasses, respectively, which describe more specific rock attributes. According to the GLiM, the Earth is covered by 64 % sediments (a third of which is carbonates), 13 % metamorphics, 7 % plutonics, and 6 % volcanics, and 10% are covered by water or ice. The high resolution of the GLiM allows observation of regional lithological distributions which often vary from the global average. The GLiM enables regional analysis of Earth surface processes at global scales.

Benthic cover map of Heron Reef derived from a high-spatial-resolution multi-spectral satellite image using object based image analysis

Providing accurate maps of coral reefs where the spatial scale and labels of the mapped features correspond to map units appropriate for examining biological and geomorphic structures and processes is a major challenge for remote sensing. The objective of this work is to assess the accuracy and relevance of the process used to derive geomorphic zone and benthic community zone maps for three western Pacific coral reefs produced from multi-scale, object-based image analysis (OBIA) of high-spatial-resolution multi-spectral images, guided by field survey data. Three Quickbird-2 multi-spectral data sets from reefs in Australia, Palau and Fiji and georeferenced field photographs were used in a multi-scale segmentation and object-based image classification to map geomorphic zones and benthic community zones. A per-pixel approach was also tested for mapping benthic community zones. Validation of the maps and comparison to past approaches indicated the multi-scale OBIA process enabled field data, operator field experience and a conceptual hierarchical model of the coral reef environment to be linked to provide output maps at geomorphic zone and benthic community scales on coral reefs. The OBIA mapping accuracies were comparable with previously published work using other methods; however, the classes mapped were matched to a predetermined set of features on the reef.

A 20 m spatial resolution seamless multisource Digital Elevation/Depth Model for Lizard Island, northern Great Barrier Reef

Hydrographers have traditionally referred to the nearshore area as the "white ribbon" area due to the challenges associated with the collection of elevation data in this highly dynamic transitional zone between terrestrial and marine environments. Accordingly, available information in this zone is typically characterised by a range of datasets from disparate sources. In this paper we propose a framework to 'fill' the white ribbon area of a coral reef system by integrating multiple elevation and bathymetric datasets acquired by a suite of remote-sensing technologies into a seamless digital elevation model (DEM). A range of datasets are integrated, including field-collected GPS elevation points, terrestrial and bathymetric LiDAR, single and multibeam bathymetry, nautical chart depths and empirically derived bathymetry estimations from optical remote sensing imagery. The proposed framework ranks data reliability internally, thereby avoiding the requirements to quantify absolute error and results in a high resolution, seamless product. Nested within this approach is an effective spatially explicit technique for improving the accuracy of bathymetry estimates derived empirically from optical satellite imagery through modelling the spatial structure of residuals. The approach was applied to data collected on and around Lizard Island in northern Australia. Collectively, the framework holds promise for filling the white ribbon zone in coastal areas characterised by similar data availability scenarios. The seamless DEM is referenced to the horizontal coordinate system MGA Zone 55 - GDA 1994, mean sea level (MSL) vertical datum and has a spatial resolution of 20 m.

Coral reef habitat maps derived from a high-spatial-resolution multi-spectral satellite image using object based image analysis

Coral reef maps at various spatial scales and extents are needed for mapping, monitoring, modelling, and management of these environments. High spatial resolution satellite imagery, pixel <10 m, integrated with field survey data and processed with various mapping approaches, can provide these maps. These approaches have been accurately applied to single reefs (10-100 km**2), covering one high spatial resolution scene from which a single thematic layer (e.g. benthic community) is mapped. This article demonstrates how a hierarchical mapping approach can be applied to coral reefs from individual reef to reef-system scales (10-1000 km**2) using object-based image classification of high spatial resolution images guided by ecological and geomorphological principles. The approach is demonstrated for three individual reefs (10-35 km**2) in Australia, Fiji, and Palau; and for three complex reef systems (300-600 km**2) one in the Solomon Islands and two in Fiji. Archived high spatial resolution images were pre-processed and mosaics were created for the reef systems. Georeferenced benthic photo transect surveys were used to acquire cover information. Field and image data were integrated using an object-based image analysis approach that resulted in a hierarchically structured classification. Objects were assigned class labels based on the dominant benthic cover type, or location-relevant ecological and geomorphological principles, or a combination thereof. This generated a hierarchical sequence of reef maps with an increasing complexity in benthic thematic information that included: 'reef', 'reef type', 'geomorphic zone', and 'benthic community'. The overall accuracy of the 'geomorphic zone' classification for each of the six study sites was 76-82% using 6-10 mapping categories. For 'benthic community' classification, the overall accuracy was 52-75% with individual reefs having 14-17 categories and reef systems 20-30 categories. We show that an object-based classification of high spatial resolution imagery, guided by field data and ecological and geomorphological principles, can produce consistent, accurate benthic maps at four hierarchical spatial scales for coral reefs of various sizes and complexities.

Habitat map of Danajon Bank, Philippines, derived from a high-spatial-resolution multi-spectral satellite image and georeferenced point intercept transect and spot-check survey field data, using an object-based image classification method

A mosaic of two WorldView-2 high resolution multispectral images (Acquisition dates: October 2010 and April 2012), in conjunction with field survey data, was used to create a habitat map of the Danajon Bank, Philippines (10°15'0'' N, 124°08'0'' E) using an object-based approach. To create the habitat map, we conducted benthic cover (seafloor) field surveys using two methods. Firstly, we undertook georeferenced point intercept transects (English et al., 1997). For ten sites we recorded habitat cover types at 1 m intervals on 10 m long transects (n= 2,070 points). Second, we conducted geo-referenced spot check surveys, by placing a viewing bucket in the water to estimate the percent cover benthic cover types (n = 2,357 points). Survey locations were chosen to cover a diverse and representative subset of habitats found in the Danajon Bank. The combination of methods was a compromise between the higher accuracy of point intercept transects and the larger sample area achievable through spot check surveys (Roelfsema and Phinn, 2008, doi:10.1117/12.804806). Object-based image analysis, using the field data as calibration data, was used to classify the image mosaic at each of the reef, geomorphic and benthic community levels. The benthic community level segregated the image into a total of 17 pure and mixed benthic classes.

High-resolution ice thickness and bed topography of a land-terminating section of the Greenland Ice Sheet, with links to GeoTIFFs and text files

We present ice thickness and bed topography maps with a high spatial resolution (250-500 m) of a land-terminating section of the Greenland Ice Sheet derived from ground-based and airborne radar surveys. The data have a total area of ~12 000 km^2 and cover the whole ablation area of the outlet glaciers of Isunnguata Sermia, Russell, Leverett, Ørkendalen and Isorlersuup up to the long-term mass balance equilibrium line altitude at ~1600 m above sea level. The bed topography shows highly variable subglacial trough systems, and the trough of Isunnguata Sermia Glacier is overdeepened and reaches an elevation of ~500 m below sea level. The ice surface is smooth and only reflects the bedrock topography in a subtle way, resulting in a highly variable ice thickness. The southern part of our study area consists of higher bed elevations compared to the northern part. The compiled data sets of ground-based and airborne radar surveys cover one of the most studied regions of the Greenland Ice Sheet and can be valuable for detailed studies of ice sheet dynamics and hydrology.

Short-term temporal and spatial variations in physical, optical and biogeochemical parameters in June of 2013 in Alfacs Bay (NW Mediterranean)

With the aim of analyzing the complex physical and biogeochemical interactions at high temporal and spatial resolution in the complex estuarine waters of Alfacs Bay, a beam attenuation-based approach was used as optical proxy of different biogeochemical variables. Thus, the dataset contains the attenuation proxies as well as laboratory results from the analysis of water samples, which were used to validate our approach. In addition, the major physical forcing in the Bay was also measured.

Retreat of top cliff of Kurungnakh Island, Lena Delta, Siberia, Russia, 2010-2014, with links to shapefiles

Thawing-induced cliff top retreat in permafrost landscapes is mainly due to thermo-erosion. Ground-ice-rich permafrost landscapes are specifically vulnerable to thermo-erosion and may show high degradation rates. Within the HGF Alliance Remote Sensing and the FP7 PAGE21 permafrost programs we investigated how SAR and optical remote sensing can contribute to the monitoring of erosion rates of ice-rich cliffs in Arctic Siberia (Lena Delta, Russia). We produced two different vector products: i) Intra-annual cliff top retreat based on TerraSAR-X (TSX) satellite data (2012-2014): High-temporal resolution time series of TSX satellite data allow the inter-annual and intra-annual monitoring of the upper cliff-line retreat also under bad weather conditions and continuous cloud coverage. This published SAR product contains the retreating upper cliff lines of a 1.5 km long part of eroding ice-rich coast of Kurungnakh Island in the central Lena Delta. The upper cliff line was mapped using a thresholding approach for images acquired in the years 2012, 2013 and 2014 for the months June (2013, 2014), July (2013, 2014), August (2012, 2013, 2014) and September (2013, 2014). The cliff top retreat vector product is called 'upper_cliff_TerraSAR-X'. While the 2014 cliff lines show a clear retreat of 2 to 3 m/month, the cliff top lines for 2012 and 2013 are not chronologically ordered. However, lines from the end of the season of a year are always close to the lines from the beginning of the next summer season, indicating low cliff retreat in winter. ii) 4-year cliff top retreat based on optical satellite data (2010-2014): Long-term cliff top retreat could be assessed with two high-spatial resolution optical satellite images (GeoEye-1, 2010-08-05 and Worldview-1, 2014-08-19). The cliff top retreat vector product is called 'upper_cliff_optical'. Results: The long-term cliff top retreat derived from optical satellite data are 35 m cliff retreat within 4 years. The higher-temporal resolution SAR data equivalently show long-term rates of 18 m within 2 years and nearly now degradation activities in winter but maximum erosion rates in summer months.The Intra-seasonal cliff top retreat lines from 2014 show a rate of 2 to 3 m per month.

In-situ Iron isotope measurements and mineral compositions from ODP Hole 206-1256D and IODP Hole 335-U1256D

In-situ Fe isotope measurements have been carried out to estimate the impact of the hydrothermal metamorphic overprint on the Fe isotopic composition of Fe-Ti-oxides and Fe-sulfides of the different lithologies of the drilled rocks from IODP Hole 1256D (eastern equatorial Pacific; 15 Ma crust formed at the East Pacific Rise). Most igneous rocks normally have a very restricted range in their 56Fe/54Fe ratio. In contrast, Fe isotope compositions of hot fluids (> 300 °C) from mid-ocean-ridge spreading centers define a narrow range that is shifted to lower delta 56Fe values by 0.2 per mil - 0.5 per mil as compared to igneous rocks. Therefore, it is expected that mineral phases that contain large amounts of Fe are especially affected by the interaction with a fluid that fractionates Fe isotopes during exsolution/precipitation of those minerals. We have used a femtosecond UV-Laser ablation system to determine mineral 56Fe/54Fe ratios of selected samples with a precision of < 0.1 per mil (2 sigma level) at micrometer-scale. We have found significant variations of the delta 56Fe (IRMM-014) values in the minerals between different samples as well as within samples and mineral grains. The overall observed scale of delta 56Fe (magnetite) in 1256D rocks ranges from - 0.12 to + 0.64 per mil, and of delta 56Fe (ilmenite) from - 0.77 to + 0.01 per mil. Pyrite in the lowermost sheeted dike section is clearly distinguishable from the other investigated lithological units, having positive delta 56Fe values between + 0.29 and + 0.56 per mil, whereas pyrite in the other samples has generally negative delta 56Fe values from - 1.10 to - 0.59 permil. One key observation is that the temperature dependent inter-mineral fractionations of Fe isotopes between magnetite and ilmenite are systematically shifted towards higher values when compared to theoretically expected values, while synthesized, well equilibrated magnetite-ilmenite pairs are compatible with the theoretical predictions. Theoretical considerations including beta-factors of different aqueous Fe-chlorides and Rayleigh-type fractionations in the presence of a hydrous, chlorine-bearing fluid can explain this observation. The disagreement between observed and theoretical equilibrium fractionation, the fact that magnetite, in contrast to ilmenite shows a slight downhole trend in the delta 56Fe values, and the observation of small scale heterogeneities within single mineral grains imply that a general re-equilibration of the magnetite-ilmenite pairs is overprinted by kinetic fractionation effects, caused by the interaction of magnetite/ilmenite with hydrothermal fluids penetrating the upper oceanic crust during cooling, or incomplete re-equilibration at low temperatures. Furthermore, the observation of significant small-scale variations in the 56Fe/54Fe ratios of single minerals in this study highlights the importance of high spatial-resolution-analyses of stable isotope ratios for further investigations.

Robust trends of landscape dynamics in the Arctic Lena Delta with temporally dense Landsat time-series stacks, with links to GeoTIFFs

Arctic permafrost landscapes are among the most vulnerable and dynamic landscapes globally, but due to their extent and remoteness most of the landscape changes remain unnoticed. In order to detect disturbances in these areas we developed an automated processing chain for the calculation and analysis of robust trends of key land surface indicators based on the full record of available Landsat TM, ETM +, and OLI data. The methodology was applied to the ~ 29,000 km**2 Lena Delta in Northeast Siberia, where robust trend parameters (slope, confidence intervals of the slope, and intercept) were calculated for Tasseled Cap Greenness, Wetness and Brightness, NDVI, and NDWI, and NDMI based on 204 Landsat scenes for the observation period between 1999 and 2014. The resulting datasets revealed regional greening trends within the Lena Delta with several localized hot-spots of change, particularly in the vicinity of the main river channels. With a 30-m spatial resolution various permafrost-thaw related processes and disturbances, such as thermokarst lake expansion and drainage, fluvial erosion, and coastal changes were detected within the Lena Delta region, many of which have not been noticed or described before. Such hotspots of permafrost change exhibit significantly different trend parameters compared to non-disturbed areas. The processed dataset, which is made freely available through the data archive PANGAEA, will be a useful resource for further process specific analysis by researchers and land managers. With the high level of automation and the use of the freely available Landsat archive data, the workflow is scalable and transferrable to other regions, which should enable the comparison of land surface changes in different permafrost affected regions and help to understand and quantify permafrost landscape dynamics.