ALBERT

Description: Establishing sustainable and responsible speleotourism development is a major challenge and involves complex activities. Adequate theoretical starting point is the application of geoethical values related to the conservation and protection of the caves to be used for touristic purposes. Positive and negative cases of human behaviors towards speleological geoheritage are discussed, in order to highlight what should be done in cave management to avoid malpractices and on what elements could be founded adequate strategies aimed at promoting sustainable speleotourism. This is important to tourism management organizations involved in the promotion of caves and in creating economic opportunities for local populations, while respecting cave ecosystems. Modern cave management must be focused on the protection of the cave ecosystems, finding ways to achieve at the same time an economic development of local communities. But this approach needs the adoption of a geoethical framework of values to be shared by all stakeholders involved so that successful cooperation can be achieved despite differences in interests and expectations. The aim of this paper is to raise the awareness about the need to apply the values of geoethics to speleotourism, stimulating new fields of discussion within the scientific and technical communities involved in studies and activities related to geotourism and geoheritage. The possibilities of developing new ways to manage caves, in order to promote a sustainable socio-economic development of local communities, have to be balanced with the protection of natural environments as much as possible. The proposed theoretical frameworks have the goal to increase the discussion on the best ways of connecting speleotourism to sustainable and responsible cave management, presenting two case studies, and pointing out potential solutions.

Description: Open access funding provided by Istituto Nazionale di Geofisica e Vulcanologia within the CRUI-CARE Agreement

Description: Published

Description: id 73

Description: 7SR AMBIENTE – Servizi e ricerca per la società

Description: 1TM. Formazione

Description: 2TM. Divulgazione Scientifica

Description: JCR Journal

Description: Large earthquakes occurring worldwide have long been recognized to be non Poisson distributed, so involving some large scale correlation mechanism, which could be internal or external to the Earth. Till now, no statistically significant correlation of the global seismicity with one of the possible mechanisms has been demonstrated yet. In this paper, we analyze 20 years of proton density and velocity data, as recorded by the SOHO satellite, and the worldwide seismicity in the corresponding period, as reported by the ISC-GEM catalogue. We found clear correlation between proton density and the occurrence of large earthquakes (M 〉 5.6), with a time shift of one day. The significance of such correlation is very high, with probability to be wrong lower than 10-5. The correlation increases with the magnitude threshold of the seismic catalogue. A tentative model explaining such a correlation is also proposed, in terms of the reverse piezoelectric effect induced by the applied electric field related to the proton density. This result opens new perspectives in seismological interpretations, as well as in earthquake forecast.

Description: Published

Description: 11495

Description: 7T. Variazioni delle caratteristiche crostali e precursori sismici

Description: JCR Journal

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: We reconstruct the composite dynamics of Mt. Vesuvius volcano in the period 2012–2019 from the study of ground deformation, seismicity, and geofluid (groundwater and fumarolic fluids) circulation and recognize complex spatio-temporal variations in these observables at medium (years) and short (months) time-scales. We interpret the observed patterns as the combined effect of structural changes affecting the volcanic edifice and variations of the dynamics of the hydrothermal system. In particular, we identify a change in the activity state of Mt. Vesuvius. After the activity reached minimum levels in 2014, the centroid of the surface manifestations migrated towards the SE. Episodic variations of co-seismic and aseismic deformation and fluid release, if analysed separately, would likely have been interpreted as pseudo-random oscillations of the background geophysical and geochemical signals. When organised in a comprehensive, multiparametric fashion, they shed light on the evolution of the volcano in 4D (x,y,z, time) space. These inferences play a crucial role in the formulation of civil protection scenarios for Mt. Vesuvius, a high risk, densely urbanized volcanic area which has never experienced unrest episodes in the modern era of instrumental volcanology.

Description: Published

Description: 965

Description: 2V. Struttura e sistema di alimentazione dei vulcani

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: A recent temperature reconstruction of global annual temperature shows Early Holocene warmth followed by a cooling trend through the Middle to Late Holocene [Marcott SA, et al., 2013, Science 339(6124):1198–1201]. This global cooling is puzzling because it is opposite from the expected and simulated global warming trend due to the retreating ice sheets and rising atmospheric greenhouse gases. Our critical reexamination of this contradiction between the reconstructed cooling and the simulated warming points to potentially significant biases in both the seasonality of the proxy reconstruction and the climate sensitivity of current climate models.

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: We present early Cretaceous to present paleobathymetric reconstructions and quantitative uncertainty estimates for the South Atlantic, offering a strong basis for studies of paleocirculation, paleoclimate and paleobiogeography. Circulation in an initially salty and anoxic ocean, restricted by the topography of the Falkland Plateau, Rio Grande Ridge and Walvis Rise, favoured deposition of thick evaporites in shallow water of the Brazilian-Angolan margins. This ceased as sea oor spreading propagated northwards, opening an equatorial gateway to shallow and intermediate circulation. This gateway, together with subsiding volcano-tectonic barriers would have played a key role in Late Cretaceous climate changes. Later deepening and widening of the South Atlantic, together with gateway opening at Drake Passage would lead, by mid-Miocene (∼15 Ma) to the establishment of modern-style thermohaline circulation.

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: Currently there is a scarcity of paleo-records related to the North Atlantic Oscillation (NAO), particularly in East-Central Europe (ECE). Here we report δ15N analysis of guano from a cave in NW Romania with the intent of reconstructing past variation in ECE hydroclimate and examine NAO impacts on winter precipitation. We argue that the δ15N values of guano indicate that the nitrogen cycle is hydrologically controlled and the δ15N values likely reflect winter precipitation related to nitrogen mineralization prior to the growing season. Drier conditions indicated by δ15N values at AD 1848–1852 and AD 1880–1930 correspond to the positive phase of the NAO. The increased frequency of negative phases of the NAO between AD 1940–1975 is contemporaneous with higher δ15N values (wetter conditions). A 4‰ decrease in δ15N values at the end of the 1970’s corresponds to a strong reduction in precipitation associated with a shift from negative to positive phase of the NAO. Using the relationship between NAO index and δ15N values in guano for the instrumental period, we reconstructed NAO-like phases back to AD 1650. Our results advocate that δ15N values of guano offer a proxy of the NAO conditions in the more distant past, helping assess its predictability.

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: The Antarctic silverfish (Pleuragramma antarctica) is a critically important forage species with a circumpolar distribution and is unique among other notothenioid species for its wholly pelagic life cycle. Previous studies have provided mixed evidence of population structure over regional and circumpolar scales. The aim of the present study was to test the recent population hypothesis for Antarctic silverfish, which emphasizes the interplay between life history and hydrography in shaping connectivity. A total of 1067 individuals were collected over 25 years from different locations on a circumpolar scale. Samples were genotyped at fifteen microsatellites to assess population differentiation and genetic structuring using clustering methods, F-statistics, and hierarchical analysis of variance. A lack of differentiation was found between locations connected by the Antarctic Slope Front Current (ASF), indicative of high levels of gene flow. However, gene flow was significantly reduced at the South Orkney Islands and the western Antarctic Peninsula where the ASF is absent. This pattern of gene flow emphasized the relevance of large-scale circulation as a mechanism for circumpolar connectivity. Chaotic genetic patchiness characterized population structure over time, with varying patterns of differentiation observed between years, accompanied by heterogeneous standard length distributions. The present study supports a more nuanced version of the genetic panmixia hypothesis that reflects physical-biological interactions over the life history.

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: With each cell division, phytoplankton create new space for primary colonization by marine bacteria. Although this surface microenvironment is available to all planktonic bacterial colonizers, we show the assembly of bacterial consortia on a cosmopolitan marine diatom to be highly specific and reproducible. While phytoplankton–bacteria interactions play fundamental roles in marine ecosystems, namely primary production and the carbon cycle, the ecological paradigm behind epiphytic microbiome assembly remains poorly understood. In a replicated and repeated primary colonization experiment, we exposed the axenic diatom Thalassiosira rotula to several complex and compositionally different bacterial inocula derived from phytoplankton species of varying degrees of relatedness to the axenic Thalassiosira host or natural seawater. This revealed a convergent assembly of diverse and compositionally different bacterial inocula, containing up to 2071 operational taxonomic units (OTUs), towards a stable and reproducible core community. Four of these OTUs already accounted for a cumulative abundance of 60%. This core community was dominated by Rhodobacteraceae (30.5%), Alteromonadaceae (27.7%), and Oceanospirillales (18.5%) which was qualitatively and quantitatively most similar to its conspecific original. These findings reject a lottery assembly model of bacterial colonization and suggest selective microhabitat filtering. This is likely due to diatom host traits such as surface properties and different levels of specialization resulting in reciprocal stable-state associations.

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: Antarctic krill (Euphausia superba) are high latitude pelagic organisms which play a key ecological role in the ecosystem of the Southern Ocean. To synchronize their daily and seasonal life-traits with their highly rhythmic environment, krill rely on the implementation of rhythmic strategies which might be regulated by a circadian clock. A recent analysis of krill circadian transcriptome revealed that their clock might be characterized by an endogenous free-running period of about 12–15 h. Using krill exposed to simulated light/dark cycles (LD) and constant darkness (DD), we investigated the circadian regulation of krill diel vertical migration (DVM) and oxygen consumption, together with daily patterns of clock gene expression in brain and eyestalk tissue. In LD, we found clear 24 h rhythms of DVM and oxygen consumption, suggesting a synchronization with photoperiod. In DD, the DVM rhythm shifted to a 12 h period, while the peak of oxygen consumption displayed a temporal advance during the subjective light phase. This suggested that in free-running conditions the periodicity of these clock-regulated output functions might reflect the shortening of the endogenous period observed at the transcriptional level. Moreover, differences in the expression patterns of clock gene in brain and eyestalk, in LD and DD, suggested the presence in krill of a multiple oscillator system. Evidence of short periodicities in krill behavior and physiology further supports the hypothesis that a short endogenous period might represent a circadian adaption to cope with extreme seasonal photoperiodic variability at high latitude.

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: High Arctic ecosystems and Indigenous livelihoods are tightly linked and exposed to climate change, yet assessing their sensitivity requires a long-term perspective. Here, we assess the vulnerability of the North Water polynya, a unique seaice ecosystem that sustains the world’s northernmost Inuit communities and several keystone Arctic species. We reconstruct mid-tolate Holocene changes in sea ice, marine primary production, and little auk colony dynamics through multi-proxy analysis of marine and lake sediment cores. Our results suggest a productive ecosystem by 4400–4200 cal yrs b2k coincident with the arrival of the first humans in Greenland. Climate forcing during the late Holocene, leading to periods of polynya instability and marine productivity decline, is strikingly coeval with the human abandonment of Greenland from c. 2200–1200 cal yrs b2k. Our long-term perspective highlights the future decline of the North Water ecosystem, due to climate warming and changing sea-ice conditions, as an important climate change risk.

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 evolution of past global ice sheets is highly uncertain. One example is the missing ice problem during the Last Glacial Maximum (LGM, 26 000-19 000 years before present) – an apparent 8-28 m discrepancy between far-field sea level indicators and modelled sea level from ice sheet reconstructions. In the absence of ice sheet reconstructions, researchers often use marine δ 18 O proxy records to infer ice volume prior to the LGM. We present a global ice sheet reconstruction for the past 80 000 years, called PaleoMIST 1.0, constructed inde- pendently of far-field sea level and δ 18 O proxy records. Our reconstruction is compatible with LGM far-field sea-level records without requiring extra ice volume, thus solving the missing ice problem. However, for Marine Isotope Stage 3 (57 000-29 000 years before present) - a pre-LGM period - our reconstruction does not match proxy-based sea level reconstructions, indicating the relationship between marine δ 18 O and sea level may be more complex than assumed.

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: Lake Baikal is inhabited by more than 300 endemic amphipod species, which are narrowly adapted to certain thermal niches due to the high interspecific competition. In contrast, the surrounding freshwater fauna is commonly represented by species with large‐scale distribution and high phenotypic thermal plasticity. Here, we investigated the thermal plasticity of the energy metabolism in two closely‐related endemic amphipod species from Lake Baikal (Eulimnogammarus verrucosus; stenothermal and Eulimnogammarus cyaneus; eurythermal) and the ubiquitous Holarctic amphipod Gammarus lacustris (eurythermal) by exposure to a summer warming scenario (6–23.6 °C; 0.8 °C d−1). In concert with routine metabolic rates, activities of key metabolic enzymes increased strongly with temperature up to 15 °C in E. verrucosus, whereupon they leveled off (except for lactate dehydrogenase). In contrast, exponential increases were seen in E. cyaneus and G. lacustris throughout the thermal trial (Q10‐values: 1.6–3.7). Cytochrome‐c‐oxidase, lactate dehydrogenase, and 3‐hydroxyacyl‐CoA dehydrogenase activities were found to be higher in G. lacustris than in E. cyaneus, especially at the highest experimental temperature (23.6 °C). Decreasing gene expression levels revealed some thermal compensation in E. cyaneus but not in G. lacustris. In all species, shifts in enzyme activities favored glycolytic energy generation in the warmth. The congruent temperature‐ dependencies of enzyme activities and routine metabolism in E. verrucosus indicate a strong feedback‐ regulation of enzymatic activities by whole organism responses. The species‐specific thermal reaction norms reflect the different ecological niches, including the spatial distribution, distinct thermal behavior such as temperature‐dependent migration, movement activity, and mating season.

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: Central Europe has experienced a severe drought almost every April for the last 14 years consecutively, driven by record high temperatures, low flows, high evapotranspiration, and high soil moisture deficit. The dynamic of this recent and recurrent mid-spring dryness is not yet understood. Here we show that the period 2007â€``2020 was characterized by a reduction of ~50% of the usual April rainfall amount over large areas in central Europe. The precipitation deficit and the record high temperatures were triggered by a multiyear recurrent high-pressure system centered over the North Sea and northern Germany and a decline in the temperature gradient between the Arctic region and the mid-latitudes, which diverted the Atlantic storm tracks northward. From a long-term perspective, the precipitation, temperature, and soil moisture anomalies observed over the last 14 years have reached the highest amplitudes over the observational record. Our study provides an in-depth analysis of the hydroclimate extremes in central Europe over the last 140 years and their atmospheric drivers, enabling us to increase our dynamical understating of long-term dry periods, which is vital to enhance forecasting and mitigation of such events.

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: Floating ice shelves are the Achilles’ heel of the Antarctic Ice Sheet. They limit Antarctica’s contribution to global sea level rise, yet they can be rapidly melted from beneath by a warming ocean. At Filchner-Ronne Ice Shelf, a decline in sea ice formation may increase basal melt rates and accelerate marine ice sheet mass loss within this century. However, the understanding of this tipping-point behavior largely relies on numerical models. Our new multi-annual observations from five hot-water drilled boreholes through Filchner-Ronne Ice Shelf show that since 2015 there has been an intensification of the density-driven ice shelf cavity-wide circulation in response to reinforced wind-driven sea ice formation in the Ronne polynya. Enhanced southerly winds over Ronne Ice Shelf coincide with westward displacements of the Amundsen Sea Low position, connecting the cavity circulation with changes in large-scale atmospheric circulation patterns as a new aspect of the atmosphere-ocean-ice shelf system.

Description: The thecosome pteropods Limacina helicina and L. retroversa are important contributors to the zooplankton community in high-latitude environments but little is known about their distribution and life cycle under polar conditions. We collected the early life stages (〈 1 mm) of the thecosome population in 2012 and 2013 at a bi-weekly to monthly resolution in fjord highly influenced by Arctic waters as well as Atlantic inflows (Adventfjorden, Svalbard, 78°N), together with environmental parameters. L. retroversa only occurred episodically, in association with the inflow of Atlantic water, with low numbers and random size distributions. This suggests that this boreal species does not fulfill its life cycle in Adventfjorden. In contrast, young specimens of L. helicina were present during the entire study. Veligers hatched in late summer/autumn and measured 0.14 mm on average. They grew with rates of 0.0006 mm day− 1 over the 10–11 months of development. Only thereafter, growth accelerated by one order of magnitude and maximal rates were reached in autumn (0.0077 mm day− 1). Our results indicate that L. helicina reaches a size of 1 mm after approximately 1.5 years in Adventfjorden. We therefore suggest that L. helicina overwinters the first year as a small juvenile and that it needs at least 2 years to reach an adult size of 5 mm in Adventfjorden. This reveals an complex and delicate aspect of the life-cycle of L. helicina and further research is needed to determine if it makes the population especially vulnerable towards climate changes.

Description: Slope failure like in the Hinlopen/Yermak Megaslide is one of the major geohazards in a changing Arctic environment. We analysed hydroacoustic and 2D high-resolution seismic data from the apparently intact continental slope immediately north of the Hinlopen/Yermak Megaslide for signs of past and future instabilities. Our new bathymetry and seismic data show clear evidence for incipient slope instability. Minor slide deposits and an internally-deformed sedimentary layer near the base of the gas hydrate stability zone imply an incomplete failure event, most probably about 30000 years ago, contemporaneous to or shortly after the Hinlopen/Yermak Megaslide. An active gas reservoir at the base of the gas hydrate stability zone demonstrate that over-pressured fluids might have played a key role in the initiation of slope failure at the studied slope, but more importantly also for the giant HYM slope failure. To date, it is not clear, if the studied slope is fully preconditioned to fail completely in future or if it might be slowly deforming and creeping at present. We detected widespread methane seepage on the adjacent shallow shelf areas not sealed by gas hydrates.

Description: Little is known about the production of fluorescent dissolved organic matter (FDOM) in the anoxic oceanic sediments. In this study, sediment pore waters were sampled from four different sites in the Chukchi-East Siberian Seas area to examine the bulk dissolved organic carbon (DOC) and their optical properties. The production of FDOM, coupled with the increase of nutrients, was observed above the sulfate-methane-transition-zone (SMTZ). The presence of FDOM was concurrent with sulfate reduction and increased alkalinity (R2 〉 0.96, p 〈 0.0001), suggesting a link to organic matter degradation. This inference was supported by the positive correlation (R2 〉 0.95, p 〈 0.0001) between the net production of FDOM and the modeled degradation rates of particulate organic carbon sulfate reduction. The production of FDOM was more pronounced in a shallow shelf site S1 with a total net production ranging from 17.9 to 62.3 RU for different FDOM components above the SMTZ depth of ca. 4.1 mbsf, which presumably underwent more accumulation of particulate organic matter than the other three deeper sites. The sediments were generally found to be the sources of CDOM and FDOM to the overlying water column, unearthing a channel of generally bio-refractory and pre-aged DOM to the oceans.

Description: The variability of sea surface temperatures (SSTs) at multidecadal and longer timescales is poorly constrained, primarily because instrumental records are short and proxy records are noisy. Through applying a new noise filtering technique to a global network of late Holocene SST proxies, we estimate SST variability between annual and millennial timescales. Filtered estimates of SST variability obtained from coral, foraminifer, and alkenone records are shown to be consistent with one another and with instrumental records in the frequency bands at which they overlap. General circulation models, however, simulate SST variability that is systematically smaller than instrumental and proxy-based estimates. Discrepancies in variability are largest at low latitudes and increase with timescale, reaching two orders of magnitude for tropical variability at millennial timescales. This result implies major deficiencies in observational estimates or model simulations, or both, and has implications for the attribution of past variations and prediction of future change.

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: Stable water isotope records from Antarctica are key for our understanding of Quaternary climate variations. However, the exact quantitative interpretation of these important climate proxy records in terms of surface temperature, ice sheet height and other climatic changes is still a matter of debate. Here we report results obtained with an atmospheric general circulation model equipped with water isotopes, run at a high-spatial horizontal resolution of one-by-one degree. Comparing different glacial maximum ice sheet reconstructions, a best model data match is achieved for the PMIP3 reconstruction. Reduced West Antarctic elevation changes between 400 and 800 m lead to further improved agreement with ice core data. Our modern and glacial climate simulations support the validity of the isotopic paleothermometer approach based on the use of present-day observations and reveal that a glacial ocean state as displayed in the GLAMAP reconstruction is suitable for capturing the observed glacial isotope changes in Antarctic ice cores.

Description: A dominant Antarctic ecological paradigm suggests that winter sea ice is generally the main feeding ground for krill larvae. Observations from our winter cruise to the southwest Atlantic sector of the Southern Ocean contradict this view and present the first evidence that the pack-ice zone is a food-poor habitat for larval development. In contrast, the more open marginal ice zone provides a more favourable food environment for high larval krill growth rates. We found that complex under-ice habitats are, however, vital for larval krill when water column productivity is limited by light, by providing structures that offer protec- tion from predators and to collect organic material released from the ice. The larvae feed on this sparse ice-associated food during the day. After sunset, they migrate into the water below the ice (upper 20 m) and drift away from the ice areas where they have previously fed. Model analyses indicate that this behaviour increases both food uptake in a patchy food environment and the likelihood of overwinter transport to areas where feeding conditions are more favourable in spring.

Description: There is a strong spatial correlation between submarine slope failures and the occurrence of gas hydrates. This has been attributed to the dynamic nature of gas hydrate systems and the potential reduction of slope stability due to bottom water warming or sea level drop. However, 30 years of research into this process found no solid supporting evidence. Here we present new reflection seismic data from the Arctic Ocean and numerical modelling results supporting a different link between hydrates and slope stability. Hydrates reduce sediment permeability and cause build-up of overpressure at the base of the gas hydrate stability zone. Resulting hydro-fracturing forms pipe structures as pathways for overpressured fluids to migrate upward. Where these pipe structures reach shallow permeable beds, this overpressure transfers laterally and destabilises the slope. This process reconciles the spatial correlation of submarine landslides and gas hydrate, and it is independent of environmental change and water depth.

Description: Arctic tundra landscapes are composed of a complex mosaic of patterned ground features, varying in soil moisture, vegetation composition, and surface hydrology over small spatial scales (10–100 m). The importance of microtopography and associated geomorphic landforms in influencing ecosystem structure and function is well founded, however, spatial data products describing local to regional scale distribution of patterned ground or polygonal tundra geomorphology are largely unavailable. Thus, our understanding of local impacts on regional scale processes (e.g., carbon dynamics) may be limited. We produced two key spatiotemporal datasets spanning the Arctic Coastal Plain of northern Alaska (~60,000 km2) to evaluate climate-geomorphological controls on arctic tundra productivity change, using (1) a novel 30m classification of polygonal tundra geomorphology and (2) decadal-trends in surface greenness using the Landsat archive (1999–2014). These datasets can be easily integrated and adapted in an array of local to regional applications such as (1) upscaling plot-level measurements (e.g., carbon/energy fluxes), (2) mapping of soils, vegetation, or permafrost, and/or (3) initializing ecosystem biogeochemistry, hydrology, and/or habitat modeling.

Description: The field of Arctic sea ice prediction on “weather time scales” is still in its infancy with little existing understanding of the limits of predictability. This is especially true for sea ice deformation along so-called Linear Kinematic Features (LKFs) including leads that are relevant for marine operations. Here the potential predictability of the sea ice pack in the wintertime Arctic up to ten days ahead is determined, exploiting the fact that sea ice-ocean models start to show skill at representing sea ice deformation at high spatial resolutions. Results are based on ensemble simulations with a high-resolution sea ice-ocean model driven by atmospheric ensemble forecasts. The predictability of LKFs as measured by different metrics drops quickly, with predictability being almost completely lost after 4–8 days. In contrast, quantities such as sea ice concentration or the location of the ice edge retain high levels of predictability throughout the full 10-day forecast period. It is argued that the rapid error growth for LKFs is mainly due to the chaotic behaviour of the atmosphere associated with the low predictability of near surface wind divergence and vorticity; initial condition uncertainty for ice thickness is found to be of minor importance as long as LKFs are initialized at the right locations.

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: Antarctic krill (Euphausia superba) is a key species in Southern Ocean ecosystem where it plays a central role in the Antarctic food web. Available information supports the existence of an endogenous timing system in krill enabling it to synchronize metabolism and behavior with an environment characterized by extreme seasonal changes in terms of day length, food availability, and surface ice extent. A screening of our transcriptome database “KrillDB” allowed us to identify the putative orthologues of 20 circadian clock components. Mapping of conserved domains and phylogenetic analyses strongly supported annotations of the identi ed sequences. Luciferase assays and co-immunoprecipitation experiments allowed us to de ne the role of the main clock components. Our ndings provide an overall picture of the molecular mechanisms underlying the functioning of the endogenous circadian clock in the Antarctic krill and shed light on their evolution throughout crustaceans speciation. Interestingly, the core clock machinery shows both mammalian and insect features that presumably contribute to an evolutionary strategy to cope with polar environment’s challenges. Moreover, despite the extreme variability characterizing the Antarctic seasonal day length, the conserved light mediated degradation of the photoreceptor EsCRY1 suggests a persisting pivotal role of light as a Zeitgeber.

Description: Arctic tundra ecosystems have experienced unprecedented change associated with climate warming over recent decades. Across the Pan-Arctic, vegetation productivity and surface greenness have trended positively over the period of satellite observation. However, since 2011 these trends have slowed considerably, showing signs of browning in many regions. It is unclear what factors are driving this change and which regions/landforms will be most sensitive to future browning. Here we provide evidence linking decadal patterns in arctic greening and browning with regional climate change and local permafrost-driven landscape heterogeneity. We analyzed the spatial variability of decadal-scale trends in surface greenness across the Arctic Coastal Plain of northern Alaska (~60,000 km²) using the Landsat archive (1999–2014), in combination with novel 30 m classifications of polygonal tundra and regional watersheds, finding landscape heterogeneity and regional climate change to be the most important factors controlling historical greenness trends. Browning was linked to increased temperature and precipitation, with the exception of young landforms (developed following lake drainage), which will likely continue to green. Spatiotemporal model forecasting suggests carbon uptake potential to be reduced in response to warmer and/or wetter climatic conditions, potentially increasing the net loss of carbon to the atmosphere, at a greater degree than previously expected.

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: Antarctic krill (Euphausia superba)—one of the most abundant animal species on Earth—exhibits a five to six year population cycle, with oscillations in biomass exceeding one order of magnitude. Previous studies have postulated that the krill cycle is induced by periodic climatological factors, but these postulated drivers neither show consistent agreement, nor are they supported by quantitative models. Here, using data analysis complemented with modelling of krill ontogeny and population dynamics, we identify intraspecific competition for food as the main driver of the krill cycle, while external climatological factors possibly modulate its phase and synchronization over large scales. Our model indicates that the cycle amplitude increases with reduction of krill loss rates. Thus, a decline of apex predators is likely to increase the oscillation amplitude, potentially destabilizing the marine food web, with drastic consequences for the entire Antarctic ecosystem.

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: 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: In the following we present a new non-invasive methodology aimed at the diagnosis of stone building materials used in historical buildings and architectural elements. This methodology consists of the integrated sequential application of in situ proximal sensing methodologies such as the 3D Terrestrial Laser Scanner for the 3D modelling of investigated objects together with laboratory and in situ non-invasive multi-techniques acoustic data, preceded by an accurate petrographical study of the investigated stone materials by optical and scanning electron microscopy. The increasing necessity to integrate different types of techniques in the safeguard of the Cultural Heritage is the result of the following two interdependent factors: 1) The diagnostic process on the building stone materials of monuments is increasingly focused on difficult targets in critical situations. In these cases, the diagnosis using only one type of non-invasive technique may not be sufficient to investigate the conservation status of the stone materials of the superficial and inner parts of the studied structures 2) Recent technological and scientific developments in the field of non-invasive diagnostic techniques for different types of materials favors and supports the acquisition, processing and interpretation of huge multidisciplinary datasets.

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: 4334

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: Understanding the temporal variation of cosmic radiation and solar activity during the Holocene is essential for studies of the solar-terrestrial relationship. Cosmic-ray produced radionuclides, such as 10Be and 14C which are stored in polar ice cores and tree rings, offer the unique opportunity to reconstruct the history of cosmic radiation and solar activity over many millennia. Although records from different archives basically agree, they also show some deviations during certain periods. So far most reconstructions were based on only one single radionuclide record, which makes detection and correction of these deviations impossible. Here we combine different 10Be ice core records from Greenland and Antarctica with the global 14C tree ring record using principal component analysis. This approach is only possible due to a new high-resolution 10Be record from Dronning Maud Land obtained within the European Project for Ice Coring in Antarctica in Antarctica. The new cosmic radiation record enables us to derive total solar irradiance, which is then used as a proxy of solar activity to identify the solar imprint in an Asian climate record. Though generally the agreement between solar forcing and Asian climate is good, there are also periods without any coherence, pointing to other forcings like volcanoes and greenhouse gases and their corresponding feedbacks. The newly derived records have the potential to improve our understanding of the solar dynamics and to quantify the solar influence on climate.

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.