Temperature, water level and bathymetry of thermokarst lakes in the continuous permafrost zone of northern Siberia - Lena River Delta, Siberia

Thermokarst lakes are typical features of the northern permafrost ecosystems, and play an important role in the thermal exchange between atmosphere and subsurface. The objective of this study is to describe the main thermal processes of the lakes and to quantify the heat exchange with the underlying sediments. The thermal regimes of five lakes located within the continuous permafrost zone of northern Siberia (Lena River Delta) were investigated using hourly water temperature and water level records covering a 3-year period (2009-2012), together with bathymetric survey data. The lakes included thermokarst lakes located on Holocene river terraces that may be connected to Lena River water during spring flooding, and a thermokarst lake located on deposits of the Pleistocene Ice Complex. Lakes were covered by ice up to 2 m thick that persisted for more than 7 months of the year, from October until about mid-June. Lake-bottom temperatures increased at the start of the ice-covered period due to upward-directed heat flux from the underlying thawed sediment. Prior to ice break-up, solar radiation effectively warmed the water beneath the ice cover and induced convective mixing. Ice break-up started at the beginning of June and lasted until the middle or end of June. Mixing occurred within the entire water column from the start of ice break-up and continued during the ice-free periods, as confirmed by the Wedderburn numbers, a quantitative measure of the balance between wind mixing and stratification that is important for describing the biogeochemical cycles of lakes. The lake thermal regime was modeled numerically using the FLake model. The model demonstrated good agreement with observations with regard to the mean lake temperature, with a good reproduction of the summer stratification during the ice-free period, but poor agreement during the ice-covered period. Modeled sensitivity to lake depth demonstrated that lakes in this climatic zone with mean depths > 5 m develop continuous stratification in summer for at least 1 month. The modeled vertical heat flux across the bottom sediment tends towards an annual mean of zero, with maximum downward fluxes of about 5 W/m**2 in summer and with heat released back into the water column at a rate of less than 1 W/m**2 during the ice-covered period. The lakes are shown to be efficient heat absorbers and effectively distribute the heat through mixing. Monthly bottom water temperatures during the ice-free period range up to 15 °C and are therefore higher than the associated monthly air or ground temperatures in the surrounding frozen permafrost landscape. The investigated lakes remain unfrozen at depth, with mean annual lake-bottom temperatures of between 2.7 and 4 °C.

Spatial distribution of key zooplankton species using continuous plankton recorder (CPR) data from the North Atlantic (2000-2009)

The continuous plankton recorder (CPR) survey is an upper layer plankton monitoring program that has regularly collected samples, at monthly intervals, in the North Atlantic and adjacent seas since 1946. Water from approximately 6 m depth enters the CPR through a small aperture at the front of the sampler and travels down a tunnel where it passes through a silk filtering mesh of 270 µm before exiting at the back of the CPR. The plankton filtered on the silk is analyzed in sections corresponding to 10 nautical miles (approx. 3 m**3 of seawater filtered) and the plankton microscopically identified (Richardson et al., 2006 and reference therein). In the present study we used the CPR data to investigate the current basin scale distribution of C. finmarchicus (C5-C6), C. helgolandicus (C5-C6), C. hyperboreus (C5-C6), Pseudocalanus spp. (C6), Oithona spp. (C1-C6), total Euphausiida, total Thecosomata and the presence/absence of Cnidaria and the Phytoplankton Colour Index (PCI). The PCI, which is a visual assessment of the greenness of the silk, is used as an indicator of the distribution of total phytoplankton biomass across the Atlantic basin (Batten et al., 2003). Monthly data collected between 2000 and 2009 were gridded using the inverse-distance interpolation method, in which the interpolated values were the nodes of a 2 degree by 2 degree grid. The resulting twelve monthly matrices were then averaged within the year and in the case of the zooplankton the data were log-transformed (i.e. log10 (x+1).

Data and Code archive for the interpolation of surface ocean carbon dioxide

We have developed a statistical gap-filling method adapted to the specific coverage and properties of observed fugacity of surface ocean CO2 (fCO2). We have used this method to interpolate the Surface Ocean CO2 Atlas (SOCAT) v2 database on a 2.5°×2.5° global grid (south of 70°N) for 1985-2011 at monthly resolution. The method combines a spatial interpolation based on a 'radius of influence' to determine nearby similar fCO2 values with temporal harmonic and cubic spline curve-fitting, and also fits long term trends and seasonal cycles. Interannual variability is established using deviations of observations from the fitted trends and seasonal cycles. An uncertainty is computed for all interpolated values based on the spatial and temporal range of the interpolation. Tests of the method using model data show that it performs as well as or better than previous regional interpolation methods, but in addition it provides a near-global and interannual coverage.