ALBERT

Description: Heinrich events are among the dominant modes of glacial climate variability. During these events, massive iceberg armadas were released by the Laurentide Ice Sheet, sailed across the Atlantic, and caused large-scale climate changes. We study these events in a fully coupled complex ice sheet–climate model with synchronous coupling between ice sheets and oceans. The ice discharges occur as internal variability of the model with a recurrence period of 5kyr, an event duration of 1–1.5kyr, and a peak discharge rate of about 50mSv, roughly consistent with reconstructions. The climate response shows a two-stage behavior, with freshwater release effects dominating the surge phase and ice-sheet elevation effects dominating in the post-surge phase. As a direct response to the freshwater discharge during the surge phase, the deepwater formation in the North Atlantic decreases and the North Atlantic deepwater cell weakens by 3.5Sv. With the reduced oceanic heat transport, the surface temperatures across the North Atlantic decrease, and the associated reduction in evaporation causes a drying in Europe. The ice discharge lowers the surface elevation in the Hudson Bay area and thus leads to increased precipitation and accelerated ice sheet regrowth in the post-surge phase. Furthermore, the jet stream widens to the north and becomes more zonal. This contributes to a weakening of the subpolar gyre, and a continued cooling over Europe even after the ice discharge. This two-stage behavior can explain previously contradicting model results and understandings of Heinrich Events.

Description: Microplastics (MPs, 〈5 mm) have been reported as emerging environmental contaminants, but reliable data are still lacking. We compared the two most promising techniques for MP analysis, namely, Raman and Fourier transform infrared (FTIR) spectroscopy, by analyzing MPs extracted from North Sea surface waters. Microplastics 〉500 μm were visually sorted and manually analyzed by μ-Raman and attenuated total reflection (ATR)-FTIR spectroscopy. Microplastics ≤500 μm were concentrated on gold-coated filters and analyzed by automated single-particle exploration coupled to μ-Raman (ASPEx-μ-Raman) and FTIR imaging (reflection mode). The number of identified MPs 〉500 μm was slightly higher for μ-Raman (+23%) than ATR-FTIR analysis. Concerning MPs ≤500 μm, ASPEx-μ-Raman quantified two-times higher MP numbers but required a four-times higher analysis time compared to FTIR imaging. Because ASPEx-μ-Raman revealed far higher MP concentrations (38–2621 particles m–3) compared to the results of previous water studies (0–559 particles m–3), the environmental concentration of MPs ≤500 μm may have been underestimated until now. This may be attributed to the exceptional increase in concentration with decreasing MP size found in this work. Our results demonstrate the need for further research to enable time-efficient routine application of ASPEx-μ-Raman for reliable MP counting down to 1 μm.

Description: The cage occupancy plays a crucial role in the thermodynamic stability of clathrate hydrates and is an important quantity for understanding the CO2–CH4 replacement phenomenon. In this work, the occupancy isotherms of pure CH4, pure CO2, and their mixture in sI and sII hydrates are studied by GCMC + MD simulations. The adsorption of CH4 and CO2 + CH4 in the sI and sII hydrates can be categorized as the one-site Langmuir type. The calculated occupancy ratio θL/θS and the abundance ratio of CO2 to CH4 vary with the temperature and pressure, which provide the prerequisite information for the prediction of CH4 recovery yield at different conditions in the CO2–CH4 gas exchange process. The phase equilibria of clathrate hydrates of pure gases and mixtures are explored and the corresponding heat of dissociation and hydration numbers are determined. The current investigation provides new perspectives to understand the mechanism behind the gas adsorption behavior of clathrate hydrates.

Description: In this study, we propose a statistical method to validate sea-level reconstructions using geological records known as sea-level indicators (SLIs). SLIs are often the only available data to retrace late-glacial relative sea level (RSL). Determining the RSL from SLI height is not straight forward, the elevation at which an SLI was found usually does not represent the past RSL. In contrast, it has to be related to past RSL by investigating sample’s type, habitat and deposition conditions. For instance, water distribution at which a specific specimen is found today can be related to the indicator's depositional height range. Furthermore, the precision of dating varies between geological samples, and, in case of radiocarbon dating, the age has to be calibrated using a non-linear calibration curve. To avoid an a-priori assumption like normal-distributed uncertainties, we define likelihood functions which take into account the indicative meaning’s available error information and calibration statistics represented by joint probabilities. For this conceptional study, we restrict ourselves to one type of indicators, shallow-water shells, which are usually considered as low-grade samples giving only a lower limit of former sea level, as the depth range in which they live spreads over several tens of meters, and does not follow a normal distribution. The presented method is aimed to serve as a strategy for glacial isostatic adjustment reconstructions, in this case for the German Paleo-Climate Modelling Initiative PalMod (https://www.palmod.de/en) and by extending it to other SLI types.

Description: Dynamic vegetation models simulate global vegetation in terms of fractional coverages 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 reconstructions are initially only available in 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 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 harmonize PFT-distributions by assigning them into nine mega-biomes that follow the definitions commonly used for vegetation reconstructions. The method works well for all state-of the art dynamic vegetation models, independent of the spatial resolution or the complexity of the models. Large biome belts (such as tropical forest) are well represented, but regionally confined biomes (warm-mixed forest, Savanna) are only partly captured. Overall, the PFT-based biomisation is able to keep up with the conventional biomisation approach of forcing biome models (here: BIOME1) with the background climate states. The new method has, however, the advantage that it allows a more direct comparison and evaluation of the vegetation distributions simulated by Earth System Models. Thereby, the new method provides a powerful tool for the evaluation of Earth System Models in general.

Description: We have developed a new module to calculate soil organic carbon (SOC) accumulation in perennially frozen ground in the land surface model JSBACH. Running this offline version of MPI-ESM we have modelled permafrost carbon accumulation and release from the Last Glacial Maximum (LGM) to the Pre-industrial (PI). Our simulated near-surface PI permafrost extent of 16.9Miokm2 is close to observational evidence. Glacial boundary conditions, especially ice sheet coverage, result in profoundly different spatial patterns of glacial permafrost extent. Deglacial warming leads to large-scale changes in soil temperatures, manifested in permafrost disappearance in southerly regions, and permafrost aggregation in formerly glaciated grid cells. In contrast to the large spatial shift in simulated permafrost occurrence, we infer an only moderate increase of total LGM permafrost area (18.3Miokm2) – together with pronounced changes in the depth of seasonal thaw. Reconstructions suggest a larger spread of glacial permafrost towards more southerly regions, but with a highly uncertain extent of non-continuous permafrost. Compared to a control simulation without describing the transport of SOC into perennially frozen ground, the implementation of our newly developed module for simulating permafrost SOC accumulation leads to a doubling of simulated LGM permafrost SOC storage (amounting to a total of ~150PgC). Despite LGM temperatures favouring a larger permafrost extent, simulated cold glacial temperatures – together with low precipitation and low CO2 levels – limit vegetation productivity and therefore prevent a larger glacial SOC build-up in our model. Changes in physical and biogeochemical boundary conditions during deglacial warming lead to an increase in mineral SOC storage towards the Holocene (168PgC at PI), which is below observational estimates (575PgC in continuous and discontinuous permafrost). Additional model experiments clarified the sensitivity of simulated SOC storage to model parameters, affecting long-term soil carbon respiration rates and simulated active layer depths. Rather than a steady increase in carbon release from the LGM to PI as a consequence of deglacial permafrost degradation, our results suggest alternating phases of soil carbon accumulation and loss as an effect of dynamic changes in permafrost extent, active layer depths, soil litter input, and heterotrophic respiration.

Description: Climate reconstructions based on proxy records recovered from marine sediments, such as alkenone records or geochemical parameters measured on foraminifera, play an important role in our understanding of the climate system. They provide information about the state of the ocean ranging back hundreds to millions of years and form the backbone of paleo-oceanography. However, there are many sources of uncertainty associated with the signal recovered from sediment archived proxies. These include seasonal or depth habitat biases in the recorded signal, a frequency dependent reduction in the amplitude of the recorded signal due to bioturbation of the sediment, aliasing of high frequency climate variation onto a nominally annual, decadal or centennial resolution signal, and additional sample processing and measurement error introduced when the proxy signal is recovered. Here we present a forward model for sediment archived proxies that jointly models the above processes, so that the magnitude of their separate and combined effects can be investigated. Applications include the interpretation and analysis of uncertainty in existing proxy records, parameter sensitivity analysis to optimize future studies, and the generation of pseudo-proxy records that can be used to test reconstruction methods. We provide examples, such as the simulation of individual foraminifera records, that demonstrate the usefulness of the forward model for paleoclimate studies. The model is implemented as a user-friendly R package, sedproxy, the use of which we hope will contribute to a better understanding of both the limitations and potential of marine sediment proxies to inform about past climate.

Description: The Argentine margin contains important sedimentological, paleontological and chemical records of regional and local tectonic evolution, sea level, climate evolution and ocean circulation since the opening of the South Atlantic in the Late Jurassic–Early Cretaceous as well as the present-day results of post-depositional chemical and biological alteration. Despite its important location, which underlies the exchange of southern- and northern-sourced water masses, the Argentine margin has not been investigated in detail using scientific drilling techniques, perhaps because the margin has the reputation of being erosional. However, a number of papers published since 2009 have reported new high-resolution and/or multichannel seismic surveys, often combined with multi-beam bathymetric data, which show the common occurrence of layered sediments and prominent sediment drifts on the Argentine and adjacent Uruguayan margins. There has also been significant progress in studying the climatic records in surficial and near-surface sediments recovered in sediment cores from the Argentine margin. Encouraged by these recent results, our 3.5-day IODP (International Ocean Discovery Program) workshop in Buenos Aires (8–11 September 2015) focused on opportunities for scientific drilling on the Atlantic margin of Argentina, which lies beneath a key portion of the global ocean conveyor belt of thermohaline circulation. Significant opportunities exist to study the tectonic evolution, paleoceanography and stratigraphy, sedimentology, and biosphere and geochemistry of this margin.

Description: The discharges from industrial processes constitute the main source of copper contamination in aqueous ecosystems. In this study we investigated the capacity of different types of biochar (derived from chicken manure, eucalyptus, corncob, olive mill and pine sawdust) to remove copper from aqueous solution in a continuous-flow system. The flow rate of the system strongly influenced the amount of copper retained. The adsorption to the corncob biochar varied from 5.51 to 3.48 mg Cu g-1 as the flux decreased from 13 to 2.5 mL min-1. The physicochemical characteristics of biochar determine the copper retention capacity and the underlying immobilization mechanisms. Biochars with high inorganic contents retain the largest amounts of copper and may be suitable for using in water treatment systems to remove heavy metals. The copper retention capacity of the biochars ranged between ~1.3 and 26 mg g-1 and varied in the following order: chicken manure 〉 olive mill 〉〉 corncob 〉 eucalyptus 〉 sawdust pine.

Description: The current generation of marine biogeochemical modules in Earth system models (ESMs) considers mainly the effect of marine biota on the carbon cycle. We propose to also implement other biologically driven mechanisms in ESMs so that more climate-relevant feedbacks are captured. We classify these mechanisms in three categories according to their functional role in the Earth system: (1) "biogeochemical pumps", which affect the carbon cycling; (2) "biological gas and particle shuttles", which affect the atmospheric composition; and (3) "biogeophysical mechanisms", which affect the thermal, optical, and mechanical properties of the ocean. To resolve mechanisms from all three classes, we find it sufficient to include five functional groups: bulk phyto- and zooplankton, calcifiers, and coastal gas and surface mat producers. We strongly suggest to account for a larger mechanism diversity in ESMs in the future to improve the quality of climate projections.

Description: The majority of methane produced in many anoxic sediments is released via ebullition. These bubbles are subject to dissolution as they rise, and dissolution rates are strongly influenced by bubble size. Current understanding of natural methane bubble size distributions is limited by the difficulty in measuring bubble sizes over wide spatial or temporal scales. Our custom optical bubble size sensors recorded bubble sizes and release timing at 8 locations in Upper Mystic Lake, MA continuously for 3 months. Bubble size distributions were spatially heterogeneous even over relatively small areas experiencing similar flux, suggesting that localized sediment conditions are important to controlling bubble size. There was no change in bubble size distributions over the 3 month sampling period, but mean bubble size was positively correlated with daily ebullition flux. Bubble data was used to verify the performance of a widely used bubble dissolution model, and the model was then used to estimate that bubble dissolution accounts for approximately 10% of methane accumulated in the hypolimnion during summer stratification, and at most 15% of the diffusive air–water–methane flux from the epilimnion.

Description: Laboratory sediment incubations and continuous ebullition monitoring over an annual cycle in the temperate Saar River, Germany confirm that impounded river zones can produce and emit methane at high rates (7 to 30 (g CH4 m–3 d–1) at 25 °C and 270 to 700 (g CH4 m–2 yr–1), respectively). Summer methane ebullition (ME) peaks were a factor of 4 to 10 times the winter minima, and sediment methane formation was dominated by the upper sediment (depths of 0.14 to 0.2 m). The key driver of the seasonal ME dynamics was temperature. An empirical model relating methane formation to temperature and sediment depth, derived from the laboratory incubations, reproduced the measured daily ebullition from winter to midsummer, although late summer and autumn simulated ME exceeded the observed ME. A possible explanation for this was substrate limitation. We recommend measurements of methanogenically available carbon sources to identify substrate limitation and help characterize variation in methane formation with depth and from site to site.