The Lena River Delta - Land cover classification of tundra environments based on Landsat 7 ETM+ data and its application for upscaling of methane emissions

The Lena River Delta, situated in Northern Siberia (72.0 - 73.8° N, 122.0 - 129.5° E), is the largest Arctic delta and covers 29,000 km**2. Since natural deltas are characterised by complex geomorphological patterns and various types of ecosystems, high spatial resolution information on the distribution and extent of the delta environments is necessary for a spatial assessment and accurate quantification of biogeochemical processes as drivers for the emission of greenhouse gases from tundra soils. In this study, the first land cover classification for the entire Lena Delta based on Landsat 7 Enhanced Thematic Mapper (ETM+) images was conducted and used for the quantification of methane emissions from the delta ecosystems on the regional scale. The applied supervised minimum distance classification was very effective with the few ancillary data that were available for training site selection. Nine land cover classes of aquatic and terrestrial ecosystems in the wetland dominated (72%) Lena Delta could be defined by this classification approach. The mean daily methane emission of the entire Lena Delta was calculated with 10.35 mg CH4/m**2/d. Taking our multi-scale approach into account we find that the methane source strength of certain tundra wetland types is lower than calculated previously on coarser scales.

Documentation of glaciers and mountain chains in the Nordvestfjord of Scorsby Sound (East Greenland, July 2016) by landscape photography during Maria S. Merian cruise MSM56

From the 12th until the 17th of July 2016, research vessel Maria S. Merian entered the Nordvestfjord of Scorsby Sound (East Greenland) as part of research cruise MSM56, "Ecological chemistry in Arctic fjords". A large variety of chemical and biological parameters of fjord and meltwater were measured during this cruise to characterize biogeochemical fluxes in arctic fjords. The photo documentation described here was a side project. It was started when we were close to the Daugaard-Jensen glacier at the end of the Nordvestfjord and realized that not many people have seen this area before and photos available for scientists are probably rare. These pictures shall help to document climate and landscape changes in a remote area of East Greenland. Pictures were taken with a Panasonic Lumix G6 equipped with either a 14-42 or 45-150 objective (zoom factor available in jpg metadata). Polarizer filters were used on both objectives. The time between taking the pictures and writing down the coordinates was maximally one minute but usually shorter. The uncertainty in position is therefore small as we were steaming slowly most of the time the pictures were taken (i.e. below 5 knots). I assume the uncertainty is in most cases below 200 m radius of the noted position. I did not check the direction I directed the camera to with a compass at the beginning. Hence, the direction that was noted is an approximation based on the navigation map and the positioning of the ship. The uncertainty was probably around +/- 40° but initially (pictures 1-17) perhaps even higher as this documentation was a spontaneous idea and it took some time to get the orientation right. It should be easy, however, to find the location of the mountains and glaciers when being on the respective positions because the mountains have a quite characteristic shape. In a later stage of this documentation, I took pictures from the bridge and used the gyros to approximate the direction the camera was pointed at. Here the uncertainty was much lower (i.e. +/- 20° or better). Directions approximated with the help of gyros have degree values in the overview table. The ship data provided in the MSM56 cruise report will contain all kinds of sensor data from Maria S. Merian sensor setup. This data can also be used to further constrain the position the pictures were taken because the exact time a photo was shot is noted in the metadata of the .jpg photo file. The shipboard clock was set on UTC. It was 57 minutes and 45 seconds behind the time in the camera. For example 12:57:45 on the camera was 12:00:00 UTC on the ship. All pictures provided here can be used for scientific purposes. In case of usage in presentations etc. please acknowledge RV Maria S. Merian (MSM56) and Lennart T. Bach as author. Please inform me and ask for reprint permission in case you want to use the pictures for scientific publications. I would like to thank all participants and the crew of Maria S. Merian Cruise 56 (MSM56, Ecological chemistry in Arctic fjords).

GPS linked photos of benthic cover transect surveys in Heron Reef

Underwater photo-transect surveys were conducted on September 23-27, 2007 at different sections of the reef flat, reef crest and reef slope in Heron Reef. This survey was done by swimming along pre-defined transect sites and taking a picture of the bottom substrate parallel to the bottom at constant vertical distance (30cm) every two to three metres. A total of 3,586 benthic photos were taken. A floating GPS setup connected to the swimmer/diver by a line enabled recording of coordinates of transect surveys. Approximation of the coordinates for each benthic photo was based on the photo timestamp and GPS coordinate time stamp, using GPS Photo Link Software. Coordinates of each photo were interpolated by finding the the gps coordinates that were logged at a set time before and after the photo was captured. The output of this process was an ArcMap point shapefile, a Google Earth KML file and a thumbnail of each benthic photo taken. The data in the ArcMap shapefile and in the Google Earth KML file consisted of the approximated coordinate of each benthic photo taken during the survey. Using the GPS Photo Link extension within the ArcMap environment, opening the ArcMap shapefile will enable thumbnail to be displayed on the associated benthic cover photo whenever hovering with the mouse over a point on the transect. By downloading the GPSPhotoLink software from the www.geospatialexperts.com, and installing it as a trial version the ArcMap exstension will be installed in the ArcMap environment.

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.

AWI Bathymetric Chart of the Gakkel Ridge (AWI BCGR) (Scale 1:150,000)

The 1 : 1,500,000 AWI Bathymetric Chart of the Gakkel Ridge (AWI BCGR) has been developed from multibeam data measured during the Arctic Mid-Ocean Ridge Expedition in 2001 (AMORE 2001, ARK-XVII/2). This expedition was conducted to investigate the Gakkel Ridge in the Arctic Ocean and was carried out by the icebreaking research vessels RV Polarstern and USCGC Healy. Polarstern is equipped with the multibeam sonar system Hydrosweep DS-2, whereas Healy carries Seabeam 2112. During the expedition an area of 8890 km length and 18 - 46 km width, situated between 82°N/8°W and 87°N/75°E, was surveyed simultaneously by both vessels. Water depths ranged from 566 to 5673 meters. Dense sea ice cover derogated the sonar measurements and decreased data quality. Data errors were corrected in an extensive post-processing. The data of two different sonar systems had to be consolidated in order to derive a high resolution bathymetry of the Gakkel Ridge. Final result was a digital terrain model (DTM) with a grid spacing of 100 meters, which was utilized for generating the map series AWI Bathymetric Chart of the Gakkel Ridge, consisting of ten map sheets.

Circumarctic wetland dataset based on ENVISAT ASAR Global Monitoring Mode with links to geotiff image and NetCDF files (2005-01 to 2011-12)

A circumpolar representative and consistent wetland map is required for a range of applications ranging from upscaling of carbon fluxes and pools to climate modelling and wildlife habitat assessment. Currently available data sets lack sufficient accuracy and/or thematic detail in many regions of the Arctic. Synthetic aperture radar (SAR) data from satellites have already been shown to be suitable for wetland mapping. Envisat Advanced SAR (ASAR) provides global medium-resolution data which are examined with particular focus on spatial wetness patterns in this study. It was found that winter minimum backscatter values as well as their differences to summer minimum values reflect vegetation physiognomy units of certain wetness regimes. Low winter backscatter values are mostly found in areas vegetated by plant communities typically for wet regions in the tundra biome, due to low roughness and low volume scattering caused by the predominant vegetation. Summer to winter difference backscatter values, which in contrast to the winter values depend almost solely on soil moisture content, show expected higher values for wet regions. While the approach using difference values would seem more reasonable in order to delineate wetness patterns considering its direct link to soil moisture, it was found that a classification of winter minimum backscatter values is more applicable in tundra regions due to its better separability into wetness classes. Previous approaches for wetland detection have investigated the impact of liquid water in the soil on backscatter conditions. In this study the absence of liquid water is utilized. Owing to a lack of comparable regional to circumpolar data with respect to thematic detail, a potential wetland map cannot directly be validated; however, one might claim the validity of such a product by comparison with vegetation maps, which hold some information on the wetness status of certain classes. It was shown that the Envisat ASAR-derived classes are related to wetland classes of conventional vegetation maps, indicating its applicability; 30% of the land area north of the treeline was identified as wetland while conventional maps recorded 1-7%.

Sea surface and the sea floor regionalization of the Southern Ocean by multivariate cluster analysis, links to ArcGIS project files

This study subdivides the Weddell Sea, Antarctica, into seafloor regions using multivariate statistical methods. These regions are categories used for comparing, contrasting and quantifying biogeochemical processes and biodiversity between ocean regions geographically but also regions under development within the scope of global change. The division obtained is characterized by the dominating components and interpreted in terms of ruling environmental conditions. The analysis uses 28 environmental variables for the sea surface, 25 variables for the seabed and 9 variables for the analysis between surface and bottom variables. The data were taken during the years 1983-2013. Some data were interpolated. The statistical errors of several interpolation methods (e.g. IDW, Indicator, Ordinary and Co-Kriging) with changing settings have been compared for the identification of the most reasonable method. The multivariate mathematical procedures used are regionalized classification via k means cluster analysis, canonical-correlation analysis and multidimensional scaling. Canonical-correlation analysis identifies the influencing factors in the different parts of the cove. Several methods for the identification of the optimum number of clusters have been tested. For the seabed 8 and 12 clusters were identified as reasonable numbers for clustering the Weddell Sea. For the sea surface the numbers 8 and 13 and for the top/bottom analysis 8 and 3 were identified, respectively. Additionally, the results of 20 clusters are presented for the three alternatives offering the first small scale environmental regionalization of the Weddell Sea. Especially the results of 12 clusters identify marine-influenced regions which can be clearly separated from those determined by the geological catchment area and the ones dominated by river discharge.

Mapped cetacean habitat suitability and richness in the Azores, links to ArcGIS files

Marine spatial planning and ecological research call for high-resolution species distribution data. However, those data are still not available for most marine large vertebrates. The dynamic nature of oceanographic processes and the wide-ranging behavior of many marine vertebrates create further difficulties, as distribution data must incorporate both the spatial and temporal dimensions. Cetaceans play an essential role in structuring and maintaining marine ecosystems and face increasing threats from human activities. The Azores holds a high diversity of cetaceans but the information about spatial and temporal patterns of distribution for this marine megafauna group in the region is still very limited. To tackle this issue, we created monthly predictive cetacean distribution maps for spring and summer months, using data collected by the Azores Fisheries Observer Programme between 2004 and 2009. We then combined the individual predictive maps to obtain species richness maps for the same period. Our results reflect a great heterogeneity in distribution among species and within species among different months. This heterogeneity reflects a contrasting influence of oceanographic processes on the distribution of cetacean species. However, some persistent areas of increased species richness could also be identified from our results. We argue that policies aimed at effectively protecting cetaceans and their habitats must include the principle of dynamic ocean management coupled with other area-based management such as marine spatial planning.

Viral, bacterial and ciliate abundance and viral-bacterial diversity in mesocosm experiments, 2007-2008, Kongsfjorden

Two mesocosm experiments, PAME-I and PAME-II were conducted in 2007 and 2008 to investigate fate of organic carbon in the arctic microbial food web. Mesocosms were nutrient fertilized initially to induce phytoplankton bloom development. In PAME-I eight units (each 700 L) formed two four point gradients of additional DOC in form of glucose (0, 0.5, 1 and 3 times Redfield ratio in terms of carbon relative to the nitrogen and phosphorus additions) (Fig. 1). All the eight units also got a daily dose of NH4+ and PO4**3- in Redfield ratio. Two gradients were set up, one with silicate addition, performed in the Arctic location Ny Ålesund, Svalbard, have previously been reported to give different food-web level responses to similar nutrient perturbations. In PAME-II all ten units (each 900 L) formed two four point gradients of additional DOC in form of glucose (0, 0.5, 1, 2 and 3 times Redfield ratio in terms of carbon relative to nitrogen and phosphorus additions). The two gradients in glucose were kept silicate replete. NH4+ was used as the DIN source in one gradient (units 1 to 5) and NO3- in the other (units 6-9). All units got a daily dose of PO4**3- in Redfield ratio. Prokaryotes and viruses were measured by flow cytometry, while ciliate abundances were counted using a Flow Cam. Viral and bacterial diversity was measured by PFGE and DGGE, respectively. In PAME-II the abundance of ciliates was lower than in PAME-I, presumably caused by higher copepod grazing. The abundances of prokaryotes and viruses were also lower in PAME-II compared to PAME-I. Further, less diversity was detected in the viral community (FCM and PFGE) in PAME-II, and no response was observed in the bacterial community structure due to addition of organic carbon.

Surface wind stress modeled for the Panama Seaway with the Community Climate System Model 3

The early Pliocene warm phase was characterized by high sea surface temperatures and a deep thermocline in the eastern equatorial Pacific. A new hypothesis suggests that the progressive closure of the Panamanian seaway contributed substantially to the termination of this zonally symmetric state in the equatorial Pacific. According to this hypothesis, intensification of the Atlantic meridional overturning circulation (AMOC) - induced by the closure of the gateway - was the principal cause of equatorial Pacific thermocline shoaling during the Pliocene. In this study, twelve Panama seaway sensitivity experiments from eight ocean/climate models of different complexity are analyzed to examine the effect of an open gateway on AMOC strength and thermocline depth. All models show an eastward Panamanian net throughflow, leading to a reduction in AMOC strength compared to the corresponding closed-Panama case. In those models that do not include a dynamic atmosphere, deepening of the equatorial Pacific thermocline appears to scale almost linearly with the throughflow-induced reduction in AMOC strength. Models with dynamic atmosphere do not follow this simple relation. There are indications that in four out of five models equatorial wind-stress anomalies amplify the tropical Pacific thermocline deepening. In summary, the models provide strong support for the hypothesized relationship between Panama closure and equatorial Pacific thermocline shoaling.