Seasonal mean air temperatures and ice phenology characteristics of Great Bear Lake and Great Slave Lake (2002-2009)

Time series of brightness temperatures (T(B)) from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) are examined to determine ice phenology variables on the two largest lakes of northern Canada: Great Bear Lake (GBL) and Great Slave Lake (GSL). T(B) measurements from the 18.7, 23.8, 36.5, and 89.0 GHz channels (H- and V- polarization) are compared to assess their potential for detecting freeze-onset/melt-onset and ice-on/ice-off dates on both lakes. The 18.7 GHz (H-pol) channel is found to be the most suitable for estimating these ice dates as well as the duration of the ice cover and ice-free seasons. A new algorithm is proposed using this channel and applied to map all ice phenology variables on GBL and GSL over seven ice seasons (2002-2009). Analysis of the spatio-temporal patterns of each variable at the pixel level reveals that: (1) both freeze-onset and ice-on dates occur on average about one week earlier on GBL than on GSL (Day of Year (DY) 318 and 333 for GBL; DY 328 and 343 for GSL); (2) the freeze-up process or freeze duration (freeze-onset to ice-on) takes a slightly longer amount of time on GBL than on GSL (about 1 week on average); (3) melt-onset and ice-off dates occur on average one week and approximately four weeks later, respectively, on GBL (DY 143 and 183 for GBL; DY 135 and 157 for GSL); (4) the break-up process or melt duration (melt-onset to ice-off) lasts on average about three weeks longer on GBL; and (5) ice cover duration estimated from each individual pixel is on average about three weeks longer on GBL compared to its more southern counterpart, GSL. A comparison of dates for several ice phenology variables derived from other satellite remote sensing products (e.g. NOAA Interactive Multisensor Snow and Ice Mapping System (IMS), QuikSCAT, and Canadian Ice Service Database) show that, despite its relatively coarse spatial resolution, AMSR-E 18.7 GHz provides a viable means for monitoring of ice phenology on large northern lakes.

Global database of surface ocean particulate organic carbon export fluxes diagnosed from the 234Th technique

Through the processes of the biological pump, carbon is exported to the deep ocean in the form of dissolved and particulate organic matter. There are several ways by which downward export fluxes can be estimated. The great attraction of the 234Th technique is that its fundamental operation allows a downward flux rate to be determined from a single water column profile of thorium coupled to an estimate of POC/234Th ratio in sinking matter. We present a database of 723 estimates of organic carbon export from the surface ocean derived from the 234Th technique. Data were collected from tables in papers published between 1985 and 2013 only. We also present sampling dates, publication dates and sampling areas. Most of the open ocean Longhurst provinces are represented by several measurements. However, the Western Pacific, the Atlantic Arctic, South Pacific and the South Indian Ocean are not well represented. There is a variety of integration depths ranging from surface to 220m. Globally the fluxes ranged from -22 to 125 mmol of C/m**2/d. We believe that this database is important for providing new global estimate of the magnitude of the biological carbon pump.

(Table S2) Biometric data, time of ice breakup and carbon isotopic composition of adipose tissue of polar bears (Ursus maritimus) between 1991-2007

Two global environmental issues, climate change and contamination by persistent organic pollutants, represent major concerns for arctic ecosystems. Yet, it is unclear how these two stressors interact in the Arctic. For instance, the influence of climate-associated changes in food web structure on exposure to pollutants within arctic ecosystems is presently unknown. Here, we report on recent changes in feeding ecology (1991-2007) in polar bears (Ursus maritimus) from the western Hudson Bay subpopulation that have resulted in increases in the tissue concentrations of several chlorinated and brominated contaminants. Differences in timing of the annual sea ice breakup explained a significant proportion of the diet variation among years. As expected from climate change predictions, this diet change was consistent with an increase in the consumed proportions of open water-associated seal species compared to ice-associated seal species in years of earlier sea ice breakup. Our results demonstrate that climate change is a modulating influence on contaminants in this polar bear subpopulation and may pose an additional and previously unidentified threat to northern ecosystems through altered exposures to contaminants.

Modelled spatial and seasonal distribution of Blue Whiting (Micromesistius poutassou) larvae in the North-East Atlantic (1951 to 2005)

Blue whiting (Micromesistius poutassou, http://www.marinespecies.org/aphia.php?p=taxdetails&id=126439) is a small mesopelagic planktivorous gadoid found throughout the North-East Atlantic. This data contains the results of a model-based analysis of larvae captured by the Continuous Plankton Recorder (CPR) during the period 1951-2005. The observations are analysed using Generalised Additive Models (GAMs) of the the spatial, seasonal and interannual variation in the occurrence of larvae. The best fitting model is chosen using the Aikaike Information Criteria (AIC). The probability of occurrence in the continous plankton recorder is then normalised and converted to a probability distribution function in space (UTM projection Zone 28) and season (day of year). The best fitting model splits the distribution into two separate spawning grounds north and south of a dividing line at 53 N. The probability distribution is therefore normalised in these two regions (ie the space-time integral over each of the two regions is 1). The modelled outputs are on a UTM Zone 28 grid: however, for convenience, the latitude ("lat") and longitude ("lon") of each of these grid points are also included as a variable in the NetCDF file. The assignment of each grid point to either the Northern or Southern component (defined here as north/south of 53 N), is also included as a further variable ("component"). Finally, the day of year ("doy") is stored as the number of days elapsed from and included January 1 (ie doy=1 on January 1) - the year is thereafter divided into 180 grid points.

Eddy covariance and environmental data of boreal bogs plant species

In boreal bogs plant species are low in number, but they differ greatly in their growth forms and photosynthetic properties. We assessed how ecosystem carbon (C) sink dynamics were affected by seasonal variations in photosynthetic rate and leaf area of different species. Photosynthetic properties (light-response parameters), leaf area development and areal cover (abundance) of the species were used to quantify species-specific net and gross photosynthesis rates (PN and PG, respectively), which were summed to express ecosystem-level PN and PG. The ecosystem-level PG was compared with a gross primary production (GPP) estimate derived from eddy covariance measurements (EC). Species areal cover rather than differences in photosynthetic properties determined the species with the highest PG of both vascular plants and Sphagna. Species-specific contributions to the ecosystem PG varied over the growing season, which in turn determined the seasonal variation in ecosystem PG. The upscaled growing-season PG estimate, 230 g C/m**2, agreed well with the GPP estimated by the EC, 243 g C/m**2. Sphagna were superior to vascular plants in ecosystem-level PG throughout the growing season but had a lower PN. PN results indicated that areal cover of the species together with their differences in photosynthetic parameters shape the ecosystem-level C balance. Species with low areal cover but high photosynthetic efficiency appear to be potentially important for the ecosystem C sink. Results imply that functional diversity may increase the stability of C sink of boreal bogs.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.