ALBERT

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

Description: Published

Description: 1681–1698

Description: 2T. Deformazione crostale attiva

Description: JCR Journal

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

Description: Copernicus

Description: Published

Description: On line

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: 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: 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: 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: 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: In order to investigate the impact of spatial resolution on the discrepancy between simulated δ18O and observed δ18O in Greenland ice cores, regional climate simulations are performed with the isotope-enabled regional climate model (RCM) COSMO_iso. For this purpose, isotope-enabled general circulation model (GCM) simulations with the ECHAM5-wiso general circulation model (GCM) under present-day conditions and the MPI-ESM-wiso GCM under mid-Holocene conditions are dynamically downscaled with COSMO_iso for the Arctic region. The capability of COSMO_iso to reproduce observed isotopic ratios in Greenland ice cores for these two periods is investigated by comparing the simulation results to measured δ18O ratios from snow pit samples, Global Network of Isotopes in Precipitation (GNIP) stations and ice cores. To our knowledge, this is the first time that a mid-Holocene isotope-enabled RCM simulation is performed for the Arctic region. Under present-day conditions, a dynamical downscaling of ECHAM5-wiso (1.1◦ × 1.1◦) with COSMO_iso to a spatial resolution of 50km improves the agreement with the measured δ18O ratios for 14 of 19 observational data sets. A further increase in the spatial resolution to 7km does not yield substantial improvements except for the coastal areas with its complex terrain. For the mid-Holocene, a fully coupled MPI-ESM-wiso time slice simulation is downscaled with COSMO_iso to a spatial resolution of 50km. In the mid-Holocene, MPI-ESM-wiso already agrees well with observations in Greenland and a downscaling with COSMO_iso does not further improve the model–data agreement. Despite this lack of improvement in model biases, the study shows that in both periods, observed δ18O values at measurement sites constitute isotope ratios which are mainly within the subgrid-scale variability of the global ECHAM5-wiso and MPI-ESM-wiso simulation results. The correct δ18O ratios are consequently not resolved in the GCM simulation results and need to be extracted by a refinement with an RCM. In this context, the RCM simulations provide a spatial δ18O distribution by which the effects of local uncertainties can be taken into account in the comparison between point measurements and model outputs. Thus, an isotope-enabled GCM–RCM model chain with realistically implemented fractionating processes constitutes a useful supplement to reconstruct regional paleo-climate conditions during the mid-Holocene in Greenland. Such model chains might also be applied to reveal the full potential of GCMs in other regions and climate periods, in which large deviations relative to observed isotope ratios are simulated.

Description: Proxy climate records are an invaluable source of information about the earth’s climate prior to the instrumental record. The temporal- and spatial-coverage of records continues to increase, however, these records of past climate are associated with significant uncertainties due to non-climate processes that influence the recorded and measured proxy values. Generally, these uncertainties are timescale-dependent and correlated in time. Accounting for structure in the errors is essential to providing realistic error estimates for smoothed or stacked records, detection of anomalies and identifying trends, but this structure is seldom accounted for. In the first of these companion articles we outlined a theoretical framework for handling proxy uncertainties by deriving the power spectrum of proxy error components from which it is possible to obtain timescale-dependent error estimates. Here in part II, we demonstrate the practical application of this theoretical framework using the example of marine sediment cores. We consider how to obtain estimates for the required parameters and give examples of the application of this approach for typical marine sediment proxy records. Our new approach of estimating and providing timescale-dependent proxy errors overcomes the limitations of simplistic single value error estimates. We aim to provide the conceptual basis for a more quantitative use of paleo-records for applications such as model-data comparison, regional and global synthesis of past climate states and data assimilation.

Description: The mid-Pliocene warm period (mPWP; ∼3.2 million years ago) is seen as the most recent time period characterized by a warm climate state, with similar to modern geography and ∼400 ppmv atmospheric CO2 concentration, and is therefore often considered an interesting analogue for near-future climate projections. Paleoenvironmental reconstructions indicate higher surface temperatures, decreasing tropical deserts, and a more humid climate in West Africa characterized by a strengthened West African Monsoon (WAM). Using model results from the second phase of the Pliocene Modelling Intercomparison Project (PlioMIP2) ensemble, we analyse changes of the WAM rainfall during the mPWP by comparing them with the control simulations for the pre-industrial period. The ensemble shows a robust increase in the summer rainfall over West Africa and the Sahara region, with an average increase of 2.5 mm/d, contrasted by a rainfall decrease over the equatorial Atlantic. An anomalous warming of the Sahara and deepening of the Saharan Heat Low, seen in 〉90 % of the models, leads to a strengthening of the WAM and an increased monsoonal flow into the continent. A similar warming of the Sahara is seen in future projections using both phase 3 and 5 of the Coupled Model Intercomparison Project (CMIP3 and CMIP5). Though previous studies of future projections indicate a west–east drying–wetting contrast over the Sahel, PlioMIP2 simulations indicate a uniform rainfall increase in that region in warm climates characterized by increasing greenhouse gas forcing. We note that this effect will further depend on the long-term response of the vegetation to the CO2 forcing.

Description: In this study, the first fully continuous monitoring of water vapour isotopic composition at Neumayer Station III, Antarctica, during the 2-year period from February 2017 to January 2019 is presented. Seasonal and synoptic-scale variations in both stable water isotopes H182O and HDO are reported, and their links to variations in key meteorological variables are analysed. In addition, the diurnal cycle of isotope variations during the summer months (December and January 2017/18 and 2018/19) has been examined. Changes in local temperature and specific humidity are the main drivers for the variability in δ18O and δD in vapour at Neumayer Station III, on both seasonal and shorter timescales. In contrast to the measured δ18O and δD variations, no seasonal cycle in the Deuterium excess signal (d) in vapour is detected. However, a rather high uncertainty in measured d values especially in austral winter limits the confidence of this finding. Overall, the d signal shows a stronger inverse correlation with specific humidity than with temperature, and this inverse correlation between d and specific humidity is stronger for the cloudy-sky conditions than for clear-sky conditions during summertime. Back-trajectory simulations performed with the FLEXPART model show that seasonal and synoptic variations in δ18O and δD in vapour coincide with changes in the main sources of water vapour transported to Neumayer Station III. In general, moisture transport pathways from the east lead to higher temperatures and more enriched δ18O values in vapour, while weather situations with southerly winds lead to lower temperatures and more depleted δ18O values. However, on several occasions, δ18O variations linked to wind direction changes were observed, which were not accompanied by a corresponding temperature change. Comparing isotopic compositions of water vapour at Neumayer Station III and snow samples taken in the vicinity of the station reveals almost identical slopes, both for the δ18O–δD relation and for the temperature–δ18O relation.

Description: Due to its dryness, the subtropical free troposphere plays a critical role in the radiative balance of the Earth's climate system. But the complex interactions of the dynamical and physical processes controlling the variability in the moisture budget of this sensitive region of the subtropical atmosphere are still not fully understood. Stable water isotopes can provide important information about several of the latter processes, namely subsidence drying, turbulent mixing, and dry and moist convective moistening. In this study, we use high-resolution simulations of the isotope-enabled version of the regional weather and climate prediction model of the Consortium for Small-Scale Modelling (COSMOiso) to investigate predominant moisture transport pathways in the Canary Islands region in the eastern subtropical North Atlantic. Comparison of the simulated isotope signals with multi-platform isotope observations (aircraft, ground- and space-based remote sensing) from a field campaign in summer 2013 shows that COSMOiso can reproduce the observed variability of stable water vapour isotopes on timescales of hours to days, thus allowing us to study the mechanisms that control the subtropical free-tropospheric humidity. Changes in isotopic signals along backward trajectories from the Canary Islands region reveal the physical processes behind the synoptic-scale isotope variability. We identify four predominant moisture transport pathways of mid-tropospheric air, each with distinct isotopic signatures: - air parcels originating from the convective boundary layer of the Saharan heat low (SHL) – these are characterised by a homogeneous isotopic composition with a particularly high δD (median mid-tropospheric δD=−122‰), which results from dry convective mixing of low-level moisture of diverse origin advected into the SHL; - air parcels originating from the free troposphere above the SHL – although experiencing the largest changes in humidity and δD during their subsidence over West Africa, these air parcels typically have lower δD values (median δD=−148‰) than air parcels originating from the boundary layer of the SHL; - air parcels originating from outside the SHL region, typically descending from tropical upper levels south of the SHL, which are often affected by moist convective injections from mesoscale convective systems in the Sahel – their isotopic composition is much less enriched in heavy isotopes (median δD=−175‰) than those from the SHL region; - air parcels subsiding from the upper-level extratropical North Atlantic – this pathway leads to the driest and most depleted conditions (median δD=−255‰) in the middle troposphere near the Canary Islands. The alternation of these transport pathways explains the observed high variability in humidity and δD on synoptic timescales to a large degree. We further show that the four different transport pathways are related to specific large-scale flow conditions. In particular, distinct differences in the location of the North African mid-level anticyclone and of extratropical Rossby wave patterns occur between the four transport pathways. Overall, this study demonstrates that the adopted Lagrangian isotope perspective enhances our understanding of air mass transport and mixing and offers a sound interpretation of the free-tropospheric variability of specific humidity and isotope composition on timescales of hours to days in contrasting atmospheric conditions over the eastern subtropical North Atlantic.

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: The concept of cloud radiative forcing (CRF) is commonly applied to quantify the impact of clouds on the surface radiative energy budget (REB). In the Arctic, specific radiative interactions between microphysical and macrophysical properties of clouds and the surface strongly modify the warming or cooling effect of clouds, complicating the estimate of CRF obtained from observations or models. Clouds tend to increase the broadband surface albedo over snow or sea ice surfaces compared to cloud-free conditions. However, this effect is not adequately considered in the derivation of CRF in the Arctic so far. Therefore, we have quantified the effects caused by surface-albedo–cloud interactions over highly reflective snow or sea ice surfaces on the CRF using radiative transfer simulations and below-cloud airborne observations above the heterogeneous springtime marginal sea ice zone (MIZ) during the Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) campaign. The impact of a modified surface albedo in the presence of clouds, as compared to cloud-free conditions, and its dependence on cloud optical thickness is found to be relevant for the estimation of the shortwave CRF. A method is proposed to consider this surface albedo effect on CRF estimates by continuously retrieving the cloud-free surface albedo from observations under cloudy conditions, using an available snow and ice albedo parameterization. Using ACLOUD data reveals that the estimated average shortwave cooling by clouds almost doubles over snow- and ice-covered surfaces (−62 W m−2 instead of −32 W m−2), if surface-albedo–cloud interactions are considered. As a result, the observed total (shortwave plus longwave) CRF shifted from a warming effect to an almost neutral one. Concerning the seasonal cycle of the surface albedo, it is demonstrated that this effect enhances shortwave cooling in periods when snow dominates the surface and potentially weakens the cooling by optically thin clouds during the summertime melting season. These findings suggest that the surface-albedo–cloud interaction should be considered in global climate models and in long-term studies to obtain a realistic estimate of the shortwave CRF to quantify the role of clouds in Arctic amplification.

Description: The evaluation of potential landslides in mountain areas is a very complex process. Currently, event understanding is scarce due to information limitations. Identifying the whole chain of events is not a straightforward task, and the impacts of mass-wasting processes depend on the conditions downstream of the origin. In this paper, we present an example that illustrates the complexities in the evaluation of the chain of events that may lead to a natural disaster. On 16 December 2017, a landslide occurred in the Yelcho mountain range (southern Chile). In that event, 7 million m3 of rocks and soil fell on the Yelcho glacier, depositing 2 million m3 on the glacier terminal, and the rest continued downstream, triggering a mudflow that hit Villa Santa Lucía in Chilean Patagonia and killing 22 people. The complex event was anticipated in the region by the National Geological and Mining Survey (Sernageomin in Spanish). However, the effects of the terrain characteristics along the run-out area were more significant than anticipated. In this work, we evaluate the conditions that enabled the mudflow that hit Villa Santa Lucía. We used the information generated by Sernageomin's professionals after the mudflow. We carried out geotechnical tests to characterize the soil. We simulated the mudflow using two hydrodynamic programs (r.avaflow and Flo-2D) that can handle the rheology of the water–soil mixture. Our results indicate that the soil is classified as volcanic pumices. This type of soil can be susceptible to the collapse of the structure when subjected to shearing (molding), flowing as a viscous liquid. From the numerical modeling, we concluded that r.avaflow performs better than Flo-2D. The mudflow was satisfactorily simulated using a water content in the mixture ranging from 30 % to 40 %, which would have required a source of about 3 million m3 of water. Coupling the simulations and the soil tests that we performed, we estimated that in the area scoured by the mudflow, there were probably around 2 800 000 m3 of water within the soil. Therefore, the conditions of the valley were crucial to enhancing the impacts of the landslide. This result is relevant because it highlights the importance of evaluating the complete chain of events to map hazards. We suggest that in future hazard mapping, geotechnical studies in combination with hydrodynamic simulation should be included, in particular when human lives are at risk.

Description: Production and reduction of nitrous oxide (N2O) by soil denitrifiers influence atmospheric concentrations of this potent greenhouse gas. Accurate projections of the net N2O flux have three key uncertainties: (1) short- vs. long-term responses to warming, (2) interactions among soil horizons, and (3) temperature responses of different steps in the denitrification pathway. We addressed these uncertainties by sampling soil from a boreal forest climate transect encompassing a 5.2 ∘C difference in the mean annual temperature and incubating the soil horizons in isolation and together at three ecologically relevant temperatures in conditions that promote denitrification. Both short-term exposure to warmer temperatures and long-term exposure to a warmer climate increased N2O emissions from organic and mineral soils; an isotopic tracer suggested that an increase in N2O production was more important than a decline in N2O reduction. Short-term warming promoted the reduction of organic horizon-derived N2O by mineral soil when these horizons were incubated together. The abundance of nirS (a precursor gene for N2O production) was not sensitive to temperature, whereas that of nosZ clade I (a gene for N2O reduction) decreased with short-term warming in both horizons and was higher from a warmer climate. These results suggest a decoupling of gene abundance and process rates in these soils that differs across horizons and timescales. In spite of these variations, our results suggest a consistent, positive response of denitrifier-mediated net N2O efflux rates to temperature across timescales in these boreal forests. Our work also highlights the importance of understanding cross-horizon N2O fluxes for developing a predictive understanding of net N2O efflux from soils.

Description: Wetlands are the largest and most uncertain natural sources of atmospheric methane (CH4). Several process-based models have been developed to quantify the magnitude and estimate spatial and temporal variations in CH4 emissions from global wetlands. Reliable models are required to estimate global wetland CH4 emissions. This study aimed to test two process-based models, CH4MODwetland and Terrestrial Ecosystem Model (TEM), against the CH4 flux measurements of marsh, swamp, peatland and coastal wetland sites across the world; specifically, model accuracy and generality were evaluated for different wetland types and in different continents, and then the global CH4 emissions from 2000 to 2010 were estimated. Both models showed similar high correlations with the observed seasonal/annual total CH4 emissions, and the regression of the observed versus computed total seasonal/annual CH4 emissions resulted in R2 values of 0.81 and 0.68 for CH4MODwetland and TEM, respectively. The CH4MODwetland produced accurate predictions for marshes, peatlands, swamps and coastal wetlands, with model efficiency (EF) values of 0.22, 0.52, 0.13 and 0.72, respectively. TEM produced good predictions for peatlands and swamps, with EF values of 0.69 and 0.74, respectively, but it could not accurately simulate marshes and coastal wetlands (EF 

Description: As knowledge about the cirrus clouds in the lower stratosphere is limited, reliable long-term measurements are needed to assess their characteristics, radiative impact and important role in upper troposphere and lower stratosphere (UTLS) chemistry. We used 6 years (2006–2012) of Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) measurements to investigate the global and seasonal distribution of stratospheric cirrus clouds and compared the MIPAS results with results derived from the latest version (V4.x) of the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) data. For the identification of stratospheric cirrus clouds, precise information on both the cloud top height (CTH) and the tropopause height is crucial. Here, we used lapse rate tropopause heights estimated from the ERA-Interim global reanalysis. Considering the uncertainties of the tropopause heights and the vertical sampling grid, we define CTHs more than 0.5 km above the tropopause as stratospheric for CALIPSO data. For MIPAS data, we took into account the coarser vertical sampling grid and the broad field of view so that we considered cirrus CTHs detected more than 0.75 km above the tropopause as stratospheric. Further sensitivity tests were conducted to rule out sampling artefacts in MIPAS data. The global distribution of stratospheric cirrus clouds was derived from night-time measurements because of the higher detection sensitivity of CALIPSO. In both data sets, MIPAS and CALIPSO, the stratospheric cirrus cloud occurrence frequencies are significantly higher in the tropics than in the extra-tropics. Tropical hotspots of stratospheric cirrus clouds associated with deep convection are located over equatorial Africa, South and Southeast Asia, the western Pacific, and South America. Stratospheric cirrus clouds were more often detected in December–February (15 %) than June–August (8 %) in the tropics (±20∘). At northern and southern middle latitudes (40–60∘), MIPAS observed about twice as many stratospheric cirrus clouds (occurrence frequencies of 4 %–5 % for MIPAS rather than about 2 % for CALIPSO). We attribute more frequent observations of stratospheric cirrus clouds with MIPAS to the higher detection sensitivity of the instrument to optically thin clouds. In contrast to the difference between daytime and night-time occurrence frequencies of stratospheric cirrus clouds by a factor of about 2 in zonal means in the tropics (4 % and 10 %, respectively) and at middle latitudes for CALIPSO data, there is little diurnal cycle in MIPAS data, in which the difference of occurrence frequencies in the tropics is about 1 percentage point in zonal mean and about 0.5 percentage point at middle latitudes. The difference between CALIPSO day and night measurements can also be attributed to their differences in detection sensitivity. Future work should focus on better understanding the origin of the stratospheric cirrus clouds and their impact on radiative forcing and climate.

Description: Air quality on our planet has been changing in particular since the industrial revolution (1750s) because of anthropogenic emissions. It is becoming increasingly important to realize air cleanliness, since clean air is as valuable a resource as clean water. A global standard for quantifying the level of air cleanliness is swiftly required, and we defined a novel concept, namely the Clean aIr Index (CII). The CII is a simple index defined by the normalization of the amount of a set of individual air pollutants. A CII value of 1 indicates completely clean air (no air pollutants), and 0 indicates the presence of air pollutants that meet the numerical environmental criteria for the normalization. In this time, the air pollutants used in the CII were taken from the Air Quality Guidelines (AQG) set by the World Health Organization (WHO), namely O3, particulate matters, NO2, and SO2. We chose Japan as a study area to evaluate CII because of the following reasons: (i) accurate validation data, as the in situ observation sites of the Atmospheric Environmental Regional Observation System (AEROS) provide highly accurate values of air pollutant amounts, and (ii) fixed numerical criteria from the Japanese Environmental Quality Standards (JEQS) as directed by the Ministry of the Environment (MOE) of Japan. We quantified air cleanliness in terms of the CII for the all 1896 municipalities in Japan and used data from Seoul and Beijing to evaluate Japanese air cleanliness. The amount of each air pollutant was calculated using a model that combined the Weather Research and Forecasting (WRF) and Community Multiscale Air Quality (CMAQ) models for 1 April 2014 to 31 March 2017. The CII values calculated by the WRF–CMAQ model and the AEROS measurements showed good agreement. The mean of the correlation coefficient for the CII values of 498 municipalities where the AEROS measurements operated was 0.66±0.05, which was higher than that of the Air Quality Index (AQI) of 0.57±0.06. The CII values averaged for the study period were 0.67, 0.52, and 0.24 in Tokyo, Seoul, and Beijing, respectively; thus, the air in Tokyo was 1.5 and 2.3 times cleaner, i.e., lower amounts of air pollutants, than the air in Seoul and Beijing, respectively. The average CII value for the all Japanese municipalities was 0.72 over the study period. The extremely clean air, CII ≈0.90, occurred in the southern remote islands of Tokyo and to the west of the Pacific coast, i.e., Kochi, Mie, and Wakayama prefectures during summer, with the transport of clean air from the ocean. We presented the top 100 clean air cities in Japan as one example of an application of CII in society. We confirmed that the CII enabled the quantitative evaluation of air cleanliness. The CII can be useful and valuable in various scenarios such as encouraging sightseeing and migration, investment and insurance business, and city planning. The CII is a simple and fair index that can be applied to all nations.

Description: Calculating a multi-model mean, a commonly used method for ensemble averaging, assumes model independence and equal model skill. Sharing of model components amongst families of models and research centres, conflated by growing ensemble size, means model independence cannot be assumed and is hard to quantify. We present a methodology to produce a weighted-model ensemble projection, accounting for model performance and model independence. Model weights are calculated by comparing model hindcasts to a selection of metrics chosen for their physical relevance to the process or phenomena of interest. This weighting methodology is applied to the Chemistry–Climate Model Initiative (CCMI) ensemble to investigate Antarctic ozone depletion and subsequent recovery. The weighted mean projects an ozone recovery to 1980 levels, by 2056 with a 95 % confidence interval (2052–2060), 4 years earlier than the most recent study. Perfect-model testing and out-of-sample testing validate the results and show a greater projective skill than a standard multi-model mean. Interestingly, the construction of a weighted mean also provides insight into model performance and dependence between the models. This weighting methodology is robust to both model and metric choices and therefore has potential applications throughout the climate and chemistry–climate modelling communities.

Description: Connections between vegetation and soil thermal dynamics are critical for estimating the vulnerability of permafrost to thaw with continued climate warming and vegetation changes. The interplay of complex biophysical processes results in a highly heterogeneous soil temperature distribution on small spatial scales. Moreover, the link between topsoil temperature and active layer thickness remains poorly constrained. Sixty-eight temperature loggers were installed at 1–3 cm depth to record the distribution of topsoil temperatures at the Trail Valley Creek study site in the northwestern Canadian Arctic. The measurements were distributed across six different vegetation types characteristic for this landscape. Two years of topsoil temperature data were analysed statistically to identify temporal and spatial characteristics and their relationship to vegetation, snow cover, and active layer thickness. The mean annual topsoil temperature varied between −3.7 and 0.1 ∘C within 0.5 km2. The observed variation can, to a large degree, be explained by variation in snow cover. Differences in snow depth are strongly related with vegetation type and show complex associations with late-summer thaw depth. While cold winter soil temperature is associated with deep active layers in the following summer for lichen and dwarf shrub tundra, we observed the opposite beneath tall shrubs and tussocks. In contrast to winter observations, summer topsoil temperature is similar below all vegetation types with an average summer topsoil temperature difference of less than 1 ∘C. Moreover, there is no significant relationship between summer soil temperature or cumulative positive degree days and active layer thickness. Altogether, our results demonstrate the high spatial variability of topsoil temperature and active layer thickness even within specific vegetation types. Given that vegetation type defines the direction of the relationship between topsoil temperature and active layer thickness in winter and summer, estimates of permafrost vulnerability based on remote sensing or model results will need to incorporate complex local feedback mechanisms of vegetation change and permafrost thaw.

Description: Evapotranspiration (ET) from tropical forests serves as a critical moisture source for regional and global climate cycles. However, the magnitude, seasonality, and interannual variability of ET in the Congo Basin remain poorly constrained due to a scarcity of direct observations, despite the Congo being the second-largest river basin in the world and containing a vast region of tropical forest. In this study, we applied a water balance model to an array of remotely sensed and in situ datasets to produce monthly, basin-wide ET estimates spanning April 2002 to November 2016. Data sources include water storage changes estimated from the Gravity Recovery and Climate Experiment (GRACE) satellites, in situ measurements of river discharge, and precipitation from several remotely sensed and gauge-based sources. An optimal precipitation dataset was determined as a weighted average of interpolated data by Nicholson et al. (2018), Climate Hazards InfraRed Precipitation with Station data version 2 (CHIRPS2) , and the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks–Climate Data Record product (PERSIANN-CDR), with the relative weights based on the error magnitudes of each dataset as determined by triple collocation. The resulting water-balance-derived ET (ETwb) features a long-term average that is consistent with previous studies (117.2±3.5 cm yr−1) but displays greater seasonal and interannual variability than seven global ET products. The seasonal cycle of ETwb generally tracks that of precipitation over the basin, with the exception that ETwb is greater in March–April–May (MAM) than in the relatively wetter September–October–November (SON) periods. This pattern appears to be driven by seasonal variations in the diffuse photosynthetically active radiation (PAR) fraction, net radiation (Rn), and soil water availability. From 2002 to 2016, Rn, PAR, and vapor-pressure deficit (VPD) all increased significantly within the Congo Basin; however, no corresponding trend occurred in ETwb. We hypothesize that the stability of ETwb over the study period despite sunnier and less humid conditions may be due to increasing atmospheric CO2 concentrations that offset the impacts of rising VPD and irradiance on stomatal water use efficiency (WUE).

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

Description: The meridional heat transport (MHT) in the North Atlantic is critically important to climate variability and the global overturning circulation. A wide range of ocean processes contribute to North Atlantic MHT, ranging from basin-scale overturning and gyre motions to mesoscale instabilities (such as eddies). However, previous analyses of “eddy” MHT in the region have mostly focused on the contributions of time-variable velocity and temperature, rather than considering the association of MHT with distinct spatial scales within the basin. In this study, a zonal spatial-scale decomposition separates large-scale from mesoscale velocity and temperature contributions to MHT, in order to characterize the physical processes driving MHT. Using this approach, we found that the mesoscale contributions to the time-mean and interannual/decadal (ID) variability of MHT in the latitude range 39–45∘ N are larger than large-scale horizontal contributions, though smaller than the overturning contributions. Considering the 40∘ N transect as a case study, large-scale ID variability is mostly generated close to the western boundary. In contrast, most ID MHT variability associated with mesoscales originates in two distinct regions: a western boundary region (70–60∘ W) associated with 1- to 4-year interannual variations and an interior region (50–35∘ W) associated with decadal variations. Surface eddy kinetic energy is not a reliable indicator of high MHT episodes, but the large-scale meridional temperature gradient is an important factor, by influencing the local temperature variance as well as the local correlation of velocity and temperature. Most of the mesoscale contribution to MHT at 40∘ N is associated with transient and propagating processes, but stationary mesoscale structures explain most of the mesoscale MHT south of the Gulf Stream separation, highlighting the differences between the temporal and spatial decomposition of meridional temperature fluxes.

Description: Accurate reference spectroscopic information for the water molecule from the microwave to the near-ultraviolet is of paramount importance in atmospheric research. A semi-empirical potential energy surface for the ground electronic state of H216O has been created by refining almost 4000 experimentally determined energy levels. These states extend into regions with large values of rotational and vibrational excitation. For all states considered in our refinement procedure, which extend to 37 000 cm−1 and J=20 (total angular momentum), the average root-mean-square deviation is approximately 0.05 cm−1. This potential energy surface offers significant improvements when compared to recent models by accurately predicting states possessing high values of J. This feature will offer significant improvements in calculated line positions for high-temperature spectra where transitions between high J states become more prominent. Combining this potential with the latest dipole moment surface for water vapour, a line list has been calculated which extends reliably to 37 000 cm−1. Obtaining reliable results in the ultraviolet is of special importance as it is a challenging spectral region for the water molecule both experimentally and theoretically. Comparisons are made against several experimental sources of cross sections in the near-ultraviolet and discrepancies are observed. In the near-ultraviolet our calculations are in agreement with recent atmospheric retrievals and the upper limit obtained using broadband spectroscopy by Wilson et al. (2016, p. 194), but they do not support recent suggestions of very strong absorption in this region.

Description: Various studies investigated the fate of evaporation and the origin of precipitation. The more recent studies among them were often carried out with the help of numerical moisture tracking. Many research questions could be answered within this context, such as dependencies of atmospheric moisture transfers between different regions, impacts of land cover changes on the hydrological cycle, sustainability-related questions, and questions regarding the seasonal and interannual variability of precipitation. In order to facilitate future applications, global datasets on the fate of evaporation and the sources of precipitation are needed. Since most studies are on a regional level and focus more on the sources of precipitation, the goal of this study is to provide a readily available global dataset on the fate of evaporation for a fine-meshed grid of source and receptor cells. The dataset was created through a global run of the numerical moisture tracking model Water Accounting Model-2layers (WAM-2layers) and focused on the fate of land evaporation. The tracking was conducted on a 1.5∘×1.5∘ grid and was based on reanalysis data from the ERA-Interim database. Climatic input data were incorporated in 3- to 6-hourly time steps and represent the time period from 2001 to 2018. Atmospheric moisture was tracked forward in time and the geographical borders of the model were located at ±79.5∘ latitude. As a result of the model run, the annual, the monthly and the interannual average fate of evaporation were determined for 8684 land grid cells (all land cells except those located within Greenland and Antarctica) and provided via source–receptor matrices. The gained dataset was complemented via an aggregation to country and basin scales in order to highlight possible usages for areas of interest larger than grid cells. This resulted in data for 265 countries and 8223 basins. Finally, five types of source–receptor matrices for average moisture transfers were chosen to build the core of the dataset: land grid cell to grid cell, country to grid cell, basin to grid cell, country to country, basin to basin. The dataset is, to our knowledge, the first ready-to-download dataset providing the overall fate of evaporation for land cells of a global fine-meshed grid in monthly resolution. At the same time, information on the sources of precipitation can be extracted from it. It could be used for investigations into average annual, seasonal, and interannual sink and source regions of atmospheric moisture from land masses for most of the regions in the world and shows various application possibilities for studying interactions between people and water, such as land cover changes or human water consumption patterns. The dataset is accessible under https://doi.org/10.1594/PANGAEA.908705 (Link et al., 2019a) and comes along with example scripts for reading and plotting the data.

Description: Optimum management of irrigated crops in regions with shallow saline groundwater requires a careful balance between application of irrigation water and upward movement of salinity from the groundwater. Few field-validated surrogate models are available to aid in the management of irrigation water under shallow groundwater conditions. The objective of this research is to develop a model that can aid in the management using a minimum of input data that are field validated. In this paper a 2-year field experiment was carried out in the Hetao irrigation district in Inner Mongolia, China, and a physically based integrated surrogate model for arid irrigated areas with shallow groundwater was developed and validated with the collected field data. The integrated model that links crop growth with available water and salinity in the vadose zone is called Evaluation of the Performance of Irrigated Crops and Soils (EPICS). EPICS recognizes that field capacity is reached when the matric potential is equal to the height above the groundwater table and thus not by a limiting hydraulic conductivity. In the field experiment, soil moisture contents and soil salt conductivity at five depths in the top 100 cm, groundwater depth, crop height, and leaf area index were measured in 2017 and 2018. The field results were used for calibration and validation of EPICS. Simulated and observed data fitted generally well during both calibration and validation. The EPICS model that can predict crop growth, soil water, groundwater depth, and soil salinity can aid in optimizing water management in irrigation districts with shallow aquifers.

Description: Recently, it has been established that interplanetary coronal mass ejections (ICMEs) can dramatically affect both trapped electron fluxes in the outer radiation belt and precipitating electron fluxes lost from the belt into the atmosphere. Precipitating electron fluxes and energies can vary over a range of timescales during these events. These variations depend on the initial energy and location of the electron population and the ICME characteristics and structures. One important factor controlling electron dynamics is the magnetic field orientation within the ejecta that is an integral part of the ICME. In this study, we examine Van Allen Probes (RBSPs) and Polar Orbiting Environmental Satellites (POESs) data to explore trapped and precipitating electron fluxes during two ICMEs. The ejecta in the selected ICMEs have magnetic cloud characteristics that exhibit the opposite sense of the rotation of the north–south magnetic field component (BZ). RBSP data are used to study trapped electron fluxes in situ, while POES data are used for electron fluxes precipitating into the upper atmosphere. The trapped and precipitating electron fluxes are qualitatively analysed to understand their variation in relation to each other and to the magnetic cloud rotation during these events. Inner magnetospheric wave activity was also estimated using RBSP and Geostationary Operational Environmental Satellite (GOES) data. In each event, the largest changes in the location and magnitude of both the trapped and precipitating electron fluxes occurred during the southward portion of the magnetic cloud. Significant changes also occurred during the end of the sheath and at the sheath–ejecta boundary for the cloud with south to north magnetic field rotation, while the ICME with north to south rotation had significant changes at the end boundary of the cloud. The sense of rotation of BZ and its profile also clearly affects the coherence of the trapped and/or precipitating flux changes, timing of variations with respect to the ICME structures, and flux magnitude of different electron populations. The differing electron responses could therefore imply partly different dominant acceleration mechanisms acting on the outer radiation belt electron populations as a result of opposite magnetic cloud rotation.

Description: In this paper, we quantify the CO2 and N2O emissions from denitrification over the Amazonian wetlands. The study concerns the entire Amazonian wetland ecosystem with a specific focus on three floodplain (FP) locations: the Branco FP, the Madeira FP and the FP alongside the Amazon River. We adapted a simple denitrification model to the case of tropical wetlands and forced it by open water surface extent products from the Soil Moisture and Ocean Salinity (SMOS) satellite. A priori model parameters were provided by in situ observations and gauging stations from the HYBAM Observatory. Our results show that the denitrification and the trace gas emissions present a strong cyclic pattern linked to the inundation processes that can be divided into three distinct phases: activation, stabilization and deactivation. We quantify the average yearly denitrification and associated emissions of CO2 and N2O over the entire watershed at 17.8 kgN ha−1 yr−1, 0.37 gC-CO2 m−2 yr−1 and 0.18 gN-N2O m−2 yr−1 respectively for the period 2011–2015. When compared to local observations, it was found that the CO2 emissions accounted for 0.01 % of the integrated ecosystem, which emphasizes the fact that minor changes to the land cover may induce strong impacts on the Amazonian carbon budget. Our results are consistent with the state of the art of global nitrogen models with a positive bias of 28 %. When compared to other wetlands in different pedoclimatic environments we found that the Amazonian wetlands have similar emissions of N2O with the Congo tropical wetlands and lower emissions than the temperate and tropical anthropogenic wetlands of the Garonne (France), the Rhine (Europe) and south-eastern Asia rice paddies. In summary our paper shows that a data-model-based approach can be successfully applied to quantify N2O and CO2 fluxes associated with denitrification over the Amazon basin. In the future, the use of higher-resolution remote sensing products from sensor fusion or new sensors like the Surface Water and Ocean Topography (SWOT) mission will permit the transposition of the approach to other large-scale watersheds in tropical environments.

Description: In the framework of the EU Copernicus programme, the European Centre for Medium-Range Weather Forecasts (ECMWF) on behalf of the Joint Research Centre (JRC) is forecasting daily fire weather indices using its medium-range ensemble prediction system. The use of weather forecasts in place of local observations can extend early warnings by up to 1–2 weeks, allowing for greater proactive coordination of resource-sharing and mobilization within and across countries. Using 1 year of pre-operational service in 2017 and the Fire Weather Index (FWI), here we assess the capability of the system globally and analyse in detail three major events in Chile, Portugal and California. The analysis shows that the skill provided by the ensemble forecast system extends to more than 10 d when compared to the use of mean climate, making a case for extending the forecast range to the sub-seasonal to seasonal timescale. However, accurate FWI prediction does not translate into accuracy in the forecast of fire activity globally. Indeed, when all fires detected in 2017 are considered, including agricultural- and human-induced burning, high FWI values only occur in 50 % of the cases and are limited to the Boreal regions. Nevertheless for very large events which were driven by weather conditions, FWI forecasts provide advance warning that could be instrumental in setting up management and containment strategies.

Description: North African dust reaches the southeastern United States every summer. Size-resolved dust mass measurements taken in Miami, Florida, indicate that more than one-half of the surface dust mass concentrations reside in particles with geometric diameters less than 2.1 µm, while vertical profiles of micropulse lidar depolarization ratios show dust reaching above 4 km during pronounced events. These observations are compared to the representation of dust in the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) aerosol reanalysis and closely related Goddard Earth Observing System model version 5 (GEOS-5) Forward Processing (FP) aerosol product, both of which assimilate satellite-derived aerosol optical depths using a similar protocol and inputs. These capture the day-to-day variability in aerosol optical depth well, in a comparison to an independent sun-photometer-derived aerosol optical depth dataset. Most of the modeled dust mass resides in diameters between 2 and 6 µm, in contrast to the measurements. Model-specified mass extinction efficiencies equate light extinction with approximately 3 times as much aerosol mass, in this size range, compared to the measured dust sizes. GEOS-5 FP surface-layer sea salt mass concentrations greatly exceed observed values, despite realistic winds and relative humidities. In combination, these observations help explain why, despite realistic total aerosol optical depths, (1) free-tropospheric model volume extinction coefficients are lower than those retrieved from the micro-pulse lidar, suggesting too-low model dust loadings in the free troposphere, and (2) model dust mass concentrations near the surface can be higher than those measured. The modeled vertical distribution of dust, when captured, is reasonable. Large, aspherical particles exceeding the modeled dust sizes are also occasionally present, but dust particles with diameters exceeding 10 µm contribute little to the measured total dust mass concentrations after such long-range transport. Remaining uncertainties warrant a further integrated assessment to confirm this study's interpretations.

Description: The impact of a hexamethyldisiloxane (HMDSO) treatment on the response of doped SnO2 sensors is investigated for acetone, carbon monoxide and hydrogen. The sensor was operated in temperature cycles based on the DSR concept (differential surface reduction). According to this concept, the rate constants for the reduction and oxidation of the surface after fast temperature changes can be evaluated and used for quantification of reducing gases as well as quantification and compensation of sensor poisoning by siloxanes, which is shown in this work. Increasing HMDSO exposure reduces the rate constants and therefore the sensitivity of the sensor more and more for all processes. On the other hand, while the rate constants for acetone and carbon monoxide are reduced nearly to zero already for short treatments, the hydrogen sensitivity remains fairly stable, which greatly increases the selectivity. During repeated HMDSO treatment the quasistatic sensitivity, i.e. equilibrium sensitivity at one point during the temperature cycle, rises at first for all gases but then drops rapidly for acetone and carbon monoxide, which can also be explained by reduced rate constants for oxygen chemisorption on the sensor surface when considering the generation of surface charge.

Description: Anaerobic nitrate-dependent Fe(II) oxidation (NDFeO) is widespread in various aquatic environments and plays a major role in iron and nitrogen redox dynamics. However, evidence for truly enzymatic, autotrophic NDFeO remains limited, with alternative explanations involving the coupling of heterotrophic denitrification with the abiotic oxidation of structurally bound or aqueous Fe(II) by reactive intermediate nitrogen (N) species (chemodenitrification). The extent to which chemodenitrification is caused (or enhanced) by ex vivo surface catalytic effects has not been directly tested to date. To determine whether the presence of either an Fe(II)-bearing mineral or dead biomass (DB) catalyses chemodenitrification, two different sets of anoxic batch experiments were conducted: 2 mM Fe(II) was added to a low-phosphate medium, resulting in the precipitation of vivianite (Fe3(PO4)2), to which 2 mM nitrite (NO2-) was later added, with or without an autoclaved cell suspension (∼1.96×108 cells mL−1) of Shewanella oneidensis MR-1. Concentrations of nitrite (NO2-), nitrous oxide (N2O), and iron (Fe2+, Fetot) were monitored over time in both set-ups to assess the impact of Fe(II) minerals and/or DB as catalysts of chemodenitrification. In addition, the natural-abundance isotope ratios of NO2- and N2O (δ15N and δ18O) were analysed to constrain the associated isotope effects. Up to 90 % of the Fe(II) was oxidized in the presence of DB, whereas only ∼65 % of the Fe(II) was oxidized under mineral-only conditions, suggesting an overall lower reactivity of the mineral-only set-up. Similarly, the average NO2- reduction rate in the mineral-only experiments (0.004±0.003 mmol L−1 d−1) was much lower than in the experiments with both mineral and DB (0.053±0.013 mmol L−1 d−1), as was N2O production (204.02±60.29 nmol L−1 d−1). The N2O yield per mole NO2- reduced was higher in the mineral-only set-ups (4 %) than in the experiments with DB (1 %), suggesting the catalysis-dependent differential formation of NO. N-NO2- isotope ratio measurements indicated a clear difference between both experimental conditions: in contrast to the marked 15N isotope enrichment during active NO2- reduction (15εNO2=+10.3 ‰) observed in the presence of DB, NO2- loss in the mineral-only experiments exhibited only a small N isotope effect (

Description: The Mediterranean (MED) Basin is a climate change hotspot that has seen drying and a pronounced increase in heatwaves over the last century. At the same time, it is experiencing increased heavy precipitation during wintertime cold spells. Understanding and quantifying the risks from compound events over the MED is paramount for present and future disaster risk reduction measures. Here, we apply a novel method to study compound events based on dynamical systems theory and analyse compound temperature and precipitation events over the MED from 1979 to 2018. The dynamical systems analysis quantifies the strength of the coupling between different atmospheric variables over the MED. Further, we consider compound warm–dry anomalies in summer and cold–wet anomalies in winter. Our results show that these warm–dry and cold–wet compound days are associated with large values of the temperature–precipitation coupling parameter of the dynamical systems analysis. This indicates that there is a strong interaction between temperature and precipitation during compound events. In winter, we find no significant trend in the coupling between temperature and precipitation. However in summer, we find a significant upward trend which is likely driven by a stronger coupling during warm and dry days. Thermodynamic processes associated with long-term MED warming can best explain the trend, which intensifies compound warm–dry events.

Description: Anhand von sedimentologischen und geländemorphologischen Untersuchungen wird die Abschmelzgeschichte des südöstlichen Chiemsee-Gletschers beschrieben. Mit dem Trockenfallen der Bad Adelholzen-Erlstätter Rinne im Verlaufe des Spätwürm entwickelt sich aus dem Abschmelzen des Eislappens in der Grabenstätter Bucht eine sich ständig tiefer legende konzentrische Abfolge von zunächst peripheren Entwässerungsrinnen, wobei die ältesten Rinnen dieser Phase bei Chieming, die jüngeren dann entsprechend weiter im Süden, in die zentripetale Richtung umschwenken. Die Entstehung des Tüttensee-Komplexes ist im Kontext dieser Entwicklung zu sehen. Er ist das Ergebnis der glazifluvialen und glazilakustrinen Sedimentation im Einflussbereich des sukzessiven Eisabbaus in der Grabenstätter Bucht in Kombination mit einer Toteisbildung im Bereich des heutigen Tüttensees. Dafür sprechen die stufenartige Abfolge der beschriebenen peripheren Abflussrinnen mit ihren immer tiefer liegenden Abflussniveaus, die Höhengleichheit von drei dieser Rinnen mit den Tüttensee-Terrassen sowie die für die jeweilige Terrassenentstehung typische glazifluviale bzw. delta-artige Sedimentstruktur und -reife. Dieses Ergebnis stellt ein Korrektiv zur Hypothese des Chiemgau-Impakts dar, wonach der Tüttensee ein Impaktkrater sein soll. Da diese nun falsifizierte Annahme vor allem im deutschsprachigen Raum von zahlreichen Medien propagiert wird, ist der folgende Artikel auf Deutsch verfasst, um einer breiten Leserschaft zugänglich zu sein.

Description: Temperature appears to be the best predictor of species composition of planktonic foraminifera communities, making it possible to use their fossil assemblages to reconstruct sea surface temperature (SST) variation in the past. However, the role of other environmental factors potentially modulating the spatial and vertical distribution of planktonic foraminifera species is poorly understood. This is especially relevant for environmental factors affecting the subsurface habitat. If such factors play a role, changes in the abundance of subsurface-dwelling species may not solely reflect SST variation. In order to constrain the effect of subsurface parameters on species composition, we here characterize the vertical distribution of living planktonic foraminifera community across an east–west transect through the subtropical South Atlantic Ocean, where SST variability was small, but the subsurface water mass structure changed dramatically. Four planktonic foraminifera communities could be identified across the top 700 m of the transect. Gyre and Agulhas Leakage surface faunas were predominantly composed of Globigerinoides ruber, Globigerinoides tenellus, Trilobatus sacculifer, Globoturborotalita rubescens, Globigerinella calida, Tenuitella iota, and Globigerinita glutinata, and these only differed in terms of relative abundances (community composition). Upwelling fauna was dominated by Neogloboquadrina pachyderma, Neogloboquadrina incompta, Globorotalia crassaformis, and Globorotalia inflata. Thermocline fauna was dominated by Tenuitella fleisheri, Globorotalia truncatulinoides, and Globorotalia scitula in the west and by G. scitula only in the east. The largest part of the standing stock was consistently found in the surface layer, but SST was not the main predictor of species composition either for the depth-integrated fauna across the stations or at individual depth layers. Instead, we identified a pattern of vertical stacking of different parameters controlling species composition, reflecting different aspects of the pelagic habitat. Whereas productivity appears to dominate in the mixed layer (0–60 m), physical properties (temperature, salinity) become important at intermediate depths and in the subsurface, a complex combination of factors including oxygen concentration is required to explain the assemblage composition. These results indicate that the seemingly straightforward relationship between assemblage composition and SST in sedimentary assemblages reflects vertically and seasonally integrated processes that are only indirectly linked to SST. It also implies that fossil assemblages of planktonic foraminifera should also contain a signature of subsurface processes, which could be used for paleoceanographic reconstructions.

Description: Light-use efficiency defines the ability of primary producers to convert sunlight energy to primary production and is computed as the ratio between the gross primary production and the intercepted photosynthetic active radiation. While this measure has been applied broadly within terrestrial ecology to investigate habitat resource-use efficiency, it remains underused within the aquatic realm. This report provides a conceptual framework to compute hourly and daily light-use efficiency using underwater O2 eddy covariance, a recent technological development that produces habitat-scale rates of primary production under unaltered in situ conditions. The analysis, tested on two benthic flux datasets, documents that hourly light-use efficiency may approach the theoretical limit of 0.125 O2 per photon under low-light conditions, but it decreases rapidly towards the middle of the day and is typically 10-fold lower on a 24 h basis. Overall, light-use efficiency provides a useful measure of habitat functioning and facilitates site comparison in time and space.

Description: Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean, where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the first Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF1) in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96∘ S, 14.35∘ W). Using observed profiles of temperature, humidity, and clouds from the field campaign alongside aerosol optical properties from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM), profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. Idealized experiments using the RRTM with and without aerosols and a range of values for the single-scattering albedo (SSA) demonstrate that shortwave (SW) heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K d−1 depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA, and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave (LW) cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single-scattering albedo assumptions in models.

Description: Information on global cropland distribution and agricultural production is critical for the world's agricultural monitoring and food security. We present datasets of cropland extent and agricultural production in a two-paper series of a cultivated planet in 2010. In the first part, we propose a new Self-adapting Statistics Allocation Model (SASAM) to develop the global map of cropland distribution. SASAM is based on the fusion of multiple existing cropland maps and multilevel statistics of the cropland area, which is independent of training samples. First, cropland area statistics are used to rank the input cropland maps, and then a scoring table is built to indicate the agreement among the input datasets. Secondly, statistics are allocated adaptively to the pixels with higher agreement scores until the cumulative cropland area is close to the statistics. The multilevel allocation results are then integrated to obtain the extent of cropland. We applied SASAM to produce a global cropland synergy map with a 500 m spatial resolution for circa 2010. The accuracy assessments show that the synergy map has higher accuracy than the input datasets and better consistency with the cropland statistics. The synergy cropland map is available via an open-data repository (https://doi.org/10.7910/DVN/ZWSFAA; Lu et al., 2020). This new cropland map has been used as an essential input to the Spatial Production Allocation Model (SPAM) for producing the global dataset of agricultural production for circa 2010, which is described in the second part of the two-paper series.

Description: The western European kinematic evolution results from the opening of the western Neotethys and the Atlantic oceans since the late Paleozoic and the Mesozoic. Geological evidence shows that the Iberian domain recorded the propagation of these two oceanic systems well and is therefore a key to significantly advancing our understanding of the regional plate reconstructions. The late-Permian–Triassic Iberian rift basins have accommodated extension, but this tectonic stage is often neglected in most plate kinematic models, leading to the overestimation of the movements between Iberia and Europe during the subsequent Mesozoic (Early Cretaceous) rift phase. By compiling existing seismic profiles and geological constraints along the North Atlantic margins, including well data over Iberia, as well as recently published kinematic and paleogeographic reconstructions, we propose a coherent kinematic model of Iberia that accounts for both the Neotethyan and Atlantic evolutions. Our model shows that the Europe–Iberia plate boundary was a domain of distributed and oblique extension made of two rift systems in the Pyrenees and in the Iberian intra-continental basins. It differs from standard models that consider left-lateral strike-slip movement localized only in the northern Pyrenees in introducing a significant strike-slip movement south of the Ebro block. At a larger scale it emphasizes the role played by the late-Permian–Triassic rift and magmatism, as well as strike-slip faulting in the evolution of the western Neotethys Ocean and their control on the development of the Atlantic rift.

Description: This paper addresses the question of how much uncertainties in CO2 fluxes over Australia can be reduced by assimilation of total-column carbon dioxide retrievals from the Orbiting Carbon Observatory-2 (OCO-2) satellite instrument. We apply a four-dimensional variational data assimilation system, based around the Community Multiscale Air Quality (CMAQ) transport-dispersion model. We ran a series of observing system simulation experiments to estimate posterior error statistics of optimized monthly-mean CO2 fluxes in Australia. Our assimilations were run with a horizontal grid resolution of 81 km using OCO-2 data for 2015. Based on four representative months, we find that the integrated flux uncertainty for Australia is reduced from 0.52 to 0.13 Pg C yr−1. Uncertainty reductions of up to 90 % were found at grid-point resolution over productive ecosystems. Our sensitivity experiments show that the choice of the correlation structure in the prior error covariance plays a large role in distributing information from the observations. We also found that biases in the observations would significantly impact the inverted fluxes and could contaminate the final results of the inversion. Biases in prior fluxes are generally removed by the inversion system. Biases in the boundary conditions have a significant impact on retrieved fluxes, but this can be mitigated by including boundary conditions in our retrieved parameters. In general, results from our idealized experiments suggest that flux inversions at this unusually fine scale will yield useful information on the carbon cycle at continental and finer scales.

Description: In this paper, we present the first multiyear time series of atmospheric ammonia (NH3) ground-based measurements in the Paris region (Créteil, 48.79∘ N, 2.44∘ E, France) retrieved with the midresolution “Observations of the Atmosphere by Solar absorption Infrared Spectroscopy” (OASIS) ground-based Fourier transform infrared solar observatory. Located in an urban region, OASIS has previously been used for monitoring air quality (tropospheric ozone and carbon monoxide) thanks to its specific column sensitivity across the whole troposphere down to the atmospheric boundary layer. A total of 4920 measurements of atmospheric total columns of ammonia have been obtained from 2009 to 2017, with uncertainties ranging from 20 % to 35 %, and have been compared with NH3 concentrations derived from the Infrared Atmospheric Sounding Interferometer (IASI). OASIS ground-based measurements show significant interannual and seasonal variabilities of atmospheric ammonia. NH3 total columns over the Paris megacity (12 million people) vary seasonally by 2 orders of magnitude from approximately 0.1×1016 molec. cm−2 in winter to 10×1016 molec. cm−2 for spring peaks, probably due to springtime spreading of fertilizers on surrounding croplands.

Description: In order to track progress towards the global climate targets, the parties that signed the Paris Climate Agreement will regularly report their anthropogenic carbon dioxide (CO2) emissions based on energy statistics and CO2 emission factors. Independent evaluation of this self-reporting system is a fast-growing research topic. Here, we study the value of satellite observations of the column CO2 concentrations to estimate CO2 anthropogenic emissions with 5 years of the Orbiting Carbon Observatory-2 (OCO-2) retrievals over and around China. With the detailed information of emission source locations and the local wind, we successfully observe CO2 plumes from 46 cities and industrial regions over China and quantify their CO2 emissions from the OCO-2 observations, which add up to a total of 1.3 Gt CO2 yr−1 that accounts for approximately 13 % of mainland China's annual emissions. The number of cities whose emissions are constrained by OCO-2 here is 3 to 10 times larger than in previous studies that only focused on large cities and power plants in different locations around the world. Our satellite-based emission estimates are broadly consistent with the independent values from China's detailed emission inventory MEIC but are more different from those of two widely used global gridded emission datasets (i.e., EDGAR and ODIAC), especially for the emission estimates for the individual cities. These results demonstrate some skill in the satellite-based emission quantification for isolated source clusters with the OCO-2, despite the sparse sampling of this instrument not designed for this purpose. This skill can be improved by future satellite missions that will have a denser spatial sampling of surface emitting areas, which will come soon in the early 2020s.

Description: We develop a new large-scale hydrological and water resources model, the Community Water Model (CWatM), which can simulate hydrology both globally and regionally at different resolutions from 30 arcmin to 30 arcsec at daily time steps. CWatM is open source in the Python programming environment and has a modular structure. It uses global, freely available data in the netCDF4 file format for reading, storage, and production of data in a compact way. CWatM includes general surface and groundwater hydrological processes but also takes into account human activities, such as water use and reservoir regulation, by calculating water demands, water use, and return flows. Reservoirs and lakes are included in the model scheme. CWatM is used in the framework of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP), which compares global model outputs. The flexible model structure allows for dynamic interaction with hydro-economic and water quality models for the assessment and evaluation of water management options. Furthermore, the novelty of CWatM is its combination of state-of-the-art hydrological modeling, modular programming, an online user manual and automatic source code documentation, global and regional assessments at different spatial resolutions, and a potential community to add to, change, and expand the open-source project. CWatM also strives to build a community learning environment which is able to freely use an open-source hydrological model and flexible coupling possibilities to other sectoral models, such as energy and agriculture.

Description: In the past decades, severe heat waves have frequently occurred in many parts of the world. These conspicuous heat waves exerted terrible influences on human health, society, the economy, agriculture, the ecosystem and so on. Based on observed daily temperatures in China, an integrated index of heat waves and extreme-temperature days was established involving the frequency, duration, intensity and scale of these events across large cities in China. Heat waves and extreme-temperature days showed an increasing trend in most regions except northwest China from 1955 to 2014. After the late 1980s, the increasing trend was more obvious than the decades before. The cities in the middle and lower reaches of the Yangtze River were threatened by the most serious heat events in the past 60 years, especially Chongqing and Changsha. Due to the subtropical monsoon climate and special terrain, Chongqing experienced the most heat events in a long period of time. In particular, there was obvious fluctuation of hot years in 31 cities, which did not continuously rise with global warming; 21 cities mainly located in the eastern and southern regions of China had an obvious rising trend; eight cities had a clear declining trend which was mainly distributed in the western and northern regions of China; and there were no extreme-temperature days in Kunming and Lhasa in the past 60 years. The study revealed an obvious differentiation of heat events for 31 cities under climate change; heat threat in most cities is increasing but declining or remaining unchanged in the other cities. The trend is likely to intensify with global warming.

Description: Understanding the behavior of halogens (Cl, Br, and I) in subduction zones is critical to constrain the geochemical cycle of these volatiles and associated trace metals, as well as to quantify the halogen fluxes to the atmosphere via volcanic degassing. Here, the partitioning of bromine between coexisting aqueous fluids and hydrous granitic melts and its speciation in slab-derived fluids have been investigated in situ up to 840 ∘C and 2.2 GPa by synchrotron x-ray fluorescence (SXRF) and x-ray absorption spectroscopy (XAS) in diamond anvil cells. The partition coefficients DBrf/m range from ∼2 to ∼15, with an average value of 6.7±3.6 (1σ) over the whole pressure–temperature (P–T) range, indicating a moderate Br enrichment in aqueous fluids, in agreement with previous work. Extended x-ray-absorption fine-structure (EXAFS) analysis further evidences a gradual evolution of Br speciation from hydrated Br ions [Br(H2O)6]− in slab dehydration fluids to more complex structures involving both Na ions and water molecules, [BrNax(H2O)y], in hydrous silicate melts and supercritical fluids released at greater depth (〉 200 km). In denser fluids (ρ 〉 1.5 g cm−3) containing 60 wt % dissolved alkali–silicates and in hydrous Na2Si2O5 melts (10 wt % H2O), Br is found to be in a “salt-like” structure involving the six nearest Na ions and several next-nearest O neighbors that are either from water molecules and/or the silicate network. Bromine (and likely chlorine and iodine) complexing with alkalis is thus an efficient mechanism for the mobilization and transport of halogens by hydrous silicate melts and silica-rich supercritical fluids. Our results suggest that both shallow dehydration fluids and deeper silicate-bearing fluids efficiently remove halogens from the slab in the sub-arc region, thus favoring an efficient transfer of halogens across subduction zones.

Description: The total Antarctic sea ice extent (SIE) experiences a distinct annual cycle, peaking in September and reaching its minimum in February. In this paper we propose a mathematical and statistical decomposition of this temporal variation in SIE. Each component is interpretable and, when combined, gives a complete picture of the variation in the sea ice. We consider timescales varying from the instantaneous and not previously defined to the multi-decadal curvilinear trend, the longest. Because our representation is daily, these timescales of variability give precise information about the timing and rates of advance and retreat of the ice and may be used to diagnose physical contributors to variability in the sea ice. We define a number of annual cycles each capturing different components of variation, especially the yearly amplitude and phase that are major contributors to SIE variation. Using daily sea ice concentration data, we show that our proposed invariant annual cycle explains 29 % more of the variation in daily SIE than the traditional method. The proposed annual cycle that incorporates amplitude and phase variation explains 77 % more variation than the traditional method. The variation in phase explains more of the variability in SIE than the amplitude. Using our methodology, we show that the anomalous decay of sea ice in 2016 was associated largely with a change of phase rather than amplitude. We show that the long term trend in Antarctic sea ice extent is strongly curvilinear and the reported positive linear trend is small and dependent strongly on a positive trend that began around 2011 and continued until 2016.