The CLIMATE data base of monthly normals for globally gridded climatic variables (version 2.1)

This database (Leemans & Cramer 1991) contains monthly averages of mean temperature, temperature range, precipitation, rain days and sunshine hours for the terrestrial surface of the globe, gridded at 0.5 degree longitude/latitude resolution. All grd-files contain the same 62483 pixels in the same order, with 30' latitude and longitude resolution. The coordinates are in degree-decimals and indicate the SW corner of each pixel. Topography is from ETOPO5 and indicates modal elevation. Data were generated from a large data base, using the partial thin-plate splining algorithm (Hutchinson & Bischof 1983). This version is widely used around the globe, notably by all groups participating in the IGBP NPP model intercomparison.

Fossil pollen record of composite sediment core TMD from Tso Moriri, analysis of modern surface pollen samples, mean annual precipitation reconstruction, and digitisation and recalibration of different discussed palaeoclimate proxy records

This paper presents a new fossil pollen record from Tso Moriri (32°54'N, 78°19'E, 4512 m a.s.l.) and seeks to reconstruct changes in mean annual precipitation (MAP) during the last 12,000 years. This high-alpine lake occupies an area of 140 km**2 in a glacial-tectonic valley in the northwestern Himalaya. The region has a cold climate, with a MAP <300 mm, and open vegetation. The hydrology is controlled by the Indian Summer Monsoon (ISM), but winter westerly-associated precipitation also affects the regional water balance. Results indicate that precipitation levels varied significantly during the Holocene. After a rapid increase in MAP, a phase of maximum humidity was reached between ca. 11 to 9.6 cal ka BP, followed by a gradual decline in MAP. This trend parallels the reduction in the Northern Hemisphere summer insolation. Comparison of different palaeoclimate proxy records reveal evidence for a stronger Holocene decrease in precipitation in the northern versus the southern parts of the ISM domain. The long-term trend of ISM weakening is overlaid with several short periods of greater dryness, which are broadly synchronous with the North Atlantic cold spells, suggesting reduced amounts of westerly-associated winter precipitation. Compared to the mid and late Holocene, it appears that westerlies had a greater influence on the western parts of the ISM domain during the early Holocene. During this period, the westerly-associated summer precipitation belt was positioned at Mediterranean latitudes and amplified the ISM-derived precipitation. The Tso Moriri pollen record and moisture reconstructions also suggest that changes in climatic conditions affected the ancient Harappan Civilisation, which flourished in the greater Indus Valley from approximately 5.2 to 3 cal ka BP. The prolonged Holocene trend towards aridity, punctuated by an interval of increased dryness (between ca. 4.5 to 4.3 cal ka BP), may have pushed the Mature Harappan urban settlements (between ca. 4.5 to 3.9 cal ka BP) to develop more efficient agricultural practices to deal with the increasingly acute water shortages. The amplified aridity associated with North Atlantic cooling between ca. 4 to 3.6 and around 3.2 cal ka BP further hindered local agriculture, possibly causing the deurbanisation that occurred from ca. 3.9 cal ka BP and eventual collapse of the Harappan Civilisation between ca. 3.5 to 3 cal ka BP.

ROCADA: Romanian daily gridded climatic dataset (1961-2013) V1.0

Daily records of nine meteorological variables covering the interval 1961-2013 were used in order to create a state-of-the-art homogenized climatic dataset over Romania at a spatial resolution of 0.1°. All meteorological stations with full data records, as well as stations with up to 30 % missing data, were used for the following variables: air pressure (150 stations); minimum, maximum, and average air temperature (150 stations); soil temperature (127 stations); precipitation (188 stations); sunshine hours (135 stations); cloud cover (104 stations); relative humidity (150 stations). For each parameter, the data series were first homogenized with the software MASH (Multiple Analysis of Series for Homogenization); then, the data series were gridded by means of the software MISH (Meteorological Interpolation based on Surface Homogenized Data).

Salinity of the Eocene Arctic Ocean from oxygen isotope analysis of fish bone carbonate

Stable isotope analysis was performed on the structural carbonate of fish bone apatite from early and early middle Eocene samples (~55 to ~45 Ma) recently recovered from the Lomonosov Ridge by Integrated Ocean Drilling Program Expedition 302 (the Arctic Coring Expedition). The d18O values of the Eocene samples ranged from -6.84 per mil to -2.96 per mil Vienna Peedee belemnite, with a mean value of -4.89 per mil, compared to 2.77 per mil for a Miocene sample in the overlying section. An average salinity of 21 to 25 per mil was calculated for the Eocene Arctic, compared to 35 per mil for the Miocene, with lower salinities during the Paleocene Eocene thermal maximum, the Azolla event at ~48.7 Ma, and a third previously unidentified event at ~47.6 Ma. At the Azolla event, where the organic carbon content of the sediment reaches a maximum, a positive d13C excursion was observed, indicating unusually high productivity in the surface waters.

A global seasonal surface ocean climatology of phytoplankton types based on CHEMTAX analysis of HPLC pigments

Much advancement has been made in recent years in field data assimilation, remote sensing and ecosystem modeling, yet our global view of phytoplankton biogeography beyond chlorophyll biomass is still a cursory taxonomic picture with vast areas of the open ocean requiring field validations. High performance liquid chromatography (HPLC) pigment data combined with inverse methods offer an advantage over many other phytoplankton quantification measures by way of providing an immediate perspective of the whole phytoplankton community in a sample as a function of chlorophyll biomass. Historically, such chemotaxonomic analysis has been conducted mainly at local spatial and temporal scales in the ocean. Here, we apply a widely tested inverse approach, CHEMTAX, to a global climatology of pigment observations from HPLC. This study marks the first systematic and objective global application of CHEMTAX, yielding a seasonal climatology comprised of ~1500 1°x1° global grid points of the major phytoplankton pigment types in the ocean characterizing cyanobacteria, haptophytes, chlorophytes, cryptophytes, dinoflagellates, and diatoms, with results validated against prior regional studies where possible. Key findings from this new global view of specific phytoplankton abundances from pigments are a) the large global proportion of marine haptophytes (comprising 32 ± 5% of total chlorophyll), whose biogeochemical functional roles are relatively unknown, and b) the contrasting spatial scales of complexity in global community structure that can be explained in part by regional oceanographic conditions. These publicly accessible results will guide future parameterizations of marine ecosystem models exploring the link between phytoplankton community structure and marine biogeochemical cycles.

Grain-size analysis and geochemistry of sediment core GeoB13602-1

Grain-size, terrigenous element and rock magnetic remanence data of Quaternary marine sediments retrieved at the NW African continental margin off Gambia (gravity core GeoB 13602-1, 13°32.71' N, 17°50.96'W) were jointly analyzed by end-member (EM) unmixing methods to distinguish and budget past terrigenous fluxes. We compare and cross-validate the identified single-parameter EM systems and develop a numerical strategy to calculate associated multi-parameter EM properties. One aeolian and two fluvial EMs were found. The aeolian EM is much coarser than the fluvial EMs and is associated with a lower goethite/hematite ratio, a higher relative concentration of magnetite and lower Al/Si and Fe/K ratios. Accumulation rates and grain sizes of the fluvial sediment appear to be primarily constrained by shore distance (i.e., sea-level fluctuations) and to a lesser extent by changes in hinterland precipitation. High dust fluxes occurred during the Last Glacial Maximum (LGM) and during Heinrich Stadials (HS) while the fluvial input remained unchanged. Our approach reveals that the LGM dust fluxes were ~7 times higher than today's. However, by far the highest dust accumulation occurred during HS 1 (~300 g m**-2 yr** -1), when dust fluxes were ~80 fold higher than today. Such numbers have not yet been reported for NW Africa, and emphasize strikingly different environmental conditions during HSs. They suggest that deflation rate and areal extent of HSs dust sources were much larger due to retreating vegetation covers. Beyond its regional and temporal scope, this study develops new, in principle, generally applicable strategies for multi-method end-member interpretation, validation and flux budgeting calibration.

Spatially explicit estimates of North Atlantic albacore tuna (Thunnus alalunga) biomass in the North Atlantic for the period 1956-2010 - Gridded data product (NetCDF)

The development of the ecosystem approach and models for the management of ocean marine resources requires easy access to standard validated datasets of historical catch data for the main exploited species, together with the model estimates achieved from these data, allowing models inter-comparison and evaluation of model skills. North Atlantic albacore tuna is exploited all year round by longline and in summer and autumn by surface fisheries and fishery statistics compiled by the International Commission for the Conservation of Atlantic Tunas (ICCAT). Catch and effort with geographical coordinates at monthly spatial resolution of 1° or 5° squares were extracted for this species with a careful definition of fisheries and data screening. Length frequencies of catch were also extracted according to the definition of fisheries for the period 1956-2010. Using these data, an application of the spatial ecosystem and population dynamics model (SEAPODYM) was developed for the North Atlantic albacore population and fisheries and provided the first spatially explicit estimate of albacore density in the North Atlantic by life stage. These densities by life stage (larval recruits, young immature fish adult mature fish and total biomass) are provided in gridded file (Netcdf) at resolution of 2° x 2° x month.

Geochemical processes and rare earth element analysis in cold-seep pore waters of Hydrate Ridge, northeast Pacific Ocean

The concentrations of rare earth elements (REEs), sulphate, hydrogen sulphide, total alkalinity, calcium, magnesium and phosphate were measured in shallow (<12 cm below seafloor) pore waters from cold-seep sediments on the northern and southern summits of Hydrate Ridge, offshore Oregon. Downward-decreasing sulphate and coevally increasing sulphide concentrations reveal sulphate reductionas dominant early diagenetic process from ~2 cm depth downwards. A strong increase of total dissolved REE concentrations is evident immediately below the sediment-water interface, which can be related to early diagenetic release of REEs into pore water resulting from the remineralization of particulate organic matter. The highest pore water REE concentrations were measured close to the sediment-water interface at ~2 cm depth. Distinct shale normalized REE patterns point to particulate organic matter and iron oxides as main REE sources in the upper ~2-cm depth interval. In general, the pore waters have shalenormalized patterns reflecting heavy REE (HREE) enrichment, which suggests preferential complexation of HREEs with carbonate ions. Below ~2 cm depth, a downward decrease in REE correlates with a decrease in pore water calcium concentrations. At this depth, the anaerobic oxidation of methane (AOM) coupled to sulphate reduction increases carbonate alkalinity through the production of bicarbonate, which results in the precipitation of carbonate minerals. It seems therefore likely that the REEs and calcium are consumed during vast AOM-induced precipitation of carbonate in shallow Hydrate Ridge sediments. The analysis of pore waters from Hydrate Ridge shed new light on early diagenetic processes at cold seeps, corroborating the great potential of REEs to identify geochemical processes and to constrain environmental conditions.

Global Lithological Map Database v1.0 (gridded to 0.5° spatial resolution)

Lithology describes the geochemical, mineralogical, and physical properties of rocks. It plays a key role in many processes at the Earth surface, especially the fluxes of matter to soils, ecosystems, rivers, and oceans. Understanding these processes at the global scale requires a high resolution description of lithology. A new high resolution global lithological map (GLiM) was assembled from existing regional geological maps translated into lithological information with the help of regional literature. The GLiM represents the rock types of the Earth surface using 1,235,400 polygons. The lithological classification consists of three levels. The first level contains 16 lithological classes comparable to previously applied definitions in global lithological maps. The additional two levels contain 12 and 14 subclasses, respectively, which describe more specific rock attributes. According to the GLiM, the Earth is covered by 64 % sediments (a third of which is carbonates), 13 % metamorphics, 7 % plutonics, and 6 % volcanics, and 10% are covered by water or ice. The high resolution of the GLiM allows observation of regional lithological distributions which often vary from the global average. The GLiM enables regional analysis of Earth surface processes at global scales.

Heat index at 2 m above ground: A globally gridded dataset based on reanalysis data from 1979-2013, links to GeoTIFFs

The increase in global mean temperatures resulting from climate change has wide reaching consequences for the earth's ecosystems and other natural systems. Many studies have been devoted to evaluating the distribution and effects of these changes. We go a step further and evaluate global changes to the heat index, a measure of temperature as perceived by humans. Heat index, which is computed from temperature and relative humidity, is more important than temperature for the health of humans and other animals. Even in cases where the heat index does not reach dangerous levels from a health perspective, it has been shown to be an important factor in worker productivity and thus in economic productivity. We compute heat index from dewpoint temperature and absolute temperature 2 m above ground from the ERA-Interim reanalysis dataset for the years 1979-2013. The data is provided aggregated to daily minima, means and maxima. Furthermore, the data is temporally aggregated to monthly and yearly values and spatially aggregated to the level of countries after being weighted by population density in order to demonstrate its usefulness for the analysis of its impact on human health and productivity. The resulting data deliver insights into the spatiotemporal development of near-ground heat index during the course of the past 3 decades. It is shown that the impact of changing heat index is unevenly distributed through space and time, affecting some areas differently than others. The likelihood of dangerous heat index events has increased globally. Also, heat index climate groups that would formerly be expected closer to the tropics have spread latitudinally to include areas closer to the poles. The data can serve in future studies as a basis for evaluating and understanding the evolution of heat index in the course of climate change, as well as its impact on human health and productivity.

(Table 1) Correlation of precipitation and winter mass balance (bw) at Storbreen, Norway

Measurements of winter balance (bw) and summer balance (bs) have been carried out at Storbreen since 1949. Here we apply a simple mass balance model to study the climate sensitivity and to reconstruct the mass balance series prior to 1949. The model is calibrated and validated with data from an automatic weather station (AWS) operating in the ablation zone of Storbreen since 2001. Regression analysis revealed that bw was best modelled using precipitation data southwest of the glacier. Results from the model compared well with reported mass balance values for the period 1949-2006, obtained correlations (r) for bw and bs varied between 0.83 and 0.87 depending on model set up. Reconstruction of the mass balance series for the period 1924/1925-1948/1949 suggested a cumulative mass deficit of c. 30 m w.e. mainly due to highly negative summer balances, but also lower bw than the average for 1949-2006. Calculated change in specific mass balance for a ±1°C change in air temperature was ±0.55 m w.e., whereas a ±10 % increase in precipitation represented a change of ± 0.20 m w.e. Model results further indicated that for a 2°C warming, the ablation season will be extended by c. 30 days and that the period of ice melt at the AWS location will increase from c. 40 to c. 80 days.

Solid phase analysis of sediment core GeoB9309-1 from the Zambezi (or Mozambique) deep-sea fan

The Zambezi deep-sea fan, the largest of its kind along the east African continental margin, is poorly studied to date, despite its potential to record marine and terrestrial climate signals in the southwest Indian Ocean. Therefore, gravity core GeoB 9309-1, retrieved from 1219 m water depth, was investigated for various geophysical (magnetic susceptibility, porosity, colour reflectance) and geochemical (pore water and sediment geochemistry, Fe and P speciation) properties. Onboard and onshore data documented a sulphate/methane transition (SMT) zone at ~ 450-530 cm sediment depth, where the simultaneous consumption of pore water sulphate and methane liberates hydrogen sulphide and bi-carbonate into the pore space. This leads to characteristic changes in the sediment and pore water chemistry, as the reduction of primary Fe (oxyhydr)oxides, the precipitation of Fe sulphides, and the mobilization of Fe (oxyhydr)oxide-bound P. These chemical processes also lead to a marked decrease in magnetic susceptibility. Below the SMT, we find a reduction of porosity, possibly due to pore space cementation by authigenic minerals. Formation of the observed geochemical, magnetic and mineralogical patterns requires a fixation of the SMT at this distinct sediment depth for a considerable time-which we calculated to be ~ 10 000 years assuming steady-state conditions-following a period of rapid upward migration towards this interval. We postulate that the worldwide sea-level rise at the last glacial/interglacial transition (~ 10 000 years B.P.) most probably caused the fixation of the SMT at its present position, through drastically reduced sediment delivery to the deep-sea fan. In addition, we report an internal redistribution of P occurring around the SMT, closely linked to the (de)coupling of sedimentary Fe and P, and leaving a characteristic pattern in the solid P record. By phosphate re-adsorption onto Fe (oxyhydr)oxides above, and formation of authigenic P minerals (e.g. vivianite) below the SMT, deep-sea fan deposits may potentially act as long-term sinks for P.

Pollen analysis and age model of sediment core GIK-15856-2

Pollen and spores from a deep-sea core located west of the Niger Delta record an uninterrupted area of lowland rain forest in West Africa from Guinea to Cameroon during the last Interglacial and the early Holocene. During other periods of the last 150 ka, a savanna corridor between the western - Guinean - and the eastern - Congolian - part of the African lowland rain forest existed. This so-called Dahomey Gap had its largest extension during Glacial Stages 6, 4, 3, and 2. Reduced surface salinity in the eastern Gulf of Guinea as recorded by dinoflagellate cysts indicates sufficient precipitation for extensive forest growth during Stages 5 and 1. The large modern extension of dry forest and savanna in West Africa cannot be solely explained by climatic factors. Mangrove expansion in and west of the Niger Delta was largest during the phases of sea-level rise of Stages 5 and 1. During Stages 6, 4, 3, and 2, shelf areas were exposed and the area of the mangrove swamps was minimal.

Stable isotope analysis and biomarkers of carbonates from the Columbia River in southwestern Washington, USA

Exotic limestone masses with silicified fossils, enclosed within deep-water marine siliciclastic sediments of the Early to Middle Miocene Astoria Formation, are exposed along the north shore of the Columbia River in southwestern Washington, USA. Samples from four localities were studied to clarify the origin and diagenesis of these limestone deposits. The bioturbated and reworked limestones contain a faunal assemblage resembling that of modern and Cenozoic deep-water methane-seeps. Five phases make up the paragenetic sequence: (1) micrite and microspar; (2) fibrous, banded and botryoidal aragonite cement, partially replaced by silica or recrystallized to calcite; (3) yellow calcite; (4) quartz replacing carbonate phases and quartz cement; and (5) equant calcite spar and pseudospar. Layers of pyrite frequently separate different carbonate phases and generations, indicating periods of corrosion. Negative d13Ccarbonate values as low as -37.6 per mill V-PDB reveal an uptake of methane-derived carbon. In other cases, d13Ccarbonate values as high as 7.1 per mill point to a residual, 13C-enriched carbon pool affected by methanogenesis. Lipid biomarkers include 13C-depleted, archaeal 2,6,10,15,19-pentamethylicosane (PMI; d13C: -128 per mill), crocetane and phytane, as well as various iso- and anteiso-carbon chains, most likely derived from sulphate-reducing bacteria. The biomarker inventory proves that the majority of the carbonates formed as a consequence of sulphate-dependent anaerobic oxidation of methane. Silicification of fossils and early diagenetic carbonate cements as well as the precipitation of quartz cement - also observed in other methane-seep limestones enclosed in sediments with abundant diatoms or radiolarians - is a consequence of a preceding increase of alkalinity due to anaerobic oxidation of methane, inducing the dissolution of silica skeletons. Once anaerobic oxidation of methane has ceased, the pH drops again and silica phases can precipitate.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.