ALBERT

Description: This study presents and discusses horizontal and vertical geodetic velocities for a low strain rate region of the south Alpine thrust front in northeastern Italy obtained by integrating GPS, interferometric synthetic aperture radar (InSAR) and leveling data. The area is characterized by the presence of subparallel, south-verging thrusts whose seismogenic potential is still poorly known. Horizontal GPS velocities show that this sector of the eastern Southern Alps is undergoing ∼1 mm a−1 of NW–SE shortening associated with the Adria–Eurasia plate convergence, but the horizontal GPS velocity gradient across the mountain front provides limited constraints on the geometry and slip rate of the several subparallel thrusts. In terms of vertical velocities, the three geodetic methods provide consistent results showing a positive velocity gradient, of ∼ 1.5 mm a−1, across the mountain front, which can hardly be explained solely by isostatic processes. We developed an interseismic dislocation model whose geometry is constrained by available subsurface geological reconstructions and instrumental seismicity. While a fraction of the measured uplift can be attributed to glacial and erosional isostatic processes, our results suggest that interseismic strain accumulation at the Montello and the Bassano–Valdobbiadene thrusts it significantly contributing to the measured uplift. The seismogenic potential of the Montello thrust turns out to be smaller than that of the Bassano–Valdobbiadene fault, whose estimated parameters (locking depth equals 9.1 km and slip rate equals 2.1 mm a−1) indicate a structure capable of potentially generating a Mw〉6.5 earthquake. These results demonstrate the importance of precise vertical ground velocity data for modeling interseismic strain accumulation in slowly deforming regions where seismological and geomorphological evidence of active tectonics is often scarce or not conclusive.

Description: Published

Description: 1681–1698

Description: 2T. Deformazione crostale attiva

Description: JCR Journal

Description: The diagnosis of the conservation state of monumental structures from constraints to the spatial distribution of their physical properties on shallow and inner materials represents one of the key objectives in the application of non-invasive techniques. In situ, CRP and 3D ultrasonic tomography can provide an effective coverage of stone materials in space and time. The intrinsic characteristics of the materials that make up a monumental structure and affect the two properties (i.e., reflectivity, longitudinal velocity) through the above methods substantially differ. Consequently, the content of their information is mainly complementary rather than redundant. In this study we present the integrated application of different non-destructive techniques i.e., Close Range Photogrammetry (CRP), and low frequency (24 KHz) ultrasonic tomography complemented by petrographycal analysis based essentially on Optical Microscopy (OM). This integrated methodology has been applied to a Carrara marble column of the Basilica of San Saturnino, in Byzantine-Proto-Romanesque style, which is part of the Paleo Christian complex of the V-VI century. This complex also includes the adjacent Christian necropolis in the square of San Cosimo in the city of Cagliari, Sardinia, Italy. The column under study is made of bare material dating back probably to the first century A.D., it was subjected to various traumas due to disassembly and transport to the site, including damage caused by the close blast of a WWII fragmentation bomb. High resolution 3D modelling of the studied artifact was computed starting from the integration of proximal sensing techniques such as CRP based on Structure from Motion (SfM), with which information about the geometrical anomalies and reflectivity of the investigated marble column surface was obtained. On the other hand, the inner parts of the studied body were successfully inspected in a non-invasive way by computing the velocity pattern of the ultrasonic signal through the investigated materials using 3D ultrasonic tomography. This technique gives information on the elastic properties of the material related with mechanical properties and a number of factors, such as presence of fractures, voids, and flaws. Extracting information on such factors from the elastic wave velocity using 3D tomography provides a non-invasive approach to analyse the property changes of the inner material of the ancient column. The integrated application of in situ CRP and ultrasonic techniques provides a full 3D high resolution model of the investigated artifact. This model enhanced by the knowledge of the petrographic characteristics of the materials, improves the diagnostic process and affords reliable information on the state of conservation of the materials used in the construction processes of the studied monumental structure. The integrated use of the non-destructive techniques described above also provides suitable data for a possible restoration and future preservation.

Description: Copernicus

Description: Published

Description: On line

Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica

Description: The composition and abundance of algal pigments provide information on characteristics of a phytoplankton community in respect to its photoacclimation, overall biomass, and taxonomic composition. Particularly, these pigments play a major role in photoprotection and in the light-driven part of photosynthesis. Most phytoplankton pigments can be measured by High Performance Liquid Chromatography (HPLC) techniques to filtered water samples. This method, like others when water samples have to be analysed in the laboratory, is time consuming and therefore only a limited number of data points can be obtained. In order to receive information on phytoplankton pigment composition with a higher temporal and spatial resolution, we have developed a method to assess pigment concentrations from continuous optical measurements. The method applies an Empirical Orthogonal Function (EOF) analysis to remote sensing reflectance data derived from ship-based hyper-spectral underwater radiometric and from multispectral satellite data (using the MERIS Polymer product developed by Steinmetz et al., 2011) measured in the Eastern Tropical Atlantic. Subsequently we developed statistically linear models with measured (collocated) pigment concentrations as the response variable and EOF loadings as predictor variables. The model results, show that surface concentrations of a suite of pigments and pigment groups can be well predicted from the ship-based reflectance measurements, even when only a multi-spectral resolution is chosen (i.e. eight bands similar to those used by MERIS). Based on the MERIS reflectance data, concentrations of total and monovinyl chlorophyll a and the groups of photoprotective and photosynthetic carotenoids can be predicted with high quality. The fitted statistical model constructed on the satellite reflectance data as input was applied to one month of MERIS Polymer data to predict the concentration of those pigment groups for the whole Eastern Tropical Atlantic area. Bootstrapping explorations of cross-validation error indicate that the method can produce reliable predictions with relatively small data sets (e.g., 〈 50 collocated values of reflectance and pigment concentration). The method allows for the derivation of time series from continuous reflectance data of various pigment groups at various regions, which can be used to study variability and change of phytoplankton composition and photo-physiology.

Description: Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of time scales. On sub-seasonal time scales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellite (12 m resolution). The high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km2) allow us to track volume losses as drivers such as the available energy change during summer in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm/day) despite the ample available energy, indicating the widespread presence of an insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm/day), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August, respectively. The observed temporal patterns are poorly related to slump characteristics like the slump area. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.

Description: A suite of oxygenated volatile organic compounds (OVOCs – acetaldehyde, acetone, propanal, butanal and butanone) were measured concurrently in the surface water and atmosphere of the South China Sea and Sulu Sea in November 2011. A strong correlation was observed between all OVOC concentrations in the surface seawater along the entire cruise track, except for acetaldehyde, suggesting similar sources and sinks in the surface ocean. Additionally, several phytoplankton groups, such as haptophytes or pelagophytes, were also correlated to all OVOCs indicating that phytoplankton may be an important source for marine OVOCs in the South China and Sulu Seas. Humic and protein like fluorescent dissolved organic matter (FDOM) components seemed to be additional precursors for butanone and acetaldehyde. The atmospheric OVOC mixing ratios were relative high compared with literature values, suggesting the coastal region of North Borneo as a local hot spot for atmospheric OVOCs. The flux of atmospheric OVOCs was largely into the ocean for all 5 gases, with a few important exceptions near the coast of Borneo. The calculated amount of OVOCs entrained into the ocean seemed to be an important source of OVOCs to the surface ocean. When the fluxes were out of the ocean, marine OVOCs were found to be enough to control the local measured OVOC distribution in the atmosphere. Based on our model calculations, at least 0.4 ppb of marine derived acetone and butanone can reach the upper troposphere, where they may have an important influence on hydrogen oxide radical formation over the western Pacific Ocean.

Description: The Extrapolar SWIFT model is a fast ozone chemistry scheme for interactive calculation of the extrapolar stratospheric ozone layer in coupled general circulation models (GCMs). In contrast to the widely used prescribed ozone, the SWIFT ozone layer interacts with the model dynamics and can respond to atmospheric variability or climatological trends. The Extrapolar SWIFT model employs a repro-modelling approach, where algebraic functions are used to approximate the numerical output of a full stratospheric chemistry and transport model (ATLAS). The full model solves a coupled chemical differential equations system with 55 initial and boundary conditions (mixing ratio of various chemical species and atmospheric parameters). Hence the rate of change of ozone over 24  h is a function of 55 variables. Using covariances between these variables, we can find linear combinations in order to reduce the parameter space to the following nine basic variables: latitude, pressure altitude, temperature, local ozone column, mixing ratio of ozone and of the ozone depleting families (Cly, Bry, NOy and HOy). We will show that these 9 variables are sufficient to characterize the rate of change of ozone. An automated procedure fits a polynomial function of fourth degree to the rate of change of ozone obtained from several simulations with the ATLAS model. One polynomial function is determined per month which yields the rate of change of ozone over 24 h. A key aspect for the robustness of the Extrapolar SWIFT model is to include a wide range of stratospheric variability in the numerical output of the ATLAS model, also covering atmospheric states that will occur in a future climate (e.g. temperature and meridional circulation changes or reduction of stratospheric chlorine loading). For validation purposes, the Extrapolar SWIFT model has been integrated into the ATLAS model replacing the full stratospheric chemistry scheme. Simulations with SWIFT in ATLAS have proven that the systematic error is small and does not accumulate during the course of a simulation. In the context of a 10 year simulation, the ozone layer, simulated by SWIFT, shows a stable annual cycle, with inter-annual variations comparable to the ATLAS model. The application of Extrapolar SWIFT requires the evaluation of polynomial functions with 30–100 terms. Nowadays, computers can calculate such polynomial functions at thousands of model grid points in seconds. SWIFT provides the desired numerical efficiency and computes the ozone layer 104 times faster than the chemistry scheme in the ATLAS CTM.

Description: Ice retreat in the eastern Eurasian Arctic is a consequence of atmospheric and oceanic processes and regional feedback mechanisms acting on the ice cover, both in winter and summer. A correct representation of these processes in numerical models is important, since it will improve predictions of sea ice anomalies along the Northeast Passage and beyond. In this study, we highlight the importance of winter ice dynamics for local summer sea ice anomalies in thickness, volume and extent. By means of airborne sea ice thickness surveys made over pack ice areas in the south-eastern Laptev Sea, we show that years of offshore-directed sea ice transport have a thinning effect on the late-winter sea ice cover. To confirm the preconditioning effect of enhanced offshore advection in late winter on the summer sea ice cover, we perform a sensitivity study using a numerical model. Results verify that the preconditioning effect plays a bigger role for the regional ice extent. Furthermore, they indicate an increase in volume export from the Laptev Sea as a consequence of enhanced offshore advection, which has far-reaching consequences for the entire Arctic sea ice mass balance. Moreover we show that ice dynamics in winter not only preconditions local summer ice extent, but also accelerate fast-ice decay.

Description: We use a comprehensive suite of partially laminated high-resolution sediment cores from the Bering Sea, covering a depth transect from 1100 m to 2700 m to study deglacial surface ocean warming patterns, associated changes in biological productivity, oxygen minimum zone dynamics and continent-ocean links through Yukon river runoff. We apply a combination of planktic and benthic isotopes, x-ray fluorescence (XRF)-derived ele- mental ratios and a multi-proxy assessment of changes in upper ocean temperatures. Severe oxygen depletions occurred during the Bølling/Allerød (B/A) and early Holocene, which is in accordance with other locations in the North Pacific, especially the Alaska margin. Detailed analysis of the timing of lamination occurrence between the different sediment cores revealed that the onset of severe anoxia at the beginning of the B/A and early Holocene is a near-synchronous event, while the disappearance of laminations is a diachronic process. The deglacial Oxygen Minimum Zone(OMZ) strengthening is mainly driven by increased export production, visible in XRF-derived elemental ratios, and corresponding high accumulation rates of biogenic components. The export production in turn is a response to rising sea surface temperatures, decreased sea ice cover and increased thermal stratification, while a major nutrient source was the eastern continental shelf, which was flooded during the deglacial global sea level rise. It is discussed controversially whether oxygenation variations in the deglacial subarctic Pacific were coupled to changes in mid-depth water chemistry, or rather a response to physical processes like deep-intermediate ocean or mixed layer warming, or stratification changes. However, knowledge of the driving forcing mechanism for OMZ strengthening is of particular importance, as these are tightly coupled to the regional marine carbon budget, e.g. via the strength and efficiency of the biological pump. Here, our laminated sediments provided the opportunity to study ocean dynamics in exceptional detail, possible on decadal to annual timescales. Due to the correlation patterns of our records to the NGRIP oxygen isotope record through layer counts we presume that (i) the presence of laminations is tightly coupled to submillennial, short-term warm phases, especially during the Bølling-Allerød (B/A), (ii) that the laminations represent annual layered sediments (varves). The latter point in conjunction with our geochemical proxies strongly supports an atmospheric teleconnection between SE Asia, the North Atlantic and the North Pacific, with observed changes in mid-depth ocean dynamics occurring on fast, nearly decadal timescales. Thus, the Bering Sea OMZ is a highly sensitive system reacting almost instantaneously to small temperature changes and therefore has the potential to influence the global carbon budget on short timescales, in particular during episodes of rapidly warming climate.

Description: The aim of the presented study was to investigate the impact on the radiation budget of a biomass-burning plume, transported from Alaska to the High Arctic region of Ny-Ålesund, Svalbard, in early July 2015. Since the mean aerosol optical depth increased by the factor of 10 above the average summer background values, this large aerosol load event is considered particularly exceptional in the last 25 years. In situ data with hygroscopic growth equations, as well as remote sensing measurements as inputs to radiative transfer models, were used, in order to estimate biases associated with (i) hygroscopicity, (ii) variability of single-scattering albedo profiles, and (iii) plane-parallel closure of the modelled atmosphere. A chemical weather model with satellite-derived biomass-burning emissions was applied to interpret the transport and transformation pathways. The provided MODTRAN radiative transfer model (RTM) simulations for the smoke event (14:00 9 July–11:30 11 July) resulted in a mean aerosol direct radiative forcing at the levels of −78.9 and −47.0 W m ^-2 at the surface and at the top of the atmosphere, respectively, for the mean value of aerosol optical depth equal to 0.64 at 550 nm. This corresponded to the average clear-sky direct radiative forcing of −43.3 W/m ^2, estimated by radiometer and model simulations at the surface. Ultimately, uncertainty associated with the plane-parallel atmosphere approximation altered results by about 2 W m^−2. Furthermore, model-derived aerosol direct radiative forcing efficiency reached on average −126 W m^−2/τ550 and −71 W^m−2/τ550 at the surface and at the top of the atmosphere, respectively. The heating rate, estimated at up to 1.8 K day^−1 inside the biomass-burning plume, implied vertical mixing with turbulent kinetic energy of 0.3 m^2s^−2

Description: A new 21.3m firn core was drilled in 2015 at a coastal Antarctic high-accumulation site in Adélie Land (66.78◦ S; 139.56◦ E, 602 m a.s.l.), named Terre Adélie 192A (TA192A). The mean isotopic values (−19.3 ‰ ± 3.1 ‰ for δ18O and 5.4 ‰±2.2 ‰ for deuterium excess) are consistent with other coastal Antarctic values. No significant isotope–temperature relationship can be evidenced at any timescale. This rules out a simple interpretation in terms of local temperature. An observed asymmetry in the δ18O seasonal cycle may be explained by the precipitation of air masses coming from the eastern and western sectors in autumn and winter, recorded in the d-excess signal showing outstanding values in austral spring versus autumn. Significant positive trends are observed in the annual d-excess record and local sea ice extent (135–145◦ E) over the period 1998–2014. However, process studies focusing on resulting isotopic compositions and particularly the deuterium excess–δ18O relationship, evidenced as a potential fingerprint of moisture origins, as well as the collection of more isotopic measurements in Adélie Land are needed for an accurate interpretation of our signals.

Description: A range of future climate scenarios are projected for high atmospheric CO2 concentrations, given uncertainties over future human actions as well as potential environmental and climatic feedbacks. The geological record offers an opportunity to understand climate system response to a range of forcings and feedbacks which operate over multiple temporal and spatial scales. Here, we examine a single interglacial during the late Pliocene (KM5c, ca. 3.205±0.01 Ma) when atmospheric CO2 exceeded pre-industrial concentrations, but were similar to today and to the lowest emission scenarios for this century. As orbital forcing and continental configurations were almost identical to today, we are able to focus on equilibrium climate system response to modern and near-future CO2. Using proxy data from 32 sites, we demonstrate that global mean sea-surface temperatures were warmer than pre-industrial values, by ∼2.3°C for the combined proxy data (foraminifera Mg∕Ca and alkenones), or by ∼3.2–3.4°C (alkenones only). Compared to the pre-industrial period, reduced meridional gradients and enhanced warming in the North Atlantic are consistently reconstructed. There is broad agreement between data and models at the global scale, with regional differences reflecting ocean circulation and/or proxy signals. An uneven distribution of proxy data in time and space does, however, add uncertainty to our anomaly calculations. The reconstructed global mean sea-surface temperature anomaly for KM5c is warmer than all but three of the PlioMIP2 model outputs, and the reconstructed North Atlantic data tend to align with the warmest KM5c model values. Our results demonstrate that even under low-CO2 emission scenarios, surface ocean warming may be expected to exceed model projections and will be accentuated in the higher latitudes.

Description: Reconstructions of global hydroclimate during the Common Era (CE; the past ∼2000 years) are important for providing context for current and future global environmental change. Stable isotope ratios in water are quantitative indicators of hydroclimate on regional to global scales, and these signals are encoded in a wide range of natural geologic archives. Here we present the Iso2k database, a global compilation of previously published datasets from a variety of natural archives that record the stable oxygen (δ18O) or hydrogen (δ2H) isotopic compositions of environmental waters, which reflect hydroclimate changes over the CE. The Iso2k database contains 759 isotope records from the terrestrial and marine realms, including glacier and ground ice (210); speleothems (68); corals, sclerosponges, and mollusks (143); wood (81); lake sediments and other terrestrial sediments (e.g., loess) (158); and marine sediments (99). Individual datasets have temporal resolutions ranging from sub-annual to centennial and include chronological data where available. A fundamental feature of the database is its comprehensive metadata, which will assist both experts and nonexperts in the interpretation of each record and in data synthesis. Key metadata fields have standardized vocabularies to facilitate comparisons across diverse archives and with climate-model-simulated fields. This is the first global-scale collection of water isotope proxy records from multiple types of geological and biological archives. It is suitable for evaluating hydroclimate processes through time and space using large-scale synthesis, model–data intercomparison and (paleo)data assimilation. The Iso2k database is available for download at https://doi.org/10.25921/57j8-vs18 (Konecky and McKay, 2020) and is also accessible via the NOAA/WDS Paleo Data landing page: https://www.ncdc.noaa.gov/paleo/study/29593 (last access: 30 July 2020).

Description: Originating from the boreal forest and often transported over large distances, driftwood characterises many Arctic coastlines. Here we present a combined assessment of radiocarbon (14C) and dendrochronological (ring width) age estimates of driftwood samples to constrain the progradation of two Holocene beach-ridge systems near the Lena Delta in the Siberian Arctic (Laptev Sea). Our data show that the 14C ages obtained on syndepositional driftwood from beach deposits yield surprisingly coherent chronologies for the coastal evolution of the field sites. The dendrochronological analysis of wood from modern driftlines revealed the origin and recent delivery of the wood from the Lena River catchments. This finding suggests that the duration transport lies within the uncertainty of state-of-the-art 14C dating and thus substantiates the validity of age indication obtained from driftwood. This observation will help to better understand changes in similar coastal environments, and to improve our knowledge about the response of coastal systems to past climate and sea-level changes.

Description: The Northeast Greenland Ice Stream (NEGIS) extends around 600 km upstream from the coast to its onset near the ice divide in interior Greenland. Several maps of surface velocity and topography of interior Greenland exist, but their accuracy is not well constrained by in situ observations. Here we present the results from a GPS mapping of surface velocity in an area located approximately 150 km from the ice divide near the East Greenland Ice-core Project (EastGRIP) deep-drilling site. A GPS strain net consisting of 63 poles was established and observed over the years 2015–2019. The strain net covers an area of 35 km by 40 km, including both shear margins. The ice flows with a uniform surface speed of approximately 55 m a^−1 within a central flow band with longitudinal and transverse strain rates on the order of 10−4 a^−1 and increasing by an order of magnitude in the shear margins. We compare the GPS results to the Arctic Digital Elevation Model and a list of satellite-derived surface velocity products in order to evaluate these products. For each velocity product, we determine the bias in and precision of the velocity compared to the GPS observations, as well as the smoothing of the velocity products needed to obtain optimal precision. The best products have a bias and a precision of ∼0.5 m a^−1. We combine the GPS results with satellite-derived products and show that organized patterns in flow and topography emerge in NEGIS when the surface velocity exceeds approximately 55 m a−1 and are related to bedrock topography.

Description: We present the Alfred Wegener Institute's contribution to the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2) wherein we employ the Community Earth System Models (COSMOS) that include a dynamic vegetation scheme. This work builds on our contribution to Phase 1 of the Pliocene Model Intercomparison Project (PlioMIP1) wherein we employed the same model without dynamic vegetation. Our input to the PlioMIP2 special issue of Climate of the Past is twofold. In an accompanying paper we compare results derived with COSMOS in the framework of PlioMIP2 and PlioMIP1. With this paper we present details of our contribution with COSMOS to PlioMIP2. We provide a description of the model and of methods employed to transfer reconstructed mid-Pliocene geography, as provided by the Pliocene Reconstruction and Synoptic Mapping Initiative Phase 4 (PRISM4), to model boundary conditions. We describe the spin-up procedure for creating the COSMOS PlioMIP2 simulation ensemble and present large-scale climate patterns of the COSMOS PlioMIP2 mid-Pliocene core simulation. Furthermore, we quantify the contribution of individual components of PRISM4 boundary conditions to characteristics of simulated mid-Pliocene climate and discuss implications for anthropogenic warming. When exposed to PRISM4 boundary conditions, COSMOS provides insight into a mid-Pliocene climate that is characterised by increased rainfall (+0.17 mm d−1) and elevated surface temperature (+3.37 ∘C) in comparison to the pre-industrial (PI). About two-thirds of the mid-Pliocene core temperature anomaly can be directly attributed to carbon dioxide that is elevated with respect to PI. The contribution of topography and ice sheets to mid-Pliocene warmth is much smaller in contrast – about one-quarter and one-eighth, respectively, and nonlinearities are negligible. The simulated mid-Pliocene climate comprises pronounced polar amplification, a reduced meridional temperature gradient, a northwards-shifted tropical rain belt, an Arctic Ocean that is nearly free of sea ice during boreal summer, and muted seasonality at Northern Hemisphere high latitudes. Simulated mid-Pliocene precipitation patterns are defined by both carbon dioxide and PRISM4 paleogeography. Our COSMOS simulations confirm long-standing characteristics of the mid-Pliocene Earth system, among these increased meridional volume transport in the Atlantic Ocean, an extended and intensified equatorial warm pool, and pronounced poleward expansion of vegetation cover. By means of a comparison of our results to a reconstruction of the sea surface temperature (SST) of the mid-Pliocene we find that COSMOS reproduces reconstructed SST best if exposed to a carbon dioxide concentration of 400 ppmv. In the Atlantic to Arctic Ocean the simulated mid-Pliocene core climate state is too cold in comparison to the SST reconstruction. The discord can be mitigated to some extent by increasing carbon dioxide that causes increased mismatch between the model and reconstruction in other regions.

Description: Earth system and climate modelling involves the simulation of processes on a wide range of scales and within and across various compartments of the Earth system. In practice, component models are often developed independently by different research groups, adapted by others to their special interests and then combined using a dedicated coupling software. This procedure not only leads to a strongly growing number of available versions of model components and coupled setups but also to model- and high-performance computing (HPC)-system-dependent ways of obtaining, configuring, building and operating them. Therefore, implementing these Earth system models (ESMs) can be challenging and extremely time consuming, especially for less experienced modellers or scientists aiming to use different ESMs as in the case of intercomparison projects. To assist researchers and modellers by reducing avoidable complexity, we developed the ESM-Tools software, which provides a standard way for downloading, configuring, compiling, running and monitoring different models on a variety of HPC systems. It should be noted that ESM-Tools is not a coupling software itself but a workflow and infrastructure management tool to provide access to increase usability of already existing components and coupled setups. As coupled ESMs are technically the more challenging tasks, we will focus on coupled setups, always implying that stand-alone models can benefit in the same way. With ESM-Tools, the user is only required to provide a short script consisting of only the experiment-specific definitions, while the software executes all the phases of a simulation in the correct order. The software, which is well documented and easy to install and use, currently supports four ocean models, three atmosphere models, two biogeochemistry models, an ice sheet model, an isostatic adjustment model, a hydrology model and a land-surface model. Compared to previous versions, ESM-Tools has lately been entirely recoded in a high-level programming language (Python) and provides researchers with an even more user-friendly interface for Earth system modelling. ESM-Tools was developed within the framework of the Advanced Earth System Model Capacity project, supported by the Helmholtz Association.

Description: We compare results obtained from modeling the mid-Pliocene warm period using the Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) with the two different modeling methodologies and sets of boundary conditions prescribed for the two phases of the Pliocene Model Intercomparison Project (PlioMIP), tagged PlioMIP1 and PlioMIP2. Here, we bridge the gap between our contributions to PlioMIP1 (Stepanek and Lohmann, 2012) and PlioMIP2 (Stepanek et al., 2020). We highlight some of the effects that differences in the chosen mid-Pliocene model setup (PlioMIP2 vs. PlioMIP1) have on the climate state as derived with COSMOS, as this information will be valuable in the framework of the model–model and model–data comparison within PlioMIP2. We evaluate the model sensitivity to improved mid-Pliocene boundary conditions using PlioMIP's core mid-Pliocene experiments for PlioMIP1 and PlioMIP2 and present further simulations in which we test model sensitivity to variations in paleogeography, orbit, and the concentration of CO2. Firstly, we highlight major changes in boundary conditions from PlioMIP1 to PlioMIP2 and also the challenges recorded from the initial effort. The results derived from our simulations show that COSMOS simulates a mid-Pliocene climate state that is 0.29°C colder in PlioMIP2 if compared to PlioMIP1 (17.82°C in PlioMIP1, 17.53°C in PlioMIP2; values based on simulated surface skin temperature). On the one hand, high-latitude warming, which is supported by proxy evidence of the mid-Pliocene, is underestimated in simulations of both PlioMIP1 and PlioMIP2. On the other hand, spatial variations in surface air temperature (SAT), sea surface temperature (SST), and the distribution of sea ice suggest improvement of simulated SAT and SST in PlioMIP2 if employing the updated paleogeography. Our PlioMIP2 mid-Pliocene simulation produces warmer SSTs in the Arctic and North Atlantic Ocean than those derived from the respective PlioMIP1 climate state. The difference in prescribed CO2 accounts for 0.5°C of temperature difference in the Arctic, leading to an ice-free summer in the PlioMIP1 simulation, and a quasi ice-free summer in PlioMIP2. Beyond the official set of PlioMIP2 simulations, we present further simulations and analyses that sample the phase space of potential alternative orbital forcings that have acted during the Pliocene and may have impacted geological records. Employing orbital forcing, which differs from that proposed for PlioMIP2 (i.e., corresponding to pre-industrial conditions) but falls into the mid-Pliocene time period targeted in PlioMIP, leads to pronounced annual and seasonal temperature variations. Our result identifies the changes in mid-Pliocene paleogeography from PRISM3 to PRISM4 as the major driver of the mid-Pliocene warmth within PlioMIP and not the minor differences in forcings.

Description: Palaeoclimate simulations improve our understanding of the climate, inform us about the performance of climate models in a different climate scenario, and help to identify robust features of the climate system. Here, we analyse Arctic warming in an ensemble of 16 simulations of the mid-Pliocene Warm Period (mPWP), derived from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The PlioMIP2 ensemble simulates Arctic (60–90 °N) annual mean surface air temperature (SAT) increases of 3.7 to 11.6 °C compared to the pre-industrial period, with a multi-model mean (MMM) increase of 7.2 °C. The Arctic warming amplification ratio relative to global SAT anomalies in the ensemble ranges from 1.8 to 3.1 (MMM is 2.3). Sea ice extent anomalies range from −3.0 to −10.4×106 km2, with a MMM anomaly of −5.6×106 km2, which constitutes a decrease of 53 % compared to the pre-industrial period. The majority (11 out of 16) of models simulate summer sea-ice-free conditions (≤1×106 km2) in their mPWP simulation. The ensemble tends to underestimate SAT in the Arctic when compared to available reconstructions, although the degree of underestimation varies strongly between the simulations. The simulations with the highest Arctic SAT anomalies tend to match the proxy dataset in its current form better. The ensemble shows some agreement with reconstructions of sea ice, particularly with regard to seasonal sea ice. Large uncertainties limit the confidence that can be placed in the findings and the compatibility of the different proxy datasets. We show that while reducing uncertainties in the reconstructions could decrease the SAT data–model discord substantially, further improvements are likely to be found in enhanced boundary conditions or model physics. Lastly, we compare the Arctic warming in the mPWP to projections of future Arctic warming and find that the PlioMIP2 ensemble simulates greater Arctic amplification than CMIP5 future climate simulations and an increase instead of a decrease in Atlantic Meridional Overturning Circulation (AMOC) strength compared to pre-industrial period. The results highlight the importance of slow feedbacks in equilibrium climate simulations, and that caution must be taken when using simulations of the mPWP as an analogue for future climate change.

Description: Antarctic geothermal heat flow (GHF) affects the temperature of the ice sheet, determining its ability to slide and internally deform, as well as the behaviour of the continental crust. However, GHF remains poorly constrained, with few and sparse local, borehole-derived estimates and large discrepancies in the magnitude and distribution of existing continent-scale estimates from geophysical models. We review the methods to estimate GHF, discussing the strengths and limitations of each approach; compile borehole and probe-derived estimates from measured temperature profiles; and recommend the following future directions. (1) Obtain more borehole-derived estimates from the subglacial bedrock and englacial temperature profiles. (2) Estimate GHF from inverse glaciological modelling, constrained by evidence for basal melting and englacial temperatures (e.g. using microwave emissivity). (3) Revise geophysically derived GHF estimates using a combination of Curie depth, seismic, and thermal isostasy models. (4) Integrate in these geophysical approaches a more accurate model of the structure and distribution of heat production elements within the crust and considering heterogeneities in the underlying mantle. (5) Continue international interdisciplinary communication and data access.

Description: Northwestern Alaska has been highly affected by changing climatic patterns with new temperature and precipitation maxima over the recent years. In particular, the Baldwin and northern Seward peninsulas are characterized by an abundance of thermokarst lakes that are highly dynamic and prone to lake drainage like many other regions at the southern margins of continuous permafrost. We used Sentinel-1 synthetic aperture radar (SAR) and Planet CubeSat optical remote sensing data to analyze recently observed widespread lake drainage. We then used synoptic weather data, climate model outputs and lake ice growth simulations to analyze potential drivers and future pathways of lake drainage in this region. Following the warmest and wettest winter on record in 2017/2018, 192 lakes were identified as having completely or partially drained by early summer 2018, which exceeded the average drainage rate by a factor of ∼ 10 and doubled the rates of the previous extreme lake drainage years of 2005 and 2006. The combination of abundant rain- and snowfall and extremely warm mean annual air temperatures (MAATs), close to 0 ∘C, may have led to the destabilization of permafrost around the lake margins. Rapid snow melt and high amounts of excess meltwater further promoted rapid lateral breaching at lake shores and consequently sudden drainage of some of the largest lakes of the study region that have likely persisted for millennia. We hypothesize that permafrost destabilization and lake drainage will accelerate and become the dominant drivers of landscape change in this region. Recent MAATs are already within the range of the predictions by the University of Alaska Fairbanks' Scenarios Network for Alaska and Arctic Planning (UAF SNAP) ensemble climate predictions in scenario RCP6.0 for 2100. With MAAT in 2019 just below 0 ∘C at the nearby Kotzebue, Alaska, climate station, permafrost aggradation in drained lake basins will become less likely after drainage, strongly decreasing the potential for freeze-locking carbon sequestered in lake sediments, signifying a prominent regime shift in ice-rich permafrost lowland regions.

Description: The ESA Earth Explorer CryoSat-2 was launched on 8 April 2010 to monitor the precise changes in the thickness of terrestrial ice sheets and marine floating ice. To do that, CryoSat orbits the planet at an altitude of around 720 km with a retrograde orbit inclination of 92∘ and a quasi repeat cycle of 369 d (30 d subcycle). To reach the mission goals, the CryoSat products have to meet the highest quality standards to date, achieved through continual improvements of the operational processing chains. The new CryoSat Ice Baseline-D, in operation since 27 May 2019, represents a major processor upgrade with respect to the previous Ice Baseline-C. Over land ice the new Baseline-D provides better results with respect to the previous baseline when comparing the data to a reference elevation model over the Austfonna ice cap region, improving the ascending and descending crossover statistics from 1.9 to 0.1 m. The improved processing of the star tracker measurements implemented in Baseline-D has led to a reduction in the standard deviation of the point-to-point comparison with the previous star tracker processing method implemented in Baseline-C from 3.8 to 3.7 m. Over sea ice, Baseline-D improves the quality of the retrieved heights inside and at the boundaries of the synthetic aperture radar interferometric (SARIn or SIN) acquisition mask, removing the negative freeboard pattern which is beneficial not only for freeboard retrieval but also for any application that exploits the phase information from SARIn Level 1B (L1B) products. In addition, scatter comparisons with the Beaufort Gyre Exploration Project (BGEP; https://www.whoi.edu/beaufortgyre, last access: October 2019) and Operation IceBridge (OIB; Kurtz et al., 2013) in situ measurements confirm the improvements in the Baseline-D freeboard product quality. Relative to OIB, the Baseline-D freeboard mean bias is reduced by about 8 cm, which roughly corresponds to a 60 % decrease with respect to Baseline-C. The BGEP data indicate a similar tendency with a mean draft bias lowered from 0.85 to −0.14 m. For the two in situ datasets, the root mean square deviation (RMSD) is also well reduced from 14 to 11 cm for OIB and by a factor of 2 for the BGEP. Observations over inland waters show a slight increase in the percentage of good observations in Baseline-D, generally around 5 %–10 % for most lakes. This paper provides an overview of the new Level 1 and Level 2 (L2) CryoSat Ice Baseline-D evolutions and related data quality assessment, based on results obtained from analyzing the 6-month Baseline-D test dataset released to CryoSat expert users prior to the final transfer to operations.

Description: The timing and intensity of snowmelt processes on sea ice are key drivers determining the seasonal sea-ice energy and mass budgets. In the Arctic, satellite passive microwave and radar observations have revealed a trend towards an earlier snowmelt onset during the last decades, which is an important aspect of Arctic amplification and sea ice decline. Around Antarctica, snowmelt on perennial ice is weak and very different than in the Arctic, with most snow surviving the summer. Here we compile time series of snowmelt-onset dates on seasonal and perennial Antarctic sea ice from 1992 to 2014/15 using active microwave observations from European Remote Sensing Satellite (ERS-1/2), Quick Scatterometer (QSCAT) and Advanced Scatterometer (ASCAT) radar scatterometers. We define two snowmelt transition stages: A weak backscatter rise indicating the initial warming and destructive metamorphism of the snowpack (pre-melt), followed by a rapid backscatter rise indicating the onset of thaw-freeze cycles (snowmelt). Results show large interannual variability with an average pre-melt onset date of 29 November and melt onset of 10 December, respectively, on perennial ice, without any significant trends over the study period, consistent with the small trends of Antarctic sea ice extent. There was a latitudinal gradient from early snowmelt onsets in mid-November in the northern Weddell Sea to late (end-December) or even absent snowmelt conditions in the southern Weddell Sea. We show that QSCAT Ku-band (13.4 GHz signal frequency) derived pre-melt and snowmelt onset dates are earlier by 20 and 18 days, respectively, than ERS and ASCAT C-band (5.6 GHz) derived dates. This offset has been considered when constructing the time series. Snowmelt onset dates from passive microwave observations (37 GHz) are later by 14 and 6 days than those from the scatterometers, respectively. Based on these characteristic differences between melt onset dates observed by different microwave wavelengths, we developed a conceptual model which illustrates how the seasonal evolution of snow temperature profiles may affect different microwave bands with different penetration depths. These suggest that future multi-frequency active/passive microwave satellite missions could be used to resolve melt processes throughout the vertical snow column of thick snow on perennial Antarctic sea ice.

Description: The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface to bottom ocean biogeochemical data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of water samples. GLODAPv2.2020 is an update of the previous version, GLODAPv2.2019. The major changes are: data from 106 more cruises added, extension of time coverage until 2019, and the inclusion of available discrete fugacity of CO2 (fCO2) values in the merged product files. GLODAPv2.2020 includes measurements from more than 1.2 million water samples from the global oceans collected on 946 cruises. The data for the 12 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) have undergone extensive quality control, especially systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but updated to WOCE exchange format and (ii) as a merged data product with adjustments applied to minimize bias. These adjustments were derived by comparing the data from the 106 new cruises with the data from the 840 quality-controlled cruises of the GLODAPv2.2019 data product. They correct for errors related to measurement, calibration, and data handling practices, while taking into account any known or likely time trends or variations in the variables evaluated. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 μmol kg−1 in dissolved inorganic carbon, 4 μmol kg−1 in total alkalinity, 0.01–0.02, depending on region, in pH, and 5 % in the halogenated transient tracers. The other variables included in the compilation, such as isotopic tracers and discrete fCO2 were not subjected to bias comparison or adjustments. The original data, their documentation and doi codes are available at the Ocean Carbon Data System of NOAA NCEI (https://www.nodc.noaa.gov/ocads/oceans/GLODAPv2_2020/, last access: 22 June 2020). This site also provides access to the merged data product, which is provided as a single global file and as four regional ones – the Arctic, Atlantic, Indian, and Pacific oceans – under https://doi.org/10.25921/2c8h-sa89 (Olsen et al., 2020). The bias corrected product files also include significant ancillary and approximated data. These were obtained by interpolation of, or calculation from, measured data. This living data update documents the GLODAPv2.2020 methods and provides a broad overview of the secondary quality control procedures and results.

Description: Recent observations of near-surface soil temperatures over the circumpolar Arctic show accelerated warming of permafrost-affected soils. The availability of a comprehensive near-surface permafrost and active layer dataset is critical to better understanding climate impacts and to constraining permafrost thermal conditions and its spatial distribution in land system models. We compiled a soil temperature dataset from 72 monitoring stations in Alaska using data collected by the US Geological Survey, the National Park Service, and the University of Alaska Fairbanks permafrost monitoring networks. The array of monitoring stations spans a large range of latitudes from 60.9 to 71.3 N and elevations from near sea level to~ 1300 m, comprising tundra and boreal forest regions. This dataset consists of monthly ground temperatures at depths up to 1 m, volumetric soil water content, snow depth, and air temperature during 1997–2016. These data have been quality controlled in collection and processing. Meanwhile, we implemented data harmonization evaluation for the processed dataset. The final product (PF-AK, v0. 1) is available at the Arctic Data Center (https://doi. org/10.18739/A2KG55).

Description: Numerous studies have addressed the possible existence of large floating ice sheets in the glacial Arctic Ocean from theoretical, modelling, or seafloor morphology perspectives. Here, we add evidence from the sediment record that support the existence of such freshwater ice caps in certain intervals, and we discuss their implications for possible non-linear and rapid behaviour of such a system in the high latitudes. We present sedimentary activities of 230Th together with 234U/238U ratios, the concentrations of manganese, sulphur and calcium in the context of lithological information and records of microfossils and their isotope composition. New analyses (PS51/038, PS72/396) and a re-analysis of existing marine sediment records (PS1533, PS1235, PS2185, PS2200, amongst others) in view of the naturally occurring radionuclide 230Thex and, where available, 10Be from the Arctic Ocean and the Nordic Seas reveal the widespread occurrence of intervals with a specific geochemical signature. The pattern of these parameters in a pan-Arctic view can best be explained when assuming the repeated presence of freshwater in frozen and liquid form across large parts of the Arctic Ocean and the Nordic Seas. Based on the sedimentary evidence and known environmental constraints at the time, we develop a glacial scenario that explains how these ice sheets, together with eustatic sea-level changes, may have affected the past oceanography of the Arctic Ocean in a fundamental way that must have led to a drastic and non-linear response to external forcing. This concept offers a possibility to explain and to some extent reconcile contrasting age models for the Late Pleistocene in the Arctic Ocean. Our view, if adopted, offers a coherent dating approach across the Arctic Ocean and the Nordic Seas, linked to events outside the Arctic.

Description: Although quantitative isotope data from speleothems has been used to evaluate isotope-enabled model simulations, currently no consensus exists regarding the most appropriate methodology through which to achieve this. A number of modelling groups will be running isotope-enabled palaeoclimate simulations in the framework of the Coupled Model Intercomparison Project Phase 6, so it is timely to evaluate different approaches to using the speleothem data for data–model comparisons. Here, we illustrate this using 456 globally distributed speleothem δ18O records from an updated version of the Speleothem Isotopes Synthesis and Analysis (SISAL) database and palaeoclimate simulations generated using the ECHAM5-wiso isotope-enabled atmospheric circulation model. We show that the SISAL records reproduce the first-order spatial patterns of isotopic variability in the modern day, strongly supporting the application of this dataset for evaluating model-derived isotope variability into the past. However, the discontinuous nature of many speleothem records complicates the process of procuring large numbers of records if data–model comparisons are made using the traditional approach of comparing anomalies between a control period and a given palaeoclimate experiment. To circumvent this issue, we illustrate techniques through which the absolute isotope values during any time period could be used for model evaluation. Specifically, we show that speleothem isotope records allow an assessment of a model’s ability to simulate spatial isotopic trends. Our analyses provide a protocol for using speleothem isotope data for model evaluation, including screening the observations to take into account the impact of speleothem mineralogy on δ18O values, the optimum period for the modern observational baseline and the selection of an appropriate time window for creating means of the isotope data for palaeo-time-slices.

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 H16O, H18O 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: Glacial isostatic adjustment (GIA) is a major source of uncertainty for ice and ocean mass balance estimates derived from satellite gravimetry. In Antarctica the gravimetric effect of cryospheric mass change and GIA are of the same order of magnitude. Inverse estimates from geodetic observations hold some promise for mass signal separation. Here, we investigate the combination of satellite gravimetry and altimetry and demonstrate that the choice of input data sets and processing methods will influence the resultant GIA inverse estimate. This includes the combination that spans the full GRACE record (April 2002–August 2016). Additionally, we show the variations that arise from combining the actual time series of the differing data sets. Using the inferred trends, we assess the spread of GIA solutions owing to (1) the choice of different degree-1 and C20 products, (2) viable candidate surface-elevation-change products derived from different altimetry missions corresponding to different time intervals, and (3) the uncertainties associated with firn process models. Decomposing the total-mass signal into the ice mass and the GIA components is strongly dependent on properly correcting for an apparent bias in regions of small signal. Here our ab initio solutions force the mean GIA and GRACE trend over the low precipitation zone of East Antarctica to be zero. Without applying this bias correction, the overall spread of total-mass change and GIA-related mass change using differing degree-1 and C20 products is 68 and 72 Gt a−1, respectively, for the same time period (March 2003–October 2009). The bias correction method collapses this spread to 6 and 5 Gt a−1, respectively. We characterize the firn process model uncertainty empirically by analysing differences between two alternative surface mass balance products. The differences propagate to a 10 Gt a−1 spread in debiased GIA-related mass change estimates. The choice of the altimetry product poses the largest uncertainty on debiased mass change estimates. The spread of debiased GIA-related mass change amounts to 15 Gt a−1 for the period from March 2003 to October 2009. We found a spread of 49 Gt a−1 comparing results for the periods April 2002–August 2016 and July 2010–August 2016. Our findings point out limitations associated with data quality, data processing, and correction for apparent biases.

Description: Landfast sea ice (fast ice) attached to Antarctic (near-)coastal elements is a critical component of the local physical and ecological systems. Through its direct coupling with the atmosphere and ocean, fast-ice properties are also a potential indicator of processes related to a changing climate. However, in situ fast-ice observations in Antarctica are extremely sparse because of logistical challenges and harsh environmental conditions. Since 2010, a monitoring program observing the seasonal evolution of fast ice in Atka Bay has been conducted as part of the Antarctic Fast Ice Network (AFIN). The bay is located on the northeastern edge of Ekström Ice Shelf in the eastern Weddell Sea, close to the German wintering station Neumayer III. A number of sampling sites have been regularly revisited each year between annual ice formation and breakup to obtain a continuous record of sea-ice and sub-ice platelet-layer thickness, as well as snow depth and freeboard across the bay. Here, we present the time series of these measurements over the last 9 years. Combining them with observations from the nearby Neumayer III meteorological observatory as well as auxiliary satellite images enables us to relate the seasonal and interannual fast-ice cycle to the factors that influence their evolution. On average, the annual consolidated fast-ice thickness at the end of the growth season is about 2 m, with a loose platelet layer of 4 m thickness beneath and 0.70 m thick snow on top. Results highlight the predominately seasonal character of the fast-ice regime in Atka Bay without a significant interannual trend in any of the observed variables over the 9-year observation period. Also, no changes are evident when comparing with sporadic measurements in the 1980s and 1990s. It is shown that strong easterly winds in the area govern the year-round snow distribution and also trigger the breakup of fast ice in the bay during summer months. Due to the substantial snow accumulation on the fast ice, a characteristic feature is frequent negative freeboard, associated flooding of the snow–ice interface, and a likely subsequent snow ice formation. The buoyant platelet layer beneath negates the snow weight to some extent, but snow thermodynamics is identified as the main driver of the energy and mass budgets for the fast-ice cover in Atka Bay. The new knowledge of the seasonal and interannual variability of fast-ice properties from the present study helps to improve our understanding of interactions between atmosphere, fast ice, ocean, and ice shelves in one of the key regions of Antarctica and calls for intensified multidisciplinary studies in this region.

Description: Permafrost temperatures are increasing in Alaska due to climate change and in some cases permafrost is thawing and degrading. In areas where degradation has already occurred the effects can be dramatic, resulting in changing ecosystems, carbon release, and damage to infrastructure. However, in many areas we lack baseline data, such as subsurface temperatures, needed to assess future changes and potential risk areas. Besides climate, the physical properties of the vegetation cover and subsurface material have a major influence on the thermal state of permafrost. These properties are often directly related to the type of ecosystem overlaying permafrost. In this paper we demonstrate that classifying the landscape into general ecotypes is an effective way to scale up permafrost thermal data collected from field monitoring sites. Additionally, we find that within some ecotypes the absence of a moss layer is indicative of the absence of near-surface permafrost. As a proof of concept, we used the ground temperature data collected from the field sites to recode an ecotype land cover map into a map of mean annual ground temperature ranges at 1 m depth based on analysis and clustering of observed thermal regimes. The map should be useful for decision making with respect to land use and understanding how the landscape might change under future climate scenarios.

Description: The Pliocene epoch has great potential to improve our understanding of the long-term climatic and environmental consequences of an atmospheric CO2 concentration near ∼400 parts per million by volume. Here we present the large-scale features of Pliocene climate as simulated by a new ensemble of climate models of varying complexity and spatial resolution based on new reconstructions of boundary conditions (the Pliocene Model Intercomparison Project Phase 2; PlioMIP2). As a global annual average, modelled surface air temperatures increase by between 1.7 and 5.2 °C relative to the pre-industrial era with a multi-model mean value of 3.2 °C. Annual mean total precipitation rates increase by 7 % (range: 2 %–13 %). On average, surface air temperature (SAT) increases by 4.3 °C over land and 2.8 °C over the oceans. There is a clear pattern of polar amplification with warming polewards of 60°N and 60°S exceeding the global mean warming by a factor of 2.3. In the Atlantic and Pacific oceans, meridional temperature gradients are reduced, while tropical zonal gradients remain largely unchanged. There is a statistically significant relationship between a model's climate response associated with a doubling in CO2 (equilibrium climate sensitivity; ECS) and its simulated Pliocene surface temperature response. The mean ensemble Earth system response to a doubling of CO2 (including ice sheet feedbacks) is 67 % greater than ECS; this is larger than the increase of 47 % obtained from the PlioMIP1 ensemble. Proxy-derived estimates of Pliocene sea surface temperatures are used to assess model estimates of ECS and give an ECS range of 2.6–4.8°C. This result is in general accord with the ECS range presented by previous Intergovernmental Panel on Climate Change (IPCC) Assessment Reports.

Description: In the last decades, changing climate conditions have had a severe impact on sea ice at the western Antarctic Peninsula (WAP), an area rapidly transforming under global warming. To study the development of spring sea ice and environmental conditions in the pre-satellite era we investigated three short marine sediment cores for their biomarker inventory with a particular focus on the sea ice proxy IPSO25 and micropaleontological proxies. The core sites are located in the Bransfield Strait in shelf to deep basin areas characterized by a complex oceanographic frontal system, coastal influence and sensitivity to large-scale atmospheric circulation patterns. We analyzed geochemical bulk parameters, biomarkers (highly branched isoprenoids, glycerol dialkyl glycerol tetraethers, sterols), and diatom abundances and diversity over the past 240 years and compared them to observational data, sedimentary and ice core climate archives, and results from numerical models. Based on biomarker results we identified four different environmental units characterized by (A) low sea ice cover and high ocean temperatures, (B) moderate sea ice cover with decreasing ocean temperatures, (C) high but variable sea ice cover during intervals of lower ocean temperatures, and (D) extended sea ice cover coincident with a rapid ocean warming. While IPSO25 concentrations correspond quite well to satellite sea ice observations for the past 40 years, we note discrepancies between the biomarker-based sea ice estimates, the long-term model output for the past 240 years, ice core records, and reconstructed atmospheric circulation patterns such as the El Niño–Southern Oscillation (ENSO) and Southern Annular Mode (SAM). We propose that the sea ice biomarker proxies IPSO25 and PIPSO25 are not linearly related to sea ice cover, and, additionally, each core site reflects specific local environmental conditions. High IPSO25 and PIPSO25 values may not be directly interpreted as referring to high spring sea ice cover because variable sea ice conditions and enhanced nutrient supply may affect the production of both the sea-ice-associated and phytoplankton-derived (open marine, pelagic) biomarker lipids. For future interpretations we recommend carefully considering individual biomarker records to distinguish between cold sea-ice-favoring and warm sea-ice-diminishing environmental conditions.

Description: The modeling of paleoclimate, using physically based tools, is increasingly seen as a strong out-of-sample test of the models that are used for the projection of future climate changes. New to the Coupled Model Intercomparison Project (CMIP6) is the Tier 1 Last Interglacial experiment for 127 000 years ago (lig127k), designed to address the climate responses to stronger orbital forcing than the midHolocene experiment, using the same state-of-the-art models as for the future and following a common experimental protocol. Here we present a first analysis of a multi-model ensemble of 17 climate models, all of which have completed the CMIP6 DECK (Diagnostic, Evaluation and Characterization of Klima) experiments. The equilibrium climate sensitivity (ECS) of these models varies from 1.8 to 5.6 ∘C. The seasonal character of the insolation anomalies results in strong summer warming over the Northern Hemisphere continents in the lig127k ensemble as compared to the CMIP6 piControl and much-reduced minimum sea ice in the Arctic. The multi-model results indicate enhanced summer monsoonal precipitation in the Northern Hemisphere and reductions in the Southern Hemisphere. These responses are greater in the lig127k than the CMIP6 midHolocene simulations as expected from the larger insolation anomalies at 127 than 6 ka. New synthesis for surface temperature and precipitation, targeted for 127 ka, have been developed for comparison to the multi-model ensemble. The lig127k model ensemble and data reconstructions are in good agreement for summer temperature anomalies over Canada, Scandinavia, and the North Atlantic and for precipitation over the Northern Hemisphere continents. The model–data comparisons and mismatches point to further study of the sensitivity of the simulations to uncertainties in the boundary conditions and of the uncertainties and sparse coverage in current proxy reconstructions. The CMIP6–Paleoclimate Modeling Intercomparison Project (PMIP4) lig127k simulations, in combination with the proxy record, improve our confidence in future projections of monsoons, surface temperature, and Arctic sea ice, thus providing a key target for model evaluation and optimization.

Description: Climate trends in the Antarctic region remain poorly characterized, owing to the brevity and scarcity of direct climate observations and the large magnitude of interannual to decadal-scale climate variability. Here, within the framework of the PAGES Antarctica2k working group, we build an enlarged database of ice core water stable isotope records from Antarctica, consisting of 112 records. We produce both unweighted and weighted isotopic (δ18O) composites and temperature reconstructions since 0 CE, binned at 5- and 10-year resolution, for seven climatically distinct regions covering the Antarctic continent. Following earlier work of the Antarctica2k working group, we also produce composites and reconstructions for the broader regions of East Antarctica, West Antarctica and the whole continent. We use three methods for our temperature reconstructions: (i) a temperature scaling based on the δ18O–temperature relationship output from an ECHAM5-wiso model simulation nudged to ERA-Interim atmospheric reanalyses from 1979 to 2013, and adjusted for the West Antarctic Ice Sheet region to borehole temperature data, (ii) a temperature scaling of the isotopic normalized anomalies to the variance of the regional reanalysis temperature and (iii) a composite-plus-scaling approach used in a previous continent-scale reconstruction of Antarctic temperature since 1 CE but applied to the new Antarctic ice core database. Our new reconstructions confirm a significant cooling trend from 0 to 1900 CE across all Antarctic regions where records extend back into the 1st millennium, with the exception of the Wilkes Land coast and Weddell Sea coast regions. Within this long-term cooling trend from 0 to 1900 CE, we find that the warmest period occurs between 300 and 1000 CE, and the coldest interval occurs from 1200 to 1900 CE. Since 1900 CE, significant warming trends are identified for the West Antarctic Ice Sheet, the Dronning Maud Land coast and the Antarctic Peninsula regions, and these trends are robust across the distribution of records that contribute to the unweighted isotopic composites and also significant in the weighted temperature reconstructions. Only for the Antarctic Peninsula is this most recent century-scale trend unusual in the context of natural variability over the last 2000 years. However, projected warming of the Antarctic continent during the 21st century may soon see significant and unusual warming develop across other parts of the Antarctic continent. The extended Antarctica2k ice core isotope database developed by this working group opens up many avenues for developing a deeper understanding of the response of Antarctic climate to natural and anthropogenic climate forcings. The first long-term quantification of regional climate in Antarctica presented herein is a basis for data–model comparison and assessments of past, present and future driving factors of Antarctic climate.

Description: The denudation history of active orogens is often interpreted in the context of modern climate gradients. Here we address the validity of this approach and ask what are the spatial and temporal variations in palaeoclimate for a latitudinally diverse range of active orogens? We do this using high-resolution (T159, ca. 80 × 80 km at the Equator) palaeoclimate simulations from the ECHAM5 global atmospheric general circulation model and a statistical cluster analysis of climate over different orogens (Andes, Himalayas, SE Alaska, Pacific NW USA). Time periods and boundary conditions considered include the Pliocene (PLIO, ∼3Ma), the Last Glacial Maximum (LGM, ∼21ka), mid-Holocene (MH, ∼6ka), and pre-industrial (PI, reference year 1850). The regional simulated climates of each orogen are described by means of cluster analyses based on the variability in precipitation, 2 m air temperature, the intra-annual amplitude of these values, and monsoonal wind speeds where appropriate. Results indicate the largest differences in the PI climate existed for the LGM and PLIO climates in the form of widespread cooling and reduced precipitation in the LGM and warming and enhanced precipitation during the PLIO. The LGM climate shows the largest deviation in annual precipitation from the PI climate and shows enhanced precipitation in the temperate Andes and coastal regions for both SE Alaska and the US Pacific Northwest. Furthermore, LGM precipitation is reduced in the western Himalayas and enhanced in the eastern Himalayas, resulting in a shift of the wettest regional climates eastward along the orogen. The cluster-analysis results also suggest more climatic variability across latitudes east of the Andes in the PLIO climate than in other time slice experiments conducted here. Taken together, these results highlight significant changes in late Cenozoic regional climatology over the last ∼3Myr. Comparison of simulated climate with proxy-based reconstructions for the MH and LGM reveal satisfactory to good performance of the model in reproducing precipitation changes, although in some cases discrepancies between neighbouring proxy observations highlight contradictions between proxy observations themselves. Finally, we document regions where the largest magnitudes of late Cenozoic changes in precipitation and temperature occur and offer the highest potential for future observational studies that quantify the impact of climate change on denudation and weathering rates.

Description: In austral spring 2016 the Antarctic region experienced anomalous sea ice retreat in all sectors, with sea ice extent in October and November 2016 being the lowest in the Southern Hemisphere over the observational period (1979–present). The extreme sea ice retreat was accompanied by widespread warming along the coastal areas as well as in the interior of the Antarctic continent. This exceptional event occurred along with a strong negative phase of the Southern Annular Mode (SAM) and the moistest and warmest spring on record, over large areas covering the Indian Ocean, the Ross Sea and the Weddell Sea. In October 2016, the positive anomalies of the totally integrated water vapor (IWV) and 2 m air temperature (T2m) over the Indian Ocean, western Pacific, Bellingshausen Sea and southern part of Ross Sea were unprecedented in the last 39 years. In October and November 2016, when the largest magnitude of negative daily sea ice concentration anomalies was observed, repeated episodes of poleward advection of warm and moist air took place. These results suggest the importance of moist and warm air intrusions into the Antarctic region as one of the main contributors to this exceptional sea ice retreat event.

Description: To resolve the mechanisms behind the major climate reorganisation which occurred between 0.9 and 1.2Ma, the recovery of a suitable 1.5 million-year-old ice core is fundamental. The quest for such an Oldest Ice core requires a number of key boundary conditions, of which the poorly known basal geothermal heat flux (GHF) is lacking. We use a transient thermodynamical 1D vertical model that solves for the rate of change of temperature in the vertical, with surface temperature and modelled GHF as boundary conditions. For each point on the ice sheet, the model is forced with variations in atmospheric conditions over the last 2Ma, and modelled ice-thickness variations. The process is repeated for a range of GHF values to determine the value of GHF that marks the limit between frozen and melting conditions over the whole ice sheet, taking into account 2Ma of climate history. These threshold values of GHF are statistically compared to existing GHF data sets. The new probabilistic GHF fields obtained for the ice sheet thus provide the missing boundary conditions in the search for Oldest Ice. High spatial resolution radar data are examined locally in the Dome Fuji and Dome C regions, as these represent the ice core community's primary drilling sites. GHF, bedrock variability, ice thickness and other essential criteria combined highlight a dozen major potential Oldest Ice sites in the vicinity of Dome Fuji and Dome C, where GHF allows for Oldest Ice.

Description: Polar ice core water isotope records are commonly used to infer past changes in Antarctic temperature, motivating an improved understanding and quantification of the temporal relationship between δ18O and temperature. This can be achieved using simulations performed by atmospheric general circulation models equipped with water stable isotopes. Here, we evaluate the skills of the high-resolution water-isotope-enabled atmospheric general circulation model ECHAM5-wiso (the European Centre Hamburg Model) nudged to European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis using simulations covering the period 1960–2013 over the Antarctic continent. We compare model outputs with field data, first with a focus on regional climate variables and second on water stable isotopes, using our updated dataset of water stable isotope measurements from precipitation, snow, and firn–ice core samples. ECHAM5-wiso simulates a large increase in temperature from 1978 to 1979, possibly caused by a discontinuity in the European Reanalyses (ERA) linked to the assimilation of remote sensing data starting in 1979. Although some model–data mismatches are observed, the (precipitation minus evaporation) outputs are found to be realistic products for surface mass balance. A warm model bias over central East Antarctica and a cold model bias over coastal regions explain first-order δ18O model biases by too strong isotopic depletion on coastal areas and underestimated depletion inland. At the second order, despite these biases, ECHAM5-wiso correctly captures the observed spatial patterns of deuterium excess. The results of model–data comparisons for the inter-annual δ18O standard deviation difer when using precipitation or ice core data. Further studies should explore the importance of deposition and post-deposition processes affecting ice core signals and not resolved in the model. These results build trust in the use of ECHAM5-wiso outputs to investigate the spatial, seasonal, and inter-annual δ18O–temperature relationships. We thus make the first Antarctica-wide synthesis of prior results. First, we show that local spatial or seasonal slopes are not a correct surrogate for inter-annual temporal slopes, leading to the conclusion that the same isotope–temperature slope cannot be applied for the climatic interpretation of Antarctic ice core for all timescales. Finally, we explore the phasing between the seasonal cycles of deuterium excess and δ18O as a source of information on changes in moisture sources affecting the δ18O–temperature relationship. The few available records and ECHAM5-wiso show different phase relationships in coastal, intermediate, and central regions. This work evaluates the use of the ECHAM5-wiso model as a tool for the investigation of water stable isotopes in Antarctic precipitation and calls for extended studies to improve our understanding of such proxies.

Description: The effect of external forcings on atmospheric circulation is debated. Due to the short observational period, the analysis of the role of external forcings is hampered, making it difficult to assess the sensitivity of atmospheric circulation to external forcings, as well as persistence of the effects. In observations, the average response to tropical volcanic eruptions is a positive North Atlantic Oscillation (NAO) during the following winter. However, past major tropical eruptions exceeding the magnitude of eruptions during the instrumental era could have had more lasting effects. Decadal NAO variability has been suggested to follow the 11-year solar cycle, and linkages have been made between grand solar minima and negative NAO. However, the solar link to NAO found by modeling studies is not unequivocally supported by reconstructions, and is not consistently present in observations for the 20th century. Here we present a reconstruction of atmospheric winter circulation for the North Atlantic region covering the period 1241–1970 CE. Based on seasonally resolved Greenland ice core records and a 1200-year-long simulation with an isotope-enabled climate model, we reconstruct sea level pressure and temperature by matching the spatiotemporal variability in the modeled isotopic composition to that of the ice cores. This method allows us to capture the primary (NAO) and secondary mode (Eastern Atlantic Pattern) of atmospheric circulation in the North Atlantic region, while, contrary to previous reconstructions, preserving the amplitude of observed year-to-year atmospheric variability. Our results show five winters of positive NAO on average following major tropical volcanic eruptions, which is more persistent than previously suggested. In response to decadal minima of solar activity we find a high-pressure anomaly over northern Europe, while a reinforced opposite response in pressure emerges with a 5-year time lag. On centennial timescales we observe a similar response of circulation as for the 5-year time-lagged response, with a high-pressure anomaly across North America and south of Greenland. This response to solar forcing is correlated to the second mode of atmospheric circulation, the Eastern Atlantic Pattern. The response could be due to an increase in blocking frequency, possibly linked to a weakening of the subpolar gyre. The long-term anomalies of temperature during solar minima shows cooling across Greenland, Iceland and western Europe, resembling the cooling pattern during the Little Ice Age (1450–1850 CE). While our results show significant correlation between solar forcing and the secondary circulation pattern on decadal (r = 0.29, p 〈 0.01) and centennial timescales (r = 0.6, p 〈 0.01), we find no consistent relationship between solar forcing and NAO. We conclude that solar and volcanic forcing impacts different modes of our reconstructed atmospheric circulation, which can aid in separating the regional effects of forcings and understanding the underlying mechanisms.

Description: North Pacific Intermediate water (NPIW) is a dominant water mass controlling ∼400-1200m depth North Pacific Ocean, meanwhile there is a cessation of North Pacific deep water (NPDW) formation in in modern observations. In contrast, paleoceanographic evidences have recorded NPDW formations during last glacial periods. This sug- gests either a rapid or gradual shutting down process of NPDW formation during the last deglaciation. Here, we use an Earth System Model to diagnose the physical and corresponding biogeochemical evolutions in the North Pacific Ocean before and after the last deglaciation, as well as potential changes during rapid climate shifts of the last deglaciation. Linked to different background climate conditions and varying Atlantic Meridional Over- turning Circulation states, we characterize the modelled NPIW and NPDW changes and builds up linkages to marine records. Our results further develop our understanding about the deglacial switch from NPDW to modern NPIW-only formation process.

Description: Permafrost presence is determined by a complex interaction of climatic, topographic, and ecological conditions operating over long time scales. In particular, vegetation and organic layer characteristics may act to protect permafrost in regions with a mean annual air temperature (MAAT) above 0°C. In this study, we document the presence of residual permafrost plateaus in the western Kenai Peninsula lowlands of south-central Alaska, a region with a MAAT of 1.5+/-1 °C (1981–2010). Continuous ground temperature measurements between 16 September 2012 and 15 September 2015, using calibrated thermistor strings, documented the presence of warm permafrost (-0.04 to -0.08 °C). Field measurements (probing) on several plateau features during the fall of 2015 showed that the depth to the permafrost table averaged 1.48m but at some locations was as shallow as 0.53 m. Late winter surveys (augering, coring, and GPR) in 2016 showed that the average seasonally frozen ground thickness was 0.45 m, overlying a talik above the permafrost table. Measured permafrost thickness ranged from 0.33 to 〉6.90 m. Manual interpretation of historic aerial photography acquired in 1950 indicates that residual permafrost plateaus covered 920 ha as mapped across portions of four wetland complexes encompassing 4810 ha. However, between 1950 and ca. 2010, permafrost plateau extent decreased by 60.0 %, with lateral feature degradation accounting for 85.0% of the reduction in area. Permafrost loss on the Kenai Peninsula is likely associated with a warming climate, wildfires that remove the protective forest and organic layer cover, groundwater flow at depth, and lateral heat transfer from wetland surface waters in the summer. Better understanding the resilience and vulnerability of ecosystem-protected permafrost is critical for mapping and predicting future permafrost extent and degradation across all permafrost regions that are currently warming. Further work should focus on reconstructing permafrost history in south-central Alaska as well as additional contemporary observations of these ecosystem-protected permafrost sites south of the regions with relatively stable permafrost.

Description: Methane plays an important role in the Earth’s atmospheric chemistry and radiative balance being the second most important greenhouse gas after carbon dioxide. Methane is released to the atmosphere by a wide number of sources, both natural and anthropogenic, with the latter being twice as large as the former (IPCC, 2007). It has recently been established that significant amounts of geological methane, produced within the Earth’s crust, are currently released naturally into the atmosphere (Etiope, 2004). Active or recent volcanic/geothermal areas represent one of these sources of geological methane. But due to the fact that methane flux measurements are laboratory intensive, very few data have been collected until now and the contribution of this source has been generally indirectly estimated (Etiope et al., 2007). The Greek territory is geodynamically very active and has many volcanic and geothermal areas. Here we report on methane flux measurements made at two volcanic/geothermal systems along the South Aegean volcanic arc: Sousaki and Nisyros. The former is an extinct volcanic area of Plio-Pleistocene age hosting nowadays a low enthalpy geothermal field. The latter is a currently quiescent active volcanic system with strong fumarolic activity due to the presence of a high enthalpy geothermal system. Both systems have gas manifestations that emit significant amounts of hydrothermal methane and display important diffuse carbon dioxide emissions from the soils. New data on methane isotopic composition and higher hydrocarbon contents point to an abiogenic origin of the hydrothermal methane in the studied systems. Measured methane flux values range from –48 to 29,000 (38 sites) and from –20 to 1100 mg/mˆ2/d (35 sites) at Sousaki and Nisyros respectively. At Sousaki measurement sites covered almost all the degassing area and the diffuse methane output can be estimated in about 20 t/a from a surface of about 10,000 mˆ2. At Nisyros measurements covered the Stephanos and Kaminakia areas, which represent only a part of the entire degassing area. The two areas show very different methane degassing pattern with latter showing much higher flux values. Methane output can be estimated in about 0.25 t/a from an area of about 30,000 mˆ2 at Stephanos and about 1 t/a from an area of about 20,000 mˆ2 at Kaminakia. The total output from the entire geothermal system of Nisyros probably should not exceed 2 t/a.

Description: Published

Description: Vienna, Austria

Description: 4.5. Studi sul degassamento naturale e sui gas petroliferi

Description: open

Description: A biomonitoring survey, above tree line level, using two endemic species (Senecio aethnensis and Rumex aethnensis) was performed on Mt. Etna, in order to evaluate the dispersion and the impact of volcanic atmospheric emissions. Samples of leaves were collected in summer 2008 from 30 sites in the upper part of the volcano (1500- 3000 m a.s.l). Acid digestion of samples was carried out with a microwave oven, and 44 elements were analyzed by using plasma spectrometry (ICP-MS and ICP-OES). The highest concentrations of all investigated elements were found in the samples collected closest to the degassing craters, and in the downwind sector, confirming that the eastern flank of Mt. Etna is the most impacted by volcanic emissions. Leaves collected along two radial transects from the active vents on the eastern flank, highlight that the levels of metals decrease one or two orders of magnitude with increasing distance from the source. This variability is higher for volatile elements (As, Bi, Cd, Cs, Pb, Sb, Tl) than for more refractory elements (Al, Ba, Sc, Si, Sr, Th, U). The two different species of plants do not show significant differences in the bioaccumulation of most of the analyzed elements, except for lanthanides, which are systematically enriched in Rumex leaves. The high concentrations of many toxic elements in the leaves allow us to consider these plants as highly tolerant species to the volcanic emissions, and suitable for biomonitoring researches in the Mt. Etna area.

Description: Published

Description: Vienna, Austria

Description: 4.4. Scenari e mitigazione del rischio ambientale

Description: open

Description: Etna volcano, Italy, hosts one of the major groundwater systems of the island of Sicily. Waters circulate within highly permeable fractured, mainly hawaiitic, volcanic rocks. Aquifers are limited downwards by the underlying impermeable sedimentary terrains. Thickness of the volcanic rocks generally does not exceed some 300 m, preventing the waters to reach great depths. This is faced by short travel times (years to tens of years) and low thermalisation of the Etnean groundwaters. Measured temperatures are, in fact, generally lower than 25 °C. But the huge annual meteoric recharge (about 0.97 kmˆ3) with a high actual infiltration coefficient (0.75) implies a great underground circulation. During their travel from the summit area to the periphery of the volcano, waters acquire magmatic heat together with volcanic gases and solutes through water-rock interaction processes. In the last 20 years the Etnean aquifers has been extensively studied. Their waters were analysed for dissolved major, minor and trace element, O, H, C, S, B, Sr and He isotopes, and dissolved gas composition. These data have been published in several articles. Here, after a summary of the obtained results, the estimation of the magmatic heat flux through the aquifer will be discussed. To calculate heat uptake during subsurface circulation, for each sampling point (spring, well or drainage gallery) the following data have been considered: flow rate, water temperature, and oxygen isotopic composition. The latter was used to calculate the mean recharge altitude through the measured local isotopic lapse rate. Mean recharge temperatures, weighted for rain amount throughout the year, were obtained from the local weather station network. Calculations were made for a representative number of sampling points (216) including all major issues and corresponding to a total water flow of about 0.315 kmˆ3/a, which is 40% of the effective meteoric recharge. Results gave a total energy output of about 140 MW/a the half of which is ascribable to only 13 sampling points. These correspond to the highest flow drainage galleries with fluxes ranging from 50 to 1000 l/s and wells with pumping rates from 70 to 250 l/s. Geographical distribution indicates that, like magmatic gas leakage, heat flow is influenced by structural features of the volcanic edifice. The major heat discharge through groundwater are all tightly connected either to the major regional tectonic systems or to the major volcanic rift zones along which the most important flank eruptions take place. But rift zones are much more important for heat upraise due to the frequent dikes injection than for gas escape because generally when dikes have been emplaced the structure is no more permeable to gases because it becomes sealed by the cooling magma.

Description: Published

Description: Vienna, Austria

Description: 1.2. TTC - Sorveglianza geochimica delle aree vulcaniche attive

Description: open

Description: Volcanoes represent an important natural source of several trace elements to the atmosphere. For some species (e.g., As, Cd, Pb and Se) they may be the main natural source and thereby strongly influencing geochemical cycles from the local to the global scale. Mount Etna is one of the most actively degassing volcanoes in the world, and it is considered to be, on the long-term average, the major atmospheric point source of many environmental harmful compounds. Their emission occurs either through continuous passive degassing from open-conduit activity or through sporadic paroxysmal eruptive activity, in the form of gases, aerosols or particulate. To estimate the environmental impact of magma-derived trace metals and their depositions processes, rainwater and snow samples were collected at Mount Etna area. Five bulk collectors have been deployed at various altitudes on the upper flanks around the summit craters of the volcano; samples were collected every two week for a period of one year and analyzed for the main chemical-physical parameters (electric conductivity and pH) and for major and trace elements concentrations. Chemical analysis of rainwater clearly shows that the volcanic contribution is always prevailing in the sampling site closest to the summit crater (about 1.5 km). In the distal sites (5.5-10 km from the summit) and downwind of the summit craters, the volcanic contribution is also detectable but often overwhelmed by anthropogenic or other natural (seawater spray, geogenic dust) contributions. Volcanic contribution may derive from both dry and wet deposition of gases and aerosols from the volcanic plume, but sometimes also from leaching of freshly emitted volcanic ashes. In fact, in our background site (7.5 km in the upwind direction) volcanic contribution has been detected only following an ash deposition event. About 30 samples of fresh snow were collected in the upper part of the volcano, during the winters 2006 and 2007 to estimate deposition processes at high altitude during cold periods. Some of the samples were collected immediately after a major explosive event from the summit craters to understand the interaction between snow and fresh erupted ash. Sulphur, Chlorine and Fluorine, are the major elements that prevailingly characterize the volcanic contribution in atmospheric precipitation on Mount Etna, but high concentrations of many trace elements are also detected in the studied samples. In particular, bulk deposition samples display high concentration of Al, Fe, Ti, Cu, As, Rb, Pb, Tl, Cd, Cr, U and Ag, in the site most exposed to the volcanic emissions: median concentration values are about two orders of magnitude higher than those measured in our background site. Also in the snow samples the volcanic signature is clearly detectable and decreases with distance from the summit craters. Some of the analysed elements display very high enrichment values with respect to the average crust and, in the closest site to the summit craters, also deposition values higher than those measured in polluted urban or industrial sites.

Description: Published

Description: Vienna, Austria

Description: 4.5. Degassamento naturale

Description: open

Description: Improving the constraints on the atmospheric fate and depletion rates of acidic compounds persistently emitted by non-erupting (quiescent) volcanoes is important for quantitatively predicting the environmental impact of volcanic gas plumes. Here, we present new experimental data coupled with modelling studies to investigate the chemical processing of acidic volcanogenic species during tropospheric dispersion. Diffusive tube samplers were deployed at Mount Etna, a very active open-conduit basaltic volcano in eastern Sicily, and Vulcano Island, a closed-conduit quiescent volcano in the Aeolian Islands (northern Sicily). Sulphur dioxide (SO2), hydrogen sulphide (H2S), hydrogen chloride (HCl) and hydrogen fluoride (HF) concentrations in the volcanic plumes (typically several minutes to a few hours old) were repeatedly determined at distances from the summit vents ranging from 0.1 to ~10 km, and under different environmental conditions. At both volcanoes, acidic gas concentrations were found to decrease exponentially with distance from the summit vents (e.g., SO2 decreases from ~10,000 μg/m3 at 0.1 km from Etna’s vents down to ~7 _μg/m3 at ~10km distance), reflecting the atmospheric dilution of the plume within the acid gas-free background troposphere. Conversely, SO2/HCl, SO2/HF, and SO2/H2S ratios in the plume showed no systematic changes with plume aging, and fit source compositions within analytical error. Assuming that SO2 losses by reaction are small during short-range atmospheric transport within quiescent (ash-free) volcanic plumes, our observations suggest that, for these short transport distances, atmospheric reactions for H2S and halogens are also negligible. The one-dimensional model MISTRA was used to simulate quantitatively the evolution of halogen and sulphur compounds in the plume of Mt. Etna. Model predictions support the hypothesis of minor HCl chemical processing during plume transport, at least in cloud-free conditions. Larger variations in the modelled SO2/HCl ratios were predicted under cloudy conditions, due to heterogeneous chlorine cycling in the aerosol phase. The modelled evolution of the SO2/H2S ratios is found to be substantially dependent on whether or not the interactions of H2S with halogens are included in the model. In the former case, H2S is assumed to be oxidized in the atmosphere mainly by OH, which results in minor chemical loss for H2S during plume aging and produces a fair match between modelled and measured SO2/H2S ratios. In the latter case, fast oxidation of H2S by Cl leads to H2S chemical lifetimes in the early plume of a few seconds, and thus SO2 to H2S ratios that increase sharply during plume transport. This disagreement between modelled and observed plume compositions suggests that more in-detail kinetic investigations are required for a proper evaluation of H2S chemical processing in volcanic plumes.

Description: Published

Description: 1441-1450

Description: 1.2. TTC - Sorveglianza geochimica delle aree vulcaniche attive

Description: 4.5. Degassamento naturale

Description: JCR Journal

Description: open

Description: After some short test surveys, during the 2004–2005 summer expedition in Antarctica, a geomagnetic French-Italian observatory was installed on the plateau (geographic coordinates: 75.1 S, 123.4 E; corrected geomagnetic coordinates: 88.9 S, 54.3 E; UT=LT−8) very close to the geomagnetic pole. In this paper we present some peculiarities of the daily variation as observed at this polar cap observatory during the years 2005 and 2006, taking into account the different Loyd seasons and different interplanetary magnetic field conditions. Some interesting results emerge from the analysis, confirming the dependence of the daily variation (and of the associated polar current systems) on the IMF Bz and By components. In particular the analysis showed that different Bz conditions correspond to different contribution to daily variation of ionospheric and field aligned currents, while particular By conditions lead to a time shift of the diurnal variation, indicating an asymmetry with respect to the noon meridian.

Description: Published

Description: 2045–2051

Description: 3.9. Fisica della magnetosfera, ionosfera e meteorologia spaziale

Description: 1.6. Osservazioni di geomagnetismo

Description: JCR Journal

Description: reserved

Description: During the 2007-2008 antarctic campaign, the Italian PNRA installed a Low Power Magnetometer within the framework of the AIMNet (Antarctic International Magnetometer Network) project, proposed and coordinated by BAS. The magnetometer is situated at Talos Dome, around 300 km geographically North-West from Mario Zucchelli Station (MZS), and approximately at the same geomagnetic latitude as MZS. In this work we present a preliminary analysis of the geomagnetic field 1-min data, and a comparison with simultaneous data from different Antarctic stations.

Description: Published

Description: Vienna, Austria

Description: 1.6. Osservazioni di geomagnetismo

Description: open

Description: Abstract. When combined, the three-dimensional imaging of different physical properties of architectural monumen- tal structures acquired through different methodologies can highlight with efficiency the characteristics of the stone building materials. In this work, we compound high res- olution Digital Color Images (DCI) and Terrestrial Laser Scanner (TLS) data for a dense 3-D reconstruction of an ancient pillar in a nineteenth century building in the town of Cagliari, Italy. The TLS technique was supported by a digital photogrammetry survey in order to obtain a natural color texturized 3-D model of the studied pillar. Geometri- cal anomaly maps showing interesting analogies were com- puted both from the 3-D model derived from the TLS ap- plication and from the high resolution 3-D model detected with the photogrammetry. Starting from the 3-D reconstruc- tion from previous techniques, an acoustic tomography in a sector of prior interest of the investigated architectural ele- ment was planned and carried out. The ultrasonic tomogra- phy proved to be an effective tool for detecting internal decay or defects, locating the position of the anomalies and estimat- ing their sizes, shapes, and characteristics in terms of elastic- mechanical properties. Finally, the combination of geophysi- cal and petrographical data sets represents a powerful method for understanding the quality of the building stone materials in the shallow and inner parts of the investigated architectural structures.

Description: Regione Autonoma della Sardegna (RAS) (Sardinian Autonomous Region), Regional Law 7th August 2007, no. 7, Promotion of scientific research and technological innovation in Sardinia (Italy).

Description: Published

Description: 57-62

Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica

Description: JCR Journal

Description: The Istituto Nazionale di Geofisica e Vulcanologia runs the Italian National Seismic Network (about 400 stations, seismometers, accelerometers and GPS antennas) and other networks at national scale for monitoring earthquakes and tsunami as a part of the National Civil Protection System coordinated by the Italian Department of Civil Protection. This work summarises the acquisition and the distribution of the data and the analysis that are carried out for seismic surveillance and tsunami alert.

Description: INGV and DPC

Description: Published

Description: 31-38

Description: 1IT. Reti di monitoraggio

Description: N/A or not JCR

Description: In this study temporal variations of coccolithophore blooms are investigated using satellite data. Eight years, from 2003 to 2010, of data of SCIAMACHY, a hyper-spectral satellite sensor on-board ENVISAT, were processed by the PhytoDOAS method to 5 monitor the biomass of coccolithophores in three selected regions. These regions are characterized by frequent occurrence of large coccolithophore blooms. The retrieval results, shown as monthly mean time-series, were compared to related satellite products, including the total surface phytoplankton, i.e., total chlorophyll-a (from GlobColour merged data) and the particulate inorganic carbon (from MODIS-Aqua). The 10 inter-annual variations of the phytoplankton bloom cycles and their maximum monthly mean values have been compared in the three selected regions to the variations of the geophysical parameters: sea-surface temperature (SST), mixed-layer depth (MLD) and surface wind speed, which are known to affect phytoplankton dynamics. For each region the anomalies and linear trends of the monitored parameters over the period of this 15 study have been computed. The patterns of total phytoplankton biomass and specific dynamics of coccolithophores chlorophyll-a in the selected regions are discussed in relation to other studies. The PhytoDOAS results are consistent with the two other ocean color products and support the reported dependencies of coccolithophore biomass’ dynamics to the compared geophysical variables. This suggests, that PhytoDOAS 20 is a valid method for retrieving coccolithophore biomass and for monitoring its bloom developments in the global oceans. Future applications of time-series studies using the PhytoDOAS data set are proposed, also using the new upcoming generations of hyper-spectral satellite sensors with improved spatial resolution.

Description: The gradual cooling of the climate during the Cenozoic has generally been attributed to a decrease in CO2 concentration in the atmosphere. The lack of transient climate models and in particular the lack of high-resolution proxy records of CO2, beyond the ice-core record prohibit however a full understanding of for example the inception of the Northern Hemisphere glaciation and mid-Pleistocene transition. Here we elaborate on an inverse modelling technique to reconstruct a continuous CO2 series over the past 20 million year (Myr), by decomposing the global deep-sea benthic d18O record into a mutually consistent temperature and sea level record, using a set of 1-D models of the major Northern and Southern Hemisphere ice sheets. We subsequently compared the modelled temperature record with ice core and proxy-derived CO2 data to create a continuous CO2 reconstruction over the past 20 Myr. Results show a gradual decline from 450 ppmv around 15 Myr ago to 225 ppmv for mean conditions of the glacial-interglacial cycles of the last 1 Myr, coinciding with a gradual cooling of the global surface temperature of 10 K. Between 13 to 3 Myr ago there is no long-term sea level variation caused by ice-volume changes. We find no evidence for a change in the long-term relation between temperature change and CO2, other than the effect following from the saturation of the absorption bands for CO2. The reconstructed CO2 record shows that the Northern Hemisphere glaciation starts once the long-term average CO2 concentration drops below 265 ppmv after a period of strong decrease in CO2. Finally, only a small long-term decline of 23 ppmv is found during the mid-Pleistocene transition, constraining theories on this major transition in the climate system. The approach is not accurate enough to revise current ideas about climate sensitivity.

Description: Permafrost presence is determined by a complex interaction of climatic, topographic, and ecological conditions operating over long time scales. In particular, vegetation and organic layer characteristics may act to protect permafrost in regions with a mean annual air temperature (MAAT) above 0 °C. In this study, we document the presence of residual permafrost plateaus on the western Kenai Peninsula lowlands of southcentral Alaska, a region with a MAAT of 1.5 ± 1 °C (1981 to 2010). Continuous ground temperature measurements between 16 September 2012 and 15 September 2015, using calibrated thermistor strings, documented the presence of warm permafrost (−0.04 to −0.08 °C). Field measurements (probing) on several plateau features during the fall of 2015 showed that the depth to the permafrost table averaged 1.48 m but was as shallow as 0.53 m. Late winter surveys (drilling, coring, and GPR) in 2016 showed that the average seasonally frozen ground thickness was 0.45 m, overlying a talik above the permafrost table. Measured permafrost thickness ranged from 0.33 to 〉 6.90 m. Manual interpretation of historic aerial photography acquired in 1950 indicates that residual permafrost plateaus covered 920 ha as mapped across portions of four wetland complexes encompassing 4810 ha. However, between 1950 and ca. 2010, permafrost plateau extent decreased by 60 %, with lateral feature degradation accounting for 85 % of the reduction in area. Permafrost loss on the Kenai Peninsula is likely associated with a warming climate, wildfires that remove the protective forest and organic layer cover, groundwater flow at depth, and lateral heat transfer from wetland surface waters in the summer. Better understanding the resilience and vulnerability of ecosystem-protected permafrost is critical for mapping and predicting future permafrost extent and degradation across all permafrost regions that are currently warming. Further work should focus on reconstructing permafrost history in southcentral Alaska as well as additional contemporary observations of these ecosystem-protected permafrost sites lying south of the regions with relatively stable permafrost.

Description: Whereas ice cores from high-accumulation sites in coastal Antarctica clearly demonstrate annual layering, it is debated whether a seasonal signal is also preserved in ice cores from lower-accumulation sites further inland and particularly on the East Antarctic Plateau. In this study, we examine 5 m of early Holocene ice from the Dome Fuji (DF) ice core at a high temporal resolution by continuous flow analysis. The ice was continuously analysed for concentrations of dust, sodium, ammonium, liquid conductivity, and water isotopic composition. Furthermore, a dielectric profiling was performed on the solid ice. In most of the analysed ice, the multi-parameter impurity data set appears to resolve the seasonal variability although the identification of annual layers is not always unambiguous. The study thus provides information on the snow accumulation process in central East Antarctica. A layer counting based on the same principles as those previously applied to the NGRIP (North Greenland Ice core Project) and the Antarctic EPICA (European Project for Ice Coring in Antarctica) Dronning Maud Land (EDML) ice cores leads to a mean annual layer thickness for the DF ice of 3.0 ± 0.3 cm that compares well to existing estimates. The measured DF section is linked to the EDML ice core through a characteristic pattern of three significant acidity peaks that are present in both cores. The corresponding section of the EDML ice core has recently been dated by annual layer counting and the number of years identified independently in the two cores agree within error estimates. We therefore conclude that, to first order, the annual signal is preserved in this section of the DF core. This case study demonstrates the feasibility of determining annually deposited strata on the central East Antarctic Plateau. It also opens the possibility of resolving annual layers in the Eemian section of Antarctic ice cores where the accumulation is estimated to have been greater than in the Holocene. © Author(s) 2015.

Description: The widely used detailed SNOWPACK model has undergone constant development over the years. A notable recent extension is the introduction of a Richards equation (RE) solver as an alternative for the bucket-type approach for describing water transport in the snow and soil layers. In addition, continuous updates of snow settling and new snow density parameterizations have changed model behavior. This study presents a detailed evaluation of model performance against a comprehensive multiyear data set from Weissfluhjoch near Davos, Switzerland. The data set is collected by automatic meteorological and snowpack measurements and manual snow profiles. During the main winter season, snow height (RMSE: 〈 4.2 cm), snow water equivalent (SWE, RMSE: 〈 40 mm w.e.), snow temperature distributions (typical deviation with measurements: 〈 1.0 °C) and snow density (typical deviation with observations: 〈 50 kg m−3) as well as their temporal evolution are well simulated in the model and the influence of the two water transport schemes is small. The RE approach reproduces internal differences over capillary barriers but fails to predict enough grain growth since the growth routines have been calibrated using the bucket scheme in the original SNOWPACK model. However, the agreement in both density and grain size is sufficient to parameterize the hydraulic properties successfully. In the melt season, a pronounced underestimation of typically 200 mm w.e. in SWE is found. The discrepancies between the simulations and the field data are generally larger than the differences between the two water transport schemes. Nevertheless, the detailed comparison of the internal snowpack structure shows that the timing of internal temperature and water dynamics is adequately and better represented with the new RE approach when compared to the conventional bucket scheme. On the contrary, the progress of the meltwater front in the snowpack as detected by radar and the temporal evolution of the vertical distribution of melt forms in manually observed snow profiles do not support this conclusion. This discrepancy suggests that the implementation of RE partly mimics preferential flow effects.

Description: The ongoing disintegration of large ice shelf parts in Antarctica raise the need for a better understanding of the physical processes that trigger critical crack growth in ice shelves. Finite elements in combination with configurational forces facilitate the analysis of single surface fractures in ice under various boundary conditions and material parameters. The principles of linear elastic fracture mechanics are applied to show the strong influence of different depth dependent functions for the density and the Young’s modulus on the stress intensity factor KI at the crack tip. Ice, for this purpose, is treated as an elastically compressible solid and the conse- quences of this choice in comparison to the predominant in- compressible approaches are discussed. The computed stress intensity factors KI for dry and water filled cracks are com- pared to critical values KIc from measurements that can be found in literature.

Description: The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (δ13Catm), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of δ13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 yr BP). The dataset is archived on the data repository PANGEA® (www.pangea.de) under 10.1594/PANGAEA.817041. The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4‰ shift to heavier values between the mean δ13Catm level in the Penultimate (~ 140 000 yr BP) and Last Glacial Maximum (~ 22 000 yr BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 yr, but with different phasing and magnitudes. Furthermore, a 5000 yr lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS5.5 (120 000 yr BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

Description: Stable carbon isotope analysis of methane (δ13C of CH4) on atmospheric samples is one key method to constrain the current and past atmospheric CH4 budget. A frequently applied measurement technique is gas chromatography (GC) isotope ratio mass spectrometry (IRMS) coupled to a combustion-preconcentration unit. This report shows that the atmospheric trace gas krypton (Kr) can severely interfere during the mass spectrometric measurement, leading to significant biases in δ13C of CH4, if krypton is not sufficiently separated during the analysis. According to our experiments, the krypton interference is likely composed of two individual effects, with the lateral tailing of the doubly charged 86Kr peak affecting the neighbouring m/z 44 and partially the m/z 45 Faraday cups. Additionally, a broad signal affecting m/z 45 and especially m/z 46 is assumed to result from scattered ions of singly charged krypton. The introduced bias in the measured isotope ratios is dependent on the chromatographic separation, the krypton-to-CH4 mixing ratio in the sample, the focusing of the mass spectrometer as well as the detector configuration and can amount to up to several per mil in δ13C. Apart from technical solutions to avoid this interference, we present correction routines to a posteriori remove the bias.

Description: Permafrost is one of the essential climate variables addressed by the Global Terrestrial Observing System (GCOS). Remote sensing data provide area-wide monitoring of e.g. surface temperatures or soil surface status (frozen or thawed state) in the Arctic and Subarctic, where ground data collection is difficult and restricted to local measurements at few monitoring sites. The task of the ESA Data User Element (DUE) Permafrost project is to build-up an Earth observation service for northern high-latitudinal permafrost applications with extensive involvement of the international permafrost research community (www.ipf.tuwien.ac.at/permafrost). The satellite-derived DUE Permafrost products are Land Surface Temperature, Surface Soil Moisture, Surface Frozen and Thawed State, Digital Elevation Model (locally as remote sensing product and circumpolar as non-remote sensing product) and Subsidence, and Land Cover. Land Surface Temperature, Surface Soil Moisture, and Surface Frozen and Thawed State will be provided for the circumpolar permafrost area north of 55° N with 25 km spatial resolution. In addition, regional products with higher spatial resolution were developed for five case study regions in different permafrost zones of the tundra and taiga (Laptev Sea [RU], Central Yakutia [RU], Western Siberia [RU], Alaska N-S transect, [US] Mackenzie River and Valley [CA]). This study shows the evaluation of two DUE Permafrost regional products, Land Surface Temperature and Surface Frozen and Thawed State, using freely available ground truth data from the Global Terrestrial Network of Permafrost (GTN-P) and monitoring data from the Russian-German Samoylov research station in the Lena River Delta (Central Siberia, RU). The GTN-P permafrost monitoring sites with their position in different permafrost zones are highly qualified for the validation of DUE Permafrost remote sensing products. Air and surface temperatures with high-temporal resolution from eleven GTN-P sites in Alaska and four sites in Siberia were used to match up LST products. Daily average GTN-P borehole- and air temperature data for three Alaskan and six Western Siberian sites were used to evaluate surface frozen and thawed. First results are promising and demonstrate the great benefit of freely available ground truth databases for remote sensing products.

Description: The Toba eruption that occurred some 74 ka ago in Sumatra, Indonesia, is among the largest volcanic events on Earth over the last 2 million years. Tephra from this eruption has been spread over vast areas in Asia, where it constitutes a major time marker close to the Marine Isotope Stage 4/5 boundary. As yet, no tephra associated with Toba has been identified in Greenland or Antarctic ice cores. Based on new accurate dating of Toba tephra and on accurately dated European stalagmites, the Toba event is known to occur between the onsets of Greenland interstadials (GI) 19 and 20. Furthermore, the existing linking of Greenland and Antarctic ice cores by gas records and by the bipolar seesaw hypothesis suggests that the Antarctic counterpart is situated between Antarctic Isotope Maxima (AIM) 19 and 20. In this work we suggest a direct synchronization of Greenland (NGRIP) and Antarctic (EDML) ice cores at the Toba eruption based on matching of a pattern of bipolar volcanic spikes. Annual layer counting between volcanic spikes in both cores allows for a unique match. We first demonstrate this bipolar matching technique at the already synchronized Laschamp geomagnetic excursion (41 ka BP) before we apply it to the suggested Toba interval. The Toba synchronization pattern covers some 2000 yr in GI-20 and AIM-19/20 and includes nine acidity peaks that are recognized in both ice cores. The suggested bipolar Toba synchronization has decadal precision. It thus allows a determination of the exact phasing of inter-hemispheric climate in a time interval of poorly constrained ice core records, and it allows for a discussion of the climatic impact of the Toba eruption in a global perspective. The bipolar linking gives no support for a long-term global cooling caused by the Toba eruption as Antarctica experiences a major warming shortly after the event. Furthermore, our bipolar match provides a way to place palaeo-environmental records other than ice cores into a precise climatic context.

Description: Here we present results of the first comprehensive study of sulphur compounds and methane in the oligotrophic tropical West Pacific Ocean. The concentrations of dimethylsuphide (DMS), dimethylsulphoniopropionate (DMSP), dimethylsulphoxide (DMSO), and methane (CH4), as well as various phytoplankton marker pigments in the surface ocean were measured along a north-south transit from Japan to Australia in October 2009. DMS (0.9 nmol l−1), dissolved DMSP (DMSPd, 1.6 nmol l−1) and particulate DMSP (DMSPp, 2 nmol l−1) concentrations were generally low, while dissolved DMSO (DMSOd, 4.4 nmol l−1) and particulate DMSO (DMSOp, 11.5 nmol l−1) concentrations were comparably enhanced. Positive correlations were found between DMSO and DMSP as well as DMSP and DMSO with chlorophyll a, which suggests a similar source for both compounds. Similar phytoplankton groups were identified as being important for the DMSO and DMSP pool, thus, the same algae taxa might produce both DMSP and DMSO. In contrast, phytoplankton seemed to play only a minor role for the DMS distribution in the western Pacific Ocean. The observed DMSPp : DMSOp ratios were very low and seem to be characteristic of oligotrophic tropical waters representing the extreme endpoint of the global DMSPp : DMSOp ratio vs. SST relationship. It is most likely that nutrient limitation and oxidative stress in the tropical West Pacific Ocean triggered enhanced DMSO production leading to an accumulation of DMSO in the sea surface. Positive correlations between DMSPd and CH4, as well as between DMSO (particulate and total) and CH4, were found along the transit. We conclude that both DMSP and DMSO serve as substrates for methanogenic bacteria in the western Pacific Ocean.

Description: The Lena Delta in Northern Siberia is one of the largest river deltas in the world. During peak discharge, after the ice melt in spring, it delivers between 60–8000 m3 s−1 of water and sediment into the Arctic Ocean. The Lena Delta and the Laptev Sea coast also constitute a continuous permafrost region. Ongoing climate change, which is particularly pronounced in the Arctic, is leading to increased rates of permafrost thaw. This has already profoundly altered the discharge rates of the Lena River. But the chemistry of the river waters which are discharged into the coastal Laptev Sea have also been hypothesized to undergo considerable compositional changes, e.g. by increasing concentrations of inorganic nutrients such as dissolved organic carbon (DOC) and methane. These physical and chemical changes will also affect the composition of the phytoplankton communities. However, before potential consequences of climate change for coastal arctic phytoplankton communities can be judged, the inherent status of the diversity and food web interactions within the delta have to be established. In 2010, as part of the AWI Lena Delta programme, the phyto- and microzooplankton community in three river channels of the delta (Trofimov, Bykov and Olenek) as well as four coastal transects were investigated to capture the typical river phytoplankton communities and the transitional zone of brackish/marine conditions. Most CTD profiles from 23 coastal stations showed very strong stratification. The only exception to this was a small, shallow and mixed area running from the outflow of Bykov channel in a northerly direction parallel to the shore. Of the five stations in this area, three had a salinity of close to zero. Two further stations had salinities of around 2 and 5 throughout the water column. In the remaining transects, on the other hand, salinities varied between 5 and 30 with depth. Phytoplankton counts from the outflow from the Lena were dominated by diatoms (Aulacoseira species) cyanobacteria (Aphanizomenon, Pseudanabaena) and chlorophytes. In contrast, in the stratified stations the plankton was mostly dominated by dinoflagellates, ciliates and nanoflagellates, with only an insignificant diatom component from the genera Chaetoceros and Thalassiosira (brackish as opposed to freshwater species). Ciliate abundance was significantly coupled with the abundance of total flagellates. A pronounced partitioning in the phytoplankton community was also discernible with depth, with a different community composition and abundance above and below the thermocline in the stratified sites. This work is a first analysis of the phytoplankton community structure in the region where Lena River discharge enters the Laptev Sea.

Description: There are a number of clear examples in the instrumental period where positive El Niño events were coincident with a severely weakened summer monsoon over India (ISM). ENSO's influence on the Indian Monsoon has therefore remained the centerpiece of various predictive schemes of ISM rainfall for over a century. The teleconnection between the monsoon and ENSO has undergone a protracted weakening since the late 1980's suggesting the strength of ENSO's influence on the monsoon may vary considerably on multidecadal timescales. The recent weakening has specifically prompted questions as to whether this shift represents a natural mode of climate variability or a fundamental change in ENSO and/or ISM dynamics due to anthropogenic warming. The brevity of empirical observations and large systematic errors in the representation of these two systems in state-of-the-art general circulation models hamper efforts to reliably assess the low frequency nature of this dynamical coupling under varying climate forcings. Here we place the 20th century ENSO-Monsoon relationship in a millennial context by assessing the phase angle between the two systems across the time spectrum using a continuous tree-ring ENSO reconstruction from North America and a speleothem oxygen isotope (δ18O) based reconstruction of the ISM. The results suggest that in the high-frequency domain (≤ 15 yr), El Niño (La Niña) events persistently lead to a weakened (strengthened) monsoon consistent with the observed relationship between the two systems during the instrumental period. However, in the low frequency domain (≥ 60 yr), periods of strong monsoon are, in general, coincident with periods of enhanced ENSO variance. This relationship is opposite to which would be predicted dynamically and leads us to conclude that ENSO is not pacing the prominent multidecadal variability that has characterized the ISM over the last millennium.

Description: Sea ice thickness information is important for sea ice modelling and ship operations. Here a method to detect the thickness of sea ice up to 50 cm during the freeze-up season based on high incidence angle observations of the Soil Moisture and Ocean Salinity (SMOS) satellite working at 1.4 GHz is suggested. By comparison of thermodynamic ice growth data with SMOS brightness temperatures, a high correlation to intensity and an anticorrelation to the difference between vertically and horizontally polarised brightness temperatures at incidence angles between 40 and 50° are found and used to develop an empirical retrieval algorithm sensitive to thin sea ice up to 50 cm thickness. The algorithm shows high correlation with ice thickness data from airborne measurements and reasonable ice thickness patterns for the Arctic freeze-up period.

Description: Following the launch of ESA's Soil Moisture and Ocean Salinity (SMOS) mission, it has been shown that brightness temperatures at a low microwave frequency of 1.4 GHz (L-band) are sensitive to sea ice properties. In the first demonstration study, sea ice thickness up to 50 cm has been derived using a semi-empirical algorithm with constant tie-points. Here, we introduce a novel iterative retrieval algorithm that is based on a thermodynamic sea ice model and a three-layer radiative transfer model, which explicitly takes variations of ice temperature and ice salinity into account. In addition, ice thickness variations within the SMOS spatial resolution are considered through a statistical thickness distribution function derived from high-resolution ice thickness measurements from NASA's Operation IceBridge campaign. This new algorithm has been used for the continuous operational production of a SMOS-based sea ice thickness data set from 2010 on. The data set is compared to and validated with estimates from assimilation systems, remote sensing data, and airborne electromagnetic sounding data. The comparisons show that the new retrieval algorithm has a considerably better agreement with the validation data and delivers a more realistic Arctic-wide ice thickness distribution than the algorithm used in the previous study (Kaleschke et al., 2012).

Description: The ice shelf caverns around Antarctica are sources of cold and fresh water which contributes to the formation of Antarctic bottom water and thus to the ventilation of the deep basins of the World Ocean. While a realistic simulation of the cavern circulation requires high resolution, because of the complicated bottom topography and ice shelf morphology, the physics of melting and freezing at the ice shelf base is relatively simple. We have developed an analytically solvable box model of the cavern thermohaline state, using the formulation of melting and freezing as in Olbers and Hellmer (2010). There is high resolution along the cavern's path of the overturning circulation whereas the cross-path resolution is fairly coarse. The circulation in the cavern is prescribed and used as a tuning parameter to constrain the solution by attempting to match observed ranges for outflow temperature and salinity at the ice shelf front as well as of the mean basal melt rate. The method, tested for six Antarctic ice shelves, can be used for a quick estimate of melt/freeze rates and the overturning rate in particular caverns, given the temperature and salinity of the inflow and the above mentioned constrains for outflow and melting. In turn, the model can also be used for testing the compatibility of remotely sensed basal mass loss with observed cavern inflow characteristics.

Description: Bio-optical measurements and sampling were carried out in the delta of the Lena River (northern Siberia, Russia) between 26 June and 4 July 2011. The aim of this study was to determine the inherent optical properties of the Lena water, i.e., absorption, attenuation, and scattering coefficients, during the period of maximum runoff. This aimed to contribute to the development of a bio-optical model for use as the basis for optical remote sensing of coastal water of the Arctic. In this context the absorption by CDOM (colored dissolved organic matter) and particles, and the concentrations of total suspended matter, phytoplankton-pigments, and carbon were measured. CDOM was found to be the most dominant parameter affecting the optical properties of the river, with an absorption coefficient of 4.5–5 m−1 at 442 nm, which was almost four times higher than total particle absorption values at visible wavelength range. The wavelenght-dependence of absorption of the different water constituents was chracterized by determining the semi logarithmic spectral slope. Mean CDOM, and detritus slopes were 0.0149 nm−1(standard deviation (stdev) = 0.0003, n = 18), and 0.0057 nm−1 (stdev = 0.0017, n = 19), respectively, values which are typical for water bodies with high concentrations of dissolved and particulate carbon. Mean chlorophyll a and total suspended matter were 1.8 mg m−3 (stdev = 0.734 n = 18) and 31.9 g m−3 (stdev = 19.94, n = 27), respectively. DOC (dissolved organic carbon) was in the range 8–10 g m−3 and the total particulate carbon (PC) in the range 0.25–1.5 g m−3. The light penetration depth (Secchi disc depth) was in the range 30–90 cm and was highly correlated with the suspended matter concentration. The period of maximum river runoff in June was chosen to obtain bio-optical data when maximum water constituents are transported into the Laptev Sea. However, we are aware that more data from other seasons and other years need to be collected to establish a general bio-optical model of the Lena water and conclusively characterize the light climate with respect to primary production.

Description: The Surface Ocean CO2 Atlas (SOCAT), an activity of the international marine carbon research community, provides access to synthesis and gridded fCO2 (fugacity of carbon dioxide) products for the surface oceans. Version 2 of SOCAT is an update of the previous release (version 1) with more data (increased from 6.3 million to 10.1 million surface water fCO2 values) and extended data coverage (from 1968–2007 to 1968–2011). The quality control criteria, while identical in both versions, have been applied more strictly in version 2 than in version 1. The SOCAT website (http://www.socat.info/) has links to quality control comments, metadata, individual data set files, and synthesis and gridded data products. Interactive online tools allow visitors to explore the richness of the data. Applications of SOCAT include process studies, quantification of the ocean carbon sink and its spatial, seasonal, year-to-year and longerterm variation, as well as initialisation or validation of ocean carbon models and coupled climate-carbon models.

Description: In a feasibility study, the potential of proxy data for the temperature and salinity during the Last Glacial Maximum (LGM, about 19 000 to 23 000 years before present) in constraining the strength of the Atlantic meridional overturning circulation (AMOC) with a general ocean circulation model was explored. The proxy data were simulated by drawing data from four different model simulations at the ocean sediment core locations of the Multiproxy Approach for the Reconstruction of the Glacial Ocean surface (MARGO) project, and perturbing these data with realistic noise estimates. The results suggest that our method has the potential to provide estimates of the past strength of the AMOC even from sparse data, but in general, paleo-sea-surface temperature data without additional prior knowledge about the ocean state during the LGM is not adequate to constrain the model. On the one hand, additional data in the deep-ocean and salinity data are shown to be highly important in estimating the LGM circulation. On the other hand, increasing the amount of surface data alone does not appear to be enough for better estimates. Finally, better initial guesses to start the state estimation procedure would greatly improve the performance of the method. Indeed, with a sufficiently good first guess, just the sea-surface temperature data from the MARGO project promise to be sufficient for reliable estimates of the strength of the AMOC.

Description: This paper describes ESM-SnowMIP, an international coordinated modelling effort to evaluate current snow schemes, including snow schemes that are included in Earth system models, in a wide variety of settings against local and global observations. The project aims to identify crucial processes and characteristics that need to be improved in snow models in the context of local- and global-scale modelling. A further objective of ESM-SnowMIP is to better quantify snow-related feedbacks in the Earth system. Although it is not part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6), ESM-SnowMIP is tightly linked to the CMIP6-endorsed Land Surface, Snow and Soil Moisture Model Intercomparison (LS3MIP).

Description: In this study, we present results obtained from modelling the mid-Pliocene warm period using the Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) with the two different sets of boundary conditions prescribed for the two phases of the Pliocene Model Intercomparison Project (PlioMIP). Boundary conditions, model forcing, and modelling methodology of the two phases of PlioMIP, tagged PlioMIP1 and PlioMIP2,differ considerably in palaeogeography, in particular with regards to the state of ocean gateways, ice-masks, vegetation and topography. Further differences between model setups as suggested for PlioMIP1 and PlioMIP2 consider updates to the concentration of atmospheric carbon dioxide (CO2), that is specified as 405 and 400 parts per million by volume (ppmv) for PlioMIP1 and PlioMIP2, respectively, as well as minor differences in the concentrations of methane (CH4) and nitrous oxide (N2O) due to changes in the protocol of the Paleoclimate Model Intercomparison Project (PMIP) from phase 3 to phase 4. With this manuscript, we bridge the gap between our contributions to PlioMIP1 (Stepanek and Lohmann, 2012) and PlioMIP2 (Stepanek et al., 2019). We highlight some of the effects that differences in the chosen Mid-Pliocene model setup (PlioMIP2 vs. PlioMIP1) have on the climate state as derived with the COSMOS, as this information will be valuable in the framework of the model-model and model-data-comparison within PlioMIP2. We evaluate the model sensitivity to improved mid-Pliocene boundary conditions using PlioMIP’s core mid-Pliocene experiments for PlioMIP1 and PlioMIP2, and present further simulations where we test model sensitivity to variations in palaeogeography, orbit and concentration of CO2. Firstly,we highlight major changes in boundary conditions from PlioMIP1 to PlioMIP2 and also the limitations recorded from the initial effort. The results derived from of our simulations show that COSMOS simulates a mid-Pliocene climate state that is 0.08 K colder in PlioMIP2, if compared to PlioMIP1. On one hand, high-latitude warming,which is supported by proxy evidence of the mid-Pliocene, is underestimated in simulations of both PlioMIP1 andPlioMIP2. On the other hand, spatial variations in surface air temperature (SAT), sea surface temperature (SST) as well as the distribution of sea ice suggest improvement of simulated SAT and SST in PlioMIP2 if employing the updated palaeogeography. The PlioMIP2 Mid-Pliocene simulation produces warmer SSTs in the Arctic and North Atlantic Ocean than derived from the respective PlioMIP1 climate state. The difference in prescribed CO2accountsfor 1.1 K of warming in the Arctic, leading to an ice-free summer in the PlioMIP1 simulation, and a quasi-ice-free summer in PlioMIP2. Furthermore, employing different orbital forcings in simulating the Mid-Pliocene lead to pronounced annual and seasonal variations, which is not accounted for by marine and terrestrial reconstruction of the time-slice.

Description: The Antarctic temperature changes over the past millennia remain more uncertain than in many other continental regions. This has several origins: (1) the number of high-resolution ice cores is small, in particular on the East Antarctic plateau and in some coastal areas in East Antarctica; (2) the short and spatially sparse instrumental records limit the calibration period for reconstructions and the assessment of the methodologies; (3) the link between isotope records from ice cores and local climate is usually complex and dependent on the spatial scales and timescales investigated. Here, we use climate model results, pseudo-proxy experiments and data assimilation experiments to assess the potential for reconstructing the Antarctic temperature over the last 2 millennia based on a new database of stable oxygen isotopes in ice cores compiled in the frame- work of Antarctica2k (Stenni et al., 2017). The well-known covariance between δ18O and temperature is reproduced in the two isotope-enabled models used (ECHAM5/MPI-OM and ECHAM5-wiso), but is generally weak over the different Antarctic regions, limiting the skill of the reconstructions. Furthermore, the strength of the link displays large variations over the past millennium, further affecting the potential skill of temperature reconstructions based on statistical methods which rely on the assumption that the last decades are a good estimate for longer temperature reconstructions. Using a data assimilation technique allows, in theory, for changes in the δ18O–temperature link through time and space to be taken into account. Pseudoproxy experiments confirm the benefits of using data assimilation methods instead of statistical methods that provide reconstructions with unrealistic variances in some Antarctic subregions. They also confirm that the relatively weak link between both variables leads to a limited potential for reconstructing temperature based on δ18O. However, the reconstruction skill is higher and more uniform among reconstruction methods when the reconstruction target is the Antarctic as a whole rather than smaller Antarctic subregions. This consistency between the methods at the large scale is also observed when reconstructing temperature based on the real δ18O regional composites of Stenni et al. (2017). In this case, temperature reconstructions based on data assimilation confirm the long-term cooling over Antarctica during the last millennium, and the later onset of anthropogenic warming compared with the simulations without data assimilation, which is especially visible in West Antarctica. Data assimilation also allows for models and direct observations to be reconciled by reproducing the east–west contrast in the recent temperature trends. This recent warming pattern is likely mostly driven by internal variability given the large spread of individual Paleoclimate Modelling Intercomparison Project (PMIP)/Coupled Model Intercomparison Project (CMIP) model realizations in simulating it. As in the pseudoproxy framework, the reconstruction methods perform differently at the subregional scale, especially in terms of the variance of the time series produced. While the potential benefits of using a data assimilation method instead of a statistical method have been highlighted in a pseudoproxy framework, the instrumental series are too short to confirm this in a realistic setup.

Description: A still open question is how equilibrium warming in response to increasing radiative forcing (equilibrium climate sensitivity S) is depending on background climate. We here bring paleo-data based evidence on the state-dependency of S by using CO2 proxy data together with model-based reconstruction of land ice albedo over the last 5 million years. We find that the land-ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcings depends on the CO2 data set used. Over the last 2 million years the combined S_[CO2,LI] from CO2 and land-ice albedo forcing is state-dependent and during interglacials at least twice as high as during glacials, thus CO2 doubling leads to an interglacial warming of 5 K. In the Pliocene data uncertainties prevents a well-supported calculation, but our analysis suggests that S_[CO2,LI] during a land-ice free northern hemisphere was smaller than during interglacials of the Pleistocene.

Description: Thermokarst lakes are typical features of the northern permafrost ecosystems, and play an important role in the thermal exchange between atmosphere and subsurface. The objective of this study is to describe the main thermal processes of the lakes and to quantify the heat exchange with the underlying sediments. The thermal regimes of five lakes located within the continuous permafrost zone of northern Siberia (Lena River Delta) were investigated using hourly water temperature and water level records covering a 3-year period (2009–2012), together with bathymetric survey data. The lakes included thermokarst lakes located on Holocene river terraces that may be connected to Lena River water during spring flooding, and a thermokarst lake located on deposits of the Pleistocene Ice Complex. Lakes were covered by ice up to 2 m thick that persisted for more than 7 months of the year, from October until about mid-June. Lake-bottom temperatures increased at the start of the ice-covered period due to upward-directed heat flux from the underlying thawed sediment. Prior to ice break-up, solar radiation effectively warmed the water beneath the ice cover and induced convective mixing. Ice break-up started at the beginning of June and lasted until the middle or end of June. Mixing occurred within the entire water column from the start of ice break-up and continued during the ice-free periods, as confirmed by the Wedderburn numbers, a quantitative measure of the balance between wind mixing and stratification that is important for describing the biogeochemical cycles of lakes. The lake thermal regime was modeled numerically using the FLake model. The model demonstrated good agreement with observations with regard to the mean lake temperature, with a good reproduction of the summer stratification during the ice-free period, but poor agreement during the ice-covered period. Modeled sensitivity to lake depth demonstrated that lakes in this climatic zone with mean depths 〉 5 m develop continuous stratification in summer for at least 1 month. The modeled vertical heat flux across the bottom sediment tends towards an annual mean of zero, with maximum downward fluxes of about 5 W m−2 in summer and with heat released back into the water column at a rate of less than 1 W m−2 during the ice-covered period. The lakes are shown to be efficient heat absorbers and effectively distribute the heat through mixing. Monthly bottom water temperatures during the ice-free period range up to 15 °C and are therefore higher than the associated monthly air or ground temperatures in the surrounding frozen permafrost landscape. The investigated lakes remain unfrozen at depth, with mean annual lake-bottom temperatures of between 2.7 and 4 °C.

Description: Ice-wedge polygons are common features of lowland tundra in the continuous permafrost zone and prone to rapid degradation through melting of ground ice. There are many interrelated processes involved in ice-wedge thermokarst and it is a major challenge to quantify their influence on the stability of the permafrost underlying the landscape. In this study we used a numerical modelling approach to investigate the degradation of ice wedges with a focus on the influence of hydrological conditions. Our study area was Samoylov Island in the Lena River delta of northern Siberia, for which we had in situ measurements to evaluate the model. The tailored version of the CryoGrid 3 land surface model was capable of simulating the changing microtopography of polygonal tundra and also regarded lateral fluxes of heat, water, and snow. We demonstrated that the approach is capable of simulating ice-wedge degradation and the associated transition from a low-centred to a high-centred polygonal microtopography. The model simulations showed ice-wedge degradation under recent climatic conditions of the study area, irrespective of hydrological conditions. However, we found that wetter conditions lead to an earlier onset of degradation and cause more rapid ground subsidence. We set our findings in correspondence to observed types of ice-wedge polygons in the study area and hypothesized on remaining discrepancies between modelled and observed ice-wedge thermokarst activity. Our quantitative approach provides a valuable complement to previous, more qualitative and conceptual, descriptions of the possible pathways of ice-wedge polygon evolution. We concluded that our study is a blueprint for investigating thermokarst landforms and marks a step forward in understanding the complex interrelationships between various processes shaping ice-rich permafrost landscapes.

Description: Warming of the Arctic led to an increase in permafrost temperatures by about 0.3 �C during the last decade. Permafrost warming is associated with increasing sediment water content, permeability, and diffusivity and could in the long term alter microbial community composition and abundance even before permafrost thaws. We studied the long-term effect (up to 2500 years) of submarine permafrost warming on microbial communities along an onshore–offshore transect on the Siberian Arctic Shelf displaying a natural temperature gradient of more than 10 �C. We analysed the in situ development of bacterial abundance and community composition through total cell counts (TCCs), quantitative PCR of bacterial gene abundance, and amplicon sequencing and correlated the microbial community data with temperature, pore water chemistry, and sediment physicochemical parameters. On timescales of centuries, permafrost warming coincided with an overall decreasing microbial abundance, whereas millennia after warming microbial abundance was similar to cold onshore permafrost. In addition, the dissolved organic carbon content of all cores was lowest in submarine permafrost after millennial-scale warming. Based on correlation analysis, TCC, unlike bacterial gene abundance, showed a significant rank-based negative correlation with increasing temperature, while bacterial gene copy numbers showed a strong negative correlation with salinity. Bacterial community composition correlated only weakly with temperature but strongly with the pore water stable isotopes �18O and �D, as well as with depth. The bacterial community showed substantial spatial variation and an overall dominance of Actinobacteria, Chloroflexi, Firmicutes, Gemmatimonadetes, and Proteobacteria, which are amongst the microbial taxa that were also found to be active in other frozen permafrost environments. We suggest that, millennia after permafrost warming by over 10 �C, microbial community composition and abundance show some indications for proliferation but mainly reflect the sedimentation history and paleoenvironment and not a direct effect through warming.

Description: The late Pleistocene Yedoma Ice Complex is an ice-rich and organic-bearing type of permafrost deposit widely distributed across Beringia and is assumed to be especially prone to deep degradation with warming temperature, which is a potential tipping point of the climate system. To better understand Yedoma formation, its local characteristics, and its regional sedimentological composition, we compiled the grain-size distributions (GSDs) of 771 samples from 23 Yedoma locations across the Arctic; samples from sites located close together were pooled to form 17 study sites. In addition, we studied 160 samples from three non-Yedoma ice-wedge polygon and floodplain sites for the comparison of Yedoma samples with Holocene depositional environments. The multimodal GSDs indicate that a variety of sediment production, transport, and depositional processes were involved in Yedoma formation. To disentangle these processes, a robust endmember modeling analysis (rEMMA) was performed. Nine robust grain-size endmembers (rEMs) characterize Yedoma deposits across Beringia. The study sites of Yedoma deposits were finally classified using cluster analysis. The resulting four clusters consisted of two to five sites that are distributed randomly across northeastern Siberia and Alaska, suggesting that the differences are associated with rather local conditions. In contrast to prior studies suggesting a largely aeolian contribution to Yedoma sedimentation, the wide range of rEMs indicates that aeolian sedimentation processes cannot explain the entire variability found in GSDs of Yedoma deposits. Instead, Yedoma sedimentation is controlled by local conditions such as source rocks and weathering processes, nearby paleotopography, and diverse sediment transport processes. Our findings support the hypothesis of a polygenetic Yedoma origin involving alluvial, fluvial, and niveo-aeolian transport; accumulation in ponding waters; and in situ frost weathering as well as postdepositional processes of solifluction, cryoturbation, and pedogenesis. The characteristic rEM composition of the Yedoma clusters will help to improve how grain-size-dependent parameters in permafrost models and soil carbon budgets are considered. Our results show the characteristic properties of ice-rich Yedoma deposits in the terrestrial Arctic. Characterizing and quantifying site-specific past depositional processes is crucial for elucidating and understanding the trajectories of this unique kind of ice-rich permafrost in a warmer future.

Description: Deformation of ice in continental sized ice sheets determines the flow behavior of ice towards the sea. Basal dislocation glide is assumed to be the dominant deformation mechanism in the creep deformation of natural ice, but non-basal glide is active as well. Knowledge of what types of deformation mechanisms are active in polar ice is critical in predicting the response of ice sheets in future warmer climates and its contribution to sea level rise, because the activity of deformation mechanisms depends critically on deformation conditions (such as temperature) as well as on the material properties (such as grain size). One of the methods to study the deformation mechanisms in natural materials is Electron Backscattered Diffraction (EBSD). We obtained ca. 50 EBSD maps of five different depths from a Greenlandic ice core (NEEM). The step size varied between 8 and 25 micron depending on the size of the deformation features. The size of the maps varied from 2000 to 10000 grid point. Indexing rates were up to 95%, partially by saving and reanalyzing the EBSP patterns. With this method we can characterize subgrain boundaries and determine the lattice rotation configurations of each individual subgrain. Combining these observations with arrangement/geometry of subgrain boundaries the dislocation types can be determined, which form these boundaries. Three main types of subgrain boundaries have been recognized in Antarctic (EDML) ice core (Weikusat et al. 2010, 2011). Here, we present the first results obtained from EBSD measurements performed on the NEEM ice core samples from the last glacial period, focusing on the relevance of dislocation activity of the possible slip systems. Preliminary results show that all three subgrain types, recognized in the EDML core, occur in the NEEM samples. In addition to the classical boundaries made up of basal dislocations, subgrain boundaries made of non-basal dislocations are also common.

Description: Ice is a common mineral at the Earth’s surface. How much ice is stored in the Greenland and Antarctic ice sheets depends on its mechanical properties. Therefore properties of ice directly impact on human society through its role in controlling sea level. The bulk behaviour of large ice masses is the result of the behaviour of the ensemble of individual ice grains. This is strongly influenced by the viscoplastic anisotropy of these grains and their lattice orientation. Numerical modelling provides a better insight into the mechanics of ice from the micro to the ice sheet scale. We present numerical simulations that predict the microstructural evolution of an aggregate of pure ice grains at different strain rates. We simulate co-axial deformation and dynamic recrystallization up to large strain using a full-field approach. The crystal plasticity code (Lebensohn et al., 2009) is used to calculate the response of a polycrystalline aggregate that deforms by purely dislocation glide, applying a Fast Fourier Transform (FFT). This code is coupled with the ELLE microstructural modelling platform to include intracrystalline recovery, as well as grain boundary migration driven by the reduction of surface and strain energies. The results show a strong effect of recrystallization on the final microstructure, producing larger and more equiaxed grains, with smooth boundaries. This effect does not significantly modify the single-maximum pattern of c-axes that are distributed at a low angle to the shortening direction. However, in experiments with significant recrystallization the a-axes rotate towards the elongation axis at the same time as the c-axes rotate towards the compression axis. If slip systems on prismatic and/or pyramidal planes are active, it is thought that a-axes gradually concentrate with depth (Miyamoto, 2005). The bulk activity of the slip systems is different depending on the relative activity of deformation versus recrystallization: the non-basal slip systems are more active at high strain in experiments with dynamic recrystallization compared to those experiments with low recrystallization activity.

Description: The ice sheets in Greenland and Antarctica contain a significant amount of air within their upper approximately thousand meters and air hydrates below. While this air is still in exchange with the atmosphere in the permeable firn, the gas is entrapped at the firn-ice transition at 60 – 120 m depth. Understanding the dominant deformation mechanisms is essential to interpret paleo-atmosphere records and to allow a more realistic model of ice sheet dynamics. Recent research shows how the presence of air bubbles can significantly influence microdynamical processes such as grain growth and grain boundary migration (Azuma et al., 2012, Roessiger et al., 2014). Therefore, numerical modelling was performed focussing on the mechanical properties of ice with air inclusions and the implications of the presence of bubbles on recrystallisation. The full-field crystal plasticity code of Lebensohn (2001), using a Fast Fourier Transform (FFT), was coupled to the 2D numerical microstructural modelling platform Elle, following the approach by Griera et al. (2013), and used to simulate dynamic recrystallization of pure ice (Montagnat et al., 2013). FFT calculates the viscoplastic response of polycrystalline and polyphase materials that deform by dislocation glide, takes into account the mechanical anisotropy of ice and calculates dislocation densities using the local gradient of the strain-rate field. Incorporating a code for polyphase grain boundary migration driven by surface and internal strain energy reduction, based on the methodology of Becker et al. (2008) and Roessiger et al. (2014), now also enables us to model deformation of ice with air bubbles. The presence of bubbles leads to an increase in strain localization, which reduces the bulk strength of the bubbly ice. In the absence of dynamic recrystallisation, air bubbles quickly collapse at low strains and spherical to elliptical bubble shapes are only maintained when recrystallisation is activated. Our modelling confirms that strain-induced grain boundary migration already occurs in the uppermost levels of ice sheets (Kipfstuhl et al. 2009, Weikusat et al. 2009).