ALBERT

Description: The climate active trace-gas carbonyl sulfide (OCS) is the most abundant sulfur gas in the atmosphere. A missing source in its atmospheric budget is currently suggested, resulting from an upward revision of the vegetation sink. Tropical oceanic emissions have been proposed to close the resulting gap in the atmospheric budget. We present a bottom-up approach including (i) new observations of OCS in surface waters of the tropical Atlantic, Pacific and Indian oceans and (ii) a further improved global box model to show that direct OCS emissions are unlikely to account for the missing source. The box model suggests an undersaturation of the surface water with respect to OCS integrated over the entire tropical ocean area and, further, global annual direct emissions of OCS well below that suggested by top-down estimates. In addition, we discuss the potential of indirect emission from CS2 and dimethylsulfide (DMS) to account for the gap in the atmospheric budget. This bottom-up estimate of oceanic emissions has implications for using OCS as a proxy for global terrestrial CO2 uptake, which is currently impeded by the inadequate quantification of atmospheric OCS sources and sinks.

Description: Much of our knowledge about future changes in precipitation relies on global (GCM) and/or regional climate models (RCM) that have resolutions which are much coarser than typical spatial scales of precipitation, particularly extremes. The major problems with these projections are both climate model biases and the gap between gridbox and point scale. Wong et al. developed a model to jointly bias correct and downscale precipitation at daily scales. This approach, however, relied on pairwise correspondence between predictor and predictand for calibration, and thus, on nudged simulations which are rarely available. Here we present an extension of this approach that separates the downscaling from the bias correction and in principle is applicable to free running GCMs/RCMs. In a first step, we bias correct RCM-simulated precipitation against gridded observations at the same scale using a parametric quantile mapping approach. To correct the whole distribution including extreme tails we apply a mixture distribution of a gamma distribution for the precipitation mass and a generalized Pareto distribution for the extreme tail. In a second step, we bridge the scale gap: we predict local variance employing a vector generalized linear gamma model (VGLM gamma) with the bias corrected time series as predictor. The VGLM gamma model is calibrated between gridded and point scale (station) observations. For evaluation we adopt the perfect predictor experimental setup of VALUE. Precipitation is in most cases improved by (parts of) our method across different European climates. The method generally performs better in summer than in winter and in winter best in the Mediterranean region with a mild winter climate and worst for continental winter climate in mid & eastern Europe or Scandinavia. A strength of this two-step method is that the best combination of bias correction and downscaling methods can be selected. This implies that the concept can be extended to a wide range of method combinations.

Description: Estuary-type circulation is a residual circulation in coastal systems with horizontal density gradients. It drives the accumulation of suspended particulate matter in coastal embayments where density gradients are sustained by some freshwater inflow from rivers. Ebenhöh et al. (Ecol Model 174(3):241–252, 2004) found that shallow water depth can explain nutrient gradients becoming established towards the coast even in the absence of river inflow. The present study follows their concept and investigates the characteristic transport of organic matter towards the coast based on idealised scenarios whereby an estuary-type circulation is maintained by surface freshwater fluxes and pronounced shoaling towards the coast. A coupled hydrodynamical and biogeochemical model is used to simulate the dynamics of nutrient gradients and to derive budgets of organic matter flux for a coastal transect. Horizontal nutrient gradients are considered only in terms of tidal asymmetries of suspended matter transport. The results show that the accumulation of organic matter near the coast is not only highly sensitive to variations in the sinking velocity of suspended matter but is also noticeably enhanced by an increase in precipitation. This scenario is comparable with North Sea conditions. By contrast, horizontal nutrient gradients would be reversed in the case of evaporation-dominated inverse estuaries (cf. reverse gradients of nutrient and organic matter concentrations). Credible coastal nutrient budget calculations are required for resolving trends in eutrophication. For tidal systems, the present results suggest that these calculations require an explicit consideration of freshwater flux and asymmetries in tidal mixing. In the present case, the nutrient budget for the vertically mixed zone also indicates carbon pumping from the shelf sea towards the coast from as far offshore as 25 km.

Description: To describe the underlying processes involved in oceanic plankton dynamics is crucial for the determination of energy and mass flux through an ecosystem and for the estimation of biogeochemical element cycling. Many planktonic ecosystem models were developed to resolve major processes so that flux estimates can be derived from numerical simulations. These results depend on the type and number of parameterizations incorporated as model equations. Furthermore, the values assigned to respective parameters specify a model's solution. Representative model results are those that can explain data; therefore, data assimilation methods are utilized to yield optimal estimates of parameter values while fitting model results to match data. Central difficulties are (1) planktonic ecosystem models are imperfect and (2) data are often too sparse to constrain all model parameters. In this review we explore how problems in parameter identification are approached in marine planktonic ecosystem modelling. We provide background information about model uncertainties and estimation methods, and how these are considered for assessing misfits between observations and model results. We explain differences in evaluating uncertainties in parameter estimation, thereby also discussing issues of parameter identifiability. Aspects of model complexity are addressed and we describe how results from cross-validation studies provide much insight in this respect. Moreover, approaches are discussed that consider time- and space-dependent parameter values. We further discuss the use of dynamical/statistical emulator approaches, and we elucidate issues of parameter identification in global biogeochemical models. Our review discloses many facets of parameter identification, as we found many commonalities between the objectives of different approaches, but scientific insight differed between studies. To learn more from results of planktonic ecosystem models we recommend finding a good balance in the level of sophistication between mechanistic modelling and statistical data assimilation treatment for parameter estimation

Description: The North Atlantic Oscillation (NAO) is the dominant mode of winter climate variability in the North Atlantic sector. The corresponding index varies on a wide range of timescales, from days and months to decades and beyond. Sub-decadal NAO variability has been well documented, but the underlying mechanism is still under discussion. Other indices of North Atlantic sector climate variability such as indices of sea surface and surface air temperature or Arctic sea ice extent also exhibit pronounced sub-decadal variability. Here, we use sea surface temperature and sea level pressure observations, and the Kiel Climate Model to investigate the dynamics of the sub-decadal NAO variability. The sub-decadal NAO variability is suggested to originate from dynamical large-scale air-sea interactions. The adjustment of the Atlantic Meridional Overturning Circulation to previous surface heat flux variability provides the memory of the coupled mode. The results stress the role of coupled feedbacks in generating sub-decadal North Atlantic sector climate variability, which is important to multiyear climate predictability in that region

Description: The Southern Ocean is a major sink for anthropogenic carbon. Yet, there is no quantitative consensus about how this sink will change when surface winds increase (as they are anticipated to do). Among the tools employed to quantify carbon uptake are global coupled ocean-circulation biogeochemical models. Because of computational limitations these models still fail to resolve potentially-important spatial scales. Instead, processes on these scales are parameterized. There is concern that deficiencies in these so-called eddy-parameterizations might imprint wrong sensitivities of projected oceanic carbon uptake. Here, we compare natural carbon uptake in the Southern Ocean simulated with contemporary eddy-parameterizations. We find that very differing parameterizations yield surprisingly similar oceanic carbon in response to strengthening winds. In contrast, we find (in an additional simulation) that the carbon uptake does differ substantially when the supply of bioavailable iron is altered within its envelope of uncertainty. We conclude that a more comprehensive understanding of bioavailable iron dynamics will substantially reduce the uncertainty of model-based projections of oceanic carbon uptake.

Description: Continental rifting and ocean basin formation can be observed at the present day in the Red Sea, which is used as the modern analogue for the formation of mid-ocean ridges. Competing theories for how spreading begins—either by quasi-instantaneous formation of a whole spreading segment or by initiation of spreading at multiple discrete “nodes” separated by thinned continental lithosphere—have been put forward based, until recently, on the observations that many seafloor features and geophysical anomalies (gravity, magnetics) along the axis of the Red Sea appeared anomalous compared to ancient and modern examples of ocean basins in other parts of the world. The latest research shows, however, that most of the differences between the Red Sea Rift (RSR) and other (ultra)slow-spreading mid-ocean ridges can be related to its relatively young age and the presence and movement of giant submarine salt flows that blanket large portions of the rift valley. In addition, the geophysical data that was previously used to support the presence of continental crust between the axial basins with outcropping oceanic crust (formerly named “spreading nodes”) can be equally well explained by processes related to the sedimentary blanketing and hydrothermal alteration. The observed spreading nodes are not separated from one another by tectonic boundaries but rather represent “windows” onto a continuous spreading axis which is locally inundated and masked by massive slumping of sediments or evaporites from the rift flanks. Volcanic and tectonic morphologies are comparable to those observed along slow and ultra-slow spreading ridges elsewhere and regional systematics of volcanic occurrences are related to variations in volcanic activity and mantle heat flow. Melt-salt interaction due to salt flows, that locally cover the active spreading segments, and the absence of large detachment faults as a result of the nearby Afar plume are unique features of the RSR. The differences and anomalies seen in the Red Sea still may be applicable to all young oceanic rifts, associated with plumes and/or evaporites, which makes the Red Sea a unique but highly relevant type example for the initiation of slow rifting and seafloor spreading and one of the most interesting targets for future ocean research.

Description: Coral reefs in the Red Sea belong to the most diverse and productive reef ecosystems worldwide, although they are exposed to strong seasonal variability, high temperature, and high salinity. These factors are considered stressful for coral reef biota and challenge reef growth in other oceans, but coral reefs in the Red Sea thrive despite these challenges. In the central Red Sea high temperatures, high salinities, and low dissolved oxygen on the one hand reflect conditions that are predicted for ‘future oceans’ under global warming. On the other hand, alkalinity and other carbonate chemistry parameters are considered favourable for coral growth. In coral reefs of the central Red Sea, temperature and salinity follow a seasonal cycle, while chlorophyll and inorganic nutrients mostly vary spatially, and dissolved oxygen and pH fluctuate on the scale of hours to days. Within these strong environmental gradients micro- and macroscopic reef communities are dynamic and demonstrate plasticity and acclimatisation potential. Epilithic biofilm communities of bacteria and algae, crucial for the recruitment of reef-builders, undergo seasonal community shifts that are mainly driven by changes in temperature, salinity, and dissolved oxygen. These variables are predicted to change with the progression of global environmental change and suggest an immediate effect of climate change on the microbial community composition of biofilms. Corals are so-called holobionts and associate with a variety of microbial organisms that fulfill important functions in coral health and productivity. For instance, coral-associated bacterial communities are more specific and less diverse than those of marine biofilms, and in many coral species in the central Red Sea they are dominated by bacteria from the genus Endozoicomonas. Generally, coral microbiomes align with ecological differences between reef sites. They are similar at sites where these corals are abundant and successful. Coral microbiomes reveal a measurable footprint of anthropogenic influence at polluted sites. Coral-associated communities of endosymbiotic dinoflagellates in central Red Sea corals are dominated by Symbiodinium from clade C. Some corals harbour the same specific symbiont with a high physiological plasticity throughout their distribution range, while others maintain a more flexible association with varying symbionts of high physiological specificity over depths, seasons, or reef locations. The coral-Symbiodinium endosymbiosis drives calcification of the coral skeleton, which is a key process that provides maintenance and formation of the reef framework. Calcification rates and reef growth are not higher than in other coral reef regions, despite the beneficial carbonate chemistry in the central Red Sea. This may be related to the comparatively high temperatures, as indicated by reduced summer calcification and long-term slowing of growth rates that correlate with ocean warming trends. Indeed, thermal limits of abundant coral species in the central Red Sea may have been exceeded, as evidenced by repeated mass bleaching events during previous years. Recent comprehensive baseline data from central Red Sea reefs allow for insight into coral reef functioning and for quantification of the impacts of environmental change in the region.

Description: The Laptev and Eastern Siberian shelves are the world’s broadest shallow shelf systems. Large Siberian rivers and coastal erosion of up to meters per summer deliver large volumes of terrestrial matter into the Arctic shelf seas. In this chapter we investigate the applicability of Ocean Colour Remote Sensing during the ice-free summer season in the Siberian Laptev Sea region. We show that the early summer river peak discharge may be traced using remote sensing in years characterized by early sea-ice retreat. In the summer time after the peak discharge, the spreading of the main Lena River plume east and north-east of the Lena River Delta into the shelf system becomes hardly traceable using optical remote sensing methods. Measurements of suspended particulate matter (SPM) and coloured dissolved organic matter (cDOM) are of the same magnitude in the coastal waters of Buor Khaya Bay as in the Lena River. Match-up analyses of in situ chlorophyll-a (Chl-a) show that standard Medium Resolution Imaging Spectrometer (MERIS) and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite-derived Chl-a is not a valid remote sensing product for the coastal waters and the inner shelf region of the Laptev Sea. All MERIS and MODIS-derived Chl-a products are overestimated by at least a factor of ten, probably due to absorption by the extraordinarily high amount of non-algal particles and cDOM in these coastal and inner-shelf waters. Instead, Ocean Colour remote sensing provides information on wide-spread resuspension over shallows and lateral advection visible in satellite-derived turbidity. Satellite Sea Surface Temperature (SST) data clearly show hydrodynamics and delineate the outflow of the Lena River for hundreds of kilometres out into the shelf seas.

Description: The size of current plankton image datasets renders manual classification virtually infeasible. The training of models for machine classification is complicated by the fact that a large number of classes consist of only a few examples. We employ the recently introduced weight imprinting technique in order to use the available training data to train accurate classifiers in absence of enough examples for some classes. The model architecture used in this work succeeds in the identification of plankton using machine learning with its unique challenges, i.e. a limited number of training examples and a severely skewed class size distribution. Weight imprinting enables a neural network to recognize small classes immediately without re-training. This permits the mining of examples for novel classes.

Description: The Alpine orogeny is well recorded onshore and offshore by tectonic inversion of the Mesozoic rift basins. Large scale linear seamounts (more than 250 km long and with up to 5 km of uplift) involving oceanic and continental lithosphere were carried on top of thrusts, such as the Gorringe seamount and the Estremadura Spur in the SouthWest and West Iberia Margin, respectively. The SouthWest Iberia Margin also recorded the westward migration of the Gibraltar Oceanic slab as the westwards propagation of the Neo-Tethys subduction. Rotation of the tectonic compression from NW-SE to WNW-ESE inPliocene times caused the development of large scale dextral wrench faults as the present day Africa-Iberia plate boundary. Neotectonics of this plate boundary caused large to mega-scale destructive earthquakes and tsunamis.

Description: Hydrothermal circulation at mid-ocean ridges and assimilation of hydrothermally altered crust or hydrothermal fluids by rising magma can be traced by measuring chlorine (Cl) excess in erupted lavas. The Red Sea Rift provides a unique opportunity to study assimilation of hydrothermally altered crust at an ultra-slow spreading ridge (maximum 1.6 cm yr−1 full spreading rate) by Cl, due to its saline seawater (40–42‰, cf. 35‰ in open ocean water), the presence of (hot) brine pools (up to 270‰ salinity and 68 °C) and the thick evaporite sequences that flank the young rift. Absolute chlorine concentrations (up to 1300 ppm) and Cl concentrations relative to minor or trace elements of similar mantle incompatibility (e.g., K, Nb) are much higher in Red Sea basalts than in basalts from average slow spreading ridges. Mantle Cl/Nb concentrations can be used to calculate the Cl-excess, above the magmatic Cl, that is present in the samples. Homogeneous within-sample Cl concentrations, high Cl/H2O, the decoupling of Cl-excess from other trace elements and its independence of the presence of highly saline seafloor brines at the site of eruption indicate that Cl is not enriched at the seafloor. Instead we find basaltic Cl-excess to be spatially closely correlated with evidence of hydrothermal activity, suggesting that deeper assimilation of hydrothermal Cl is the dominant Cl-enrichment process. A proximity of samples to both evaporite outcrops and bathymetric signs of volcanism on the seafloor enhance Cl-excess in basalts. The basaltic Cl-excess can be used as a tracer together with new bathymetric maps as well as indications of hydrothermal venting (hot brine pools, metalliferous Hydrothermal circulation at mid-ocean ridges and assimilation of hydrothermally altered crust or hydrothermal fluids by rising magma can be traced by measuring chlorine (Cl) excess in erupted lavas. The Red Sea Rift provides a unique opportunity to study assimilation of hydrothermally altered crust at an ultra-slow spreading ridge (maximum 1.6 cm yr−1 full spreading rate) by Cl, due to its saline seawater (40–42‰, cf. 35‰ in open ocean water), the presence of (hot) brine pools (up to 270‰ salinity and 68 °C) and the thick evaporite sequences that flank the young rift. Absolute chlorine concentrations (up to 1300 ppm) and Cl concentrations relative to minor or trace elements of similar mantle incompatibility (e.g., K, Nb) are much higher in Red Sea basalts than in basalts from average slow spreading ridges. Mantle Cl/Nb concentrations can be used to calculate the Cl-excess, above the magmatic Cl, that is present in the samples. Homogeneous within-sample Cl concentrations, high Cl/H2O, the decoupling of Cl-excess from other trace elements and its independence of the presence of highly saline seafloor brines at the site of eruption indicate that Cl is not enriched at the seafloor. Instead we find basaltic Cl-excess to be spatially closely correlated with evidence of hydrothermal activity, suggesting that deeper assimilation of hydrothermal Cl is the dominant Cl-enrichment process. A proximity of samples to both evaporite outcrops and bathymetric signs of volcanism on the seafloor enhance Cl-excess in basalts. The basaltic Cl-excess can be used as a tracer together with new bathymetric maps as well as indications of hydrothermal venting (hot brine pools, metalliferous sediments) to predict where hydrothermal venting or now inactive hydrothermal vent fields can be expected. Sites of particular interest for future hydrothermal research are the Mabahiss Deep, the Thetis-HadarbaHatiba Deeps and Shagara-Aswad-Erba Deeps (especially their large axial domes), and Poseidon Deep. Older hydrothermal vent fields may be present at the Nereus and Suakin Deeps. These sites significantly increase the potential of hydrothermal vent field prospection in the Red Sea.

Description: The deep-sea brines of the Red Sea are unusual extreme environments and form characteristically steep gradients across the brine-seawater interfaces. Due to their unusual nature and unique combination of physical-chemical conditions these interfaces provide an interesting source of new findings in the fields of geochemistry, geology, microbiology, biotechnology, virology, and general biology. The current chapter summarizes recent and new results in the study of geochemistry and life at the interfaces of brine-filled deeps of the Red Sea.

Description: S. 1 - Gene Expression Profiling in Fish Toxicology: A Review / Girish Kumar and Nancy D. Denslow S. 39 - Review of the Ecotoxicological Properties of the Methylenedianiline Substances / T. Schupp, H. Allmendinger, B.T.A. Bossuyt, B. Hidding, B. Tury, and R.J. West S. 73 - Cadmium Bioavailability, Uptake, Toxicity and Detoxification in Soil-Plant System / Muhammad Shahid, Camille Dumat, Sana Khalid, Nabeel Khan Niazi, and Paula M.C. Antunes S. 139 - Radionuclides: Accumulation and Transport in Plants / D.K. Gupta, S. Chatterjee, S. Datta, A.V. Voronina, and C. Walther

Description: α2-macroglobulins are broad-spectrum endopeptidase inhibitors, which have to date been characterised from metazoans (vertebrates and invertebrates) and Gram-negative bacteria. Their structural and biochemical properties reveal two related modes of action: the “Venus flytrap” and the “snap-trap” mechanisms. In both cases, peptidases trigger a massive conformational rearrangement of α2-macroglobulin after cutting in a highly flexible bait region, which results in their entrapment. In some homologs, a second action takes place that involves a highly reactive β-cysteinyl-γ-glutamyl thioester bond, which covalently binds cleaving peptidases and thus contributes to the further stabilization of the enzyme:inhibitor complex. Trapped peptidases are still active, but have restricted access to their substrates due to steric hindrance. In this way, the human α2-macroglobulin homolog regulates proteolysis in complex biological processes, such as nutrition, signalling, and tissue remodelling, but also defends the host organism against attacks by external toxins and other virulence factors during infection and envenomation. In parallel, it participates in several other biological functions by modifying the activity of cytokines and regulating hormones, growth factors, lipid factors and other proteins, which has a great impact on physiology. Likewise, bacterial α2-macroglobulins may participate in defence by protecting cell wall components from attacking peptidases, or in host-pathogen interactions through recognition of host peptidases and/or antimicrobial peptides. α2-macroglobulins are more widespread than initially thought and exert multifunctional roles in both eukaryotes and prokaryotes, therefore, their on-going study is essential.

Description: The distribution of the main water masses in the Atlantic Ocean are investigated with the Optimal Multi-Parameter (OMP) method. The properties of the main water masses in the Atlantic Ocean are described in a companion article; here these definitions are used to map out the general distribution of those water masses. Six key properties, including conservative (potential temperature and salinity) and non-conservative (oxygen, silicate, phosphate and nitrate), are incorporated into the OMP analysis to determine the contribution of the water masses in the Atlantic Ocean based on the GLODAP v2 observational data. To facilitate the analysis the Atlantic Ocean is divided into four vertical layers based on potential density. Due to the high seasonal variability in the mixed layer, this layer is excluded from the analysis. Central waters are the main water masses in the upper/central layer, generally featuring high potential temperature and salinity and low nutrient concentrations and are easily distinguished from the intermediate water masses. In the intermediate layer, the Antarctic Intermediate Water (AAIW) from the south can be detected to ~30°N, whereas the Subarctic Intermediate Water (SAIW), having similarly low salinity to the AAIW flows from the north. Mediterranean Overflow Water (MOW) flows from the Strait of Gibraltar as a high salinity water. NADW dominates the deep and overflow layer both in the North and South Atlantic. In the bottom layer, AABW is the only natural water mass with high silicate signature spreading from the Antarctic to the North Atlantic. Due to the change of water mass properties, in this work we renamed to North East Antarctic Bottom Water NEABW north of the equator. Similarly, the distributions of Labrador Sea Water (LSW), Iceland Scotland Overflow Water (ISOW), and Denmark Strait Overflow Water (DSOW) forms upper and lower portion of NADW, respectively roughly south of the Grand Banks between ~50 and 66°N. In the far south the distributions of Circumpolar Deep Water (CDW) and Weddell Sea Bottom Water (WSBW) are of significance to understand the formation of the AABW.

Description: Organic falls can form nutrient-rich, ephemeral hotspots of productivity and biodiversity at the deep-sea floor, especially in food-poor abyssal plains. We report here the first wood falls and second carcass fall recorded from the Clarion-Clipperton Zone in the tropical eastern Pacific Ocean, an area that could be mined for polymetallic nodules in the future. A small cetacean fall in the mobile-scavenger stage likely recently arrived on the seafloor was observed, whereas most of the wood falls were highly degraded. There were multiple species in attendance at the wood falls including organic-fall specialists such as Xylophagaidae molluscs. Many of the taxa attending the carcass fall were known mobile scavengers that regularly attend bait parcels in the Pacific Ocean. These results further confirm that wood falls can occur at large distances (〉1450 km) from major land masses, providing an adequate supply of wood to the abyssal seafloor for colonization by wood-boring molluscs and associated fauna. Organic falls may be regionally abundant and are likely to influence species and habitat diversity in the abyssal areas of the Clarion-Clipperton Zone.

Description: This paper addresses visual navigation of autonomous underwater vehicles (AUVs) with and without a given map, where the latter is called Simultaneous Localization and Mapping (SLAM). We summarize the challenges and opportunities in underwater environments that make visual navigation different from land navigation and also briefly survey the current state-of-the-art in this area. Then as a position paper we argue why many of these challenges could be met by a proper modeling of uncertainties in the SLAM representation. This would in particular allow the SLAM algorithm to thoroughly handle the ambiguity between “I see the same feature again.”, “I see a different but similar looking feature.” and “The environment has changed and the feature moved.”.

Description: The northeast Atlantic encompasses archetypal examples of volcanic rifted margins. Twenty-five years after the last ODP (Ocean Drilling Program) leg on these volcanic margins, the reasons for excess melting are still disputed with at least three competing hypotheses being discussed. We are proposing a new drilling campaign that will constrain the timing, rates of volcanism, and vertical movements of rifted margins. This will allow us to parameterise geodynamic models that can distinguish between the hypotheses. Furthermore, the drilling-derived data will help us to understand the role of breakup magmatism as a potential driver for the Palaeocene–Eocene thermal maximum (PETM) and its influence on the oceanographic circulation in the earliest phase of the northeast Atlantic Ocean formation. Tackling these questions with a new drilling campaign in the northeast Atlantic region will advance our understanding of the long-term interactions between tectonics, volcanism, oceanography, and climate and the functioning of subpolar northern ecosystems and climate during intervals of extreme warmth.

Description: The El Niño Southern Oscillation (ENSO) with its warm (El Niño) and cold (La Niña) phase has strong impacts on marine ecosystems off Peru. This influence extends from changes in nutrient availability to productivity and oxygen levels. While several studies have demonstrated the influence of ENSO events on biological productivity, less is known about their impact on oxygen concentrations. In situ observations along the Peruvian and Chilean coast have shown a strong water column oxygenation during the 1997/1998 strong El Niño event. These observations suggest a deepening of the oxygen minimum zone (OMZ) along the continental shelf. However, due to reduced spatial coverage of the existing in situ observations, no studies have yet demonstrated the OMZ response to El Niño events in the whole Eastern Tropical South Pacific (ETSP). Furthermore, most studies have focused on El Niño events. Much less attention was given to the oxygen dynamics under La Niña influence. Here, we provide a comprehensive analysis of the ENSO influence on OMZ dynamics. Interannual variability of the OMZ during the period 1990–2010 is derived from a regional coupled physical-biogeochemical model forced with realistic atmospheric and lateral boundary conditions. Our results show a reduction of the vertical extent and a deepening of suboxic waters (SW) during the El Niño phase. During the La Niña phase, there is a vertical expansion of SW. These fluctuations in OMZ extent are due to changes in oxygen supply into its core depth mainly from lateral margins. During the El Niño phase, the enhanced lateral oxygen supply from the subtropics is the main reason for the reduction of SW in both coastal and offshore regions. During the La Niña phase, the oxygenated subtropical waters are blocked by the poleward transport along the southern margin of the OMZ. Consequently, oxygen concentrations within the OMZ are reduced and suboxic conditions expand during La Niña. The detailed analysis of transport pathways presented here provides new insights into how ENSO variability affects the oxygen-sensitive marine biogeochemistry of the ETSP.

Description: Oxygen minimum zones (OMZs) in the open ocean occur below the surface in regions of weak ventilation and high biological productivity. Very low levels of dissolved oxygen affect marine life and alter biogeochemical cycles. One of the most intense but least understood OMZs in the world is located in the Arabian Sea in a depth range between 300 to 1000 m. Within the last decades observations suggest a decreasing oxygen trend. Thus, an improved understanding of the crucial processes is necessary for a reliable assessment of the future development of the Arabian Sea OMZ. This study uses a combination of observational data as well as reanalysis velocity fields from the ocean model Hycom (Hybrid Coordinate Ocean Model) to explore the ventilation dynamics of the Arabian Sea OMZ. Our results show that the OMZ features a strong seasonal cycle with regional differences that is correlated with the monsoon system: In the eastern basin, the OMZ is strongest during the winter monsoon with a core thickness of 1000 m depth and oxygen values of less than 5 µmol/kg. Ventilation during that phase is dominated by Persian Gulf water, that clockwise circles the perimeter of the basin and enters the OMZ from the north. During the summer monsoon ventilation from the southeast leads to higher oxygen values indicating a reverse flow along the Indian coast in the intermediate layer compared to the southeastward surface currents. The seasonal cycle in the western basin has the same seasonality as the one in the eastern basin with a core thickness of 900 m during the winter monsoon. The oxygen supply during the summer monsoon is weaker compared to the eastern basin and correlates with the ventilation of Persian Gulf (Red Sea) water during the summer monsoon (autumn inter-monsoon) phase. As the interior exchange between the eastern and western basin is weak, the more pronounced OMZ in the eastern basin is explained by prolonged ventilation time scales. For the eastern (western) basin Persian Gulf water needs 2–3 (1–2) years and Red Sea water 7–8 (3–4) years to ventilate the OMZ.

Description: The application of acoustic techniques, such as multibeam echosounders, has permitted the identification of Maltese submarine landscapes and landforms that were progressively inundated during the postglacial sea-level rise. Remarkably, geomorphological features due to fluvial, gravity-induced and karst processes that took place under former subaerial conditions can be clearly recognised on the present seafloor around the Maltese archipelago, and they were only slightly modified by sea action during the postglacial transgression phases. The analysis of the submerged landforms described in this chapter is crucial for understanding the evolution of the Maltese Islands during the last ca. 20,000 years.

Description: The study presents the results of continuous measurements of carbon dioxide concentration in the atmospheric surface layer at Tiksi and Cape Baranov Arctic stations over the period of August 2010–May 2017 and over the whole 2016, respectively. The amplitude of diurnal variations in the CO2 concentration in Tiksi from June to September is 1.1 ± 1.3, 2.4 ± 2.0, 4.1 ± 2.3, and 2.0 ± 2.4 ppm. Diurnal variations in CO2 at Cape Baranov station are absent. The observed seasonal variations in the CO2 concentration are compared with the data of the MBL empirical model for the marine atmospheric boundary layer of the Arctic region. In 2016, the difference between the observed and model concentrations at Tiksi and Cape Baranov stations amounted to 1.7 and 0.5 ppm, respectively, in winter and −3.0 and −1.9 ppm, respectively, in summer. It is shown that wildfires in Siberia caused a long synchronous increase in the CO2 concentration by 20 ppm in Tiksi and by 15 ppm at Cape Baranov station.

Description: The water of Atlantic origin (Atlantic water) is of special importance for the formation of the hydrological regime of the Arctic Ocean and provides a “building material” for the main water masses in the surface and intermediate layers. Atlantic water are structurally included to the Arctic Transpolar System (ATS); it is a multicomponent physical object whose state defines the role of the Arctic Ocean in the planetary climate. The recent advances in the Atlantic water research are discussed, in particular, the role of Atlantic water in the Arctic sea ice reduction. Particular attention is paid to the possible activation of feedbacks in ATS that may result in an accelerated Arctic sea ice loss recorded after 2007.

Description: Three most common Arctic foraminiferal species of the family Cassidulinidae from the Laptev Sea sediment cores were analyzed for paleoenvironmental reconstructions. This paleontological study is focused on morphological characteristics of apertural elements and the wall structure. The species Cassidulina neoteretis is considered an indicator of the Atlantic influence in contrast to the Arctic species Cassidulina reniforme and Islandiella norcrossi.

Description: The ecology and taxonomy of Islandiella norcrossi (Cushman, 1933), a typical representative of the Arctic calcareous benthic foraminifers, are analyzed based on microfossils from the Late Pleistocene-Holocene deposits of the Laptev Sea. This species is distributed on the shelf and continental slope of the Arctic seas, and indicates normal marine conditions with a seasonal ice cover and cold Arctic water masses. The species is also associated with a high seasonal productivity in sea-ice marginal zone. Tests are abundant in the Laptev Sea sediments, and show some variability in the size and shape of chambers, which was possibly caused by harsh environmental conditions of the Arctic (limited period of growth and reproduction); aberrant forms with an additional aperture being also present. The population contains tests of macrospherical and microspherical generations, as well as juveniles.

Description: Atmospheric processes, ranging from microscale turbulence to severe storms on the synoptic scale, impact the continuous ground motion of the earth and have the potential to induce strong broad-band noise in seismological recordings. We designed a target-oriented experiment to quantify the influence of wind on ground motion velocity in the Dead Sea valley. For the period from March 2014 to February 2015, a seismological array, consisting of 15 three-component short-period and broad-band stations, was operated near Madaba, Jordan, complemented by one meteorological tower providing synchronized, continuous three-component measurements of wind speed. Results reveal a pronounced, predominantly linear increase of the logarithmic power of ground motion velocity with rising mean horizontal wind speed at all recording stations. Measurements in rough, mountainous terrain further identify a strong dependency of wind-induced noise on surface characteristics, such as topography and, therefore, demonstrate the necessity to consider wind direction as well. To assess the noise level of seismological recordings with respect to a dynamically changing wind field, we develop a methodology to account for the dependency of power spectral density of ground motion velocity on wind speed and wind direction for long, statistically significant periods. We further introduce the quantitative measure of the ground motion susceptibility to estimate the vulnerability of seismological recordings to the presence of wind.

Description: Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

Description: A set of ensemble integrations from the Coupled Model Intercomparison Project phase 5, with historical forcing plus RCP4.5 scenario, are used to explore if state-of-the-art climate models are able to simulate previously reported linkages between sea-ice concentration (SIC) anomalies over the eastern Arctic, namely in the Greenland–Barents–Kara Seas, and lagged atmospheric circulation that projects on the North Atlantic Oscillation (NAO)/Arctic Oscillation (AO). The study is focused on variability around the long-term trends, so that all anomalies are detrended prior to analysis; the period of study is 1979–2013. The model linkages are detected by applying maximum covariance analysis. As also found in observational data, all the models considered here show a statistically significant link with sea-ice reduction over the eastern Arctic followed by a negative NAO/AO-like pattern. If the simulated relationship is found at a lag of one month, the results suggest that a stratospheric pathway could be at play as the driving mechanism; in observations this is preferentially shown for SIC in November. The interference of a wave-like anomaly over Eurasia, accompanying SIC changes, with the climatological wave pattern appears to be key in setting the mediating role of the stratosphere. On the other hand, if the simulated relationship is found at a lag of two months, the results suggest that tropospheric dynamics are dominant, presumably due to transient eddy feedback; in observations this is preferentially shown for SIC in December. The results shown here and previous evidence from atmosphere-only experiments emphasize that there could be a detectable influence of eastern Arctic SIC variability on mid-latitude atmospheric circulation anomalies. Even if the mechanisms are robust among the models, the timing of the simulated linkages strongly depends on the model and does not generally mimic the observational ones. This implies that the atmospheric sensitivity to sea-ice changes largely depends on the mean-flow and parameterizations, which could lead to misleading conclusions elsewhere if a multi-model ensemble-mean approach is adopted. It might also represent an important source of uncertainty in climate prediction and projection. Modelling efforts are hence further required to improve representation of the background atmospheric circulation and reduce biases, in order to attain more accurate covariability.

Description: Erstmals trägt dieses nationale Assessment den Forschungsstand zum Klimawandel umfassend für alle Themenbereiche und gesellschaftlichen Sektoren zusammen. Womit müssen wir in Deutschland rechnen, welche Auswirkungen werden die Klimaveränderungen auf Wirtschaft und Gesellschaft haben, und wie können wir uns wappnen? 126 Autoren aus ganz Deutschland äußern sich zu Themen wie bereits beobachtete und zukünftige Veränderungen, Wetterkatastrophen und deren Folgen, den Projektionen für die Zukunft, den Risiken sowie möglichen Anpassungsstrategien. Die Autoren stellen in verständlicher Sprache den aktuellen Forschungsstand dar und veranschaulichen die wichtigsten Gedanken in Grafiken und Tabellen. Alle Texte wurden mehrfach wissenschaftlich begutachtet. Klimawandel in Deutschland ist die erste Gesamtschau zu dem Themenkomplex, benennt offene Fragestellungen und liefert eine Grundlage für Entscheidungen im Zusammenhang mit dem Klimawandel. Deutschland reiht sich damit ein in die Liste von Ländern wie die Vereinigten Staaten, Österreich und Großbritannien, in denen derartige Berichte bereits vorliegen. Die Herausgeber: Prof. Dr. Guy Brasseur und Prof. Dr. Daniela Jacob, ehemaliger Direktor und gegenwärtige Direktorin des Climate Center Germany/Helmholtz-Zentrum Geesthacht, und Susanne Schuck-Zöller, die auch die Projektleitung übernommen hat, werden bei diesem Buch unterstützt von einem Editorial Board, dem elf herausragende Wissenschaftler aus den wichtigsten Klimaforschungseinrichtungen in Deutschland angehören

Description: Dome ports act as spherical windows in underwater housings through which a camera can observe objects in the water. As compared to flat glass interfaces, they do not limit the field of view, and they do not cause refraction of light observed by a pinhole camera positioned exactly in the center of the dome. Mechanically adjusting a real lens to this position is a challenging task, in particular for those integrated in deep sea housings. In this contribution a mechanical adjustment procedure based on straight line observations above and below water is proposed that allows for accurate alignments. Additionally, we show a chessboard-based method employing an underwater/above-water image pair to estimate potentially remaining offsets from the dome center to allow refraction correction in photogrammetric applications. Besides providing intuition about the severity of refraction in certain settings, we demonstrate the methods on real data for acrylic and glass domes in the water.

Description: We present here the first results, for the preindustrial and mid-Holocene climatological periods, of the newly developed isotope-enhanced version of the fully coupled Earth system model MPI-ESM, called hereafter MPI-ESM-wiso. The water stable isotopes H162O, H182O and HDO have been implemented into all components of the coupled model setup. The mid-Holocene provides the opportunity to evaluate the model response to changes in the seasonal and latitudinal distribution of insolation induced by different orbital forcing conditions. The results of our equilibrium simulations allow us to evaluate the performance of the isotopic model in simulating the spatial and temporal variations of water isotopes in the different compartments of the hydrological system for warm climates. For the preindustrial climate, MPI-ESM-wiso reproduces very well the observed spatial distribution of the isotopic content in precipitation linked to the spatial variations in temperature and precipitation rate. We also find a good model–data agreement with the observed distribution of isotopic composition in surface seawater but a bias with the presence of surface seawater that is too 18O-depleted in the Arctic Ocean. All these results are improved compared to the previous model version ECHAM5/MPIOM. The spatial relationships of water isotopic composition with temperature, precipitation rate and salinity are consistent with observational data. For the preindustrial climate, the interannual relationships of water isotopes with temperature and salinity are globally lower than the spatial ones, consistent with previous studies. Simulated results under mid-Holocene conditions are in fair agreement with the isotopic measurements from ice cores and continental speleothems. MPI-ESM-wiso simulates a decrease in the isotopic composition of precipitation from North Africa to the Tibetan Plateau via India due to the enhanced monsoons during the mid-Holocene. Over Greenland, our simulation indicates a higher isotopic composition of precipitation linked to higher summer temperature and a reduction in sea ice, shown by positive isotope–temperature gradient. For the Antarctic continent, the model simulates lower isotopic values over the East Antarctic plateau, linked to the lower temperatures during the mid-Holocene period, while similar or higher isotopic values are modeled over the rest of the continent. While variations of isotopic contents in precipitation over West Antarctica between mid-Holocene and preindustrial periods are partly controlled by changes in temperature, the transport of relatively 18O-rich water vapor near the coast to the western ice core sites could play a role in the final isotopic composition. So, more caution has to be taken about the reconstruction of past temperature variations during warm periods over this area. The coupling of such a model with an ice sheet model or the use of a zoomed grid centered on this region could help to better describe the role of the water vapor transport and sea ice around West Antarctica. The reconstruction of past salinity through isotopic content in sea surface waters can be complicated for regions with strong ocean dynamics, variations in sea ice regimes or significant changes in freshwater budget, giving an extremely variable relationship between the isotopic content and salinity of ocean surface waters over small spatial scales. These complicating factors demonstrate the complexity of interpreting water isotopes as past climate signals of warm periods like the mid-Holocene. A systematic isotope model intercomparison study for further insights on the model dependency of these results would be beneficial.

Description: Probabilistic spatial reconstructions of past climate states are valuable to quantitatively study the climate system under different forcing conditions because they combine the information contained in a proxy synthesis into a comprehensible product. Unfortunately, they are subject to a complex uncertainty structure due to complicated proxy–climate relations and sparse data, which makes interpolation between samples difficult. Bayesian hierarchical models feature promising properties to handle these issues, like the possibility to include multiple sources of information and to quantify uncertainties in a statistically rigorous way. We present a Bayesian framework that combines a network of pollen and macrofossil samples with a spatial prior distribution estimated from a multi-model ensemble of climate simulations. The use of climate simulation output aims at a physically reasonable spatial interpolation of proxy data on a regional scale. To transfer the pollen data into (local) climate information, we invert a forward version of the probabilistic indicator taxa model. The Bayesian inference is performed using Markov chain Monte Carlo methods following a Metropolis-within-Gibbs strategy. Different ways to incorporate the climate simulations into the Bayesian framework are compared using identical twin and cross-validation experiments. Then, we reconstruct the mean temperature of the warmest and mean temperature of the coldest month during the mid-Holocene in Europe using a published pollen and macrofossil synthesis in combination with the Paleoclimate Modelling Intercomparison Project Phase III mid-Holocene ensemble. The output of our Bayesian model is a spatially distributed probability distribution that facilitates quantitative analyses that account for uncertainties.

Description: We discuss the results of measurements in the region of Cape Baranov (the Severnaya Zemlya archipelago) of the set of physicochemical characteristics of atmospheric aerosol: aerosol optical depth, aerosol and black carbon concentrations, elemental and ion compositions of aerosol, organic and elemental carbon contents in aerosol, as well as the isotopic composition of carbon in the aerosol and snow samples. It is shown that the average values of most aerosol characteristics, measured in April–June 2018, are a little lower than in the Arctic settlement Barentsburg (Spitsbergen archipelago) and several-fold smaller than in the south of Western Siberia in the same period.

Description: Plastic has become indispensable for human life. When plastic debris is discarded into waterways, these items can interact with organisms. Of particular concern are microscopic plastic particles (microplastics) which are subject to ingestion by several taxa. This review summarizes the results of cutting-edge research about the interactions between a range of aquatic species and microplastics, including effects on biota physiology and secondary ingestion. Uptake pathways via digestive or ventilatory systems are discussed, including (1) the physical penetration of microplastic particles into cellular structures, (2) leaching of chemical additives or adsorbed persistent organic pollutants (POPs), and (3) consequences of bacterial or viral microbiota contamination associated with microplastic ingestion. Following uptake, a number of individual-level effects have been observed, including reduction of feeding activities, reduced growth and reproduction through cellular modifications, and oxidative stress. Microplastic-associated effects on marine biota have become increasingly investigated with growing concerns regarding human health through trophic transfer. We argue that research on the cellular interactions with microplastics provide an understanding of their impact to the organisms’ fitness and, therefore, its ability to sustain their functional role in the ecosystem. The review summarizes information from 236 scientific publications. Of those, only 4.6% extrapolate their research of microplastic intake on individual species to the impact on ecosystem functioning. We emphasize the need for risk evaluation from organismal effects to an ecosystem level to effectively evaluate the effect of microplastic pollution on marine environments. Further studies are encouraged to investigate sublethal effects in the context of environmentally relevant microplastic pollution conditions.

Description: Large calderas are among the Earth's major volcanic features. They are associated with large magma reservoirs and elevated geothermal gradients. Caldera-forming eruptions result from the withdrawal and collapse of the magma chambers and produce large-volume pyroclastic deposits and later-stage deformation related to post-caldera resurgence and volcanism. Unrest episodes are not always followed by an eruption; however, every eruption is preceded by unrest. The Campi Flegrei caldera (CFc), located along the eastern Tyrrhenian coastline in southern Italy, is close to the densely populated area of Naples. It is one of the most dangerous volcanoes on Earth and represents a key example of an active, resurgent caldera. It has been traditionally interpreted as a nested caldera formed by collapses during the 100–200 km3 Campanian Ignimbrite (CI) eruption at ∼39 ka and the 40 km3 eruption of the Neapolitan Yellow Tuff (NYT) at ∼15 ka. Recent studies have suggested that the CI may instead have been fed by a fissure eruption from the Campanian Plain, north of Campi Flegrei. A MagellanPlus workshop was held in Naples, Italy, on 25–28 February 2017 to explore the potential of the CFc as target for an amphibious drilling project within the International Ocean Discovery Program (IODP) and the International Continental Drilling Program (ICDP). It was agreed that Campi Flegrei is an ideal site to investigate the mechanisms of caldera formation and associated post-caldera dynamics and to analyze the still poorly understood interplay between hydrothermal and magmatic processes. A coordinated onshore–offshore drilling strategy has been developed to reconstruct the structure and evolution of Campi Flegrei and to investigate volcanic precursors by examining (a) the succession of volcanic and hydrothermal products and related processes, (b) the inner structure of the caldera resurgence, (c) the physical, chemical, and biological characteristics of the hydrothermal system and offshore sediments, and (d) the geological expression of the phreatic and hydromagmatic eruptions, hydrothermal degassing, sedimentary structures, and other records of these phenomena. The deployment of a multiparametric in situ monitoring system at depth will enable near-real-time tracking of changes in the magma reservoir and hydrothermal system.

Description: The international endeavour to retrieve a continuous ice core, which spans the middle Pleistocene climate transition ca. 1.2–0.9 Myr ago, encompasses a multitude of field and model-based pre-site surveys. We expand on the current efforts to locate a suitable drilling site for the oldest Antarctic ice core by means of 3-D continental ice-sheet modelling. To this end, we present an ensemble of ice-sheet simulations spanning the last 2 Myr, employing transient boundary conditions derived from climate modelling and climate proxy records. We discuss the imprint of changing climate conditions, sea level and geothermal heat flux on the ice thickness, and basal conditions around previously identified sites with continuous records of old ice. Our modelling results show a range of configurational ice-sheet changes across the middle Pleistocene transition, suggesting a potential shift of the West Antarctic Ice Sheet to a marine-based configuration. Despite the middle Pleistocene climate reorganisation and associated ice-dynamic changes, we identify several regions conducive to conditions maintaining 1.5 Myr (million years) old ice, particularly around Dome Fuji, Dome C and Ridge B, which is in agreement with previous studies. This finding strengthens the notion that continuous records with such old ice do exist in previously identified regions, while we are also providing a dynamic continental ice-sheet context.

Description: While nonlinear stochastic partial differential equations arise naturally in spatiotemporal modeling, inference for such systems often faces two major challenges: sparse noisy data and ill-posedness of the inverse problem of parameter estimation. To overcome the challenges, we introduce a strongly regularized posterior by normalizing the likelihood and by imposing physical constraints through priors of the parameters and states. We investigate joint parameter-state estimation by the regularized posterior in a physically motivated nonlinear stochastic energy balance model (SEBM) for paleoclimate reconstruction. The high-dimensional posterior is sampled by a particle Gibbs sampler that combines a Markov chain Monte Carlo (MCMC) method with an optimal particle filter exploiting the structure of the SEBM. In tests using either Gaussian or uniform priors based on the physical range of parameters, the regularized posteriors overcome the ill-posedness and lead to samples within physical ranges, quantifying the uncertainty in estimation. Due to the ill-posedness and the regularization, the posterior of parameters presents a relatively large uncertainty, and consequently, the maximum of the posterior, which is the minimizer in a variational approach, can have a large variation. In contrast, the posterior of states generally concentrates near the truth, substantially filtering out observation noise and reducing uncertainty in the unconstrained SEBM.

Description: Stable carbon isotopes of sediment organic matter (δ13COM) are widely applied in paleoenvironmental studies. Interpretations of δ13COM, however, remain challenging and factors that influence δ13COM may not apply across all lakes. Common explanations for stratigraphic shifts in δ13COM include changes in lake productivity or changes in inputs of allochthonous OM. We investigated the influence of different oxygen conditions (oxic versus anoxic) on the δ13COM values in the sediments of Lake Tiefer See. We analysed (1) a long sediment core from the deepest part of the lake, (2) two short, sediment–water interface cores from shallower water depths, and (3) OM in the water column, i.e. from sediment traps. Fresh OM throughout the entire water column showed a relatively constant δ13COM value of approximately − 30.5‰. Similar values, about − 31‰, were obtained for well-varved sediments in both the long and short, sediment–water interface cores. In contrast, δ13COM values from non-varved sediments in all cores were significantly less negative (− 29‰). The δ13COM values in the sediment–water interface cores from different water depths differ for sediments of the same age, if oxygen conditions at the time of deposition were different at these sites, as suggested by the state of varve preservation. Sediments deposited from AD 1924 to 1980 at 62 m water depth are varved and exhibit δ13COM values around − 31‰, whereas sediments of the same age in the core from 35 m water depth are not varved and show less negative δ13COM values of about − 29‰. The relation between varve occurrence and δ13COM values suggests that δ13COM is associated with oxygen conditions because varve preservation depends on hypolimnetic anoxia. A mechanism that likely influences δ13COM is selective degradation of OM under oxic conditions, such that organic components with more negative δ13COM are preferably decomposed, leading to less negative δ13COM values in the remaining, undegraded OM pool. Greater decomposition of OM in non-varved sediments is supported by lower TOC concentrations in these deposits (~ 5%) compared to well-varved sediments (~ 15%). Even in lakes that display small variations in productivity and terrestrial OM input through time, large spatial and temporal differences in hypolimnetic oxygen concentrations may be an important factor controlling sediment δ13COM.

Description: We present a software toolbox that allows the efficient collection, management and maintenance of larger paleoceanographic data sets. The program combines a graphical user interface (GUI) with a simple document-based database and functionality for visualization, stratigraphy (visual alignment and radiocarbon calibration), age modelling and efficient ensemble time-series generation to create larger homogenous data compilations. Data can be imported from Excel or text files, are stored locally in netCDF format and can be easily exchanged between collaborating scientists. Within a data collection, data can be imported either to proxy-specific sub-collections or to a multi-proxy (“miscellaneous”) sub-collection that allows the import and management of any downcore data. A single age model is shared among all proxies of a core within a collection. The stand-alone software can be used with Windows and macOS and does not require web access. Installers of the current version for both Windows 10 and macOS including the C++ code can be downloaded from https://www.marum.de/Stefan-Mulitza/PaleoDataView.html (last access: 5 December 2019) along with a detailed user guide.

Description: Combining ocean general circulation models with proxy data via data assimilation is a means to obtain estimates of past ocean states that are consistent with model physics as well as with proxy data. The climate during the Last Glacial Maximum (LGM, 19–23 ka) was substantially different from today. Even though boundary conditions are comparatively well known, the large-scale patterns of the ocean circulation during this time remain uncertain. Previous efforts to combine ocean models with proxy data have shown dissimilar results regarding the state of the ocean, in particular of the Atlantic Meridional Overturning Circulation. Here, we present a new LGM ocean state estimate that extents previous estimates by using global benthic as well as planktic data on the oxygen isotopic composition of calcite. It is further constrained by global seasonal and annual sea surface temperature (SST) reconstructions. The estimate shows an Atlantic Ocean that is similar to the Late Holocene Atlantic Ocean but with a reduced formation of Antarctic Bottom Water, in contrast to results of previous studies. The results indicate that SST and oxygen isotopic data alone do not require the presence of a shallower North Atlantic Deep Water and a more extensive Antarctic Bottom Water, and highlight the need for more proxy data of different types to obtain reliable ocean state estimates. Additional adjoint sensitivity experiments reveal that data from the deep North Atlantic and from the global deep Southern Ocean are most important to constrain the Atlantic Meridional Overturning Circulation.

Description: In den tropischen und subtropischen Meeren existieren in mittleren Tiefen riesige sauerstoffarme Zonen. Im Zuge des Klimawandels dehnen sie sich immer stärker aus. Auch in Küstenregionen entstehen durch Stickstoffbelastung aus der Landwirtschaft lebensfeindliche Zonen ohne Sauerstoff – mit verheerenden Folgen für das marine Ökosystem.

Description: The repercussions of surface ocean currents for the near-surface wind and the air-sea momentum flux are investigated in two versions of a global climate model with eddying ocean. The focus is on the effect of mesoscale ocean current features at scales of less than 150 km, by considering high-pass filtered, monthly-mean model output fields. We find a clear signature of a mesoscale oceanic imprint in the wind fields over the energetic areas of the oceans, particularly along the extensions of the western boundary currents and the Antarctic Circumpolar Current. These areas are characterized by a positive correlation between mesoscale perturbations in the curl of the surface currents and the wind curl. The coupling coefficients are spatially non-uniform and show a pronounced seasonal cycle. The positive feedback of mesoscale current features on the near-surface wind acts in opposition to their damping effect on the wind stress. A tentative incorporation of this feedback in the surface stress formulation of an eddy-permitting global ocean-only model leads to a gain in the kinetic energy of up to 10 %, suggesting a fundamental shortcoming of present ocean model configurations.

Description: The invasive round goby has established a viable population within 9 years of its first introduction to Lithuanian coastal waters (SE Baltic Sea). During its expansion phase, abundances increased 23-fold, which led to the near complete eradication of its main prey, the blue mussel, at 〈 20 m depth. The round goby population showed a stabilizing trend after blue mussel biomass was depleted; however, their abundance has not declined. The round goby feeds efficiently on newly settled mollusks, causing a severe constraint for blue mussel recovery. Changes in blue mussel availability and size structure induced a dietary shift in wintering long-tailed duck towards fish prey. An energetically dense food source sustains a good body condition in long-tailed ducks, however the change in trophic position (from 3.1 to 4.3 trophic level) suggests the potential for a reduction in their carrying capacity. Results from this study also show that coastal habitats with low and unpredictable population dynamics of blue mussel become less attractive wintering sites for long-tailed duck in the Baltic Sea. We also document a cascading effect of invasive species in the food web.

Description: The species composition of many groups of marine plankton appears well predicted by sea surface temperature (SST). Consequently, fossil plankton assemblages have been widely used to reconstruct past SST. Most applications of this approach make use of the highest possible taxonomic resolution. However, not all species are sensitive to temperature, and their distribution may be governed by other parameters. There are thus reasons to question the merit of including information about all species, both for transfer function performance and for its effect on reconstructions. Here we investigate the effect of species selection on planktonic foraminifera transfer functions. We assess species importance for transfer function models using a random forest technique and evaluate the performance of models with an increasing number of species. Irrespective of using models that use the entire training set (weighted averaging) or models that use only a subset of the training set (modern analogue technique), we find that the majority of foraminifera species does not carry useful information for temperature reconstruction. Less than one-third of the species in the training set is required to provide a temperature estimate with a prediction error comparable to a transfer function that uses all species in the training set. However, species selection matters for paleotemperature estimates. We find that transfer function models with a different number of species but with the same error may yield different reconstructions of sea surface temperature when applied to the same fossil assemblages. This ambiguity in the reconstructions implies that fossil assemblage change reflects a combination of temperature and other environmental factors. The contribution of the additional factors is site and time specific, indicating ecological and geological complexity in the formation of the sedimentary assemblages. The possibility of obtaining multiple different reconstructions from a single sediment record presents a previously unrecognized source of uncertainty for sea surface temperature estimates based on planktonic foraminifera assemblages. This uncertainty can be evaluated by determining the sensitivity of the reconstructions to species pruning.

Description: We introduce a new parallelizable numerical multiscale method for advection-dominated problems as they often occur in engineering and geosciences. State of the art multiscale simulation methods work well in situations in which stationary and elliptic scenarios prevail but are prone to fail when the model involves dominant lower order terms which is common in applications. We suggest to overcome the associated difficulties through a reconstruction of subgrid variations into a modified basis by solving many independent (local) inverse problems that are constructed in a semi-Lagrangian step. Globally the method looks like a Eulerian method with multiscale stabilized basis. The method is extensible to other types of Galerkin methods, higher dimensions, nonlinear problems and can potentially work with real data. We provide examples inspired by tracer transport in climate systems in one and two dimensions and numerically compare our method to standard methods.

Description: Natural polysaccharides of marine origin are gaining interest in biomedical applications. Seaweeds are most abundant source of polysaccharides, as alginates, agar and agarose as well as Carrageenans. Even cellulose and amylose have been extracted from the macroalgae. Chitin and chitosan are derived from the exoskeleton of marine crustaceans. Interdisciplinary fields involving various science and technology aspects such as cell sciences, biomaterials, medical sciences and engineering are referred to as tissue engineering, which is an upcoming new field intended to replace biological functions in human body. Tissue engineered scaffolds and artificial organs developed by such technique has replace injured parts in human body. Technological advancements have made it possible to obtain active ingredient in marine organisms by controlling the growth and isolation conditions. Present review has focused on progress in discovering and producing new applications of marine polysaccharides in biomedical and tissue engineering.

Description: Paleoclimate data assimilation (DA) is a promising technique to systematically combine the information from climate model simulations and proxy records. Here, we investigate the assimilation of tree-ring-width (TRW) chronologies into an atmospheric global climate model using ensemble Kalman filter (EnKF) techniques and a process-based tree-growth forward model as an observation operator. Our results, within a perfect-model experiment setting, indicate that the "online DA" approach did not outperform the "off-line" one, despite its considerable additional implementation complexity. On the other hand, it was observed that the nonlinear response of tree growth to surface temperature and soil moisture does deteriorate the operation of the time-averaged EnKF methodology. Moreover, for the first time we show that this skill loss appears significantly sensitive to the structure of the growth rate function, used to represent the principle of limiting factors (PLF) within the forward model. In general, our experiments showed that the error reduction achieved by assimilating pseudo-TRW chronologies is modulated by the magnitude of the yearly internal variability in the model. This result might help the dendrochronology community to optimize their sampling efforts.

Description: Glacial cycles of the late Quaternary are controlled by the asymmetrically varying mass balance of continental ice sheets in the Northern Hemisphere. Surface mass balance is governed by processes of ablation and accumulation. Here two ablation schemes, the positive-degree-day (PDD) method and the surface energy balance (SEB) approach, are compared in transient simulations of the last glacial cycle with the Earth system model of intermediate complexity CLIMBER-2. The standard version of the CLIMBER-2 model incorporates the SEB approach and simulates ice volume variations in reasonable agreement with paleoclimate reconstructions during the entire last glacial cycle. Using results from the standard CLIMBER-2 model version, we simulated ablation with the PDD method in offline mode by applying different combinations of three empirical parameters of the PDD scheme. We found that none of the parameter combinations allow us to simulate a surface mass balance of the American and European ice sheets that is similar to that obtained with the standard SEB method. The use of constant values for the empirical PDD parameters led either to too much ablation during the first phase of the last glacial cycle or too little ablation during the final phase. We then substituted the standard SEB scheme in CLIMBER-2 with the PDD scheme and performed a suite of fully interactive (online) simulations of the last glacial cycle with different combinations of PDD parameters. The results of these simulations confirmed the results of the offline simulations: no combination of PDD parameters realistically simulates the evolution of the ice sheets during the entire glacial cycle. The use of constant parameter values in the online simulations leads either to a buildup of too much ice volume at the end of glacial cycle or too little ice volume at the beginning. Even when the model correctly simulates global ice volume at the last glacial maximum (21 ka), it is unable to simulate complete deglaciation during the Holocene. According to our simulations, the SEB approach proves superior for simulations of glacial cycles.

Description: Using new high-resolution 10Be measurements in the NGRIP, EDML and Vostok ice cores, together with previously published data from EDC, we present an improved synchronization between Greenland and Antarctic ice cores during the Laschamp geomagnetic excursion  ∼ 41kyr ago. We estimate the precision of this synchronization to be ±20 years, an order of magnitude better than previous work. We discuss the implications of this new synchronization for making improved estimates of the depth difference between ice and enclosed gas of the same age (Δdepth), difference between age of ice and enclosed gas at the same depth (Δage) in the EDC and EDML ice cores, spectral properties of the 10Be profiles and phasing between Dansgaard–Oeschger-10 (in NGRIP) and AIM-10 (in EDML and EDC).

Description: The Last Glacial Maximum (LGM, 21 000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to the present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional (Tier 2) sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalize on the relative abundance of paleoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Description: Global warming has led to a strong deterioration of the Arctic sea ice cover. Ice thickness, age and coverage have been strongly declining in recent years. Brine channels that form in sea ice when seawater freezes represent a unique habitat for bacteria, algae, proto- and small metazoans. We hypothesized that the loss of multi-year ice and the more prevalent formation of first-year ice even in central regions of the Arctic will lead to changes in the Arctic sea ice meiofauna community composition. We therefore analysed the sea ice meiofauna community composition of three different ice types sampled in summer and autumn 2007. Young, thin ice of few cm thickness was typified by taxa of pelagic origin or with good swimming abilities (ciliates, pelagic foraminifera, rotifers and platyhelminthes). Harpacticoid copepods and nematodes with poor swimming abilities were prevalent in older, thicker (〉0.5 m) first- and multi-year ice. Brash ice—which was likely a mix of older broken ice, slush and pancake ice—was characterized by a high abundance of platyhelminthes and rotifers. An experimental analysis of colonization efficiencies of artificial thin ice also revealed that species with poor swimming ability are less successful to colonize newly forming thin ice. We conclude that observed and predicted changes in the ice formation regime will likely result in changes in the composition of Arctic sea ice communities. We predict negative effects particularly for species with low dispersal capacities like harpacticoid copepods and endemic nematodes, as these are less successful in colonizing newly forming thin ice.

Description: The shallow ice approximation (SIA) is commonly used in ice-sheet models to simplify the force balance equations within the ice. However, the SIA cannot adequately reproduce the dynamics of the fast flowing ice streams usually found at the margins of ice sheets. To overcome this limitation, recent studies have introduced heuristic hybrid combinations of the SIA and the shelfy stream approximation. Here, we implement four different hybrid schemes into a model of the Antarctic Ice Sheet in order to compare their performance under present-day conditions. For each scheme, the model is calibrated using an iterative technique to infer the spatial variability in basal sliding parameters. Model results are validated against topographic and velocity data. Our analysis shows that the iterative technique compensates for the differences between the schemes, producing similar ice-sheet configurations through quantitatively different results of the sliding coefficient calibration. Despite this we observe a robust agreement in the reconstructed patterns of basal sliding parameters. We exchange the calibrated sliding parameter distributions between the schemes to demonstrate that the results of the model calibration cannot be straightforwardly transferred to models based on different approximations of ice dynamics. However, easily adaptable calibration techniques for the potential distribution of basal sliding coefficients can be implemented into ice models to overcome such incompatibility, as shown in this study

Description: In spite of significant progress in paleoclimate reconstructions and modelling of different aspects of the past glacial cycles, the mechanisms which transform regional and seasonal variations in solar insolation into long-term and global-scale glacial–interglacial cycles are still not fully understood – in particular, in relation to CO2 variability. Here using the Earth system model of intermediate complexity CLIMBER-2 we performed simulations of the co-evolution of climate, ice sheets, and carbon cycle over the last 400000 years using the orbital forcing as the only external forcing. The model simulates temporal dynamics of CO2, global ice volume, and other climate system characteristics in good agreement with paleoclimate reconstructions. These results provide strong support for the idea that long and strongly asymmetric glacial cycles of the late Quaternary represent a direct but strongly nonlinear response of the Northern Hemisphere ice sheets to orbital forcing. This response is strongly amplified and globalised by the carbon cycle feedbacks. Using simulations performed with the model in different configurations, we also analyse the role of individual processes and sensitivity to the choice of model parameters. While many features of simulated glacial cycles are rather robust, some details of CO2 evolution, especially during glacial terminations, are sensitive to the choice of model parameters. Specifically, we found two major regimes of CO2 changes during terminations: in the first one, when the recovery of the Atlantic meridional overturning circulation (AMOC) occurs only at the end of the termination, a pronounced overshoot in CO2 concentration occurs at the beginning of the interglacial and CO2 remains almost constant during the interglacial or even declines towards the end, resembling Eemian CO2 dynamics. However, if the recovery of the AMOC occurs in the middle of the glacial termination, CO2 concentration continues to rise during the interglacial, similar to the Holocene. We also discuss the potential contribution of the brine rejection mechanism for the CO2 and carbon isotopes in the atmosphere and the ocean during the past glacial termination.

Description: The composition of planktonic foraminiferal (PF) calcite is routinely used to reconstruct climate variability. However, PF ecology leaves a large imprint on the proxy signal: seasonal and vertical habitats of PF species vary spatially, causing variable offsets from annual mean surface conditions recorded by sedimentary assemblages. PF seasonality changes with temperature in a way that minimises the environmental change that individual species experience and it is not unlikely that changes in depth habitat also result from such habitat tracking. While this behaviour could lead to an underestimation of spatial or temporal trends as well as of variability in proxy records, most palaeoceanographic studies are (implicitly) based on the assumption of a constant habitat. Up to now, the effect of habitat tracking on foraminifera proxy records has not yet been formally quantified on a global scale. Here we attempt to characterise this effect on the amplitude of environmental change recorded in sedimentary PF using core top δ18O data from six species. We find that the offset from mean annual near-surface δ18O values varies with temperature, with PF δ18O indicating warmer than mean conditions in colder waters (on average by −0.1 ‰ (equivalent to 0.4 °C) per °C), thus providing a first-order quantification of the degree of underestimation due to habitat tracking. We use an empirical model to estimate the contribution of seasonality to the observed difference between PF and annual mean δ18O and use the residual Δδ18O to assess trends in calcification depth. Our analysis indicates that given an observation-based model parametrisation calcification depth increases with temperature in all species and sensitivity analysis suggests that a temperature-related seasonal habitat adjustment is essential to explain the observed isotope signal. Habitat tracking can thus lead to a significant reduction in the amplitude of recorded environmental change. However, we show that this behaviour is predictable. This allows accounting for habitat tracking, enabling more meaningful reconstructions and improved data–model comparison.

Description: Large-scale atmospheric patterns are examined on orbital timescales using a climate model which explicitly resolves the atmosphere–ocean–sea ice dynamics. It is shown that, in contrast to boreal summer where the climate mainly follows the local radiative forcing, the boreal winter climate is strongly determined by modulation of circulation modes linked to the Arctic Oscillation/North Atlantic Oscillation (AO/NAO) and the Northern/Southern Annular Modes. We find that during a positive phase of the AO/NAO the convection in the tropical Pacific is below normal. The related atmospheric circulation provides an atmospheric bridge for the precessional forcing inducing a non-uniform temperature anomalies with large amplitudes over the continents. We argue that this is important for mechanisms responsible for multi-millennial climate variability and glacial inception.

Description: Coccolithophores are calcifying phytoplankton and major contributors to both the organic and inorganic oceanic carbon pumps. Their export fluxes, species composition, and seasonal patterns were determined in two sediment trap moorings (M4 at 12° N, 49° W and M2 at 14° N, 37° W) collecting settling particles synchronously from October 2012 to November 2013 at 1200 m of water depth in the open equatorial North Atlantic. The two trap locations showed a similar seasonal pattern in total coccolith export fluxes and a predominantly tropical coccolithophore settling assemblage. Species fluxes were dominated throughout the year by lower photic zone (LPZ) taxa (Florisphaera profunda, Gladiolithus flabellatus) but also included upper photic zone (UPZ) taxa (Umbellosphaera spp., Rhabdosphaera spp., Umbilicosphaera spp., Helicosphaera spp.). The LPZ flora was most abundant during fall 2012, whereas the UPZ flora was more important during summer. In spite of these similarities, the western part of the study area produced persistently higher fluxes, averaging 241×107 ± 76×107 coccoliths m−2 d−1 at station M4 compared to only 66×107 ± 31×107 coccoliths m−2 d−1 at station M2. Higher fluxes at M4 were mainly produced by the LPZ species, favoured by the westward deepening of the thermocline and nutricline. Still, most UPZ species also contributed to higher fluxes, reflecting enhanced productivity in the western equatorial North Atlantic. Such was the case of two marked flux peaks of the more opportunistic species Gephyrocapsa muellerae and Emiliania huxleyi in January and April 2013 at M4, indicating a fast response to the nutrient enrichment of the UPZ, probably by wind-forced mixing. Later, increased fluxes of G. oceanica and E. huxleyi in October–November 2013 coincided with the occurrence of Amazon-River-affected surface waters. Since the spring and fall events of 2013 were also accompanied by two dust flux peaks, we propose a scenario in which atmospheric dust also provided fertilizing nutrients to this area. Enhanced surface buoyancy associated with the river plume indicates that the Amazon acted not only as a nutrient source, but also as a surface density retainer for nutrients supplied from the atmosphere. Nevertheless, lower total coccolith fluxes during these events compared to the maxima recorded in November 2012 and July 2013 indicate that transient productivity by opportunistic species was less important than "background" tropical productivity in the equatorial North Atlantic. This study illustrates how two apparently similar sites in the tropical open ocean actually differ greatly in ecological and oceanographic terms. The results presented here provide valuable insights into the processes governing the ecological dynamics and the downward export of coccolithophores in the tropical North Atlantic.

Description: Little is known about the climate evolution on the Kamchatka Peninsula during the last glacial–interglacial transition as existing climate records do not reach beyond 12 ka BP. In this study, a summer temperature record for the past 20 kyr is presented. Branched glycerol dialkyl glycerol tetraethers, terrigenous biomarkers suitable for continental air temperature reconstructions, were analyzed in a sediment core from the western continental margin off Kamchatka in the marginal northwest Pacific (NW Pacific). The record suggests that summer temperatures on Kamchatka during the Last Glacial Maximum (LGM) equaled modern temperatures. We suggest that strong southerly winds associated with a pronounced North Pacific High pressure system over the subarctic NW Pacific accounted for the warm conditions. A comparison with an Earth system model reveals discrepancies between model and proxy-based reconstructions for the LGM temperature and atmospheric circulation in the NW Pacific realm. The deglacial temperature development is characterized by abrupt millennial-scale temperature oscillations. The Bølling–Allerød warm phase and the Younger Dryas cold spell are pronounced events, suggesting a connection to North Atlantic climate variability.

Description: The oceans absorb about a quarter of the annually produced anthropogenic atmospheric carbon dioxide (CO2), resulting in a decrease in surface water pH, a process termed ocean acidification (OA). Surprisingly little is known about how OA affects the physiology of heterotrophic bacteria or the coupling of heterotrophic bacteria to phytoplankton when nutrients are limited. Previous experiments were, for the most part, undertaken during productive phases or following nutrient additions designed to stimulate algal blooms. Therefore, we performed an in situ large-volume mesocosm ( ∼  55 m3) experiment in the Baltic Sea by simulating different fugacities of CO2 (fCO2) extending from present to future conditions. The study was conducted in July–August after the nominal spring bloom, in order to maintain low-nutrient conditions throughout the experiment. This resulted in phytoplankton communities dominated by small-sized functional groups (picophytoplankton). There was no consistent fCO2-induced effect on bacterial protein production (BPP), cell-specific BPP (csBPP) or biovolumes (BVs) of either free-living (FL) or particle-associated (PA) heterotrophic bacteria, when considered as individual components (univariate analyses). Permutational Multivariate Analysis of Variance (PERMANOVA) revealed a significant effect of the fCO2 treatment on entire assemblages of dissolved and particulate nutrients, metabolic parameters and the bacteria–phytoplankton community. However, distance-based linear modelling only identified fCO2 as a factor explaining the variability observed amongst the microbial community composition, but not for explaining variability within the metabolic parameters. This suggests that fCO2 impacts on microbial metabolic parameters occurred indirectly through varying physicochemical parameters and microbial species composition. Cluster analyses examining the co-occurrence of different functional groups of bacteria and phytoplankton further revealed a separation of the four fCO2-treated mesocosms from both control mesocosms, indicating that complex trophic interactions might be altered in a future acidified ocean. Possible consequences for nutrient cycling and carbon export are still largely unknown, in particular in a nutrient-limited ocean.

Description: The large variety of atmospheric circulation systems affecting the eastern Asian climate is reflected by the complex Asian vegetation distribution. Particularly in the transition zones of these circulation systems, vegetation is supposed to be very sensitive to climate change. Since proxy records are scarce, hitherto a mechanistic understanding of the past spatio-temporal climate–vegetation relationship is lacking. To assess the Holocene vegetation change and to obtain an ensemble of potential mid-Holocene biome distributions for eastern Asia, we forced the diagnostic biome model BIOME4 with climate anomalies of different transient Holocene climate simulations performed in coupled atmosphere–ocean(–vegetation) models. The simulated biome changes are compared with pollen-based biome records for different key regions. In all simulations, substantial biome shifts during the last 6000 years are confined to the high northern latitudes and the monsoon–westerly wind transition zone, but the temporal evolution and amplitude of change strongly depend on the climate forcing. Large parts of the southern tundra are replaced by taiga during the mid-Holocene due to a warmer growing season and the boreal treeline in northern Asia is shifted northward by approx. 4° in the ensemble mean, ranging from 1.5 to 6° in the individual simulations, respectively. This simulated treeline shift is in agreement with pollen-based reconstructions from northern Siberia. The desert fraction in the transition zone is reduced by 21 % during the mid-Holocene compared to pre-industrial due to enhanced precipitation. The desert–steppe margin is shifted westward by 5° (1–9° in the individual simulations). The forest biomes are expanded north-westward by 2°, ranging from 0 to 4° in the single simulations. These results corroborate pollen-based reconstructions indicating an extended forest area in north-central China during the mid-Holocene. According to the model, the forest-to-non-forest and steppe-to-desert changes in the climate transition zones are spatially not uniform and not linear since the mid-Holocene.

Description: The Regional Antarctic ice and Global Ocean (RAnGO) model has been developed to study the interaction between the world ocean and the Antarctic ice sheet. The coupled model is based on a global implementation of the Finite Element Sea-ice Ocean Model (FESOM) with a mesh refinement in the Southern Ocean, particularly in its marginal seas and in the sub-ice-shelf cavities. The cryosphere is represented by a regional setup of the ice flow model RIMBAY comprising the Filchner–Ronne Ice Shelf and the grounded ice in its catchment area up to the ice divides. At the base of the RIMBAY ice shelf, melt rates from FESOM's ice-shelf component are supplied. RIMBAY returns ice thickness and the position of the grounding line. The ocean model uses a pre-computed mesh to allow for an easy adjustment of the model domain to a varying cavity geometry. RAnGO simulations with a 20th-century climate forcing yield realistic basal melt rates and a quasi-stable grounding line position close to the presently observed state. In a centennial-scale warm-water-inflow scenario, the model suggests a substantial thinning of the ice shelf and a local retreat of the grounding line. The potentially negative feedback from ice-shelf thinning through a rising in situ freezing temperature is more than outweighed by the increasing water column thickness in the deepest parts of the cavity. Compared to a control simulation with fixed ice-shelf geometry, the coupled model thus yields a slightly stronger increase in ice-shelf basal melt rates.

Description: Two interglacial epochs are included in the suite of Paleoclimate Modeling Intercomparison Project (PMIP4) simulations in the Coupled Model Intercomparison Project (CMIP6). The experimental protocols for simulations of the mid-Holocene (midHolocene, 6000 years before present) and the Last Interglacial (lig127k, 127 000 years before present) are described here. These equilibrium simulations are designed to examine the impact of changes in orbital forcing at times when atmospheric greenhouse gas levels were similar to those of the preindustrial period and the continental configurations were almost identical to modern ones. These simulations test our understanding of the interplay between radiative forcing and atmospheric circulation, and the connections among large-scale and regional climate changes giving rise to phenomena such as land–sea contrast and high-latitude amplification in temperature changes, and responses of the monsoons, as compared to today. They also provide an opportunity, through carefully designed additional sensitivity experiments, to quantify the strength of atmosphere, ocean, cryosphere, and land-surface feedbacks. Sensitivity experiments are proposed to investigate the role of freshwater forcing in triggering abrupt climate changes within interglacial epochs. These feedback experiments naturally lead to a focus on climate evolution during interglacial periods, which will be examined through transient experiments. Analyses of the sensitivity simulations will also focus on interactions between extratropical and tropical circulation, and the relationship between changes in mean climate state and climate variability on annual to multi-decadal timescales. The comparative abundance of paleoenvironmental data and of quantitative climate reconstructions for the Holocene and Last Interglacial make these two epochs ideal candidates for systematic evaluation of model performance, and such comparisons will shed new light on the importance of external feedbacks (e.g., vegetation, dust) and the ability of state-of-the-art models to simulate climate changes realistically.

Description: Polycystina (~400–800 living species and several thousand extinct forms) and Phaeodaria (~400–500 living species) are exclusively marine, open-ocean planktonic protists, most of which possess elaborate siliceous skeletons. The cytoplasm is divided into an internal part (endoplasm) separated from the external, more vacuolated one (ectoplasm) by a perforated membrane – the central capsule. The Polycystina protrude long and slender cytoplasmic projections (axopodia) supported internally by a rigid central rod (axoneme); while the Phaeodria have anetwork ofperipheral finely interconnectedpseudopodia.Afew Polycystina are colonial, but most, as well as all Phaeodaria, are solitary, around 40 μm to almost 2 mm in size. Most polycystine species peak in abundance between 0 and 100 m, whereasphaeodarianstendtolivedeeper,oftenbelow300m.Polycystineshavea rich fossil record dating from the Cambrian and are important for stratigraphic, paleoecologic, and evolutionary studies. The world-wide biogeography and diversity of radiolarians is chiefly governed by water temperature. Radiolarian prey includes bacteria, algae, protozoa, and microinvertebrates. Many surfacedwelling species of Polycystina possess symbiotic algae and photosynthetic cyanobacteria that provide nourishment to the host. Some colonial radiolaria reproduce by binary fission of the central capsules. Sexual reproduction of polycystines or Phaeodaria has not been confirmed, but the release of motile swarmers, likely gametes, has been widely documented. In species with a radial symmetry (Spumellaria) shell-growth is centrifugal, whereas in the Nassellaria the internal cephalic elements and the cephalis appear first. Individual longevity is estimated to range between 2 and 3 weeks and 1–2 months.

Description: We present the first results of the implementation of stable water isotopes in the Massachusetts Institute of Technology general circulation model (MITgcm). The model is forced with the isotopic content of precipitation and water vapor from an atmospheric general circulation model (NCAR IsoCAM), while the fractionation during evaporation is treated explicitly in the MITgcm. Results of the equilibrium simulation under pre-industrial conditions are compared to observational data and measurements of plankton tow records (the oxygen isotopic composition of planktic foraminiferal calcite). The broad patterns and magnitude of the stable water isotopes in annual mean seawater are well captured in the model, both at the sea surface as well as in the deep ocean. However, the surface water in the Arctic Ocean is not depleted enough, due to the absence of highly depleted precipitation and snowfall. A model–data mismatch is also recognizable in the isotopic composition of the seawater–salinity relationship in midlatitudes that is mainly caused by the coarse grid resolution. Deep-ocean characteristics of the vertical water mass distribution in the Atlantic Ocean closely resemble observational data. The reconstructed δ18Oc at the sea surface shows a good agreement with measurements. However, the model–data fit is weaker when individual species are considered and deviations are most likely attributable to the habitat depth of the foraminifera. Overall, the newly developed stable water isotope package opens wide prospects for long-term simulations in a paleoclimatic context.

Description: Here we outline how to identify hydrogenase enzymes from metagenomic libraries through an activity-based screening approach. A metagenomic fosmid library is constructed in E. coli and the fosmids are transferred into a hydrogenase deletion mutant of Shewanella oneidensis (ΔhyaB) via triparental mating. If a fosmid exhibits hydrogen uptake activity, S. oneidensis' phenotype is restored and hydrogenase activity is indicated by a color change of the medium from yellow to colorless. This new method enables screening of 48 metagenomic fosmid clones in parallel.

Description: Wetlands are one of the most significant natural sources of methane (CH4) to the atmosphere. They emit CH4 because decomposition of soil organic matter in waterlogged anoxic conditions produces CH4, in addition to carbon dioxide (CO2). Production of CH4 and how much of it escapes to the atmosphere depend on a multitude of environmental drivers. Models simulating the processes leading to CH4 emissions are thus needed for upscaling observations to estimate present CH4 emissions and for producing scenarios of future atmospheric CH4 concentrations. Aiming at a CH4 model that can be added to models describing peatland carbon cycling, we composed a model called HIMMELI that describes CH4 build-up in and emissions from peatland soils. It is not a full peatland carbon cycle model but it requires the rate of anoxic soil respiration as input. Driven by soil temperature, leaf area index (LAI) of aerenchymatous peatland vegetation, and water table depth (WTD), it simulates the concentrations and transport of CH4, CO2, and oxygen (O2) in a layered one-dimensional peat column. Here, we present the HIMMELI model structure and results of tests on the model sensitivity to the input data and to the description of the peat column (peat depth and layer thickness), and demonstrate that HIMMELI outputs realistic fluxes by comparing modeled and measured fluxes at two peatland sites. As HIMMELI describes only the CH4-related processes, not the full carbon cycle, our analysis revealed mechanisms and dependencies that may remain hidden when testing CH4 models connected to complete peatland carbon models, which is usually the case. Our results indicated that (1) the model is flexible and robust and thus suitable for different environments; (2) the simulated CH4 emissions largely depend on the prescribed rate of anoxic respiration; (3) the sensitivity of the total CH4 emission to other input variables is mainly mediated via the concentrations of dissolved gases, in particular, the O2 concentrations that affect the CH4 production and oxidation rates; (4) with given input respiration, the peat column description does not significantly affect the simulated CH4 emissions in this model version.

Description: Mesocosm experiments on phytoplankton dynamics under high CO2 concentrations mimic the response of marine primary producers to future ocean acidification. However, potential acidification effects can be hindered by the high standard deviation typically found in the replicates of the same CO2 treatment level. In experiments with multiple unresolved factors and a sub-optimal number of replicates, post-processing statistical inference tools might fail to detect an effect that is present. We propose that in such cases, data-based model analyses might be suitable tools to unearth potential responses to the treatment and identify the uncertainties that could produce the observed variability. As test cases, we used data from two independent mesocosm experiments. Both experiments showed high standard deviations and, according to statistical inference tools, biomass appeared insensitive to changing CO2 conditions. Conversely, our simulations showed earlier and more intense phytoplankton blooms in modeled replicates at high CO2 concentrations and suggested that uncertainties in average cell size, phytoplankton biomass losses, and initial nutrient concentration potentially outweigh acidification effects by triggering strong variability during the bloom phase. We also estimated the thresholds below which uncertainties do not escalate to high variability. This information might help in designing future mesocosm experiments and interpreting controversial results on the effect of acidification or other pressures on ecosystem functions

Description: Following the recent Global Carbon Project (GCP) synthesis of the decadal methane (CH4) budget over 2000–2012 (Saunois et al., 2016), we analyse here the same dataset with a focus on quasi-decadal and inter-annual variability in CH4 emissions. The GCP dataset integrates results from top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models (including process-based models for estimating land surface emissions and atmospheric chemistry), inventories of anthropogenic emissions, and data-driven approaches. The annual global methane emissions from top-down studies, which by construction match the observed methane growth rate within their uncertainties, all show an increase in total methane emissions over the period 2000–2012, but this increase is not linear over the 13 years. Despite differences between individual studies, the mean emission anomaly of the top-down ensemble shows no significant trend in total methane emissions over the period 2000–2006, during the plateau of atmospheric methane mole fractions, and also over the period 2008–2012, during the renewed atmospheric methane increase. However, the top-down ensemble mean produces an emission shift between 2006 and 2008, leading to 22 [16–32] Tg CH4 yr−1 higher methane emissions over the period 2008–2012 compared to 2002–2006. This emission increase mostly originated from the tropics, with a smaller contribution from mid-latitudes and no significant change from boreal regions. The regional contributions remain uncertain in top-down studies. Tropical South America and South and East Asia seem to contribute the most to the emission increase in the tropics. However, these two regions have only limited atmospheric measurements and remain therefore poorly constrained. The sectorial partitioning of this emission increase between the periods 2002–2006 and 2008–2012 differs from one atmospheric inversion study to another. However, all top-down studies suggest smaller changes in fossil fuel emissions (from oil, gas, and coal industries) compared to the mean of the bottom-up inventories included in this study. This difference is partly driven by a smaller emission change in China from the top-down studies compared to the estimate in the Emission Database for Global Atmospheric Research (EDGARv4.2) inventory, which should be revised to smaller values in a near future. We apply isotopic signatures to the emission changes estimated for individual studies based on five emission sectors and find that for six individual top-down studies (out of eight) the average isotopic signature of the emission changes is not consistent with the observed change in atmospheric 13CH4. However, the partitioning in emission change derived from the ensemble mean is consistent with this isotopic constraint. At the global scale, the top-down ensemble mean suggests that the dominant contribution to the resumed atmospheric CH4 growth after 2006 comes from microbial sources (more from agriculture and waste sectors than from natural wetlands), with an uncertain but smaller contribution from fossil CH4 emissions. In addition, a decrease in biomass burning emissions (in agreement with the biomass burning emission databases) makes the balance of sources consistent with atmospheric 13CH4 observations. In most of the top-down studies included here, OH concentrations are considered constant over the years (seasonal variations but without any inter-annual variability). As a result, the methane loss (in particular through OH oxidation) varies mainly through the change in methane concentrations and not its oxidants. For these reasons, changes in the methane loss could not be properly investigated in this study, although it may play a significant role in the recent atmospheric methane changes as briefly discussed at the end of the paper.

Description: Continuous records of the atmospheric greenhouse gases (GHGs) CO2, CH4, and N2O are necessary input data for transient climate simulations, and their associated radiative forcing represents important components in analyses of climate sensitivity and feedbacks. Since the available data from ice cores are discontinuous and partly ambiguous, a well-documented decision process during data compilation followed by some interpolating post-processing is necessary to obtain those desired time series. Here, we document our best possible data compilation of published ice core records and recent measurements on firn air and atmospheric samples spanning the interval from the penultimate glacial maximum ( ∼  156 kyr BP) to the beginning of the year 2016 CE. We use the most recent age scales for the ice core data and apply a smoothing spline method to translate the discrete and irregularly spaced data points into continuous time series. These splines are then used to compute the radiative forcing for each GHG using well-established, simple formulations. We compile only a Southern Hemisphere record of CH4 and discuss how much larger a Northern Hemisphere or global CH4 record might have been due to its interpolar difference. The uncertainties of the individual data points are considered in the spline procedure. Based on the given data resolution, time-dependent cutoff periods of the spline, defining the degree of smoothing, are prescribed, ranging from 5000 years for the less resolved older parts of the records to 4 years for the densely sampled recent years. The computed splines seamlessly describe the GHG evolution on orbital and millennial timescales for glacial and glacial–interglacial variations and on centennial and decadal timescales for anthropogenic times. Data connected with this paper, including raw data and final splines, are available at doi:10.1594/PANGAEA.871273.

Description: Basalt weathering is one of many relevant processes balancing the global carbon cycle via land–ocean alkalinity fluxes. The CO2 consumption by weathering can be calculated using alkalinity and is often scaled with runoff and/or temperature. Here, it is tested if the surface age distribution of a volcanic system derived by geological maps is a useful proxy for changes in alkalinity production with time. A linear relationship between temperature normalized alkalinity fluxes and the Holocene area fraction of a volcanic field was identified using information from 33 basalt volcanic fields, with an r2=0.93. This relationship is interpreted as an aging function and suggests that fluxes from Holocene areas are ∼10 times higher than those from old inactive volcanic fields. However, the cause for the decrease with time is probably a combination of effects, including a decrease in alkalinity production from material in the shallow critical zone as well as a decline in hydrothermal activity and magmatic CO2 contribution. The addition of fresh reactive material on top of the critical zone has an effect in young active volcanic settings which should be accounted for, too. A comparison with global models suggests that global alkalinity fluxes considering Holocene basalt areas are ∼60 % higher than the average from these models imply. The contribution of Holocene areas to the global basalt alkalinity fluxes is today however only ∼5 %, because identified, mapped Holocene basalt areas cover only ∼1 % of the existing basalt areas. The large trap basalt proportion on the global basalt areas today reduces the relevance of the aging effect. However, the aging effect might be a relevant process during periods of globally intensive volcanic activity, which remains to be tested.

Description: Dynamic vegetation models simulate global vegetation in terms of fractional coverage of a few plant functional types (PFTs). Although these models often share the same concept, they differ with respect to the number and kind of PFTs, complicating the comparability of simulated vegetation distributions. Pollen-based vegetation reconstructions are initially only available in the form of time series of individual taxa that are not distinguished in the models. Thus, to evaluate simulated vegetation distributions, the modelling results and pollen-based vegetation reconstructions have to be converted into a comparable format. The classical approach is the method of biomisation, but hitherto PFT-based biomisation methods were only available for individual models. We introduce and evaluate a simple, universally applicable technique to harmonise PFT distributions by assigning them into nine mega-biomes, using only assumptions on the minimum PFT cover fractions and few bioclimatic constraints (based on the 2 m temperature). These constraints mainly follow the limitation rules used in the classical biome models (here BIOME4). We test the method for six state-of-the-art dynamic vegetation models that are included in Earth system models based on pre-industrial, mid-Holocene and Last Glacial Maximum simulations. The method works well, independent of the spatial resolution or the complexity of the models. Large biome belts (such as tropical forest) are generally better represented than regionally confined biomes (warm–temperate forest, savanna). The comparison with biome distributions inferred via the classical biomisation approach of forcing biome models (here BIOME1) with the simulated climate states shows that the PFT-based biomisation is even able to keep up with the classical method. However, as the new method considers the PFT distributions actually calculated by the Earth system models, it allows for a direct comparison and evaluation of simulated vegetation distributions which the classical method cannot do. Thereby, the new method provides a powerful tool for the evaluation of Earth system models in general.

Description: The sedimentary stable nitrogen isotope compositions of bulk organic matter (δ15Nbulk) and silicon isotope composition of diatoms (δ30SiBSi) both mainly reflect the degree of past nutrient utilization by primary producers. However, in ocean areas where anoxic and suboxic conditions prevail, the δ15Nbulk signal ultimately recorded within the sediments is also influenced by water column denitrification causing an increase in the subsurface δ15N signature of dissolved nitrate (δ15NO3−) upwelled to the surface. Such conditions are found in the oxygen minimum zone off Peru, where at present an increase in subsurface δ15NO3− from North to South along the shelf is observed due to ongoing denitrification within the pole-ward flowing subsurface waters, while the δ30Si signature of silicic acid (δ30Si(OH)4) at the same time remains unchanged. Here, we present three new δ30SiBSi records between 11° S and 15° S and compare these to previously published δ30SiBSi and δ15Nbulk records from Peru covering the past 600 years. We present a new approach to calculate past subsurface δ15NO3− signatures based on the correlation of δ30SiBSi and δ15Nbulk signatures at a latitudinal resolution for different time periods. Our results show source water δ15NO3− compositions during the last 200 years, the Current Warm Period (CWP) and during short-term arid events prior to that, which are close to modern values increasing southward from 7 to 10 ‰ (between 11° S and 15° S). In contrast, humid conditions during the Little Ice Age (LIA) reflect consistently low δ15NO3− values between 6 and 7.5‰. Furthermore, we are able to relate the short-term variability in both isotope compositions to changes in the ratio of nutrients (NO3− : Si(OH)4) taken up by different dominating phytoplankton groups (diatoms and non-siliceous phytoplankton) under the variable climatic conditions of the past 600 years.

Description: Differences in habitat and diet between species are often associated with morphological differences. Habitat and trophic adaptation have therefore been proposed as important drivers of speciation and adaptive radiation. Importantly, habitat and diet shifts likely impose changes in exposure to different parasites and infection risk. As strong selective agents influencing survival and mate choice, parasites might play an important role in host diversification. We explore this possibility for the adaptive radiation of Lake Tanganyika (LT) cichlids. We first compare metazoan macroparasites infection levels between cichlid tribes. We then describe the cichlids’ genetic diversity at the major histocompatibility complex (MHC), which plays a key role in vertebrate immunity. Finally, we evaluate to what extent trophic ecology and morphology explain variation in infection levels and MHC, accounting for phylogenetic relationships. We show that different cichlid tribes in LT feature partially non-overlapping parasite communities and partially non-overlapping MHC diversity. While morphology explained 15% of the variation in mean parasite abundance, trophic ecology accounted for 16% and 22% of the MHC variation at the nucleotide and at the amino acid level, respectively. Parasitism and immunogenetic adaptation may thus add additional dimensions to the LT cichlid radiation.

Description: The open ocean is a major source of nitrous oxide (N2O), an atmospheric trace gas attributable to global warming and ozone depletion. Intense sea-to-air N2O fluxes occur in major oceanic upwelling regions such as the eastern tropical South Pacific (ETSP). The ETSP is influenced by the El Niño–Southern Oscillation that leads to inter-annual variations in physical, chemical, and biological properties in the water column. In October 2015, a strong El Niño event was developing in the ETSP; we conduct field observations to investigate (1) the N2O production pathways and associated biogeochemical properties and (2) the effects of El Niño on water column N2O distributions and fluxes using data from previous non-El Niño years. Analysis of N2O natural abundance isotopomers suggested that nitrification and partial denitrification (nitrate and nitrite reduction to N2O) were occurring in the near-surface waters; indicating that both pathways contributed to N2O effluxes. Higher-than-normal sea surface temperatures were associated with a deepening of the oxycline and the oxygen minimum layer. Within the shelf region, surface N2O supersaturation was nearly an order of magnitude lower than that of non-El Niño years. Therefore, a significant reduction of N2O efflux (75 %–95 %) in the ETSP occurred during the 2015 El Niño. At both offshore and coastal stations, the N2O concentration profiles during El Niño showed moderate N2O concentration gradients, and the peak N2O concentrations occurred at deeper depths during El Niño years; this was likely the result of suppressed upwelling retaining N2O in subsurface waters. At multiple stations, water-column inventories of N2O within the top 1000 m were up to 160 % higher than those measured in non-El Niño years, indicating that subsurface N2O during El Niño could be a reservoir for intense N2O effluxes when normal upwelling is resumed after El Niño.

Description: Common problems in state-of-the-art climate models are a cold sea surface temperature (SST) bias in the equatorial Pacific and the underestimation of the two most important atmospheric feedbacks operating in the El Niño/Southern Oscillation (ENSO): the positive, i.e. amplifying wind-SST feedback and the negative, i.e. damping heat flux-SST feedback. To a large extent, the underestimation of those feedbacks can be explained by the cold equatorial SST bias, which shifts the rising branch of the Pacific Walker Circulation (PWC) too far to the west by up to 30°, resulting in an erroneous convective response during ENSO events. Based on simulations from the Kiel Climate Model (KCM) and the 5th phase of Coupled Model Intercomparison Project (CMIP5), we investigate how well ENSO dynamics are simulated in case of underestimated ENSO atmospheric feedbacks (EAF), with a special focus on ocean–atmosphere coupling over the equatorial Pacific. While models featuring realistic atmospheric feedbacks simulate ENSO dynamics close to observations, models with underestimated EAF exhibit fundamental biases in ENSO dynamics. In models with too weak feedbacks, ENSO is not predominantly wind-driven as observed; instead ENSO is driven significantly by a positive shortwave radiation feedback. Thus, although these models simulate ENSO, which in terms of simple indices is consistent with observations, it originates from very different dynamics. A too weak oceanic forcing on the SST via the positive thermocline, the Ekman and the zonal advection feedback is compensated by weaker atmospheric heat flux damping. The latter is mainly caused by a biased shortwave-SST feedback that erroneously is positive in most climate models. In the most biased models, the shortwave-SST feedback contributes to the SST anomaly growth to a similar degree as the ocean circulation. Our results suggest that a broad continuum of ENSO dynamics can exist in climate models and explain why climate models with less than a half of the observed EAF strength can still depict realistic ENSO amplitude.

Description: Diseases increasingly threaten aquaculture of kelps and other seaweeds. At the same time, protection concepts that are based upon application of biocides are usually not applicable, as such compounds would be rapidly diluted in the sea, causing ecological damage. An alternative concept could be the application of immune stimulants to prevent and control diseases in farmed seaweeds. We here present a pilot study that investigated the effects of oligoalginate elicitation on juvenile and adult sporophytes of Saccharina japonica cultivated in China and on adult sporophytes of Saccharina latissima cultivated in Germany. In two consecutive years, treatment with oligoalginate clearly reduced the detachment of S. japonica juveniles from their substrate curtains during the nursery stage in greenhouse ponds. Oligoalginate elicitation also decreased the density of endobionts and the number of bacterial cells on sporophytes of S. latissima that were cultivated on sea-based rafts. However, the treatment increased the susceptibility of kelp adults to settlement of epibionts (barnacles in Germany and filamentous algal epiphytes in China). In addition, oligoalginate elicitation accelerated the aging of S. japonica adults. Based upon these findings, oligoalginate elicitation could be a feasible way to provide “environmentally friendly” protection of kelp juveniles in nurseries. The same treatment causes not only beneficial, but also unwanted effects in adult kelp sporophytes. Therefore, it is not recommended as a treatment after the juvenile stage is completed. Future tests with other elicitors and other cultivated seaweed species may allow for the development of more feasible applications of targeted defense elicitation in seaweed aquaculture.

Description: The early Eocene (56 to 48 million years ago) is inferred to have been the most recent time that Earth's atmospheric CO2 concentrations exceeded 1000 ppm. Global mean temperatures were also substantially warmer than present day. As such, study of early Eocene climate provides insight into how a super-warm Earth system behaves and offers an opportunity to 10 evaluate climate models under conditions of high greenhouse gas forcing. The Deep Time Model Intercomparison Project (DeepMIP) is a systematic model-model and model-data intercomparison of three early Paleogene time slices: latest Paleocene, Paleocene-Eocene thermal maximum and early Eocene climatic optimum. A previous article outlined the model experimental design for climate model simulations. In this article, we outline the methodologies to be used for the compilation and analysis of climate proxy data, primarily proxies for temperature and CO2. This paper establishes the protocols for a concerted and 15 coordinated effort to compile the climate proxy records across a wide geographic range. The resulting climate "atlas" will be used to constrain and evaluate climate models for the three selected time intervals, and provide insights into the mechanisms that control these warm climate states. We provide version 0.1 of this database, in anticipation that this will be expanded in subsequent publications.

Description: The present study details the effects of basin-scale hydrographic characteristics of the Red Sea on the macroecology of Chaetognatha, a major plankton component in the pelagic realm. The hydrographic attributes and circulation of the Red Sea as a result of its limited connection with the northern Indian Ocean make it a unique ecohydrographic region in the world ocean. Here, we aimed to identify the prime determinants governing the community structure and vertical distribution of the Cheatognatha in this ecologically significant world ocean basin. The intrusion of Gulf of Aden Water influenced the Chaetognatha community composition in the south, whereas the overturning circulation altered their vertical distribution in the north. The existence of hypoxic waters (〈 100 µmol kg−1) at mid-depth also influenced their vertical distribution. The detailed evaluation of the responses of the different life stages of Chaetognatha revealed an increased susceptibility of adult individuals to hypoxic waters compared to immature stages. Higher oxygen demands of the adults for the egg and sperm production might have prevented them from inhabiting the oxygen-deficient mid-depth zones. The carbon and nitrogen content of the Copepoda and Chaetognatha communities and the quantification of the predation impact of Chaetognatha on Copepoda based on the feeding rate helped in corroborating the significant trophic link between these two prey–predator taxa. The observed influences of physical and chemical attributes on the distribution of Chaetognatha can be used as a model example for the role of the hydrography on the zooplankton community of the Red Sea.

Description: As coastal areas become increasingly vulnerable to climate change, the study of nearshore sediment textures along the littoral cell of the Medjerda delta in the Gulf of Tunis, southern Mediterranean coast can provide valuable information (i) on the origin (continental or marine) of the sediment, (ii) its transport direction, and (iii) constitutes an important tool in the assessment of coastal sensitivity. A total of 120 sediments samples underwent grain size analysis and statistic parameters have been calculated. These allowed the identification of five different Sedimentary Types (ST). Accordingly, using grain size indexes (i.e. Mz, SKI and Ku), Sediment Trend Analysis (STA) modeling tools were applied to define the seasonal sediment transport pathways throughout the nearshore of the Medjerda sedimentary cell. Results show that grain size distribution (GSD) and STA model pathways are determined by cross-shore geomorphology, location of the sediment-cell, seasonal incident wave and local terrestrial supply. The appearance in an atypical seabed location of the finer (Mo = 0.1 mm) and the coarser STs (Mo = 0.8 mm) can be indicative of human influence since the coarser particles are usually retained by dam structures. Moreover, the bimodality and the increased distribution of mud are also related to the seasonal incident wave winnowing of the historic deltaic plain submerged by the relative rise in sea level. The evolution of the sediment pattern towards a greater proportion of very fine grains indicates a deficit of sediment supply, particularly of the coarser grains, and demonstrates the coastal vulnerability of the Gulf of Tunis due to anthropic effects.

Description: Sea surface salinity is one of the most important parameters to reconstruct in paleoclimatology, reflecting amongst others the hydrological cycle, paleo-density, ice volume, and regional and global circulation of water masses. Recent culture studies and a Red Sea field study revealed a significant positive relation between salinity and Na incorporation within benthic and planktonic foraminiferal shells. However, these studies reported varying partitioning of Na between and within the same species. The latter could be associated with ontogenetic variations, most likely spine loss. Varying Na concentrations were observed in different parts of foraminiferal shells, with especially spines and regions close to the primary organic sheet being enriched in Na. In this study, we unravel the Na composition of different components of the planktonic foraminiferal shell wall using Electron Probe Micro Analysis (EPMA) and solution-ICP-MS. A model is presented to interpret EPMA data for spines and spine bases to quantitatively assess differences in composition and contribution to whole shell Na/Ca signals. The same model can also be applied to other spatial inhomogeneities observed in foraminiferal shell chemistry, like elemental (e.g. Mg, Na, S) banding and/or hotspots. The relative contribution of shell calcite, organic linings, spines and spine bases to whole shell Na chemistry is considered quantitatively. This study shows that whereas the high Na areas may be susceptible to taphonomy, the Na chemistry of the shell itself seems relatively robust. Comparing both shell and spine Na/Ca values with salinity shows that shell chemistry records salinity, albeit with a very modest slope.

Description: A broad variety of materials of biological origin have been successfully used in recent decades for the removal of pollutants from waters. These biosorbents include natural polymers that play a key role for adsorption. It is therefore critical to understand the physicochemical properties of the chemical groups of these biopolymers. The acid–base properties of biomass are affected by pH, ionic strength and medium composition. Nevertheless, these parameters are not always considered during biosorption studies. According to the literature, less than 3% of biosorption reports include studies on proton binding. Moreover, in 60% of these papers, there is key experimental information missing such as the calibration of the electrodes employed for potentiometric titrations. We consider therefore that there is an important need for reviewing the role of proton binding in biosorption studies. This review outlines the major advances on data interpretation and modelling of proton binding on biosorbents. In addition, we discuss issues concerning the acid–base properties of biosorbents.

Description: Hydrothermal fluids on the modern seafloor are important carriers of base and precious metals in a wide range of volcanic and tectonic settings. The concentrations and distribution, especially of gold and silver, in associated seafloor massive sulfide (SMS) deposits are strongly influenced by variable source rocks, fluid chemistry, and precipitation mechanisms. Compositional data of 130 SMS deposits around the world show a large range of gold and silver grades, in part reflecting strong buffering of the hydrothermal fluids by their host rocks. Geochemical reaction-path modeling shows that in most cases the investigated hydrothermal fluids are undersaturated with gold and silver, and solubilities can be orders of magnitude higher than the Au and Ag concentrations measured in the corresponding fluids. Precipitation of gold during conductive cooling of mid-ocean ridge black smoker (MOR) fluids occurs at low temperatures but can be very rapid, with 〉 90% of the gold deposited in the first 25 °C of cooling below ~ 150 °C. The result is a Zn–Au polymetallic assemblage with Au and Ag deposited at the same time together with Pb and sulfosalts. In ultramafic-dominated (UM) systems, the strongly reduced hydrothermal fluids promote the deposition of gold at higher temperatures and explain the correlation between gold and copper in these deposits. In this case, the lower stability of the AuHS° complex at low ƒO2 (buffered by fayalite, magnetite, and quartz) results in gold deposition at 〉 250 °C with early bornite and chalcopyrite and before sphalerite and silver, producing a high-temperature Cu–Au assemblage. In sediment-hosted (SED) systems, the much higher pH stabilizes Au(HS)2− and keeps gold in solution to very low temperatures, after the precipitation of chalcopyrite, sphalerite, and galena, resulting in Au-poor polymetallic sulfides and very late-stage deposition of gold, commonly with amorphous silica. In arc-related (ARC) systems, gold deposition occurs at somewhat higher temperatures than in the MOR case, in part because the fluids start with higher gold concentrations. This can be explained by probable direct magmatic contributions, and the high ƒO2 of the fluids, which promotes the solubility of gold at the source. During cooling, gold precipitates at about 160 °C with sphalerite, tennantite, silver, and galena, resulting in an Au-rich polymetallic sulfide assemblage. The mixing of hydrothermal fluids with seawater generally causes oxidation and eventually a decrease in the pH at a mixing ratio of 1:1, causing an initial increase in the solubility of gold and silver. This can delay gold deposition from aqueous species to very low temperatures. These complex systematics make prediction of Au and Ag grades difficult. However, important new data are coming to light on the actual concentrations of the precious metals in hydrothermal fluids. In particular, the input of magmatic volatiles and leaching of pre-existing gold can lead to significant increases in the Au and Ag concentrations of the venting fluids and earlier deposition. In several cases, it appears that at least part of the gold load is present as nanoparticles in suspension, allowing bulk gold concentrations that may be far in excess of liquid saturation. Boiling at the seafloor is now widely observed, even at great water depths close to the critical point of seawater. Model calculations of phase separation during boiling show the competing effects on gold solubility of H2, H2S, and CO2 partitioning into the vapor, which can result in highly variable gold-to-base metal ratios in the deposits. Flashing of the vent fluids into steam at high temperatures is also commonly observed and can lead to spectacular Au grades, with a strong Cu–Au association in the deepest and hottest vents.

Description: The ecological approach to comparative cognition emphasizes that the ecological and social environment are important predictors of cognitive performance. We used this approach to test whether differences in habitat use and social behavior in the facultative Caribbean cleaning goby Elacatinus prochilos predict differences in learning performance in two discriminatory two-choice tasks. This species has two behavioral ecotypes: one that frequently engages in cleaning interactions and inhabits corals in male–female pairs (cleaning gobies) and another that rarely engages in cleaning interactions and inhabits barrel sponges in large groups (sponge-dwellers). We predicted that cleaning gobies would outperform sponge-dwellers in a pattern-cued task, which consisted of identifying the pattern on a plate that consistently provided food, while sponge-dwellers would outperform cleaning gobies in a spatial task, which consisted of identifying the location of the plate. Contrary to our predictions, there was no difference in performance between the two ecotypes. Most of the gobies performed poorly in the pattern-cued task and well in the spatial task. A possible explanation for these results is that the association of a pattern with positive and negative reinforcement may not be a pre-requisite for engaging in cleaning interactions, while spatial skills might be equally required in both ecotypes. Alternatively, the two ecotypes can flexibly adjust to new feeding conditions, which would explain their similar performance in the spatial task. Further research should investigate which aspects of E. prochilos’ social and ecological environment might impose challenges that require spatial cognition and whether individuals can flexibly adjust to new habitats and feeding conditions.

Description: Precipitation in California is modulated by variability in the tropical Pacific associated with El Nino/Southern Oscillation (ENSO): more rainfall is expected during El Nino episodes, and reduced rainfall during La Nina. It has been suggested that besides the shape and location of the sea surface temperature ( SST) anomaly this remote connection depends on the strength and location of the atmospheric convection response in the tropical Pacific. Here we show in a perturbed physics ensemble of the Kiel Climate Model and CMIP5 models that due to a cold equatorial SST bias many climate models are in a La Nina-like mean state, resulting in a too westward position of the rising branch of the Pacific Walker Circulation. This in turn results in a convective response along the equator during ENSO events that is too far west in comparison to observations. This effect of the equatorial cold SST bias is not restricted to the tropics, moreover it leads to a too westward SLP response in the North Pacific and too westward precipitation response that does not reach California. Further we show that climate models with a reduced equatorial cold SST bias have a more realistic representation of the spatial asymmetry of the teleconnections between El Nino and La Nina.

Description: Due to its remoteness, the deep-sea floor remains an understudied ecosystem of our planet. The patchiness of existing data sets makes it difficult to draw conclusions about processes that apply to a wider area. In our study we show how different settings and processes determine sediment heterogeneity on small spatial scales. We sampled solid phase and porewater from the upper 10 m of an approximately 7.4×13 km2 area in the Peru Basin, in the southeastern equatorial Pacific Ocean, at 4100 m water depth. Samples were analyzed for trace metals, including rare earth elements and yttrium (REY), as well as for particulate organic carbon (POC), CaCO3, and nitrate. The analyses revealed the surprisingly high spatial small-scale heterogeneity of the deep-sea sediment composition. While some cores have the typical green layer from Fe(II) in the clay minerals, this layer is missing in other cores, i.e., showing a tan color associated with more Fe(III) in the clay minerals. This is due to varying organic carbon contents: nitrate is depleted at 2–3 m depth in cores with higher total organic carbon contents but is present throughout cores with lower POC contents, thus inhibiting the Fe(III)-to-Fe(II) reduction pathway in organic matter degradation. REY show shale-normalized (SN) patterns similar to seawater, with a relative enrichment of heavy REY over light REY, positive LaSN anomaly, negative CeSN anomaly, and positive YSN anomaly and correlate with the Fe-rich clay layer and, in some cores, also correlate with P. We therefore propose that Fe-rich clay minerals, such as nontronite, as well as phosphates, are the REY-controlling phases in these sediments. Variability is also seen in dissolved Mn and Co concentrations between sites and within cores, which might be due to dissolving nodules in the suboxic sediment, as well as in concentration peaks of U, Mo, As, V, and Cu in two cores, which might be related to deposition of different material at lower-lying areas or precipitation due to shifting redox boundaries.

Description: Global climatologies of the seawater CO2 chemistry variables are necessary to assess the marine carbon cycle in depth. The climatologies should adequately capture seasonal variability to properly address ocean acidification and similar issues related to the carbon cycle. Total alkalinity (A(T)) is one variable of the seawater CO2 chemistry system involved in ocean acidification and frequently measured. We used the Global Ocean Data Analysis Project version 2.2019 (GLODAPv2) to extract relationships among the drivers of the A(T) variability and A(T) concentration using a neural network (NNGv2) to generate a monthly climatology. The GLODAPv2 quality-controlled dataset used was modeled by the NNGv2 with a root-mean-squared error (RMSE) of 5.3 mu mol kg(-1). Validation tests with independent datasets revealed the good generalization of the network. Data from five ocean time-series stations showed an acceptable RMSE range of 3-6.2 mu mol kg(-1). Successful modeling of the monthly A(T) variability in the time series suggests that the NNGv2 is a good candidate to generate a monthly climatology. The climatological fields of A(T) were obtained passing through the NNGv2 the World Ocean Atlas 2013 (WOA13) monthly climatologies of temperature, salinity, and oxygen and the computed climatologies of nutrients from the previous ones with a neural network. The spatiotemporal resolution is set by WOA13: 1 degrees x 1 degrees in the horizontal, 102 depth levels (0-5500 m) in the vertical and monthly (0-1500 m) to annual (1550-5500 m) temporal resolution. The product is distributed through the data repository of the Spanish National Research Council (CSIC; https://doi.org/10.20350/digitalCSIC/8644, Broullon et al., 2019).

Description: The flow (flux) of climate-critical gases, such as carbon dioxide (CO2), between the ocean and the atmosphere is a fundamental component of our climate and an important driver of the biogeochemical systems within the oceans. Therefore, the accurate calculation of these air–sea gas fluxes is critical if we are to monitor the oceans and assess the impact that these gases are having on Earth's climate and ecosystems. FluxEngine is an open-source software toolbox that allows users to easily perform calculations of air–sea gas fluxes from model, in situ, and Earth observation data. The original development and verification of the toolbox was described in a previous publication. The toolbox has now been considerably updated to allow for its use as a Python library, to enable simplified installation, to ensure verification of its installation, to enable the handling of multiple sparingly soluble gases, and to enable the greatly expanded functionality for supporting in situ dataset analyses. This new functionality for supporting in situ analyses includes user-defined grids, time periods and projections, the ability to reanalyse in situ CO2 data to a common temperature dataset, and the ability to easily calculate gas fluxes using in situ data from drifting buoys, fixed moorings, and research cruises. Here we describe these new capabilities and demonstrate their application through illustrative case studies. The first case study demonstrates the workflow for accurately calculating CO2 fluxes using in situ data from four research cruises from the Surface Ocean CO2 ATlas (SOCAT) database. The second case study calculates air–sea CO2 fluxes using in situ data from a fixed monitoring station in the Baltic Sea. The third case study focuses on nitrous oxide (N2O) and, through a user-defined gas transfer parameterisation, identifies that biological surfactants in the North Atlantic could suppress individual N2O sea–air gas fluxes by up to 13 %. The fourth and final case study illustrates how a dissipation-based gas transfer parameterisation can be implemented and used. The updated version of the toolbox (version 3) and all documentation is now freely available.

Description: Climate reconstructions are means to extract the signal from uncertain paleo-observations, so-called proxies. It is essential to evaluate these reconstructions to understand and quantify their uncertainties. Similarly, comparing climate simulations and proxies requires approaches to bridge the temporal and spatial differences between both and to address their specific uncertainties. One way to achieve these two goals is so-called pseudoproxies. These are surrogate proxy records within the virtual reality of a climate simulation. They in turn depend on an understanding of the uncertainties of the real proxies including the noise characteristics disturbing the original environmental signal. Common pseudoproxy approaches so far concentrate on data with high temporal resolution over the last approximately 2000 years. Here we provide a simple but flexible noise model for potentially low-resolution sedimentary climate proxies for temperature on millennial timescales, the code for calculating a set of pseudoproxies from a simulation, and one example of pseudoproxies. The noise model considers the influence of other environmental variables, a dependence on the climate state, a bias due to changing seasonality, modifications of the archive (for example bioturbation), potential sampling variability, and a measurement error. Model, code, and data allow us to develop new ways of comparing simulation data with proxies on long timescales. Code and data are available at https://doi.org/10.17605/OSF.IO/ZBEHX (Bothe et al., 2018).

Description: The spatial pattern of the first mode of interannual variability associated with the East Asian summer monsoon (EASM), obtained from a multivariate Empirical Orthogonal Functions (MV-EOF) analysis, corresponds to the Pacific–Japan (PJ) pattern and is referred to as the PJ-mode. The present study investigates the interannual variation of the PJ-mode from the perspective of the intraseasonal timescale. In particular, the impact of the Madden–Julian oscillation (MJO) on the interannual variation of the PJ-mode is investigated. The results show that the MJO has a significant influence on the interannual variation of the PJ-mode mainly in the lower troposphere (850 hPa) and that the former accounts for approximately 11% of the amplitude of the latter. The major part of the contribution comes from a change in frequency of the different phases of the MJO, especially that of MJO phase 6. This suggests that intraseasonal variation of the convection anomalies over the tropical eastern Indian and western Pacific Oceans plays an important role in the interannual variation of the PJ-mode. In addition, MJO phase 7 also contributes to the interannual variability of the PJ-mode, in this case induced by both the change in frequency and the change in circulation anomalies associated with MJO phase 7.

Description: Atmospheric deposition is an important source of micronutrients to the ocean, but atmospheric deposition fluxes remain poorly constrained in most ocean regions due to the limited number of field observations of wet and dry atmospheric inputs. Here we present the distribution of dissolved aluminium (dAl), as a tracer of atmospheric inputs, in surface waters of the Atlantic Ocean along GEOTRACES sections GA01, GA06, GA08, and GA10. We used the surface mixed layer concentrations of dAl to calculate atmospheric deposition fluxes using a simple steady state model. We have optimized the aerosol Al fractional solubility, dAl residence time within the surface mixed layer and depth of the surface mixed layer for each separate cruise to calculate the atmospheric deposition fluxes. We calculated the lowest deposition fluxes of 0.15 ± 0.1 and 0.27 ± 0.13 g m−2 yr−1 for the South and North Atlantic Ocean (〉 40° S and 〉 40° N), respectively, and highest fluxes of 2.67 ± 1.96 and 3.82 ± 2.72 g m−2 yr−1 for the South East Atlantic and tropical Atlantic Ocean, respectively. Overall, our estimations are comparable to atmospheric dust deposition model estimates and reported field-based atmospheric deposition estimates. We note that our estimates diverge from atmospheric dust deposition model flux estimates in regions influenced by riverine Al inputs and in upwelling regions. As dAl is a key trace element in the GEOTRACES Programme, the approach presented in this study allows calculations of atmospheric deposition fluxes at high spatial resolution for remote ocean regions.

Description: This book includes invited contributions presenting the latest research on the oceanography and environment of the Red Sea. In addition to covering topics relevant to research in the region and providing insights into marine science for non-experts, it is also of interest to those involved in the management of coastal zones and encourages further research on the Red Sea

Description: The nature of the Ionian Sea crust has been the subject of scientific debate for more than 30 years, mainly because seismic imaging of the deep crust and upper mantle of the Ionian Abyssal Plain (IAP) has not been conclusive to date. The IAP is sandwiched between the Calabrian and Hellenic subduction zones in the central Mediterranean. To univocally confirm the proposed oceanic nature of the IAP crust as a remnant of the Tethys ocean and to confute its interpretation as a strongly thinned part of the African continental crust, a NE-SW oriented 131 km long seismic refraction and wide-angle reflection profile consisting of eight ocean bottom seismometers and hydrophones was acquired in 2014. A P-wave velocity model developed from travel time forward modelling is refined by gravimetric data and synthetic modelling of the seismic data. A roughly 6km thick crust with velocities ranging from 5.1km/s to 7.2km/s, top to bottom, can be traced throughout the IAP. In the vicinity of the Medina Seamounts at the southern IAP boundary, the crust thickens to about 9km and seismic velocities decrease to 6.8km/s at the crust-mantle boundary. The seismic velocity distribution and depth of the crust-mantle boundary in the IAP document its oceanic nature, and support the interpretation of the IAP as a remnant of the Tethys oceanic lithosphere formed during the Permian and Triassic period.

Description: The nature of the warm climates of the Cretaceous has been enigmatic since the first numerical climate models were run in the late 1970s. Quantitative simulations of the paleoclimate have consistently failed to agree with information from plant and animal fossils and climate sensitive sediments. The ‘cold continental interior paradox’ (first described by DeConto et al. in Barrera E, Johnson C (eds) Evolution of the Cretaceous Ocean/climate system, vol 332. Geological Society of America Special Paper, Boulder, pp 391–406, 1999), has been an enigma, with extensive continental interiors, especially in northeast Asia, modeled as below freezing in spite of plant and other evidence to the contrary. We reconsider the paleoelevations of specific areas, particularly along the northeastern Siberian continental margin, where paleofloras indeed indicate higher temperatures than suggested by current climate models. Evidence for significant masses of ice on land during even the otherwise warmest times of the Cretaceous is solved by reinterpretation of the δ18O record of fossil plankton. The signal interpreted as an increase in ice volume on land is the same as the signal for an increase in the volume of groundwater reservoirs on land. The problem of a warm Arctic, where fossil floras indicate that they never experienced freezing conditions in winter, could not be solved by numerical simulations using higher CO2 equivalent greenhouse gas concentrations. We propose a solution by assuming that paleoelevations were less than today and that there were much more extensive wetlands (lakes, meandering rivers, swamps, bogs) on the continents than previously assumed. Using ~ 8 × CO2 equivalent greenhouse gas concentrations and assuming 50–75% water surfaces providing water vapor as a supplementary greenhouse gas on the continents reduces the meridional temperature gradients. Under these conditions the equatorial to polar region temperature gradients produce conditions compatible with fossil and sedimentological evidence.

Description: Climate model components involve both high-dimensional input and output fields. It is desirable to efficiently generate spatio-temporal outputs of these models for applications in integrated assessment modelling or to assess the statistical relationship between such sets of inputs and outputs, for example, uncertainty analysis. However, the need for efficiency often compromises the fidelity of output through the use of low complexity models. Here, we develop a technique which combines statistical emulation with a dimensionality reduction technique to emulate a wide range of outputs from an atmospheric general circulation model, PLASIM, as functions of the boundary forcing prescribed by the ocean component of a lower complexity climate model, GENIE-1. Although accurate and detailed spatial information on atmospheric variables such as precipitation and wind speed is well beyond the capability of GENIE-1’s energy-moisture balance model of the atmosphere, this study demonstrates that the output of this model is useful in predicting PLASIM’s spatio-temporal fields through multi-level emulation. Meaningful information from the fast model, GENIE-1 was extracted by utilising the correlation between variables of the same type in the two models and between variables of different types in PLASIM. We present here the construction and validation of several PLASIM variable emulators and discuss their potential use in developing a hybrid model with statistical components.

Description: We describe the development of the “Paleoclimate PLASIM-GENIE emulator” PALEO-PGEM and its application to derive a downscaled high-resolution spatiotemporal description of the climate of the last five million years. The 5-million-year time frame is interesting for a range of paleo-environmental questions, not least because it encompasses the evolution of humans. However, the choice of time-frame was primarily pragmatic; tectonic changes can be neglected to first order, so that it is reasonable to consider climate forcing restricted to the Earth's orbital configuration, ice-sheet state and the concentration of atmosphere CO2. The approach uses the Gaussian process emulation of the singular value decomposition of boundary-condition ensembles of the intermediate complexity atmosphere-ocean GCM PLASIM-GENIE. Spatial fields of bioclimatic variables of surface air temperature (warmest and coolest seasons) and precipitation (wettest and driest seasons) are emulated at 1,000 year intervals, driven by time-series of scalar boundary-condition forcing (CO2, orbit and ice-volume), and assuming the climate is in quasi-equilibrium. Paleoclimate anomalies at climate model resolution are interpolated onto the observed modern climatology to produce a high-resolution spatiotemporal paleoclimate reconstruction of the Pliocene-Pleistocene.