Secchi depth maps for the Baltic Sea, 2003-2012, with links to GeoTIFs

Secchi depth is a measure of water transparency. In the Baltic Sea region, Secchi depth maps are used to assess eutrophication and as input for habitat models. Due to their spatial and temporal coverage, satellite data would be the most suitable data source for such maps. But the Baltic Sea's optical properties are so different from the open ocean that globally calibrated standard models suffer from large errors. Regional predictive models that take the Baltic Sea's special optical properties into account are thus needed. This paper tests how accurately generalized linear models (GLMs) and generalized additive models (GAMs) with MODIS/Aqua and auxiliary data as inputs can predict Secchi depth at a regional scale. It uses cross-validation to test the prediction accuracy of hundreds of GAMs and GLMs with up to 5 input variables. A GAM with 3 input variables (chlorophyll a, remote sensing reflectance at 678 nm, and long-term mean salinity) made the most accurate predictions. Tested against field observations not used for model selection and calibration, the best model's mean absolute error (MAE) for daily predictions was 1.07 m (22%), more than 50% lower than for other publicly available Baltic Sea Secchi depth maps. The MAE for predicting monthly averages was 0.86 m (15%). Thus, the proposed model selection process was able to find a regional model with good prediction accuracy. It could be useful to find predictive models for environmental variables other than Secchi depth, using data from other satellite sensors, and for other regions where non-standard remote sensing models are needed for prediction and mapping. Annual and monthly mean Secchi depth maps for 2003-2012 come with this paper as Supplementary materials.

Two historical maps from nineteenth-century Palestine, with links to digitized maps in shapefile format

Reconstructing past landscapes from historical maps requires quantifying the accuracy and completeness of these sources. The accuracy and completeness of two historical maps of the same period covering the same area in Israel were examined: the 1:63,360 British Palestine Exploration Fund map (1871-1877) and the 1:100,000 French Levés en Galilée (LG) map (1870). These maps cover the mountainous area of the Galilee (northern Israel), a region with significant natural and topographical diversity, and a long history of human presence. Land-cover features from both maps, as well as the contours drawn on the LG map, were digitized. The overall correspondence between land-cover features shown on both maps was 59% and we found that the geo-referencing method employed (transformation type and source of control points) did not significantly affect these correspondence measures. Both maps show that in the 1870s, 35% of the Galilee was covered by Mediterranean maquis, with less than 8% of the area used for permanent agricultural cropland (e.g., plantations). This article presents how the reliability of the maps was assessed by using two spatial historical sources, and how land-cover classes that were mapped with lower certainty and completeness are identified. Some of the causes that led to observed differences between the maps, including mapping scale, time of year, and the interests of the surveyors, are also identified.

Maps and acoustic profiles of free gas distribution, Western Baltic Sea and Kattgat-Skagerrak area

A shallow gas depth-contour map covering the Skagerrak-western Baltic Sea region has been constructed using a relatively dense grid of existing shallow seismic lines. The digital map is stored as an ESRI shape file in order to facilitate comparison with other data from the region. Free gas usually occurs in mud and sandy mud but is observed only when sediment thickness exceeds a certain threshold value, depending on the water depth of the area in question. Gassy sediments exist at all water depths from approx. 20 m in the coastal waters of the Kattegat to 360 m in the Skagerrak. In spite of the large difference in water depths, the depth of free gas below seabed varies only little within the region, indicating a relatively fast movement of methane in the gas phase towards the seabed compared to the rate of diffusion of dissolved methane. Seeps of old microbial methane occur in the northern Kattegat where a relatively thin cover of sandy sediments exists over shallow, glacially deformed Pleistocene marine sediments. Previous estimates of total methane escape from the area may be correct but the extrapolation of local methane seepage rate data to much larger areas on the continental shelf is probably not justified. Preliminary data on porewater chemistry were compared with the free gas depth contours in the Aarhus Bay area, which occasionally suffers from oxygen deficiency, in order to examine if acoustic gas mapping may be used for monitoring the condition of the bay.

Magnetic anomaly map of the Weddell Sea Region, Antarctic (Scale 1:2 500 000)

This paper describes a 1 : 2 500 000 scale aeromagnetic anomaly map produced by the joint efforts of VNIIOkeangeologia, Polar Marine Geological Research Expedition (PMGRE) and the Alfred Wegener Institute for Polar and Marine Research (AWl) for the Weddell Sea region covering 1 850 000 km' of West Antarctica. Extensive regional magnetic survey flights with line-spacing of about 20 km and 5 km were carried out by the PMGRE between 1977 and 1989. In course of these investigations the PMGRE flew 9 surveys with flight-line spacing of 20 km and 6 surveys with flight-line spacing of 5 km mainly over the mountain areas of southern Palmer Land, western Dronning Maud Land, Coats Land and Pensacola Mountains, over the Ronne lee Shelf and the Filchner Ice Shelf and the central part of the Weddell Sea. More than 215 000 line-kilometers of total field aeromagnetic data have been acquired by using an Ilyushin Il-14 ski-equipped aircraft. Survey operations were centered on the field base stations Druzhnaya-1, -2, and -3, from which the majority of the Weddell Sea region network was completed. The composite map of the Weddell Sea region is prepared in colour, showing magnetic anomaly contours at intervals of 50-100 nT with supplemental contours at an interval of 25 nT in low gradient areas, on a polar stereographic projection. The compiled colour magnetic anomaly map of the Weddell Sea region demonstrates that features of large areal extent, such as geologic provinces, fold-belts, ancient eratonic fragments and other regional structural features can be readily delineated. The map allows a comparison of regional magnetic features with similar-scale geological structures on geological and geophysical maps. It also provides a database for the future production of the ''Digital Magnetic Anomaly Map of Antarctica'' in the framework of the Scientific Committee on Antarctic Research/International Association of Geomagnetism and Aeronomy (SCAR/IAGA) compilation.

Gravity mapping in the southern Weddell Sea region, Antarctica (Scale 1:2 000 000)

New maps of free-air and the Bouguer gravity anomalies on the Weddell Sea sector (70-81° S, 6-75° W) of Antarctica are presented. These maps are based on the first computer compilation of available gravity data collected by ''Sevmorgeologia'' in 1976-89 in the southern Weddell Sea and adjacent coasts of western Dronning Maud Land (WDML) and Coats Land. The accomplished gravity studies comprise airborne observations with a line spacing of about 20 km and conventional measurements at over-the-ice points, which were spaced at 10-30 km and supplemented by seismic soundings. Hence, anomalies on the maps represent mainly large-scale and deep crustal features. The dominant feature in free-air gravity map is a large dipolar gravity anomaly stretching along the continental margin. Following the major grain of seabed morphology this shelf-edge/slope anomaly (SESA) is clearly divided into three segments characterized by diverse anomaly amplitudes, wavelengths and trends. They are associated with continental margins of different geotectonic provinces of Antarctica surrounding the Weddell Sea. Apparent distinctions in the SESA signatures are interpreted as the gravity expression of tectonic, deep crustal structure segmentation of the continental margin. The prominent gravity highs (100-140 mGal) of the shelf edge anomaly mapped along WDML are assumed to represent high-density mantle injections intruded into the middle/lower crust during initial rifting of continental breakup. Enlarged wavelengths and diminished amplitudes of the gravity anomaly westwards, along the Weddell Sea embayment (WSE) margin, reflect a widening of the continental slope and a significant increase in thickness of underlying sediment strata. Low amplitude, negative free-air anomalies in the Filchner-Ronne Ice Shelves (FRIS) contrast sharply with the dominating positive anomalies offshore. This indicates a greater sedimentary thickness of the basin in this area. Crustal response to the enlarged sediment load is impressed in mostly positive features of the Bouguer gravity field observed here. Two pronounced positive Bouguer anomalies of 50-70 mGal and an average widths of 200 km dominate the Weddell Sea embayment margins towards the Antarctic Peninsula and the East Antarctic craton. They correlate well with very deep seabed troughs (> 1000 m below sea level). The gravity highs are most likely caused by a shallow upper mantle underneath graben-rift structures evolved at the margins of the WSE basin. A regional zone (> 100 km in width) of the prominent Bouguer and free-air negative anomalies (-40 to -60 mGal) adjacent Coats Land to the north of the ice shelf edge may indicate the presence of the thick old cratonic crust far offshore beneath the Weddell Sea Embayment.

Daily pan-Arctic sea-ice lead maps for 2003-2015, with links to maps in netCDF format

The presence of sea-ice leads represents a key feature of the Arctic sea ice cover. Leads promote the flux of sensible and latent heat from the ocean to the cold winter atmosphere and are thereby crucial for air-sea-ice-ocean interactions. We here apply a binary segmentation procedure to identify leads from MODIS thermal infrared imagery on a daily time scale. The method separates identified leads into two uncertainty categories, with the high uncertainty being attributed to artifacts that arise from warm signatures of unrecognized clouds. Based on the obtained lead detections, we compute quasi-daily pan-Arctic lead maps for the months of January to April, 2003-2015. Our results highlight the marginal ice zone in the Fram Strait and Barents Sea as the primary region for lead activity. The spatial distribution of the average pan-Arctic lead frequencies reveals, moreover, distinct patterns of predominant fracture zones in the Beaufort Sea and along the shelf-breaks, mainly in the Siberian sector of the Arctic Ocean as well as the well-known polynya and fast-ice locations. Additionally, a substantial inter-annual variability of lead occurrences in the Arctic is indicated.

Geospatial analysis of groundwater depth for manual well drilling suitability in the Zinder Region, Niger

Manual and low-tech well drilling techniques have potential to assist in reaching the United Nations' millennium development goal for water in sub-Saharan Africa. This study used publicly available geospatial data in a regression tree analysis to predict groundwater depth in the Zinder region of Niger to identify suitable areas for manual well drilling. Regression trees were developed and tested on a database for 3681 wells in the Zinder region. A tree with 17 terminal leaves provided a range of ground water depth estimates that were appropriate for manual drilling, though much of the tree's complexity was associated with depths that were beyond manual methods. A natural log transformation of groundwater depth was tested to see if rescaling dataset variance would result in finer distinctions for regions of shallow groundwater. The RMSE for a log-transformed tree with only 10 terminal leaves was almost half that of the untransformed 17 leaf tree for groundwater depths less than 10 m. This analysis indicated important groundwater relationships for commonly available maps of geology, soils, elevation, and enhanced vegetation index from the MODIS satellite imaging system.

Probability sampling protocol of classification maps from spaceborne/airborne image

To deliver sample estimates provided with the necessary probability foundation to permit generalization from the sample data subset to the whole target population being sampled, probability sampling strategies are required to satisfy three necessary not sufficient conditions: (i) All inclusion probabilities be greater than zero in the target population to be sampled. If some sampling units have an inclusion probability of zero, then a map accuracy assessment does not represent the entire target region depicted in the map to be assessed. (ii) The inclusion probabilities must be: (a) knowable for nonsampled units and (b) known for those units selected in the sample: since the inclusion probability determines the weight attached to each sampling unit in the accuracy estimation formulas, if the inclusion probabilities are unknown, so are the estimation weights. This original work presents a novel (to the best of these authors' knowledge, the first) probability sampling protocol for quality assessment and comparison of thematic maps generated from spaceborne/airborne Very High Resolution (VHR) images, where: (I) an original Categorical Variable Pair Similarity Index (CVPSI, proposed in two different formulations) is estimated as a fuzzy degree of match between a reference and a test semantic vocabulary, which may not coincide, and (II) both symbolic pixel-based thematic quality indicators (TQIs) and sub-symbolic object-based spatial quality indicators (SQIs) are estimated with a degree of uncertainty in measurement in compliance with the well-known Quality Assurance Framework for Earth Observation (QA4EO) guidelines. Like a decision-tree, any protocol (guidelines for best practice) comprises a set of rules, equivalent to structural knowledge, and an order of presentation of the rule set, known as procedural knowledge. The combination of these two levels of knowledge makes an original protocol worth more than the sum of its parts. The several degrees of novelty of the proposed probability sampling protocol are highlighted in this paper, at the levels of understanding of both structural and procedural knowledge, in comparison with related multi-disciplinary works selected from the existing literature. In the experimental session the proposed protocol is tested for accuracy validation of preliminary classification maps automatically generated by the Satellite Image Automatic MapperT (SIAMT) software product from two WorldView-2 images and one QuickBird-2 image provided by DigitalGlobe for testing purposes. In these experiments, collected TQIs and SQIs are statistically valid, statistically significant, consistent across maps and in agreement with theoretical expectations, visual (qualitative) evidence and quantitative quality indexes of operativeness (OQIs) claimed for SIAMT by related papers. As a subsidiary conclusion, the statistically consistent and statistically significant accuracy validation of the SIAMT pre-classification maps proposed in this contribution, together with OQIs claimed for SIAMT by related works, make the operational (automatic, accurate, near real-time, robust, scalable) SIAMT software product eligible for opening up new inter-disciplinary research and market opportunities in accordance with the visionary goal of the Global Earth Observation System of Systems (GEOSS) initiative and the QA4EO international guidelines.

Inventory of thermo-erosional valleys and streams in three ice-rich permafrost lowlands adjacent to the Laptev Sea, with links to maps in ArcGIS format

This dataset provides an inventory of thermo-erosional landforms and streams in three lowland areas underlain by ice-rich permafrost of the Yedoma-type Ice Complex at the Siberian Laptev Sea coast. It consists of two shapefiles per study region: one shapefile for the digitized thermo-erosional landforms and streams, one for the study area extent. Thermo-erosional landforms were manually digitized from topographic maps and satellite data as line features and subsequently analyzed in a Geographic Information System (GIS) using ArcGIS 10.0. The mapping included in particular thermo-erosional gullies and valleys as well as streams and rivers, since development of all of these features potentially involved thermo-erosional processes. For the Cape Mamontov Klyk site, data from Grosse et al. [2006], which had been digitized from 1:100000 topographic map sheets, were clipped to the Ice Complex extent of Cape Mamontov Klyk, which excludes the hill range in the southwest with outcropping bedrock and rocky slope debris, coastal barrens, and a large sandy floodplain area in the southeast. The mapped features (streams, intermittent streams) were then visually compared with panchromatic Landsat-7 ETM+ satellite data (4 August 2000, 15 m spatial resolution) and panchromatic Hexagon data (14 July 1975, 10 m spatial resolution). Smaller valleys and gullies not captured in the maps were subsequently digitized from the satellite data. The criterion for the mapping of linear features as thermo-erosional valleys and gullies was their clear incision into the surface with visible slopes. Thermo-erosional features of the Lena Delta site were mapped on the basis of a Landsat-7 ETM+ image mosaic (2000 and 2001, 30 m ground resolution) [Schneider et al., 2009] and a Hexagon satellite image mosaic (1975, 10 m ground resolution) [G. Grosse, unpublished data] of the Lena River Delta within the extent of the Lena Delta Ice Complex [Morgenstern et al., 2011]. For the Buor Khaya Peninsula, data from Arcos [2012], which had been digitized based on RapidEye satellite data (8 August 2010, 6.5 m ground resolution), were completed for smaller thermo-erosional features using the same RapidEye scene as a mapping basis. The spatial resolution, acquisition date, time of the day, and viewing geometry of the satellite data used may have influenced the identification of thermo-erosional landforms in the images. For Cape Mamontov Klyk and the Lena Delta, thermo-erosional features were digitized using both Hexagon and Landsat data; Hexagon provided higher resolution and Landsat provided the modern extent of features. Allowance of up to decameters was made for the lateral expansion of features between Hexagon and Landsat acquisitions (between 1975 and 2000).

(Table 1) Sample locations and clay mineral data of sediment surface samples from the Barents Sea region

A new surface sediment sample set gained in the western Barents Sea by the MAREANO program has been analysed for basic clay mineral assemblages. Distribution maps including additional samples from earlier German research cruises to and off Svalbard are compiled. Some trends in the clay mineral assemblages are related to the sub-Barents Sea geology because the Quaternary sediment cover is rather thin. Additionally, land masses like Svalbard and northern Scandinavia dominate the clay mineral signal with their erosional products. Dense bottom water, very often of brine origin, that flows within deep troughs, such as the Storfjorden or Bear Island Trough, transport the clay mineral signal from their origin to the Norwegian-Greenland Sea.