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.