EOT11a - Global Empirical Ocean Tide model from multi-mission satellite altimetry, with links to model results

EOT11a is a global (E)mpirical (O)cean (T)ide model derived in 2011 by residual analysis of multi-mission satellite (a)ltimeter data. EOT11a includes amplitudes and phases of the main astronomical tides M2, S2, N2, K2, 2N2, O1, K1, P2, and Q1, the non-linear constituent M4, the long period tides Mm and Mf, and the radiational tide S1. Ocean tides as well as loading tides are provided. EOT11a was computed by means of residual tidal analysis of multi-mission altimeter data from TOPEX/Poseidon, ERS-2, ENVISAT, and Jason-1/2, as far as acquired between September 1992 and April 2010. The resolution of 7.5'x7.5' is identical with FES2004 which was used as reference model for the residual tide analysis. The development of EOT11a was funded by the Deutsche Forschungsgemeinschaft (DFG) under grant BO1228/6-2.

Holocene carbon cycle dynamics, links to model files

We are investigating the late Holocene rise in CO2 by performing four experiments with the climate-carbon-cycle model CLIMBER2-LPJ. Apart from the deep sea sediments, important carbon cycle processes considered are carbon uptake or release by the vegetation, carbon uptake by peatlands, and CO 2 release due to shallow water sedimentation of CaCO3. Ice core data of atmospheric CO2 between 8 ka BP and preindustrial climate can only be reproduced if CO2 outgassing due to shallow water sedimentation of CaCO3 is considered. In this case the model displays an increase of nearly 20 ppmv CO2 between 8 ka BP and present day. Model configurations that do not contain this forcing show a slight decrease in atmospheric CO2. We can therefore explain the late Holocene rise in CO2 by invoking natural forcing factors only, and anthropogenic forcing is not required to understand preindustrial CO2 dynamics.

High resolution bathymetric compilation for Potter Cove, WAP, Antarctica, with links to data in ArcGIS format

The bathymetry raster with a resolution of 5 m x 5 m was processed from unpublished single beam data from the Argentine Antarctica Institute (IAA, 2010) and multibeam data from the United Kingdom Hydrographic Office (UKHO, 2012) with a cell size of 5 m x 5 m. A coastline digitized from a satellite image (DigitalGlobe, 2014) supplemented the interpolation process. The 'Topo to Raster' tool in ArcMap 10.3 was used to merge the three data sets, while the coastline represented the 0-m-contour to the interpolation process ('contour type option').

Ocean temperature response to idealized Gleissberg and de Vries solar cycles in a comprehensive climate model

The ~90-year Gleissberg and ~200-year de Vries cycles have been identified as two distinctive quasi-periodic components of Holocene solar activity. Evidence exists for the impact of such multi-decadal to centennial-scale variability in total solar irradiance (TSI) on climate, but concerning the ocean, this evidence is mainly restricted to the surface response. Here we use a comprehensive global climate model to study the impact of idealized solar forcing, representing the Gleissberg and de Vries cycles, on global ocean potential temperature at different depth levels, after a recent proxy record indicates a signal of TSI anomalies in the northeastern Atlantic at mid-depth. Potential impacts of TSI anomalies on deeper oceanic levels are climatically relevant due to their possible effect on ocean circulation by altering water mass characteristics. Simulated solar anomalies are shown to penetrate the ocean down to at least deep-water levels. Despite the fact that the two forcing periods differ only by a factor of ~2, the spatial pattern of response is significantly distinctive between the experiments, suggesting different mechanisms for solar signal propagation. These are related to advection by North Atlantic Deep Water flow (200-year forcing), and barotropic adjustment in the South Atlantic in response to a latitudinal shift of the westerly wind belt (90-year forcing).

Direct visualization of solute locations in laboratory ice samples, links to videos

Many important chemical reactions occur in polar snow, where solutes may be present in several reservoirs, including at the air-ice interface and in liquid-like regions within the ice matrix. Some recent laboratory studies suggest chemical reaction rates may differ in these two reservoirs. While investigations have examined where solutes are found in natural snow and ice, similar research has not identified solute locations in laboratory samples, nor the possible factors controlling solute segregation. To address this, we examined solute locations in ice samples prepared from either aqueous cesium chloride (CsCl) or Rose Bengal solutions that were frozen using several different methods. Samples frozen in a laboratory freezer had the largest liquid-like inclusions and air bubbles, while samples frozen in a custom freeze chamber had somewhat smaller air bubbles and inclusions; in contrast, samples frozen in liquid nitrogen showed much smaller concentrated inclusions and air bubbles, only slightly larger than the resolution limit of our images (~2 µm). Freezing solutions in plastic versus glass vials had significant impacts on the sample structure, perhaps because the poor heat conductivity of plastic vials changes how heat is removed from the sample as it cools. Similarly, the choice of solute had a significant impact on sample structure, with Rose Bengal solutions yielding smaller inclusions and air bubbles compared to CsCl solutions frozen using the same method. Additional experiments using higher-resolution imaging of an ice sample show that CsCl moves in a thermal gradient, supporting the idea that the solutes in ice are present in liquid-like regions. Our work shows that the structure of laboratory ice samples, including the location of solutes, is sensitive to freezing method, sample container, and solute characteristics, requiring careful experimental design and interpretation of results.

Wetland maps including open water extent dynamics based on ENVISAT ASAR WS for Siberia, 2007 and 2008, links to GeoTIFFs

Wetlands store large amounts of carbon, and depending on their status and type, they release specific amounts of methane gas to the atmosphere. The connection between wetland type and methane emission has been investigated in various studies and utilized in climate change monitoring and modelling. For improved estimation of methane emissions, land surface models require information such as the wetland fraction and its dynamics over large areas. Existing datasets of wetland dynamics present the total amount of wetland (fraction) for each model grid cell, but do not discriminate the different wetland types like permanent lakes, periodically inundated areas or peatlands. Wetland types differently influence methane fluxes and thus their contribution to the total wetland fraction should be quantified. Especially wetlands of permafrost regions are expected to have a strong impact on future climate due to soil thawing. In this study ENIVSAT ASAR Wide Swath data was tested for operational monitoring of the distribution of areas with a long-term SW near 1 (hSW) in northern Russia (SW = degree of saturation with water, 1 = saturated), which is a specific characteristic of peatlands. For the whole northern Russia, areas with hSW were delineated and discriminated from dynamic and open water bodies for the years 2007 and 2008. The area identified with this method amounts to approximately 300,000 km**2 in northern Siberia in 2007. It overlaps with zones of high carbon storage. Comparison with a range of related datasets (static and dynamic) showed that hSW represents not only peatlands but also temporary wetlands associated with post-forest fire conditions in permafrost regions. Annual long-term monitoring of change in boreal and tundra environments is possible with the presented approach. Sentinel-1, the successor of ENVISAT ASAR, will provide data that may allow continuous monitoring of these wetland dynamics in the future complementing global observations of wetland fraction.

High-resolution ice thickness and bed topography of a land-terminating section of the Greenland Ice Sheet, with links to GeoTIFFs and text files

We present ice thickness and bed topography maps with a high spatial resolution (250-500 m) of a land-terminating section of the Greenland Ice Sheet derived from ground-based and airborne radar surveys. The data have a total area of ~12 000 km^2 and cover the whole ablation area of the outlet glaciers of Isunnguata Sermia, Russell, Leverett, Ørkendalen and Isorlersuup up to the long-term mass balance equilibrium line altitude at ~1600 m above sea level. The bed topography shows highly variable subglacial trough systems, and the trough of Isunnguata Sermia Glacier is overdeepened and reaches an elevation of ~500 m below sea level. The ice surface is smooth and only reflects the bedrock topography in a subtle way, resulting in a highly variable ice thickness. The southern part of our study area consists of higher bed elevations compared to the northern part. The compiled data sets of ground-based and airborne radar surveys cover one of the most studied regions of the Greenland Ice Sheet and can be valuable for detailed studies of ice sheet dynamics and hydrology.

Near-bottom zooplankton aggregations in Kongsfjorden with link to images from the optical zooplankton sensor MOKI

Near-bottom zooplankton communities have rarely been studied despite numerous reports of high zooplankton concentrations, probably due to methodological constraints. In Kongsfjorden, Svalbard, the near-bottom layer was studied for the first time by combining daytime deployments of a remotely operated vehicle (ROV), the optical zooplankton sensor moored on-sight key species investigation (MOKI), and Tucker trawl sampling. ROV data from the fjord entrance and the inner fjord showed high near-bottom abundances of euphausiids with a mean concentration of 17.3 ± 3.5 n/100 m**3. With the MOKI system, we observed varying numbers of euphausiids, amphipods, chaetognaths, and copepods on the seafloor at six stations. Light-induced zooplankton swarms reached densities in the order of 90,000 (euphausiids), 120,000 (amphipods), and 470,000 ind/m**3 (chaetognaths), whereas older copepodids of Calanus hyperboreus and C. glacialis did not respond to light. They were abundant at the seafloor and 5 m above and showed maximum abundance of 65,000 ind/m**3. Tucker trawl data provided an overview of the seasonal vertical distribution of euphausiids. The most abundant species Thysanoessa inermis reached near-bottom concentrations of 270 ind/m**3. Regional distribution was neither related to depth nor to location in the fjord. The taxa observed were all part of the pelagic community. Our observations suggest the presence of near-bottom macrozooplankton also in other regions and challenge the current view of bentho-pelagic coupling. Neglecting this community may cause severe underestimates of the stock of elagic zooplankton, especially predatory species, which link secondary production with higher trophic levels.

Modelled global ocean temperature from 6 ka to present

We compare the ocean temperature evolution of the Holocene as simulated by climate models and reconstructed from marine temperature proxies. This site provides informations about the Holocene temperature trends as simulated by the models. We use transient simulations from a coupled atmosphere-ocean general circulation model, as well as an ensemble of time slice simulations from the Paleoclimate Modelling Intercomparison Project. The general pattern of sea surface temperature (SST) in the models shows a high latitude cooling and a low latitude warming. The proxy dataset comprises a global compilation of marine alkenone- and Mg/Ca-derived SST estimates. Independently of the choice of the climate model, we observe significant mismatches between modelled and estimated SST amplitudes in the trends for the last 6000 years. Alkenone-based SST records show a similar pattern as the simulated annual mean SSTs, but the simulated SST trends underestimate the alkenone-based SST trends by a factor of two to five. For Mg/Ca, no significant relationship between model simulations and proxy reconstructions can be detected. We tested if such discrepancies can be caused by too simplistic interpretations of the proxy data. We tested different seasons and depths in the model to compare the proxy data trends, and can reconcile only part of the mismatches on a regional scale. We therefore considered the additional environmental factor changes in the planktonic organisms' habitat depth and a time-shift in the recording season to diagnose whether invoking those environmental factors can help reconciling the proxy records and the model simulations. We find that invoking shifts in the living season and habitat depth can remove some of the model-data discrepancies in SST trends. Regardless whether such adjustments in the environmental parameters during the Holocene are realistic, they indicate that when modeled temperature trends are set up to allow drastic shifts in the ecological behavior of planktonic organisms, they do not capture the full range of reconstructed SST trends. Our findings indicate that climate model and reconstructed temperature trends are to a large degree only qualitatively comparable, thus providing a challenge for the interpretation of proxy data as well as the models' sensitivity to orbital forcing.