FESOM Arctic Ocean sea ice concentration and thickness 1995-2014, links to movies in mp4 format

Northern Hemisphere sea ice from a Finite-Element Sea-Ice Ocean Model (FESOM) 4.5 km resolution simulation carried out by researchers from Alfred Wegener Institute (AWI), Germany. Concentration is shown with color; thickness is shown with shading. A global 1 degree mesh is used, with the "Arctic Ocean" locally refined to 4.5 km. South of CAA and Fram Strait the resolution is not refined in this simulation. The animation indicates that the 4.5 km model resolution helps to represent the small scale sea ice features, although much higher resolution is required to fully resolve the ice leads. The animation is created by Michael Böttinger from DKRZ (https://www.dkrz.de).

Discarge compilation from The Global River Discharge (RivDIS) Project

The Global River Discharge (RivDIS) data set contains monthly discharge measurements for 1018 stations located throughout the world. The period of record varies widely from station to station, with a mean of 21.5 years. These data were digitized from published UNESCO archives by Charles Voromarty, Balaze Fekete, and B.A. Tucker of the Complex Systems Research Center (CSRC) at the University of New Hampshire. River discharge is typically measured through the use of a rating curve that relates local water level height to discharge. This rating curve is used to estimate discharge from the observed water level. The rating curves are periodically rechecked and recalibrated through on-site measurement of discharge and river stage.

Land Surface Temperature continuous time series from SSM/I passive microwave radiometers, over northern latitudes (Latitudes > 45°N), during summer snow free period (2000-2011), with links to NetCDF files

Remote sensing instruments are key players to map land surface temperature (LST) at large temporal and spatial scales. In this paper, we present how we combine passive microwave and thermal infrared data to estimate LST during summer snow-free periods over northern high latitudes. The methodology is based on the SSM/I-SSMIS 37 GHz measurements at both vertical and horizontal polarizations on a 25 km × 25 km grid size. LST is retrieved from brightness temperatures introducing an empirical linear relationship between emissivities at both polarizations as described in Royer and Poirier (2010). This relationship is calibrated at pixel scale, using cloud-free independent LST data from MODIS instruments. The SSM/I-SSMIS and MODIS data are synchronized by fitting a diurnal cycle model built on skin temperature reanalysis provided by the European Centre for Medium-Range Weather Forecasts (ECMWF). The resulting temperature dataset is provided at 25 km scale and at an hourly time step during the ten-year analysis period (2000-2011). This new product was locally evaluated at five experimental sites of the EU-PAGE21 project against air temperature measurements and meteorological model reanalysis, and compared to the MODIS LST product at both local and circumpolar scale. The results giving a mean RMSE of the order of 2.2 K demonstrate the usefulness of the microwave product, which is unaffected by clouds as opposed to thermal infrared products and offers a better resolution compared to model reanalysis.

Downscaled land surface temperatures obtained by a particle smoothing approach, with links to NetCDF files

This dataset contains continuous time series of land surface temperature (LST) at spatial resolution of 300m around the 12 experimental sites of the PAGE21 project (grant agreement number 282700, funded by the EC seventh Framework Program theme FP7-ENV-2011). This dataset was produced from hourly LST time series at 25km scale, retrieved from SSM/I data (André et al., 2015, doi:10.1016/j.rse.2015.01.028) and downscaled to 300m using a dynamic model and a particle smoothing approach. This methodology is based on two main assumptions. First, LST spatial variability is mostly explained by land cover and soil hydric state. Second, LST is unique for a land cover class within the low resolution pixel. Given these hypotheses, this variable can be estimated using a land cover map and a physically based land surface model constrained with observations using a data assimilation process. This methodology described in Mechri et al. (2014, doi:10.1002/2013JD020354) was applied to the ORCHIDEE land surface model (Krinner et al., 2005, doi:10.1029/2003GB002199) to estimate prior values of each land cover class provided by the ESA CCI-Land Cover product (Bontemps et al., 2013) at 300m resolution . The assimilation process (particle smoother) consists in simulating ensemble of LST time series for each land cover class and for a large number of parameter sets. For each parameter set, the resulting temperatures are aggregated considering the grid fraction of each land cover and compared to the coarse observations. Miniminizing the distance between the aggregated model solutions and the observations allow us to select the simulated LST and the corresponding parameter sets which fit the observations most closely. The retained parameter sets are then duplicated and randomly perturbed before simulating the next time window. At the end, the most likely LST of each land cover class are estimated and used to reconstruct LST maps at 300m resolution using ESA CCI-Land Cover. The resulting temperature maps on which ice pixels were masked, are provided at daily time step during the nine-year analysis period (2000-2009).

On the Southern Ocean CO2 uptake and the role of the biological carbon pump in the 21st century, links to supplementary material

We use a suite of eight ocean biogeochemical/ecological general circulation models from the MAREMIP and CMIP5 archives to explore the relative roles of changes in winds (positive trend of Southern Annular Mode, SAM) and in warming- and freshening-driven trends of upper ocean stratification in altering export production and CO2 uptake in the Southern Ocean at the end of the 21st century. The investigated models simulate a broad range of responses to climate change, with no agreement ona dominance of either the SAM or the warming signal south of 44° S. In the southernmost zone, i.e., south of 58° S, they concur on an increase of biological export production, while between 44 and 58° S the models lack consensus on the sign of change in export. Yet, in both regions, the models show an enhanced CO2 uptake during spring and summer. This is due to a larger CO 2 (aq) drawdown by the same amount of summer export production at a higher Revelle factor at the end of the 21st century. This strongly increases the importance of the biological carbon pump in the entire Southern Ocean. In the temperate zone, between 30 and 44° S all models show a predominance of the warming signal and a nutrient-driven reduction of export production. As a consequence, the share of the regions south of 44° S to the total uptake of the Southern Ocean south of 30° S is projected to increase at the end of the 21st century from 47 to 66% with a commensurable decrease to the north. Despite this major reorganization of the meridional distribution of the major regions of uptake, the total uptake increases largely in line with the rising atmospheric CO2. Simulations with the MITgcm-REcoM2 model show that this is mostly driven by the strong increase of atmospheric CO2, with the climate-driven changes of natural CO2 exchange offsetting that trend only to a limited degree (~10%) and with negligible impact of climate effects on anthropogenic CO2 uptake when integrated over a full annual cycle south of 30° S.

Circumpolar Landscape Units, links to GeoTIFFs

Requirements for space based monitoring of permafrost features had been already defined within the IGOS Cryosphere Theme Report at the start of the IPY in 2007 (IGOS, 2007). The WMO Polar Space Task Group (PSTG, http://www.wmo.int/pages/prog/sat/pstg_en.php) identified the need to review the requirements for permafrost monitoring and to update these requirements in 2013. Relevant surveys with focus on satellite data are already available from the ESA DUE Permafrost User requirements survey (2009), the United States National Research Council (2014) and the ESA - CliC - IPA - GTN -P workshop in February 2014. These reports have been reviewed and specific needs discussed within the community and a white paper submitted to the WMO PSTG. Acquisition requirements for monitoring of especially terrain changes (incl. rock glaciers and coastal erosion) and lakes (extent, ice properties etc.) with respect to current satellite missions have been specified. About 50 locations ('cold spots') where permafrost (Arctic and Antarctic) in situ monitoring has been taking place for many years or where field stations are currently established have been identified. These sites have been proposed to the WMO Polar Space Task Group as focus areas for future monitoring by high resolution satellite data. The specifications of these sites including meta-data on site instrumentation have been published as supplement to the white paper (Bartsch et al. 2014, doi:10.1594/PANGAEA.847003). The representativity of the 'cold spots' around the arctic has been in the following assessed based on a landscape units product which has been developed as part of the FP7 project PAGE21. The ESA DUE Permafrost service has been utilized to produce a pan-arctic database (25km, 2000-2014) comprising Mean Annual Surface Temperature, Annual and summer Amplitude of Surface Temperature, Mean Summer (July-August) Surface Temperature. Surface status (frozen/unfrozen) related products have been also derived from the ESA DUE Permafrost service. This includes the length of unfrozen period, first unfrozen day and first frozen day. In addition, SAR (ENVISAT ASAR GM) statistics as well as topographic parameters have been considered. The circumpolar datasets have been assessed for their redundancy in information content. 12 distinct units could be derived. The landscape units reveal similarities between North Slope Alaska and the region from the Yamal Peninsula to the Yenisei estuary. Northern Canada is characterized by the same landscape units like western Siberia. North-eastern Canada shows similarities to the Laptev coast region. This paper presents the result of this assessment and formulates recommendations for extensions of the in situ monitoring networks and categorizes the sites by satellite data requirements (specifically Sentinels) with respect to the landscape type and related processes.

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.

Projections for Land Use Change in Brazil:2000-2050, links to files in ArcGIS shapefile format

This data set contains the inputs and the results of the REDD+ Policy Assessment Centre project (REDD-PAC) project (http://www.redd-pac.org), developed by a consortium of research institutes (IIASA, INPE, IPEA, UNEP-WCMC), supported by Germany's International Climate Initiative. Taking a new land use map of Brazil for 2000 as input, the research team used the global economic model GLOBIOM to project land use changes in Brazil up to 2050. Model projections show that Brazil has the potential to balance its goals of protecting the environment and becoming a major global producer of food and biofuels. The model results were taken into account by Brazilian decision-makers when developing the country's intended nationally determined contribution (INDC).

Total sediment thickness grid of the Southern Pacific Ocean off West Antarctica, links to NetCDF files

Paleotopographic models of the West Antarctic margin, which are essential for robust simulations of paleoclimate scenarios, lack information on sediment thickness and geodynamic conditions, resulting in large uncertainties. A new total sediment thickness grid spanning the Ross Sea-Amundsen Sea-Bellingshausen Sea basins is presented and is based on all the available seismic reflection, borehole, and gravity modeling data offshore West Antarctica. This grid was combined with NGDC's global 5 arc minute grid of ocean sediment thickness (Whittaker et al., 2013, doi:10.1002/ggge.20181) and extends the NGDC grid further to the south. Sediment thickness along the West Antarctic margin tends to be 3-4 km larger than previously assumed. The sediment volume in the Bellingshausen, Amundsen, and Ross Sea basins amounts to 3.61, 3.58, and 2.78 million km³, respectively. The residual basement topography of the South Pacific has been revised and the new data show an asymmetric trend over the Pacific-Antarctic Ridge. Values are anomalously high south of the spreading ridge and in the Ross Sea area, where the topography seems to be affected by persistent mantle processes. In contrast, the basement topography offshore Marie Byrd Land cannot be attributed to dynamic topography, but rather to crustal thickening due to intraplate volcanism. Present-day dynamic topography models disagree with the presented revised basement topography of the South Pacific, rendering paleotopographic reconstructions with such a limited dataset still fairly uncertain.

Size-partitioned phytoplankton carbon concentrations retrieved from ocean color data, links to data in netCDF format

Owing to their important roles in biogeochemical cycles, phytoplankton functional types (PFTs) have been the aim of an increasing number of ocean color algorithms. Yet, none of the existing methods are based on phytoplankton carbon (C) biomass, which is a fundamental biogeochemical and ecological variable and the "unit of accounting" in Earth system models. We present a novel bio-optical algorithm to retrieve size-partitioned phytoplankton carbon from ocean color satellite data. The algorithm is based on existing methods to estimate particle volume from a power-law particle size distribution (PSD). Volume is converted to carbon concentrations using a compilation of allometric relationships. We quantify absolute and fractional biomass in three PFTs based on size - picophytoplankton (0.5-2 µm in diameter), nanophytoplankton (2-20 µm) and microphytoplankton (20-50 µm). The mean spatial distributions of total phytoplankton C biomass and individual PFTs, derived from global SeaWiFS monthly ocean color data, are consistent with current understanding of oceanic ecosystems, i.e., oligotrophic regions are characterized by low biomass and dominance of picoplankton, whereas eutrophic regions have high biomass to which nanoplankton and microplankton contribute relatively larger fractions. Global climatological, spatially integrated phytoplankton carbon biomass standing stock estimates using our PSD-based approach yield - 0.25 Gt of C, consistent with analogous estimates from two other ocean color algorithms and several state-of-the-art Earth system models. Satisfactory in situ closure observed between PSD and POC measurements lends support to the theoretical basis of the PSD-based algorithm. Uncertainty budget analyses indicate that absolute carbon concentration uncertainties are driven by the PSD parameter No which determines particle number concentration to first order, while uncertainties in PFTs' fractional contributions to total C biomass are mostly due to the allometric coefficients. The C algorithm presented here, which is not empirically constrained a priori, partitions biomass in size classes and introduces improvement over the assumptions of the other approaches. However, the range of phytoplankton C biomass spatial variability globally is larger than estimated by any other models considered here, which suggests an empirical correction to the No parameter is needed, based on PSD validation statistics. These corrected absolute carbon biomass concentrations validate well against in situ POC observations.