ALBERT

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

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

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

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

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

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

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

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

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

Description: INGV and DPC

Description: Published

Description: 31-38

Description: 1IT. Reti di monitoraggio

Description: N/A or not JCR

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: 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: Dynamical processes during the formation phase of the Arctic stratospheric vortex in autumn (from September to December) can introduce considerable interannual variability in the amount of ozone that is incorporated into the vortex. Chemistry in autumn tends to remove part of this variability because ozone relaxes towards equilibrium. As a quantitative measure of how important dynamical variability during vortex formation is for the winter ozone abundances above the Arctic we analyze which fraction of an ozone anomaly induced during vortex formation persists until early winter (3 January). The work is based on the Lagrangian Chemistry Transport Model ATLAS. In a case study, model runs for the winter 1999–2000 are used to assess the fate of an ozone anomaly artificially introduced during the vortex formation phase on 16 September. The runs provide information about the persistence of the induced ozone anomaly as a function of time, potential temperature and latitude. The induced ozone anomaly survives longer inside the polar vortex compared to outside the vortex. Half of the initial perturbation survives until 3 January at 540 K inside the polar vortex, with a rapid fall off towards higher levels, mainly due to NOx induced chemistry. Above 750 K the signal falls to values below 0.5%. Hence, dynamically induced ozone variability from the early vortex formation phase cannot significantly contribute to early winter variability above 750 K. At lower levels increasingly larger fractions of the initial perturbation survive, reaching 90% at 450 K. In this vertical range dynamical processes during the vortex formation phase are crucial for the ozone abundance in early winter.

Description: We present first results from a coupled model setup, consisting of the state-of-the-art ice sheet model RIMBAY (Revised Ice Model Based on frAnk pattYn), and the community earth system model COSMOS. We show that special care has to be provided in order to ensure physical distributions of the forcings as well as numeric stability of the involved models. We demonstrate that a suitable statistical downscaling is crucial for ice sheet stability, especially for southern Greenland where surface temperatures are close to the melting point. The downscaling of net snow accumulation is based on an empirical relationship between surface slope and rainfall. The simulated ice sheet does not show dramatic loss of ice volume for pre-industrial conditions and is comparable with present-day ice orography. A sensitivity study with high CO2 level is used to demonstrate the effects of dynamic ice sheets onto climate compared to the standard setup with prescribed ice sheets.

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: 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: Sea urchins as broadcasting spawners, release their gametes into open water for fertilization, thus being particularly vulnerable to ocean acidification. In this study, we assessed the effects of different pH scenarios on fertilization success of Strongylocen- 5 trotus droebachiensis, collected at Spitsbergen, Arctic. We achieved acidification by bubbling CO2 into filtered seawater using partial pressures (pCO2) of 180, 380, 980, 1400 and 3000 μatm. Untreated filtered seawater was used as control. We recorded fertilization rates and diagnosed morphological aberrations after post-fertilization periods of 1 h and 3 h under different exposure conditions in experiments with and without 10 pre-incubation of the eggs prior to fertilization. In parallel, we conducted measurements of intracellular pH changes using BCECF/AM in unfertilized eggs exposed to a range of acidified seawater. We observed increasing rates of polyspermy in relation to higher seawater pCO2, which might be due to failures in the formation of the fertilization envelope. In addition, our experiments showed anomalies in fertilized eggs: incomplete 15 lifting-off of the fertilization envelope and blebs of the hyaline layer. Other drastic malformations consisted of constriction, extrusion, vacuolization or degeneration (observed as a gradient from the cortex to the central region of the cell) of the egg cytoplasm, and irregular cell divisions until 2- to 4-cell stages. The intracellular pH (pHi) decreased significantly from 1400 μatm on. All results indicate a decreasing fertilization success 20 at CO2 concentrations from 1400 μatm upwards. Exposure time to low pH might be a threatening factor for the cellular buffer capacity, viability, and development after fertilization.

Description: 〈jats:p〉Abstract. The spatial and temporal variability of a low-centred polygon on the eastern floodplain area of the lower Anabar River (72.070° N, 113.921° E, northern Yakutia, Siberia) has been investigated using a multi-method approach. The present-day vegetation in each square metre was analysed revealing a community of Larix shrubby Betula and Salix on the polygon rim, a dominance of Carex and Andromeda polifolia in the rim-to-pond transition zone, and a predominantly monospecific Scorpidium scorpioides coverage within the pond. The TOC content, TOC/TN ratio, grain-size, vascular plant macrofossils, moss remains, diatoms, and pollen were analysed for two vertical sections and a sediment core from a transect across the polygon. Radiocarbon dating indicates that the formation of the polygon started at least 1500 yr ago; the general positions of the pond and rim have not changed since that time. Two types of pond vegetation were identified, indicating two contrasting development stages of the polygon. The first was a well-established moss association dominated by submerged or floating Scorpidium scorpioides and/or Drepanocladus spp. and overgrown by epiphytic diatoms such as Tabellaria flocculosa and Eunotia taxa. This stage coincides temporally with a period in which the polygon was only drained by lateral subsurface water flow, as indicated by mixed grain sizes. A different moss association occurred during times of repeated river flooding (indicated by homogeneous medium-grained sand that probably accumulated during the annual spring snow melt), characterized by an abundance of Meesia triquetra and a dominance of benthic diatoms (e.g. Navicula vulpina), indicative of a relatively high pH and a high tolerance of disturbance. A comparison of the local polygon vegetation (inferred from moss and macrofossil spectra) with the regional vegetation (inferred from pollen spectra) indicated that the moss association with Scorpidium scorpioides became established during relatively favourable climatic conditions while the association dominated by Meesia triquetra occurred during periods of harsh climatic conditions. Our study revealed a strong riverine influence (in addition to climatic influences) on polygon development and the type of peat accumulated. 〈/jats:p〉

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.

Description: Source attribution and process analysis for atmospheric mercury in eastern China simulated by CMAQ-Hg Atmospheric Chemistry and Physics, 15, 8767-8779, 2015 Author(s): J. Zhu, T. Wang, J. Bieser, and V. Matthias The contribution from different emission sources and atmospheric processes to gaseous elemental mercury (GEM), gaseous oxidized mercury (GOM), particulate bound mercury (PBM) and mercury deposition in eastern China were quantified using the Community Multi-scale Air Quality (CMAQ-Hg) modeling system run with a nested domain. Natural sources (NAT) and six categories of anthropogenic mercury sources (ANTH) including cement production (CEM), domestic life (DOM), industrial boilers (IND), metal production (MET), coal-fired power plants (PP) and traffic (TRA) were considered for source apportionment. NAT were responsible for 36.6 % of annual averaged GEM concentration, which was regarded as the most important source for GEM in spite of obvious seasonal variation. Among ANTH, the influence of MET and PP on GEM were most evident especially in winter. ANTH dominated the variations of GOM and PBM concentrations with contributions of 86.7 and 79.1 %, respectively. Among ANTH, IND were the largest contributor for GOM (57.5 %) and PBM (34.4 %) so that most mercury deposition came from IND. The effect of mercury emitted from out of China was indicated by a 〉 30 % contribution to GEM concentration and wet deposition. The contributions from nine processes – consisting of emissions (EMIS), gas-phase chemical production/loss (CHEM), horizontal advection (HADV), vertical advection (ZADV), horizontal advection (HDIF), vertical diffusion (VDIF), dry deposition (DDEP), cloud processes (CLDS) and aerosol processes (AERO) – were calculated for process analysis with their comparison in urban and non-urban regions of the Yangtze River delta (YRD). EMIS and VDIF affected surface GEM and PBM concentrations most and tended to compensate each other all the time in both urban and non-urban areas. However, DDEP was the most important removal process for GOM with 7.3 and 2.9 ng m −3 reduced in the surface of urban and non-urban areas, respectively, in 1 day. The diurnal profile variation of processes revealed the transportation of GOM from urban area to non-urban areas and the importance of CHEM/AERO in higher altitudes which partly caused diffusion of GOM downwards to non-urban areas. Most of the anthropogenic mercury was transported and diffused away from urban areas by HADV and VDIF and increased mercury concentrations in non-urban areas by HADV. Natural emissions only influenced CHEM and AERO more significantly than anthropogenic. Local emissions in the YRD contributed 8.5 % more to GEM and ~ 30 % more to GOM and PBM in urban areas compared to non-urban areas.

Description: Seasonal variation of secondary organic aerosol tracers in Central Tibetan Plateau Atmospheric Chemistry and Physics, 15, 8781-8793, 2015 Author(s): R.-Q. Shen, X. Ding, Q.-F. He, Z.-Y. Cong, Q.-Q. Yu, and X.-M. Wang Secondary organic aerosol (SOA) affects the earth's radiation balance and global climate. High-elevation areas are sensitive to global climate change. However, at present, SOA origins and seasonal variations are understudied in remote high-elevation areas. In this study, particulate samples were collected from July 2012 to July 2013 at the remote Nam Co (NC) site, Central Tibetan Plateau and analyzed for SOA tracers from biogenic (isoprene, monoterpenes and β-caryophyllene) and anthropogenic (aromatics) precursors. Among these compounds, isoprene SOA (SOA I ) tracers represented the majority (26.6 ± 44.2 ng m −3 ), followed by monoterpene SOA (SOA M ) tracers (0.97 ± 0.57 ng m −3 ), aromatic SOA (SOA A ) tracer (2,3-dihydroxy-4-oxopentanoic acid, DHOPA, 0.25 ± 0.18 ng m −3 ) and β-caryophyllene SOA tracer (β-caryophyllenic acid, 0.09 ± 0.10 ng m −3 ). SOA I tracers exhibited high concentrations in the summer and low levels in the winter. The similar temperature dependence of SOA I tracers and isoprene emission suggested that the seasonal variation of SOA I tracers at the NC site was mainly influenced by the isoprene emission. The ratio of high-NO x to low-NO x products of SOA I (2-methylglyceric acid to 2-methyltetrols) was highest in the winter and lowest in the summer, due to the influence of temperature and relative humidity. The seasonal variation of SOA M tracers was impacted by monoterpenes emission and gas-particle partitioning. During the summer to the fall, temperature effect on partitioning was the dominant process influencing SOA M tracers' variation; while the temperature effect on emission was the dominant process influencing SOA M tracers' variation during the winter to the spring. SOA M tracer levels did not elevate with increased temperature in the summer, probably resulting from the counteraction of temperature effects on emission and partitioning. The concentrations of DHOPA were 1–2 orders of magnitude lower than those reported in the urban regions of the world. Due to the transport of air pollutants from the adjacent Bangladesh and northeastern India, DHOPA presented relatively higher levels in the summer. In the winter when air masses mainly came from northwestern India, mass fractions of DHOPA in total tracers increased, although its concentrations declined. The SOA-tracer method was applied to estimate secondary organic carbon (SOC) from these four precursors. The annual average of SOC was 0.22 ± 0.29 μgC m −3 , with the biogenic SOC (sum of isoprene, monoterpenes and β-caryophyllene) accounting for 75 %. In the summer, isoprene was the major precursor with its SOC contributions of 81 %. In the winter when the emission of biogenic precursors largely dropped, the contributions of aromatic SOC increased. Our study implies that anthropogenic pollutants emitted in the Indian subcontinent could be transported to the TP and have an impact on SOC over the remote NC.

Description: Influence of biomass burning on CCN number and hygroscopicity during summertime in the eastern Mediterranean Atmospheric Chemistry and Physics Discussions, 15, 21539-21582, 2015 Author(s): A. Bougiatioti, S. Bezantakos, I. Stavroulas, N. Kalivitis, P. Kokkalis, G. Biskos, N. Mihalopoulos, A. Papayannis, and A. Nenes This study investigates the CCN activity and hygroscopic properties of particles influenced by biomass burning in the eastern Mediterranean. Air masses sampled were subject to a range of atmospheric processing (several hours up to 3 days). Values of the hygroscopicity parameter, κ, were derived from cloud condensation nuclei (CCN) measurements and a Hygroscopic Tandem Differential Mobility Analyzer (HTDMA). An Aerosol Chemical Speciation Monitor (ACSM) was also used to determine the chemical composition and mass concentration of non-refractory components of the submicron aerosol fraction. During fire events, the increased organic content (and lower inorganic fraction) of the aerosol decreases the hygroscopicity parameter, κ, for all particle sizes. The reason, however, for this decrease was not the same for all size modes; smaller particle sizes appeared to be richer in less hygroscopic, less CCN-active components due to coagulation processes while larger particles become less hygroscopic during the biomass burning events due to condensation of less hygroscopic gaseous compounds. In addition, smaller particles exhibited considerable chemical dispersion (where hygroscopicity varied up to 100 % for particles of same size); larger particles, however, exhibited considerably less dispersion owing to the effects of aging and retained high levels of CCN activity. These conclusions are further supported by the observed mixing state determined by the HTDMA measurements. ACSM measurements indicate that the bulk composition reflects the hygroscopicity and chemical nature of the largest particles and a large fraction of the CCN concentrations sampled. Based on Positive Matrix Factorization (PMF) analysis of the organic ACSM spectra, CCN concentrations follow a similar trend with the BBOA component, with enhancements of CCN in biomass burning plumes ranging between 65 and 150 %, for supersaturations ranging between 0.2 and 0.7 %. Using multilinear regression, we determine the hygroscopicity of the prime organic aerosol components (BBOA, OOA-BB and OOA); it is found that the total organic hygroscopicity is very close to the inferred hygroscopicity of the oxygenated organic aerosol components. Finally, the transformation of freshly-emitted biomass burning (BBOA) to more oxidized organic aerosol (OOA-BB) can result in a two-fold increase of the inferred organic hygroscopicity. Almost 10 % of the total aerosol hygroscopicity is related to the two biomass burning components (BBOA and OOA-BB), which in turn contribute almost 35 % to the fine-particle organic water of the aerosol. This is important as organic water can contribute to the atmospheric chemistry and the direct radiative forcing.

Description: Measurement-based direct radiative effect by brown carbon over Indo-Gangetic Plain Atmospheric Chemistry and Physics Discussions, 15, 21583-21606, 2015 Author(s): A. Arola, G. L. Schuster, M. R. A. Pitkänen, O. Dubovik, H. Kokkola, A. V. Lindfors, T. Mielonen, T. Raatikainen, S. Romakkaniemi, S. N. Tripathi, and H. Lihavainen The importance of light absorbing organic aerosols, often called brown carbon (BrC), has become evident in recent years. However, there are relatively few measurement-based estimates for the direct radiative effect of BrC so far. In those earlier studies, the AErosol RObotic NETwork (AERONET) measured Aerosol Absorption Optical Depth (AAOD) and Absorption Angstrom Exponent (AAE) have been exploited. However, these two pieces of information are clearly not sufficient to separate properly carbonaceous aerosols from dust, while imaginary indices of refraction would contain more and better justified information for this purpose. This is first time that the direct radiative effect (DRE) of BrC is estimated by exploiting the AERONET-retrieved imaginary indices. We estimated it for four sites in Indo-Gangetic Plain (IGP), Karachi, Lahore, Kanpur and Gandhi College. We found a distinct seasonality, which was generally similar among all the sites, but with slightly different strengths. The monthly warming effect up to 0.5 W m -2 takes place during spring season. On the other hand, BrC results in overall cooling effect in the winter season, which can reach levels close to −1W m -2 . We then estimated similarly also DRE of black carbon and total aerosol, in order to assess the relative significance of BrC radiative effect in the radiative effects of other components. Even though BrC impact seems minor in this context, we demonstrated that it is not insignificant and moreover that it is crucial to perform spectrally resolved radiative transfer calculations to obtain good estimates for DRE of BrC.

Description: Contrail life cycle and properties from 1 year of MSG/SEVIRI rapid-scan images Atmospheric Chemistry and Physics, 15, 8739-8749, 2015 Author(s): M. Vázquez-Navarro, H. Mannstein, and S. Kox The automatic contrail tracking algorithm (ACTA) – developed to automatically follow contrails as they age, drift and spread – enables the study of a large number of contrails and the evolution of contrail properties with time. In this paper we present a year's worth of tracked contrails, from August 2008 to July 2009 in order to derive statistically significant mean values. The tracking is performed using the 5 min rapid-scan mode of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on board the Meteosat Second Generation (MSG) satellites. The detection is based on the high spatial resolution of the images provided by the Moderate Resolution Imaging Spectroradiometer on board the Terra satellite (Terra/MODIS), where a contrail detection algorithm (CDA) is applied. The results show the satellite-derived average lifetimes of contrails and contrail-cirrus along with the probability density function (PDF) of other geometric characteristics such as mean coverage, distribution and width. In combination with specifically developed algorithms (RRUMS; Rapid Retrieval of Upwelling irradiance from MSG/SEVIRI and COCS (Cirrus Optical properties derived from CALIOP and SEVIRI), explained below) it is possible to derive the radiative forcing (RF), energy forcing (EF), optical thickness (τ) and altitude of the tracked contrails. Mean values here retrieved are duration, 1 h; length, 130 km; width, 8 km; altitude, 11.7 km; optical thickness, 0.34. Radiative forcing and energy forcing are shown for land/water backgrounds in day/night situations.

Description: Accounting for the effects of sastrugi in the CERES clear-sky Antarctic shortwave angular distribution models Atmospheric Measurement Techniques, 8, 3163-3175, 2015 Author(s): J. Corbett and W. Su The Cloud and the Earth's Radiant Energy System (CERES) instruments on NASA's Terra, Aqua and Soumi NPP satellites are used to provide a long-term measurement of Earth's energy budget. To accomplish this, the radiances measured by the instruments must be inverted to fluxes by the use of a scene-type-dependent angular distribution model (ADM). For permanent snow scenes over Antarctica, shortwave (SW) ADMs are created by compositing radiance measurements over the full viewing zenith and azimuth range. However, the presence of small-scale wind blown roughness features called sastrugi cause the BRDF (bidirectional reflectance distribution function) of the snow to vary significantly based upon the solar azimuth angle and location. This can result in monthly regional biases between −12 and 7.5 Wm −2 in the inverted TOA (top-of-atmosphere) SW flux. The bias is assessed by comparing the CERES shortwave fluxes derived from nadir observations with those from all viewing zenith angles, as the sastrugi affect fluxes inverted from the oblique viewing angles more than for the nadir viewing angles. In this paper we further describe the clear-sky Antarctic ADMs from Su et al. (2015). These ADMs account for the sastrugi effect by using measurements from the Multi-Angle Imaging Spectro-Radiometer (MISR) instrument to derive statistical relationships between radiance from different viewing angles. We show here that these ADMs reduce the bias and artifacts in the CERES SW flux caused by sastrugi, both locally and Antarctic-wide. The regional monthly biases from sastrugi are reduced to between −5 and 7 Wm −2 , and the monthly-mean biases over Antarctica are reduced by up to 0.64 Wm −2 , a decrease of 74 %. These improved ADMs are used as part of the Edition 4 CERES SSF (Single Scanner Footprint) data.

Description: Space-borne observation of methane from atmospheric infrared sounder version 6: validation and implications for data analysis Atmospheric Measurement Techniques Discussions, 8, 8563-8597, 2015 Author(s): X. Xiong, F. Weng, Q. Liu, and E. Olsen Atmospheric Methane (CH 4 ) is generated as a standard product in recent version of the hyperspectral Atmospheric Infrared Sounder (AIRS-V6) aboard NASA's Aqua satellite at the NASA Goddard Earth Sciences Data and Information Services Center (NASA/GES/DISC). Significant improvements in AIRS-V6 was expected but without a thorough validation. This paper first introduced the improvements of CH 4 retrieval in AIRS-V6 and some characterizations, then presented the results of validation using ~ 1000 aircraft profiles from several campaigns spread over a couple of years and in different regions. It was found the mean biases of AIRS CH 4 at layers 343–441 and 441–575 hPa are −0.76 and −0.05 % and the RMS errors are 1.56 and 1.16 %, respectively. Further analysis demonstrates that the errors in the spring and in the high northern latitudes are larger than in other seasons or regions. The error is correlated with Degree of Freedoms (DOFs), particularly in the tropics or in the summer, and cloud amount, suggesting that the "observed" spatiotemporal variation of CH 4 by AIRS is imbedded with some artificial impact from the retrieval sensitivity in addition to its variation in reality, so the variation of information content in the retrievals needs to be taken into account in data analysis of the retrieval products. Some additional filtering (i.e. rejection of profiles with obvious oscillation as well as those deviating greatly from the norm) for quality control is recommended for the users to better utilize AIRS-V6 CH 4 , and their implementation in the future versions of the AIRS retrieval algorithm is under consideration.

Description: Impact of aerosols on the OMI tropospheric NO 2 retrievals over industrialized regions: how accurate is the aerosol correction of cloud-free scenes via a simple cloud model? Atmospheric Measurement Techniques Discussions, 8, 8385-8437, 2015 Author(s): J. Chimot, T. Vlemmix, J. P. Veefkind, J. F. de Haan, and P. F. Levelt The Ozone Monitoring Instrument (OMI) instrument has provided daily global measurements of tropospheric NO 2 for more than a decade. Numerous studies have drawn attention to the complexities related to measurements of tropospheric NO 2 in the presence of aerosols. Fine particles affect the OMI spectral measurements and the length of the average light path followed by the photons. However, they are not explicitly taken into account in the current OMI tropospheric NO 2 retrieval chain. Instead, the operational OMI O 2 -O 2 cloud retrieval algorithm is applied both to cloudy scenes and to cloud free scenes with aerosols present. This paper describes in detail the complex interplay between the spectral effects of aerosols, the OMI O 2 -O 2 cloud retrieval algorithm and the impact on the accuracy of the tropospheric NO 2 retrievals through the computed Air Mass Factor (AMF) over cloud-free scenes. Collocated OMI NO 2 and MODIS Aqua aerosol products are analysed over East China, in industrialized area. In addition, aerosol effects on the tropospheric NO 2 AMF and the retrieval of OMI cloud parameters are simulated. Both the observation-based and the simulation-based approach demonstrate that the retrieved cloud fraction linearly increases with increasing Aerosol Optical Thickness (AOT), but the magnitude of this increase depends on the aerosol properties and surface albedo. This increase is induced by the additional scattering effects of aerosols which enhance the scene brightness. The decreasing effective cloud pressure with increasing AOT represents primarily the absorbing effects of aerosols. The study cases show that the actual aerosol correction based on the implemented OMI cloud model results in biases between −20 and −40 % for the DOMINO tropospheric NO 2 product in cases of high aerosol pollution (AOT ≥ 0.6) and elevated particles. On the contrary, when aerosols are relatively close to the surface or mixed with NO 2 , aerosol correction based on the cloud model results in overestimation of the DOMINO tropospheric NO 2 product, between 10 and 20 %. These numbers are in line with comparison studies between ground-based and OMI tropospheric NO 2 measurements under conditions with high aerosol pollution and elevated particles. This highlights the need to implement an improved aerosol correction in the computation of tropospheric NO 2 AMFs.

Description: Sensitivity of thermal infrared sounders to the chemical and micro-physical properties of UTLS secondary sulphate aerosols Atmospheric Measurement Techniques Discussions, 8, 8439-8481, 2015 Author(s): P. Sellitto and B. Legras Monitoring upper tropospheric-lower stratospheric (UTLS) secondary sulphate aerosols and their chemical and micro-physical properties from satellite nadir observations is crucial to better understand their formation and evolution processes and then to estimate their impact to the UTLS chemistry, and on regional and global radiative balance. Here we present a study aimed at the evaluation of the sensitivity of thermal infrared (TIR) satellite nadir observations to the chemical composition and the size distribution of idealized UTLS sulphate aerosol layers. The extinction properties of sulphuric acid/water droplets, for different sulphuric acid mixing ratios and temperatures, are systematically analysed. The extinction coefficients are derived by means of a Mie code, using refractive indexes taken from the GEISA (Gestion et Étude des Informations Spectroscopiques Atmosphériques: Management and Study of Spectroscopic Information) spectroscopic database and log-normal size distributions with different effective radii and number concentrations. IASI (Infrared Atmospheric Sounding Interferometer) pseudo-observations are generated using forward radiative transfer calculations performed with the 4A (Automatized Atmospheric Absorption Atlas) radiative transfer model, to estimate the impact of the extinction of idealized aerosol layers, at typical UTLS conditions, on the brightness temperature spectra observed by this satellite instrument. We found a marked and typical spectral signature of these aerosol layers between 700 and 1200 cm −1 , due to the absorption bands of the sulphate and bi-sulphate ions and the undissociated sulphuric acid, with the main absorption peaks at 1170 and 905 cm −1 . The dependence of the aerosol spectral signature to the sulphuric acid mixing ratio, and effective number concentration and radius, as well as the role of interferring parameters like the ozone, sulphur dioxide, carbon dioxide and ash absorption, and temperature and water vapour profile uncertainties, are analyzed and critically discussed.

Description: Global cloud top height retrieval using SCIAMACHY limb spectra: model studies and first results Atmospheric Measurement Techniques Discussions, 8, 8295-8352, 2015 Author(s): K.-U. Eichmann, L. Lelli, C. von Savigny, H. Sembhi, and J. P. Burrows Cloud top heights (CTH) were retrieved for the period 1 January 2003 to 7 April 2012 using height-resolved limb spectra measured with the Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) on board ENVISAT (ENVIronmental SATellite). In this study, we tested the sensitivity of the colour index method used in the retrieval code SCODA (SCIAMACHY Cloud Detection Algorithm) and the accuracy of the retrieved CTHs in comparison to other methods. Sensitivity studies using the radiative transfer model SCIATRAN showed that the method is capable of generally detecting cloud tops down to about 5 km and very thin cirrus clouds even up to the tropopause. Volcanic particles can also be detected that occasionally reach the lower stratosphere. Low clouds at 2–3 km can only be retrieved under very clean atmospheric conditions, as light scattering of aerosols interferes with the cloud retrieval. Upper tropospheric ice clouds are detectable for cloud optical depths down to about τ N = 0.005, which is in the subvisual range. The detection sensitivity decreases towards the surface. An optical thickness of roughly 0.1 was the lower detection limit for water cloud top heights at 5 km. This value is much lower than thresholds reported for the passive cloud detection in nadir viewing direction. Comparisons with SCIAMACHY nadir cloud top heights, calculated with the Semi-Analytical CloUd Retrieval Algorithm (SACURA), showed a good agreement in the global cloud field distribution. But only opaque clouds (τ N 〉 5) are detectable with the nadir passive retrieval technique in the UV-visible and infrared wavelength range. So due to the frequent occurrence of thin and sub-visual cirrus clouds in the tropics, large cloud top height deviations were detected between both viewing geometries. Also the land/sea contrast seen in nadir retrievals was not detected in limb mode. Co-located cloud top height measurements of the limb viewing Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on ENVISAT for the period from January 2008 to March 2012 were compared, showing good agreement to within 1 km, which is smaller than the vertical field of view of both instruments. Lower stratospheric aerosols from volcanic eruptions occasionally interfered with the cloud retrieval and inhibited detection of tropospheric clouds. Examples of the impact of these events are shown for the volcanoes Kasatochi in August 2008, Sarychev Peak in June 2009, and Nabro in June 2010. Long-lasting aerosol layers were detected after these events in the Northern Hemisphere down to the tropics. Particle top heights up to about 22 km were retrieved in 2009, when the enhanced lower stratospheric aerosol layer persisted for about 7 months. Up to about 82 % of the Northern hemispheric lower stratosphere between 30° and 70° was covered by scattering particles in August 2009 and nearly half in October 2008.

Description: Operation of the Airmodus A11 nano Condensation Nucleus Counter at various inlet pressures, various operation temperatures and design of a new inlet system Atmospheric Measurement Techniques Discussions, 8, 8483-8508, 2015 Author(s): J. Kangasluoma, A. Franchin, J. Duplissy, L. Ahonen, F. Korhonen, M. Attoui, J. Mikkilä, K. Lehtipalo, J. Vanhanen, M. Kulmala, and T. Petäjä Measuring sub-3 nm particles outside of controlled laboratory conditions is a challenging task, as many of the instruments are operated at their limits and are subjected to changing ambient conditions. In this study, we advance the current understanding on the operation of Airmodus A11 nano Condensation Nucleus Counter (nCNC), which consists of a A10 Particle Size Magnifier (PSM) and A20 condensation particle counter (CPC). We explore the effect of the inlet line pressure on the measured particle concentration. We identify two different regions inside the instrument where supersaturation of working fluid can take place. We show the possibility of varying the cut-off of the instrument from 1 to 6 nm, a wider size range than the one usually covered by the PSM. We also present a new inlet system, which allows automated measurements of the background, minimizes the diffusion losses in the sampling line and is equipped with an electrostatic filter to remove ions. Finally, our view of the guidelines for optimal use of the Airmodus nCNC are provided.

Description: Looking beyond stratification: a model-based analysis of the biological drivers of oxygen depletion in the North Sea Biogeosciences Discussions, 12, 12543-12610, 2015 Author(s): F. Große, N. Greenwood, M. Kreus, H. J. Lenhart, D. Machoczek, J. Pätsch, L. A. Salt, and H. Thomas The problem of low oxygen conditions, often referred to as hypoxia, occurs regularly in the North Sea, a temperate European shelf sea. Stratification represents a major process regulating the seasonal dynamics of bottom oxygen. However, lowest oxygen conditions in the North Sea do not occur in the regions of strongest stratification. This suggests that stratification is an important prerequisite for hypoxia, but that the complex interaction between hydrodynamics and the biological processes drives its development. In this study we use the ecosystem model HAMSOM-ECOHAM5 to provide a general characteristic of the different North Sea oxygen regimes, and to quantify the impact of the different physical and biological factors driving the oxygen dynamics below the thermocline and in the bottom layer. We show that the North Sea can be subdivided into three different regimes in terms of oxygen dynamics: (1) a highly productive, non-stratified coastal regime, (2) a productive, seasonally stratified regime with a small sub-thermocline volume, and (3) a productive, seasonally stratified regime with a large sub-thermocline volume, with regime 2 being highly susceptible to hypoxic conditions. Our analysis of the different processes driving the oxygen development reveals that inter-annual variations in the oxygen conditions are caused by variations in primary production, while spatial differences can be attributed to differences in stratification and water depth. In addition, we show that benthic bacteria represent the main oxygen consumers in the bottom layer, consistently accounting for more than 50 % of the overall consumption. By providing these valuable insights, we show that ecosystem models can be a useful tool for the interpretation of observations and the estimation of the impact of anthropogenic drivers on the North Sea oxygen conditions.

Description: Impact of the oceanic geothermal heat flux on a glacial ocean state Climate of the Past Discussions, 11, 3597-3624, 2015 Author(s): M. Ballarotta, F. Roquet, S. Falahat, Q. Zhang, and G. Madec The oceanic geothermal heating (OGH) has a significant impact on the present-day ocean state, but its role during glacial periods, when the ocean circulation and stratification were different from those of today, remains poorly known. In the present study, we analyzed the response of the glacial ocean to OGH, by comparing ocean simulations of the Last Glacial Maximum (LGM, ∼ 21 ka ago) including or not geothermal heating. We found that applying the OGH warmed the Antarctic Bottom Waters (AABW) by ∼ 0.4 °C and increased the abyssal circulation by 15 to 30 % north of 30° S in the deep Pacific and Atlantic basins. The geothermally heated deep waters were then advected toward the Southern Ocean where they upwelled to the surface due to the Ekman transport. The extra heat transport towards Antarctica acted to reduce the amount of sea ice contributing to the freshening of the whole AABW overturning cell. The global amount of salt being conserved, this bottom freshening induced a salinification of the North Atlantic and North Pacific surface and intermediate waters, contributing to the deepening of the North Atlantic Deep Water. This indirect mechanism is responsible for the largest observed warming, found in the North Atlantic deep western boundary current between 2000 and 3000 m (up to 2 °C). The characteristic time scale of the ocean response to the OGH corresponds to an advective time scale (associated with the overturning of the AABW cell) rather than a diffusive time scale. The OGH might facilitate the transition from a glacial to an inter-glacial state but its effect on the deep stratification seems insufficient to drive alone an abrupt climate change.

Description: Woodland Survey of Great Britain 1971–2001 Earth System Science Data, 7, 203-214, 2015 Author(s): C. M. Wood, S. M. Smart, and R. G. H. Bunce The Woodland Survey of Great Britain is a unique data set, consisting of a detailed range of ecological measurements at a national scale, covering a time span of 30 years. A set of 103 woods spread across Britain were first surveyed in 1971, which were again surveyed in 2000–2003 (for convenience referred to subsequently as the "2001 survey"). Standardised methods of describing the trees, shrubs, ground flora, soils and general habitats present were used for both sets of surveys. The sample of 1648 plots spread through 103 woodland sites located across Britain makes it probably the most extensive quantitative ecological woodland survey undertaken in Britain; it is also notable for the range of sites that have been revisited after such a long interval. The data set provides a unique opportunity to explore the effects of a range of potential drivers of woodland change that operated between 1971 and 2001. The data set is available in four discrete parts, which have been assigned the following DOIs: doi:10.5285/4d93f9ac-68e3-49cf-8a41-4d02a7ead81a (Kirby et al., 2013b), doi:10.5285/d6409d40-58fe-4fa7-b7c8-71a105b965b4 (Kirby et al., 2013d), doi:10.5285/fb1e474d-456b-42a9-9a10-a02c35af10d2 (Kirby et al., 2013c), doi:10.5285/2d023ce9-6dbe-4b4f-a0cd-34768e1455ae (Kirby et al., 2013a).

Description: Estimating greenhouse gas emissions from travel – a GIS-based study Geographica Helvetica, 70, 185-192, 2015 Author(s): S. Kuonen Conferences, meetings and congresses are an important part of today's economic and scientific world. But the environmental impact, especially from greenhouse gas emissions associated with travel, can be extensive. Anthropogenic greenhouse gas (GHG) emissions account for the warming of the atmosphere and oceans. This study draws on the need to quantify and reduce greenhouse gas emissions associated with travel activities and aims to give suggestions for organizers and participants on possible ways to reduce greenhouse gas emissions, demonstrated on the example of the European Geography Association (EGEA) Annual Congress 2013 in Wasilkow, Poland. The lack of a comprehensive methodology for the estimation of greenhouse gas emissions from travel led to an outline of a methodology that uses geographic information systems (GIS) to calculate travel distances. The calculation of travel distances in GIS is adapted from actual transportation infrastructure, derived from the open-source platform OpenStreetMap. The methodology also aims to assess the possibilities to reduce GHG emissions by choosing different means of transportation and a more central conference location. The results of the participants of the EGEA congress, who shared their travel data for this study, show that the total travel distance adds up to 238 000 km, with average travel distance of 2429 km per participant. The travel activities of the participants in the study result in total GHG emissions of 39 300 kg CO 2 -eq including both outward and return trip. On average a participant caused GHG emissions of 401 kg CO 2 -eq. In addition, the analysis of the travel data showed differences in travel behaviour depending on the distance between conference site and point of origin. The findings on travel behaviour have then been used to give an estimation of total greenhouse gas emissions from travel for all participants of the conference, which result in a total amount of 79 711 kg CO 2 -eq. The potential for reducing greenhouse gas emissions by substituting short flights with train rides and car rides with bus and train rides is limited. Only 6 % of greenhouse gas emissions could be saved by applying these measures. Further considerable savings could only be made by substituting longer flights (32.6 %) or choosing a more central conference location (26.3 %).

Description: Singular vector based targeted observations of chemical constituents: description and first application of the EURAD-IM-SVA Geoscientific Model Development Discussions, 8, 6267-6307, 2015 Author(s): N. Goris and H. Elbern Measurements of the large dimensional chemical state of the atmosphere provide only sparse snapshots of the state of the system due to their typically insufficient temporal and spatial density. In order to optimize the measurement configurations despite those limitations, the present work describes the identification of sensitive states of the chemical system as optimal target areas for adaptive observations. For this purpose, the technique of singular vector analysis (SVA), which has been proved effective for targeted observations in numerical weather predication, is implemented into the chemical transport model EURAD-IM (EURopean Air pollution and Dispersion – Inverse Model) yielding the EURAD-IM-SVA. Besides initial values, emissions are investigated as critical simulation controlling targeting variables. For both variants, singular vectors are applied to determine the optimal placement for observations and moreover to quantify which chemical compounds have to be observed with preference. Based on measurements of the airship based ZEPTER-2 campaign, the EURAD-IM-SVA has been evaluated by conducting a comprehensive set of model runs involving different initial states and simulation lengths. Since the considered cases are restricted in terms of considered chemical compounds and selected areas, they allow for a retracing of the results and a confirmation of their correctness. Our analysis shows that the optimal placement for observations of chemical species is not entirely determined by mere transport and mixing processes. Rather, a combination of initial chemical concentrations, chemical conversions, and meteorological processes determine the influence of chemical compounds and regions. We furthermore demonstrate that the optimal placement of observations of emission strengths is highly dependent on the location of emission sources and that the benefit of including emissions as target variables outperforms the value of initial value optimisation with growing simulation length. The obtained results confirm the benefit of considering both initial values and emission strengths as target variables and of applying the EURAD-IM-SVA for measurement decision guidance with respect to chemical compounds.

Description: A global scale mechanistic model of the photosynthetic capacity Geoscientific Model Development Discussions, 8, 6217-6266, 2015 Author(s): A. A. Ali, C. Xu, A. Rogers, R. A. Fisher, S. D. Wullschleger, N. G. McDowell, E. C. Massoud, J. A. Vrugt, J. D. Muss, J. B. Fisher, P. B. Reich, and C. J. Wilson Although plant photosynthetic capacity as determined by the maximum carboxylation rate (i.e., V c, max25 ) and the maximum electron transport rate (i.e., J max25 ) at a reference temperature (generally 25 °C) is known to vary substantially in space and time in response to environmental conditions, it is typically parameterized in Earth system models (ESMs) with tabulated values associated to plant functional types. In this study, we developed a mechanistic model of leaf utilization of nitrogen for assimilation (LUNA V1.0) to predict the photosynthetic capacity at the global scale under different environmental conditions, based on the optimization of nitrogen allocated among light capture, electron transport, carboxylation, and respiration. The LUNA model was able to reasonably well capture the observed patterns of photosynthetic capacity in view that it explained approximately 55 % of the variation in observed V c, max25 and 65 % of the variation in observed J max25 across the globe. Our model simulations under current and future climate conditions indicated that V c, max25 could be most affected in high-latitude regions under a warming climate and that ESMs using a fixed V c, max25 or J max25 by plant functional types were likely to substantially overestimate future global photosynthesis.

Description: Uncertainties in calculating precipitation climatology in East Asia Hydrology and Earth System Sciences Discussions, 12, 7765-7783, 2015 Author(s): J. Kim and S. K. Park This study examines the uncertainty in calculating the fundamental climatological characteristics of precipitation in the East Asia region from multiple fine-resolution gridded analysis datasets based on in-situ rain gauge observations. Five observation-based gridded precipitation datasets are used to derive the long-term means, standard deviations in lieu of interannual variability and linear trends over the 28-year period from 1980 to 2007. Both the annual and summer (June–July–August) mean precipitation is examined. The agreement amongst these precipitation datasets are examined using multiple metrics including the signal-to-noise ratio (SNR) defined as the ratio between long-term means and the corresponding standard deviations, and Taylor diagrams which allows examinations of the pattern correlation, the standard deviation, and the centered root mean square error. It is found that the five gauge-based precipitation analysis datasets agree well in the long-term mean and interannual variability in most of the East Asia region including eastern China, Manchuria, South Korea, and Japan, which are densely populated and have fairly high density observation networks. The regions of large inter-dataset variations include Tibetan Plateau, Mongolia, northern Indo-China, and North Korea. The regions of large uncertainties are typically lightly populated and are characterized by severe terrain and/or extreme high elevations. Unlike the long-term mean and interannual variability, agreements between datasets in the linear trend is weak, both for the annual and summer mean values. In most of the East Asia region, the SNR for the linear trend is below 0.5, i.e., the inter-dataset variability exceeds the multi-data ensemble mean. The uncertainty in the spatial distribution of long-term means among these datasets occurs both in the spatial pattern and variability, but the uncertainty for the interannual variability and time trend is much larger in the variability than in the pattern correlation. Thus, care must be taken in using long-term trends calculated from gridded precipitation analysis data for climate studies over the East Asia region.

Description: Qualitative soil moisture assessment in semi-arid Africa – the role of experience and training on inter-rater reliability Hydrology and Earth System Sciences, 19, 3505-3516, 2015 Author(s): M. Rinderer, H. C. Komakech, D. Müller, G. L. B. Wiesenberg, and J. Seibert Soil and water management is particularly relevant in semi-arid regions to enhance agricultural productivity. During periods of water scarcity, soil moisture differences are important indicators of the soil water deficit and are traditionally used for allocating water resources among farmers of a village community. Here we present a simple, inexpensive soil wetness classification scheme based on qualitative indicators which one can see or touch on the soil surface. It incorporates the local farmers' knowledge on the best soil moisture conditions for seeding and brick making in the semi-arid environment of the study site near Arusha, Tanzania. The scheme was tested twice in 2014 with farmers, students and experts (April: 40 persons, June: 25 persons) for inter-rater reliability, bias of individuals and functional relation between qualitative and quantitative soil moisture values. During the test in April farmers assigned the same wetness class in 46 % of all cases, while students and experts agreed on about 60 % of all cases. Students who had been trained in how to apply the method gained higher inter-rater reliability than their colleagues with only a basic introduction. When repeating the test in June, participants were given improved instructions, organized in small subgroups, which resulted in a higher inter-rater reliability among farmers. In 66 % of all classifications, farmers assigned the same wetness class and the spread of class assignments was smaller. This study demonstrates that a wetness classification scheme based on qualitative indicators is a robust tool and can be applied successfully regardless of experience in crop growing and education level when an in-depth introduction and training is provided. The use of a simple and clear layout of the assessment form is important for reliable wetness class assignments.