Surface velocity fields, digital elevation models and ice front positions derived from multi-mission SAR remote sensing data at Sjögren Inlet glaciers, Antarctic Peninsula

Substantial retreat or disintegration of numerous ice shelves have been observed on the Antarctic Peninsula. The ice shelf in the Prince Gustav Channel retreated gradually since the late 1980's and broke-up in 1995. Tributary glaciers reacted with speed-up, surface lowering and increased ice discharge, consequently contributing to sea level rise. We present a detailed long-term study (1993-2014) on the dynamic response of Sjögren Inlet glaciers to the disintegration of Prince Gustav Ice Shelf. We analyzed various remote sensing datasets to observe the reactions of the glaciers to the loss of the buttressing ice shelf. A strong increase in ice surface velocities was observed with maximum flow speeds reaching 2.82±0.48 m/d in 2007 and 1.50±0.32 m/d in 2004 at Sjögren and Boydell glaciers respectively. Subsequently, the flow velocities decelerated, however in late 2014, we still measured about two times the values of our first measurements in 1996. The tributary glaciers retreated 61.7±3.1 km² behind the former grounding line of the ice shelf. In regions below 1000 m a.s.l., a mean surface lowering of -68±10 m (-3.1 m/a) was observed in the period 1993-2014. The lowering rate decreased to -2.2 m/a in recent years. Based on the surface lowering rates, geodetic mass balances of the glaciers were derived for different time steps. High mass loss rate of -1.21±0.36 Gt/a was found in the earliest period (1993-2001). Due to the dynamic adjustments of the glaciers to the new boundary conditions the ice mass loss reduced to -0.59±0.11 Gt/a in the period 2012-2014, resulting in an average mass loss rate of -0.89±0.16 Gt/a (1993-2014). Including the retreat of the ice front and grounding line, a total mass change of -38.5±7.7 Gt and a contribution to sea level rise of 0.061±0.013 mm were computed. Analysis of the ice flux revealed that available bedrock elevation estimates at Sjögren Inlet are too shallow and are the major uncertainty in ice flux computations. This temporally dense time series analysis of Sjögren Inlet glaciers shows that the adjustments of tributary glaciers to ice shelf disintegration are still going on and provides detailed information of the changes in glacier dynamics.

Downwelling spectral irradiance as measured in different depths, remote sensing reflectance and phytoplankton pigment concentration during POLARSTERN cruises ANT-XXIV/1, ANT-XXIV/4, ANT-XXVI/4 and Maria S. Merian cruise MSM18/3

The composition and abundance of algal pigments provide information on phytoplankton community characteristics such as photoacclimation, overall biomass and taxonomic composition. In particular, pigments play a major role in photoprotection and in the light-driven part of photosynthesis. Most phytoplankton pigments can be measured by high-performance liquid chromatography (HPLC) techniques applied to filtered water samples. This method, as well as other laboratory analyses, is time consuming and therefore limits the number of samples that can be processed in a given time. In order to receive information on phytoplankton pigment composition with a higher temporal and spatial resolution, we have developed a method to assess pigment concentrations from continuous optical measurements. The method applies an empirical orthogonal function (EOF) analysis to remote-sensing reflectance data derived from ship-based hyperspectral underwater radiometry and from multispectral satellite data (using the Medium Resolution Imaging Spectrometer - MERIS - Polymer product developed by Steinmetz et al., 2011, doi:10.1364/OE.19.009783) measured in the Atlantic Ocean. Subsequently we developed multiple linear regression models with measured (collocated) pigment concentrations as the response variable and EOF loadings as predictor variables. The model results show that surface concentrations of a suite of pigments and pigment groups can be well predicted from the ship-based reflectance measurements, even when only a multispectral resolution is chosen (i.e., eight bands, similar to those used by MERIS). Based on the MERIS reflectance data, concentrations of total and monovinyl chlorophyll a and the groups of photoprotective and photosynthetic carotenoids can be predicted with high quality. As a demonstration of the utility of the approach, the fitted model based on satellite reflectance data as input was applied to 1 month of MERIS Polymer data to predict the concentration of those pigment groups for the whole eastern tropical Atlantic area. Bootstrapping explorations of cross-validation error indicate that the method can produce reliable predictions with relatively small data sets (e.g., < 50 collocated values of reflectance and pigment concentration). The method allows for the derivation of time series from continuous reflectance data of various pigment groups at various regions, which can be used to study variability and change of phytoplankton composition and photophysiology.

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.

Arctic Siberia: Modern pollen dataset and pollen spectra changes over the last 7,000 years from a lacustrine sediment core (southern Taymyr region)

Palaeoecological investigations in the larch forest-tundra ecotone in northern Siberia have the potential to reveal Holocene environmental variations, which likely have consequences for global climate change because of the strong high-latitude feedback mechanisms. A sediment core, collected from a small lake (radius ~100 m), was used to reconstruct the development of the lake and its catchment as well as vegetation and summer temperatures over the last 7100 calibrated years. A multi-proxy approach was taken including pollen and sedimentological analyses. Our data indicate a gradual replacement of open larch forests by tundra with scattered single trees as found today in the vicinity of the lake. An overall trend of cooling summer temperature from a ~2 °C warmer-than-present mid-Holocene summer temperatures until the establishment of modern conditions around 3000 years ago is reconstructed based on a regional pollen-climate transfer function. The inference of regional vegetation changes was compared to local changes in the lake's catchment. An initial small water depression occurred from 7100 to 6500 cal years BP. Afterwards, a small lake formed and deepened, probably due to thermokarst processes. Although the general trends of local and regional environmental change match, the lake catchment changes show higher variability. Furthermore, changes in the lake catchment slightly precede those in the regional vegetation. Both proxies highlight that marked environmental changes occurred in the Siberian forest-tundra ecotone over the course of the Holocene.

Satellite-based daily lake areas in the Yangtze Basin, China, 2000-2011

The Yangtze River Basin downstream of China's Three Gorges Dam (TGD) (thereafter referred to as "downstream" basin) hosts the largest cluster of freshwater lakes in East Asia. These lakes are crucial water stocks to local biophysical environments and socioeconomic development. Existing studies document that individual lakes in this region have recently experienced dramatic changes under the context of enduring meteorological drought, continuous population growth, and extensive water regulation since TGD's initial impoundment (i.e., June, 2003). However, spatial and temporal patterns of lake dynamics across the complete downstream Yangtze basin remain poorly characterized. Using daily MODIS imagery and an advanced thematic mapping scheme, this study presents a comprehensive monitoring of area dynamics in the downstream lake system at a 10-day temporal resolution during 2000-2011. The studied lakes constitute ~76% (~11,400 km**2) of the total downstream lake area, including the entire +70 major lakes larger than 20 km**2. The results reveal a decadal net decline in lake inundation area across the downstream Yangtze Basin, with a cumulative decrease of 849 km**2 or 7.4% from 2000 to 2011. Despite an excessive precipitation anomaly in the year 2010, the decreasing trend was tested significant in all seasons. The most substantial decrease in the post-TGD period appears in fall (1.1%/yr), which intriguingly coincides with the TGD water storage season. Regional lake dynamics exhibit contrasting spatial patterns, manifested as evident decrease and increase of aggregated lake areas respectively within and beyond the Yangtze Plain. This contrast suggests a marked vulnerability of lakes in the Yangtze Plain, to not only local meteorological variability but also intensified human water regulations from both the upstream Yangtze main stem (e.g., the TGD) and tributaries (e.g., lakes/reservoirs beyond the Yangtze Plain). The produced lake mapping result and derived lake area dynamics across the downstream Yangtze Basin provides a crucial monitoring basis for continuous investigations of changing mechanisms in the Yangtze lake system.

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.

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.