ALBERT

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: Published

Description: 57-62

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

Description: JCR Journal

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

Description: INGV and DPC

Description: Published

Description: 31-38

Description: 1IT. Reti di monitoraggio

Description: N/A or not JCR

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: Ice cores are the only climate archives incorporating paleo-atmosphere as individual gas inclusions, enabling the extraction and analysis of the contained gasses. A firm understanding of the processes involved is mandatory for a reliable interpretation of the gas records. One prominent process is the transition from air bubbles to crystalline air hydrates, which is known to influence, at least temporarily, the gas mixing ratios by diffusion and fractionation. This transition is still not understood completely and the existing theories do not explain the large diversity of observed hydrate morphologies. Raman tomographic measurements using the AWI cryo-Raman system provide 3D reconstructions of air hydrate morphologies. The results show complex growth structures that emphasize the importance of crystallography, microstructure and ice rheology for the hydrate formation process. Accurate hydrate volumes can be calculated from the 3D objects, improving the estimates of total gas contents.

Description: The Asian monsoon system is an important tipping element in Earth's climate with a large impact on human societies in the past and present. In light of the potentially severe impacts of present and future anthropogenic climate change on Asian hydrology, it is vital to understand the forcing mechanisms of past climatic regime shifts in the Asian monsoon domain. Here we use novel recurrence network analysis techniques for detecting episodes with pronounced non-linear changes in Holocene Asian monsoon dynamics recorded in speleothems from caves distributed throughout the major branches of the Asian monsoon system. A newly developed multi-proxy methodology explicitly considers dating uncertainties with the COPRA (COnstructing Proxy Records from Age models) approach and allows for detection of continental-scale regime shifts in the complexity of monsoon dynamics. Several epochs are characterised by non-linear regime shifts in Asian monsoon variability, including the periods around 8.5–7.9, 5.7–5.0, 4.1–3.7, and 3.0–2.4 ka BP. The timing of these regime shifts is consistent with known episodes of Holocene rapid climate change (RCC) and high-latitude Bond events. Additionally, we observe a previously rarely reported non-linear regime shift around 7.3 ka BP, a timing that matches the typical 1.0–1.5 ky return intervals of Bond events. A detailed review of previously suggested links between Holocene climatic changes in the Asian monsoon domain and the archaeological record indicates that, in addition to previously considered longer-term changes in mean monsoon intensity and other climatic parameters, regime shifts in monsoon complexity might have played an important role as drivers of migration, pronounced cultural changes, and the collapse of ancient human societies.

Description: Thermokarst lakes are important emitters of methane, a potent greenhouse gas. However, accurate estimation of methane flux from thermokarst lakes is difficult due to their remoteness and observational challenges associated with the heterogeneous nature of ebullition (bubbling). We used multi-temporal high-resolution (9–11 cm) aerial images of an interior Alaskan thermokarst lake, Goldstream Lake, acquired 2 and 4 days following freeze-up in 2011 and 2012, respectively, to characterize methane ebullition seeps and to estimate whole-lake ebullition. Bubbles impeded by the lake ice sheet form distinct white patches as a function of bubbling rate vs. time as ice thickens. Our aerial imagery thus captured in a single snapshot the ebullition events that occurred before the image acquisition. Image analysis showed that low-flux A- and B-type seeps are associated with low brightness patches and are statistically distinct from high-flux C-type and Hotspot seeps associated with high brightness patches. Mean whole-lake ebullition based on optical image analysis in combination with bubble-trap flux measurements was estimated to be 174 ± 28 and 216 ± 33 mL gas m−2 d−1 for the years 2011 and 2012, respectively. A large number of seeps demonstrated spatio-temporal stability over our two-year study period. A strong inverse exponential relationship (R2 ≥ 0.79) was found between percent surface area of lake ice covered with bubble patches and distance from the active thermokarst lake margin. Our study shows that optical remote sensing is a powerful tool to map ebullition seeps on lake ice, to identify their relative strength of ebullition and to assess their spatio-temporal variability.

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. This update of GLODAPv2, v2.2019, adds data from 116 cruises to the previous version, extending its coverage in time from 2013 to 2017, while also adding some data from prior years. GLODAPv2.2019 includes measurements from more than 1.1 million water samples from the global oceans collected on 840 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 116 new cruises with the data from the 724 quality-controlled cruises of the GLODAPv2 data product. They correct for errors related to measurement, calibration, and data handling practices, taking into account any known or likely time trends or variations. 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 in pH, and 5 % in the halogenated transient tracers. The compilation also includes data for several other variables, such as isotopic tracers. These were not subjected to bias comparison or adjustments. The original data, their documentation and DOI codes are available in the Ocean Carbon Data System of NOAA NCEI (https://www.nodc.noaa.gov/ocads/oceans/GLODAPv2_2019/, last access: 17 September 2019). 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/xnme-wr20 (Olsen et al., 2019). The product files also include significant ancillary and approximated data. These were obtained by interpolation of, or calculation from, measured data. This paper documents the GLODAPv2.2019 methods and provides a broad overview of the secondary quality control procedures and results.

Description: We present a Lagrangian convective transport scheme developed for global chemistry and transport models, which considers the variable residence time that an air parcel spends in convection. This is particularly important for accurately simulating the tropospheric chemistry of short-lived species, e.g., for determining the time available for heterogeneous chemical processes on the surface of cloud droplets. In current Lagrangian convective transport schemes air parcels are stochastically redistributed within a fixed time step according to estimated probabilities for convective entrainment as well as the altitude of detrainment. We introduce a new scheme that extends this approach by modeling the variable time that an air parcel spends in convection by estimating vertical updraft velocities. Vertical updraft velocities are obtained by combining convective mass fluxes from meteorological analysis data with a parameterization of convective area fraction profiles. We implement two different parameterizations: a parameterization using an observed constant convective area fraction profile and a parameterization that uses randomly drawn profiles to allow for variability. Our scheme is driven by convective mass fluxes and detrainment rates that originate from an external convective parameterization, which can be obtained from meteorological analysis data or from general circulation models. We study the effect of allowing for a variable time that an air parcel spends in convection by performing simulations in which our scheme is implemented into the trajectory module of the ATLAS chemistry and transport model and is driven by the ECMWF ERA-Interim reanalysis data. In particular, we show that the redistribution of air parcels in our scheme conserves the vertical mass distribution and that the scheme is able to reproduce the convective mass fluxes and detrainment rates of ERA-Interim. We further show that the estimated vertical updraft velocities of our scheme are able to reproduce wind profiler measurements performed in Darwin, Australia, for velocities larger than 0.6 m s−1. SO2 is used as an example to show that there is a significant effect on species mixing ratios when modeling the time spent in convective updrafts compared to a redistribution of air parcels in a fixed time step. Furthermore, we perform long-time global trajectory simulations of radon-222 and compare with aircraft measurements of radon activity.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Annual Review of Marine Science 9 (2017): 173-203, doi:10.1146/annurev-marine-010816-060733.

Description: The events that followed the Tohoku earthquake and tsunami on March 11, 2011, included the loss of power and overheating at the Fukushima Daiichi nuclear power plants, which led to extensive releases of radioactive gases, volatiles, and liquids, particularly to the coastal ocean. The fate of these radionuclides depends in large part on their oceanic geochemistry, physical processes, and biological uptake. Whereas radioactivity on land can be resampled and its distribution mapped, releases to the marine environment are harder to characterize owing to variability in ocean currents and the general challenges of sampling at sea. Five years later, it is appropriate to review what happened in terms of the sources, transport, and fate of these radionuclides in the ocean. In addition to the oceanic behavior of these contaminants, this review considers the potential health effects and societal impacts.

Description: K.B. was supported in part by the Gordon and Betty Moore Foundation and the Deerbrook Charitable Trust. P.M. was supported in part by the Generalitat de Catalunya through MERS (grant 2014 SGR 1356), the European Commission 7th Framework COMET-FRAME project (grant agreement 604974), and the Ministerio de Economía y Competitividad of Spain (project CTM2011-15152-E). S.C. was supported in part by the French program Investissement d'Avenir run by the National Research Agency (AMORAD project, grant ANR-11-RSNR-0002). D.O. was supported in part by the Center for Environmental Radioactivity (NFR Centers of Excellence grant 223268/F50). J.N.S. was supported in part by the Marine Environmental Observation, Prediction, and Response Network.

Description: Coeval changes in atmospheric CO2 and 14C contents during the last deglaciation are often attributed to ocean circulation changes that released carbon stored in the deep ocean during the Last Glacial Maximum (LGM). Work is being done to generate records that allow for the identification of the exact mechanisms leading to the accumulation and release of carbon from the oceanic reservoir, but these mechanisms are still the subject of debate. Here we present foraminifera 14C data from five cores in a transect across the Chilean continental margin between ~540 and ~3,100 m depth spanning the last 20,000 years. Our data reveal that during the LGM, waters at ~2,000 m were 50% to 80% more depleted in Δ14C than waters at ~1,500 m when compared to modern values, consistent with the hypothesis of a glacial deep ocean carbon reservoir that was isolated from the atmosphere. During the deglaciation, our intermediate water records reveal homogenization in the Δ14C values between ~800 and ~1,500 m from ~16.5–14.5 ka cal BP to ~14–12 ka cal BP, which we interpret as deeper penetration of Antarctic Intermediate Water. While many questions still remain, this process could aid the ventilation of the deep ocean at the beginning of the deglaciation, contributing to the observed ~40 ppm rise in atmospheric pCO2.

Description: Holocene permafrost from ice wedge polygons in the vicinity of large seabird breeding colonies in the Thule District, NW Greenland, was drilled to explore the relation between permafrost aggradation and seabird presence. The latter is reliant on the presence of the North Water Polynya (NOW) in the northern Baffin Bay. The onset of peat accumulation associated with the arrival of little auks (Alle alle) in a breeding colony at Annikitisoq, north of Cape York, is radiocarbon-dated to 4400 cal BP. A thick-billed murre (Uria lomvia) colony on Appat (Saunders Island) in the mouth of the Wolstenholme Fjord started 5650 cal BP. Both species provide marine-derived nutrients (MDNs) that fertilize vegetation and promote peat growth. The geochemical signature of organic matter left by the birds is traceable in the frozen Holocene peat. The peat accumulation rates at both sites are highest after the onset, decrease over time, and were about 2-times faster at the little auk site than at the thick-billed murre site. High accumulation rates induce shorter periods of organic matter (OM) decomposition before it enters the perennially frozen state. This is seen in comparably high C=N ratios and less depleted 13C, pointing to a lower degree of OM decomposition at the little auk site, while the opposite pattern can be discerned at the thick-billed murre site. Peat accumulation rates correspond to 15N trends, where decreasing accumulation led to increasing depletion in 15N as seen in the little-auk-related data. In contrast, the more decomposed OM of the thick-billed murre site shows almost stable 15N. Late Holocene wedge ice fed by cold season precipitation was studied at the little auk site and provides the first stable-water isotopic record from Greenland with mean 18O of 8:00:8, mean D of 36:25:7, mean d excess of 7:70:7, and a 18O-D slope of 7.27, which is close to those of the modern Thule meteoric water line. The syngenetic ice wedge polygon development is mirrored in testacean records of the little auk site and delineates polygon low-center, dry-out, and polygon-high-center stages. The syngenetic permafrost formation directly depending on peat growth (controlled by bird activity) falls within the period of neoglacial cooling and the establishment of the NOW, thus indirectly following the Holocene climate trends.

Description: The Arctic Ocean is known to be contaminated by various persistent organic pollutants (POPs). The Fram Strait, the only deepwater passage to the Arctic Ocean (from the Atlantic Ocean), represents an unquantified gateway for POPs fluxes into and out of the Arctic. Polyethylene passive samplers were deployed in vertical profiles in the Fram Strait and in air and surface water in the Canadian Archipelago to determine the concentrations, profiles, and mass fluxes of dissolved polychlorinated biphenyls (PCBs) and organochlorine pesticides. In the Fram Strait, higher concentrations of ΣPCBs (1.3–3.6 pg/L) and dichlorodiphenyltrichloroethanes (ΣDDTs, 5.2–9.1 pg/L) were observed in the deepwater masses (below 1,000 m), similar to nutrient-like vertical profiles. There was net southward transport of hexachlorobenzene and hexachlorocyclohexanes (ΣHCHs) of 0.70 and 14 Mg/year but a net northward transport of ΣPCBs at 0.16 Mg/year through the Fram Strait.

Description: Coastal erosion and flooding transform terrestrial landscapes into marine environments. In the Arctic, these processes inundate terrestrial permafrost with seawater and create submarine permafrost. Permafrost begins to warm under marine conditions, which can destabilize the sea floor and may release greenhouse gases. We report on the transition of terrestrial to submarine permafrost at a site where the timing of inundation can be inferred from the rate of coastline retreat. On Muostakh Island in the central Laptev Sea, East Siberia, changes in annual coastline position have been measured for decades and vary highly spatially. We hypothesize that these rates are inversely related to the inclination of the upper surface of submarine ice-bonded permafrost (IBP) based on the consequent duration of inundation with increasing distance from the shoreline. We compared rapidly eroding and stable coastal sections of Muostakh Island and find permafrost-table inclinations, determined using direct current resistivity, of 1 and 5 %, respectively. Determinations of submarine IBP depth from a drilling transect in the early 1980s were compared to resistivity profiles from 2011. Based on borehole observations, the thickness of unfrozen sediment overlying the IBP increased from 0 to 14m below sea level with increasing distance from the shoreline. The geoelectrical profiles showed thickening of the unfrozen sediment overlying ice-bonded permafrost over the 28 years since drilling took place. We use geoelectrical estimates of IBP depth to estimate permafrost degradation rates since inundation. Degradation rates decreased from over 0.4ma-1 following inundation to around 0.1ma-1 at the latest after 60 to 110 years and remained constant at this level as the duration of inundation increased to 250 years. We suggest that long-term rates are lower than these values, as the depth to the IBP increases and thermal and porewater solute concentration gradients over depth decrease. For the study region, recent increases in coastal erosion rate and changes in benthic temperature and salinity regimes are expected to affect the depth to submarine permafrost, leading to coastal regions with shallower IBP.