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.

Phosphorus dynamics in waters of the Pacific Ocean

Tagged phosphorus was used to measure principal indices of mineral phosphorus variations in the euphotic zone of the East Pacific, i.e. total rate of uptake of phosphate phosphorus by microplankton (A_t), fraction consumed by phytoplankton (A_p/A_t), and turnover time (T). A_t reached its greatest values (150-280 ng/l/hour) in the upwelling zone of the Peru traverse, where development of phytoplankton was induced by upwelling. In other areas of this traverse values were 40-80 ng/l/hour in surface layers. In less productive waters on two other profiles (off Central America and California), values were lower, between 20 and 40 ng/l. On the vertical profile maxima of A_t were found at the upper boundary of the thermocline. Turnover time of PO4 phosphorus (T) in zones of phytoplankton abundance was very short, between 1.5 and 4 days. At most other stations it was 10-40 days, increasing to 100-200 days or longer at the lower boundary of the euphotic zone. In areas of phytoplankton abundance it accounted for 60-80% of total uptake of PO4 phosphorus. But in zones of elevated bacterial abundance, A_p/A_t fell to 20-40%. Data indicating lack of correlation between PO4 phosphorus and productivity are presented. It is emphasized that the above measures of PO4 phosphorus dynamics can be used for obtaining measures of functional condition and successional phase of marine plankton communities.

Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

Modelling calving front dynamics using a level-set method: application to Jakobshavn Isbræ, West Greenland, links model results in MATLAB format

Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates. We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales.

Dynamics of halocarbons in coastal surface waters during short term mesocosm experiments

The aim of the present study was to evaluate the influence of different light quality, especially ultraviolet radiation (UVR), on the dynamics of volatile halogenated organic compounds (VHOCs) at the sea surface. Short term experiments were conducted with floating gas-tight mesocosms of different optical qualities. Six halocarbons (CH3I, CHCl3, CH2Br2, CH2ClI, CHBr3 and CH2I2), known to be produced by phytoplankton, together with a variety of biological and environmental variables were measured in the coastal southern Baltic Sea and in the Raunefjord (North Sea). These experiments showed that ambient levels of UVR have no significant influence on VHOC dynamics in the natural systems. We attribute it to the low radiation doses that phytoplankton cells receive in a normal turbulent surface mixed layer. The VHOC concentrations were influenced by their production and removal processes, but they were not correlated with biological or environmental parameters investigated. Diatoms were most likely the dominant biogenic source of VHOCs in the Baltic Sea experiment, whereas in the Raunefjord experiment macroalgae probably contributed strongly to the production of VHOCs. The variable stable carbon isotope signatures (d13C values) of bromoform (CHBr3) also indicate that different autotrophic organisms were responsible for CHBr3 production in the two coastal environments. In the Raunefjord, despite strong daily variations in CHBr3 concentration, the carbon isotopic ratio was fairly stable with a mean value of -26 per mil. During the declining spring phytoplankton bloom in the Baltic Sea, the d13C values of CHBr3 were enriched in 13C and showed noticeable diurnal changes (-12 per mil ± 4). These results show that isotope signature analysis is a useful tool to study both the origin and dynamics of VHOCs in natural systems.