Phytoplakton pigment concentrations during POLARSTERN cruise ANT-XXVIII/3 and SONNE cruise SO218

Diatoms are the major marine primary producers on the global scale and, recently, several methods have been developed to retrieve their abundance or dominance from satellite remote sensing data. In this work, we highlight the importance of the Southern Ocean (SO) in developing a global algorithm for diatom using an Abundance Based Approach (ABA). A large global in situ data set of phytoplankton pigments was compiled, particularly with more samples collected in the SO. We revised the ABA to take account of the information on the penetration depth (Zpd) and to improve the relationship between diatoms and total chlorophyll-a (TChla). The results showed that there is a distinct relationship between diatoms and TChla in the SO, and a new global model (ABAZpd) improved the estimation of diatoms abundance by 28% in the SO compared with the original ABA model. In addition, we developed a regional model for the SO which further improved the retrieval of diatoms by 17% compared with the global ABAZpd model. As a result, we found that diatom may be more abundant in the SO than previously thought. Linear trend analysis of diatom abundance using the regional model for the SO showed that there are statistically significant trends, both increasing and decreasing, in diatom abundance over the past eleven years in the region.

Tab. 1: Ar40/Ar39 results from Kuyul terrane, northeastern Russia

The Kamchatka Peninsula of northeastern Russia is located along the northwestern margin of the Bering Sea and consists of zones of complexly deformed accreted terranes. Along the northern portion of the peninsula, progressing from then orthwestem Bering Sea inland the Olyutorskiy, Ukelayat, and Koryak superterranes area acreted to the Okhotsk-Chukotsk volcanic-plutonic bell in northern-most Kamchatka. A sedimentary sequence of Albian to Maastrichtian age overlap terranes and units of the Koryak superterrane and constrains their accretion time with this region of the North America plate. Ophiolite complexes, widespread within the Koryak superterrane, are associated with serpentinite melanges and some of the ophiolite terranes include large portions of weakly serpentinized hyperbasites, layered gabbro, sheeted dikes, and pillow basalts outcropping as internally coherent blocks within a sheared melange matrix. Interpretation of magnetic anomalies allow the correlation of the Ukelayat with the West Kamchatka and Sredinny Range superterranes. The Olyutorskiy composite terrane may be correlated with the central and southern Kamchatka Peninsula Litke, Eastern Ranges and Vetlov composite terranes. The most "out-board" of the central and southern Kamchatka Peninsula terranes is the Kronotsky composite terrane, weil exposed along the Kamchatka, Kronotsky and Shipunsky Capes. Using regional geological constraints, paleomagnetism, and plate kinematic models for the Pacific basin a regional model can be proposed in which accretion of the Koryak composite terrane to the North America plate occurs during the Campanian-Maastrichtian, followed by the accretion of the Olyutorskiy composite terrane in the Middle Eocene, and the Late Oligocene-Early Miocene collision of the Kronotsky composite terrane. A revised age estimate of a key overlapping sedirnentary sequence of the Koryak superterrane, calibrated with new Ar40/Ar39 data, supports its Late Cretaceous accretion age.

Secchi depth maps for the Baltic Sea, 2003-2012, with links to GeoTIFs

Secchi depth is a measure of water transparency. In the Baltic Sea region, Secchi depth maps are used to assess eutrophication and as input for habitat models. Due to their spatial and temporal coverage, satellite data would be the most suitable data source for such maps. But the Baltic Sea's optical properties are so different from the open ocean that globally calibrated standard models suffer from large errors. Regional predictive models that take the Baltic Sea's special optical properties into account are thus needed. This paper tests how accurately generalized linear models (GLMs) and generalized additive models (GAMs) with MODIS/Aqua and auxiliary data as inputs can predict Secchi depth at a regional scale. It uses cross-validation to test the prediction accuracy of hundreds of GAMs and GLMs with up to 5 input variables. A GAM with 3 input variables (chlorophyll a, remote sensing reflectance at 678 nm, and long-term mean salinity) made the most accurate predictions. Tested against field observations not used for model selection and calibration, the best model's mean absolute error (MAE) for daily predictions was 1.07 m (22%), more than 50% lower than for other publicly available Baltic Sea Secchi depth maps. The MAE for predicting monthly averages was 0.86 m (15%). Thus, the proposed model selection process was able to find a regional model with good prediction accuracy. It could be useful to find predictive models for environmental variables other than Secchi depth, using data from other satellite sensors, and for other regions where non-standard remote sensing models are needed for prediction and mapping. Annual and monthly mean Secchi depth maps for 2003-2012 come with this paper as Supplementary materials.

Marine Isotope Stage (MIS) 11 results of model data with Community Climate System Model version 3 including the Dynamic Global Vegetation Model (CCSM3-DGVM) in NetCDF format at global and regional scale

The climate of Marine Isotope Stage (MIS) 11, the interglacial roughly 400,000 years ago, is investigated for four time slices, 416, 410, 400, and 394 ka. The overall picture is that MIS 11 was a relatively warm interglacial in comparison to preindustrial, with Northern Hemisphere (NH) summer temperatures early in MIS 11 (416-410 ka) warmer than preindustrial, though winters were cooler. Later in MIS 11, especially around 400 ka, conditions were cooler in the NH summer, mainly in the high latitudes. Climate changes simulated by the models were mainly driven by insolation changes, with the exception of two local feedbacks that amplify climate changes. Here, the NH high latitudes, where reductions in sea ice cover lead to a winter warming early in MIS 11, as well as the tropics, where monsoon changes lead to stronger climate variations than one would expect on the basis of latitudinal mean insolation change alone, are especially prominent. The results support a northward expansion of trees at the expense of grasses in the high northern latitudes early during MIS 11, especially in northern Asia and North America.

Regional geoid model for the area of the subglacial Lake Vostok

We present a regional geoid model for the area of Lake Vostok, Antarctica, from a combination of local airborne gravity, ice-surface and ice-thickness data and a lake bathymetry model. The topography data are used for residual terrain modelling (RTM) in a remove-compute-restore approach together with the GOCE satellite model GOCO03S. The disturbing potential at the Earth's surface, i.e. the quasigeoid, is predicted by least-squares collocation (LSC) and subsequently converted to geoid heights. Compared to GOCO03S our regional solution provides an additional short-wavelength signal of up to 1.48 m, or 0.56 m standard deviation, respectively. More details can be found in Schwabe et. al (2014).

Regional ocean model for the Mediterranean, links to model files

A regional ocean general circulation model of the Mediterranean is used to study the climate of the Last Glacial Maximum. The atmospheric forcing for these simulations has been derived from simulations with an atmospheric general circulation model, which in turn was forced with surface conditions from a coarse resolution earth system model. The model is successful in reproducing the general patterns of reconstructed sea surface temperature anomalies with the strongest cooling in summer in the northwestern Mediterranean and weak cooling in the Levantine, although the model underestimates the extent of the summer cooling in the western Mediterranean. However, there is a strong vertical gradient associated with this pattern of summer cooling, which makes the comparison with reconstructions complicated. The exchange with the Atlantic is decreased to roughly one half of its present value, which can be explained by the shallower Strait of Gibraltar as a consequence of lower global sea level. This reduced exchange causes a strong increase of salinity in the Mediterranean in spite of reduced net evaporation.

ETOPO1 Global Relief Model converted to PanMap layer format

ETOPO1 is a 1 arc-minute global relief model of Earth's surface that integrates land topography and ocean bathymetry. It was built from numerous global and regional data sets. Data were converted to the PanMap layer format in 14 contour lines from 500 to 7000 meter in steps of 500 m. The link provides a zip-archive (1.1 MB) with *.lay files. The PanMap Mini-GIS software is published at doi:10.1594/PANGAEA.104840.

Three time-slice simulations of the Eemian (eem), the mid-Holocene (hol), and a pre-industrial (pri) control run with the coupled atmosphere-ocean-model ECHAM5/MPI-OM

Orbital forcing does not only exert direct insolation effects, but also alters climate indirectly through feedback mechanisms that modify atmosphere and ocean dynamics and meridional heat and moisture transfers. We investigate the regional effects of these changes by detailed analysis of atmosphere and ocean circulation and heat transports in a coupled atmosphere-ocean-sea ice-biosphere general circulation model (ECHAM5/JSBACH/MPI-OM). We perform long term quasi equilibrium simulations under pre-industrial, mid-Holocene (6000 years before present - yBP), and Eemian (125 000 yBP) orbital boundary conditions. Compared to pre-industrial climate, Eemian and Holocene temperatures show generally warmer conditions at higher and cooler conditions at lower latitudes. Changes in sea-ice cover, ocean heat transports, and atmospheric circulation patterns lead to pronounced regional heterogeneity. Over Europe, the warming is most pronounced over the north-eastern part in accordance with recent reconstructions for the Holocene. We attribute this warming to enhanced ocean circulation in the Nordic Seas and enhanced ocean-atmosphere heat flux over the Barents Shelf in conduction with retreat of sea ice and intensified winter storm tracks over northern Europe.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Trace metal concentration and age model of ODP Hole 165-1002C sediments

Titanium and iron concentration data from the anoxic Cariaco Basin, off the Venezuelan coast, can be used to infer variations in the hydrological cycle over northern South America during the past 14,000 years with subdecadal resolution. Following a dry Younger Dryas, a period of increased precipitation and riverine discharge occurred during the Holocene 'thermal maximum'. Since ~5400 years ago, a trend toward drier conditions is evident from the data, with high-amplitude fluctuations and precipitation minima during the time interval 3800 to 2800 years ago and during the 'Little Ice Age'. These regional changes in precipitation are best explained by shifts in the mean latitude of the Atlantic Intertropical Convergence Zone (ITCZ), potentially driven by Pacific-based climate variability. The Cariaco Basin record exhibits strong correlations with climate records from distant regions, including the high-latitude Northern Hemisphere, providing evidence for global teleconnections among regional climates.

Proxy data and model simulations of sediment core GeoB7920-2

The evolution of the northwest African hydrological balance throughout the Pleistocene epoch influenced the migration of prehistoric humans**1. The hydrological balance is also thought to be important to global teleconnection mechanisms during Dansgaard-Oeschger and Heinrich events**2. However, most high-resolution African climate records do not span the millennial-scale climate changes of the last glacial-interglacial cycle**1, 3, 4, 5, or lack an accurate chronology**6. Here, we use grain-size analyses of siliciclastic marine sediments from off the coast of Mauritania to reconstruct changes in northwest African humidity over the past 120,000 years. We compare this reconstruction to simulations of palaeo-humidity from a coupled atmosphere-ocean-vegetation model. These records are in good agreement, and indicate the reoccurrence of precession-forced humid periods during the last interglacial period similar to the Holocene African Humid Period. We suggest that millennial-scale arid events are associated with a reduction of the North Atlantic meridional overturning circulation and that millennial-scale humid events are linked to a regional increase of winter rainfall over the coastal regions of northwest Africa.

World distribution of land cover changes, model in NetCDF format

The role of Pre- and Protohistoric anthropogenic land cover changes needs to be quantified i) to establish a baseline for comparison with current human impact on the environment and ii) to separate it from naturally occurring changes in our environment. Results are presented from the simple, adaptation-driven, spatially explicit Global Land Use and technological Evolution Simulator (GLUES) for pre-Bronze age demographic, technological and economic change. Using scaling parameters from the History Database of the Global Environment as well as GLUES-simulated population density and subsistence style, the land requirement for growing crops is estimated. The intrusion of cropland into potentially forested areas is translated into carbon loss due to deforestation with the dynamic global vegetation model VECODE. The land demand in important Prehistoric growth areas - converted from mostly forested areas - led to large-scale regional (country size) deforestation of up to 11% of the potential forest. In total, 29 Gt carbon were lost from global forests between 10 000 BC and 2000 BC and were replaced by crops; this value is consistent with other estimates of Prehistoric deforestation. The generation of realistic (agri-)cultural development trajectories at a regional resolution is a major strength of GLUES. Most of the pre-Bronze age deforestation is simulated in a broad farming belt from Central Europe via India to China. Regional carbon loss is, e.g., 5 Gt in Europe and the Mediterranean, 6 Gt on the Indian subcontinent, 18 Gt in East and Southeast Asia, or 2.3 Gt in subsaharan Africa.