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

Description: Aerosol–cloud interactions are the largest source of uncertainty in quantifying anthropogenic radiative forcing. The large uncertainty is, in part, due to the difficulty of predicting cloud microphysical parameters, such as the cloud droplet number concentration (Nd). Even though rigorous first-principle approaches exist to calculate Nd, the cloud and aerosol research community also relies on empirical approaches such as relating Nd to aerosol mass concentration. Here we analyze relationships between Nd and cloud water chemical composition, in addition to the effect of environmental factors on the degree of the relationships. Warm, marine, stratocumulus clouds off the California coast were sampled throughout four summer campaigns between 2011 and 2016. A total of 385 cloud water samples were collected and analyzed for 80 chemical species. Single- and multispecies log–log linear regressions were performed to predict Nd using chemical composition. Single-species regressions reveal that the species that best predicts Nd is total sulfate (Radj2=0.40). Multispecies regressions reveal that adding more species does not necessarily produce a better model, as six or more species yield regressions that are statistically insignificant. A commonality among the multispecies regressions that produce the highest correlation with Nd was that most included sulfate (either total or non-sea-salt), an ocean emissions tracer (such as sodium), and an organic tracer (such as oxalate). Binning the data according to turbulence, smoke influence, and in-cloud height allowed for examination of the effect of these environmental factors on the composition–Nd correlation. Accounting for turbulence, quantified as the standard deviation of vertical wind speed, showed that the correlation between Nd with both total sulfate and sodium increased at higher turbulence conditions, consistent with turbulence promoting the mixing between ocean surface and cloud base. Considering the influence of smoke significantly improved the correlation with Nd for two biomass burning tracer species in the study region, specifically oxalate and iron. When binning by in-cloud height, non-sea-salt sulfate and sodium correlated best with Nd at cloud top, whereas iron and oxalate correlated best with Nd at cloud base.

Description: Temperature, H2O, and O3 profiles, as well as CO2, N2O, CH4, chlorofluorocarbon-12 (CFC-12), and sea surface temperature (SST) scalar anomalies are computed using a clear subset of AIRS observations over ocean for the first 16 years of NASA's Earth-Observing Satellite (EOS) Aqua Atmospheric Infrared Sounder (AIRS) operation. The AIRS Level-1c radiances are averaged over 16 d and 40 equal-area zonal bins and then converted to brightness temperature anomalies. Geophysical anomalies are retrieved from the brightness temperature anomalies using a relatively standard optimal estimation approach. The CO2, N2O, CH4, and CFC-12 anomalies are derived by applying a vertically uniform multiplicative shift to each gas in order to obtain an estimate for the gas mixing ratio. The minor-gas anomalies are compared to the National Oceanic and Atmospheric Administration (NOAA) Earth System Research Laboratory (ESRL) in situ values and used to estimate the radiometric stability of the AIRS radiances. Similarly, the retrieved SST anomalies are compared to the SST values used in the ERA-Interim reanalysis and to NOAA's Optimum Interpolation SST (OISST) product. These intercomparisons strongly suggest that many AIRS channels are stable to better than 0.02 to 0.03 K per decade, well below climate trend levels, indicating that the AIRS blackbody is not drifting. However, detailed examination of the anomaly retrieval residuals (observed – computed) shows various small unphysical shifts that correspond to AIRS hardware events (shutdowns, etc.). Some examples are given highlighting how the AIRS radiance stability could be improved, especially for channels sensitive to N2O and CH4. The AIRS shortwave channels exhibit larger drifts that make them unsuitable for climate trending, and they are avoided in this work. The AIRS Level 2 surface temperature retrievals only use shortwave channels. We summarize how these shortwave drifts impacts recently published comparisons of AIRS surface temperature trends to other surface climatologies.

Description: The CloudRoots field experiment was designed to obtain a comprehensive observational dataset that includes soil, plant, and atmospheric variables to investigate the interaction between a heterogeneous land surface and its overlying atmospheric boundary layer at the sub-hourly and sub-kilometre scale. Our findings demonstrate the need to include measurements at leaf level to better understand the relations between stomatal aperture and evapotranspiration (ET) during the growing season at the diurnal scale. Based on these observations, we obtain accurate parameters for the mechanistic representation of photosynthesis and stomatal aperture. Once the new parameters are implemented, the model reproduces the stomatal leaf conductance and the leaf-level photosynthesis satisfactorily. At the canopy scale, we find a consistent diurnal pattern on the contributions of plant transpiration and soil evaporation using different measurement techniques. From highly resolved vertical profile measurements of carbon dioxide (CO2) and other state variables, we infer a profile of the CO2 assimilation in the canopy with non-linear variations with height. Observations taken with a laser scintillometer allow us to quantify the non-steadiness of the surface turbulent fluxes during the rapid changes driven by perturbation of photosynthetically active radiation by cloud flecks. More specifically, we find 2 min delays between the cloud radiation perturbation and ET. To study the relevance of advection and surface heterogeneity for the land–atmosphere interaction, we employ a coupled surface–atmospheric conceptual model that integrates the surface and upper-air observations made at different scales from leaf to the landscape. At the landscape scale, we calculate a composite sensible heat flux by weighting measured fluxes with two different land use categories, which is consistent with the diurnal evolution of the boundary layer depth. Using sun-induced fluorescence measurements, we also quantify the spatial variability of ET and find large variations at the sub-kilometre scale around the CloudRoots site. Our study shows that throughout the entire growing season, the wide variations in stomatal opening and photosynthesis lead to large diurnal variations of plant transpiration at the leaf, plant, canopy, and landscape scales. Integrating different advanced instrumental techniques with modelling also enables us to determine variations of ET that depend on the scale where the measurement were taken and on the plant growing stage.

Description: Most available verification metrics for ensemble forecasts focus on univariate quantities. That is, they assess whether the ensemble provides an adequate representation of the forecast uncertainty about the quantity of interest at a particular location and time. For spatially indexed ensemble forecasts, however, it is also important that forecast fields reproduce the spatial structure of the observed field and represent the uncertainty about spatial properties such as the size of the area for which heavy precipitation, high winds, critical fire weather conditions, etc., are expected. In this article we study the properties of the fraction of threshold exceedance (FTE) histogram, a new diagnostic tool designed for spatially indexed ensemble forecast fields. Defined as the fraction of grid points where a prescribed threshold is exceeded, the FTE is calculated for the verification field and separately for each ensemble member. It yields a projection of a – possibly high-dimensional – multivariate quantity onto a univariate quantity that can be studied with standard tools like verification rank histograms. This projection is appealing since it reflects a spatial property that is intuitive and directly relevant in applications, though it is not obvious whether the FTE is sufficiently sensitive to misrepresentation of spatial structure in the ensemble. In a comprehensive simulation study we find that departures from uniformity of the FTE histograms can indeed be related to forecast ensembles with biased spatial variability and that these histograms detect shortcomings in the spatial structure of ensemble forecast fields that are not obvious by eye. For demonstration, FTE histograms are applied in the context of spatially downscaled ensemble precipitation forecast fields from NOAA's Global Ensemble Forecast System.

Description: In this paper, the performance of three machine-learning methods for predicting short-term evolution and for reproducing the long-term statistics of a multiscale spatiotemporal Lorenz 96 system is examined. The methods are an echo state network (ESN, which is a type of reservoir computing; hereafter RC–ESN), a deep feed-forward artificial neural network (ANN), and a recurrent neural network (RNN) with long short-term memory (LSTM; hereafter RNN–LSTM). This Lorenz 96 system has three tiers of nonlinearly interacting variables representing slow/large-scale (X), intermediate (Y), and fast/small-scale (Z) processes. For training or testing, only X is available; Y and Z are never known or used. We show that RC–ESN substantially outperforms ANN and RNN–LSTM for short-term predictions, e.g., accurately forecasting the chaotic trajectories for hundreds of numerical solver's time steps equivalent to several Lyapunov timescales. The RNN–LSTM outperforms ANN, and both methods show some prediction skills too. Furthermore, even after losing the trajectory, data predicted by RC–ESN and RNN–LSTM have probability density functions (pdf's) that closely match the true pdf – even at the tails. The pdf of the data predicted using ANN, however, deviates from the true pdf. Implications, caveats, and applications to data-driven and data-assisted surrogate modeling of complex nonlinear dynamical systems, such as weather and climate, are discussed.

Description: In the western Mediterranean Sea, Levantine intermediate waters (LIW), which circulate below the surface productive zone, progressively accumulate nutrients along their pathway from the Tyrrhenian Sea to the Algerian Basin. This study addresses the role played by diffusion in the nutrient enrichment of the LIW, a process particularly relevant inside step-layer structures extending down to deep waters – structures known as thermohaline staircases. Profiling float observations confirmed that staircases develop over epicentral regions confined in large-scale circulation features and maintained by saltier LIW inflows on the periphery. Thanks to a high profiling frequency over the 4-year period 2013–2017, float observations reveal the temporal continuity of the layering patterns encountered during the cruise PEACETIME and document the evolution of layer properties by about +0.06 ∘C in temperature and +0.02 in salinity. In the Algerian Basin, the analysis of in situ lateral density ratios untangled double-diffusive convection as a driver of thermohaline changes inside epicentral regions and isopycnal diffusion as a driver of heat and salt exchanges with the surrounding sources. In the Tyrrhenian Sea, the nitrate flux across thermohaline staircases, as opposed to the downward salt flux, contributes up to 25 % of the total nitrate pool supplied to the LIW by vertical transfer. Overall, however, the nutrient enrichment of the LIW is driven mostly by other sources, coastal or atmospheric, as well as by inputs advected from the Algerian Basin.

Description: Geologic dating methods for the most part do not directly measure ages. Instead, interpreting a geochemical observation as a geologically useful parameter – an age or a rate – requires an interpretive middle layer of calculations and supporting data sets. These are the subject of active research and evolve rapidly, so any synoptic analysis requires repeated recalculation of large numbers of ages from a growing data set of raw observations, using a constantly improving calculation method. Many important applications of geochronology involve regional or global analyses of large and growing data sets, so this characteristic is an obstacle to progress in these applications. This paper describes the ICE-D (Informal Cosmogenic-Nuclide Exposure-age Database) database project, a prototype computational infrastructure for dealing with this obstacle in one geochronological application – cosmogenic-nuclide exposure dating – that aims to enable visualization or analysis of diverse data sets by making middle-layer calculations dynamic and transparent to the user. An important aspect of this concept is that it is designed as a forward-looking research tool rather than a backward-looking archive: only observational data (which do not become obsolete) are stored, and derived data (which become obsolete as soon as the middle-layer calculations are improved) are not stored but instead calculated dynamically at the time data are needed by an analysis application. This minimizes “lock-in” effects associated with archiving derived results subject to rapid obsolescence and allows assimilation of both new observational data and improvements to middle-layer calculations without creating additional overhead at the level of the analysis application.

Description: The invasion of aquifer microbial communities by aboveground microorganisms, a phenomenon known as community coalescence, is likely to be exacerbated in groundwaters fed by stormwater infiltration systems (SISs). Here, the incidence of this increased connectivity with upslope soils and impermeabilized surfaces was assessed through a meta-analysis of 16S rRNA gene libraries. Specifically, DNA sequences encoding 16S rRNA V5-V6 regions from free-living and attached aquifer bacteria (i.e., water and biofilm samples) were analysed upstream and downstream of a SIS and compared with those from bacterial communities from watershed runoffs and surface sediments from the SIS detention and infiltration basins. Significant bacterial transfers were inferred by the SourceTracker Bayesian approach, with 23 % to 57 % of the aquifer bacterial biofilms being composed of taxa from aboveground sediments and urban runoffs. Sediments from the detention basin were found more significant contributors of taxa involved in the buildup of these biofilms than soils from the infiltration basin. Inferred taxa among the coalesced biofilm community were predicted to be high in hydrocarbon degraders such as Sphingobium and Nocardia. The 16S rRNA-based bacterial community structure of the downstream-SIS aquifer waters showed lower coalescence with aboveground taxa (8 % to 38 %) than those of biofilms and higher numbers of taxa predicted to be involved in the N and S cycles. A DNA marker named tpm enabled the tracking of bacterial species from 24 genera including Pseudomonas, Aeromonas and Xanthomonas, among these communities. Several tpm sequence types were found to be shared between the aboveground and aquifer samples. Reads related to Pseudomonas were allocated to 50 species, of which 16 were found in the aquifer samples. Several of these aquifer species were found to be involved in denitrification but also hydrocarbon degradation (P. aeruginosa, P. putida and P. fluorescens). Some tpm sequence types allocated to P. umsongensis and P. chengduensis were found to be enriched among the tpm-harbouring bacteria, respectively, of the aquifer biofilms and waters. Reads related to Aeromonas were allocated to 11 species, but only those from A. caviae were recovered aboveground and in the aquifer samples. Some tpm sequence types of the X. axonopodis phytopathogen  were recorded in higher proportions among the tpm-harbouring bacteria of the aquifer waters than in the aboveground samples. A significant coalescence of microbial communities from an urban watershed with those of an aquifer was thus observed, and recent aquifer biofilms were found to be significantly colonized by runoff-opportunistic taxa able to use urban C sources from aboveground compartments.

Description: Landslides are the main source of sediment in most mountain ranges. Rivers then act as conveyor belts, evacuating landslide-derived sediment. Sediment dynamics are known to influence landscape evolution through interactions among landslide sediment delivery, fluvial transport and river incision into bedrock. Sediment delivery and its interaction with river incision therefore control the pace of landscape evolution and mediate relationships among tectonics, climate and erosion. Numerical landscape evolution models (LEMs) are well suited to study the interactions among these surface processes. They enable evaluation of a range of hypotheses at varying temporal and spatial scales. While many models have been used to study the dynamic interplay between tectonics, erosion and climate, the role of interactions between landslide-derived sediment and river incision has received much less attention. Here, we present HyLands, a hybrid landscape evolution model integrated within the TopoToolbox Landscape Evolution Model (TTLEM) framework. The hybrid nature of the model lies in its capacity to simulate both erosion and deposition at any place in the landscape due to fluvial bedrock incision, sediment transport, and rapid, stochastic mass wasting through landsliding. Fluvial sediment transport and bedrock incision are calculated using the recently developed Stream Power with Alluvium Conservation and Entrainment (SPACE) model. Therefore, rivers can dynamically transition from detachment-limited to transport-limited and from bedrock to bedrock–alluvial to fully alluviated states. Erosion and sediment production by landsliding are calculated using a Mohr–Coulomb stability analysis, while landslide-derived sediment is routed and deposited using a multiple-flow-direction, nonlinear deposition method. We describe and evaluate the HyLands 1.0 model using analytical solutions and observations. We first illustrate the functionality of HyLands to capture river dynamics ranging from detachment-limited to transport-limited conditions. Second, we apply the model to a portion of the Namche Barwa massif in eastern Tibet and compare simulated and observed landslide magnitude–frequency and area–volume scaling relationships. Finally, we illustrate the relevance of explicitly simulating landsliding and sediment dynamics over longer timescales for landscape evolution in general and river dynamics in particular. With HyLands we provide a new tool to understand both the long- and short-term coupling between stochastic hillslope processes, river incision and source-to-sink sediment dynamics.

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: Chemistry plays an indispensable role in investigations of the atmosphere; however, many climate models either ignore or greatly simplify atmospheric chemistry, limiting both their accuracy and their scope. We present the development and evaluation of the online global atmospheric chemical model BCC-GEOS-Chem v1.0, coupling the GEOS-Chem chemical transport model (CTM) as an atmospheric chemistry component in the Beijing Climate Center atmospheric general circulation model (BCC-AGCM). The GEOS-Chem atmospheric chemistry component includes detailed tropospheric HOx–NOx–volatile organic compounds–ozone–bromine–aerosol chemistry and online dry and wet deposition schemes. We then demonstrate the new capabilities of BCC-GEOS-Chem v1.0 relative to the base BCC-AGCM model through a 3-year (2012–2014) simulation with anthropogenic emissions from the Community Emissions Data System (CEDS) used in the Coupled Model Intercomparison Project Phase 6 (CMIP6). The model captures well the spatial distributions and seasonal variations in tropospheric ozone, with seasonal mean biases of 0.4–2.2 ppbv at 700–400 hPa compared to satellite observations and within 10 ppbv at the surface to 500 hPa compared to global ozonesonde observations. The model has larger high-ozone biases over the tropics which we attribute to an overestimate of ozone chemical production. It underestimates ozone in the upper troposphere which is likely due either to the use of a simplified stratospheric ozone scheme or to biases in estimated stratosphere–troposphere exchange dynamics. The model diagnoses the global tropospheric ozone burden, OH concentration, and methane chemical lifetime to be 336 Tg, 1.16×106 molecule cm−3, and 8.3 years, respectively, which is consistent with recent multimodel assessments. The spatiotemporal distributions of NO2, CO, SO2, CH2O, and aerosol optical depth are generally in agreement with satellite observations. The development of BCC-GEOS-Chem v1.0 represents an important step for the development of fully coupled earth system models (ESMs) in China.

Description: We report the results of a multiproxy study that combines structural analysis of a fracture–stylolite network and isotopic characterization of calcite vein cements and/or fault coating. Together with new paleopiezometric and radiometric constraints on burial evolution and deformation timing, these results provide a first-order picture of the regional fluid systems and pathways that were present during the main stages of contraction in the Tuscan Nappe and Umbria–Marche Apennine Ridge (northern Apennines). We reconstruct four steps of deformation at the scale of the belt: burial-related stylolitization, Apenninic-related layer-parallel shortening with a contraction trending NE–SW, local extension related to folding, and late-stage fold tightening under a contraction still striking NE–SW. We combine the paleopiezometric inversion of the roughness of sedimentary stylolites – that constrains the range of burial depth of strata prior to layer-parallel shortening – with burial models and U–Pb absolute dating of fault coatings in order to determine the timing of development of mesostructures. In the western part of the ridge, layer-parallel shortening started in Langhian time (∼15 Ma), and then folding started at Tortonian time (∼8 Ma); late-stage fold tightening started by the early Pliocene (∼5 Ma) and likely lasted until recent/modern extension occurred (∼3 Ma onward). The textural and geochemical (δ18O, δ13C, Δ47CO2 and 87Sr∕86Sr) study of calcite vein cements and fault coatings reveals that most of the fluids involved in the belt during deformation either are local or flowed laterally from the same reservoir. However, the western edge of the ridge recorded pulses of eastward migration of hydrothermal fluids (〉140 ∘C), driven by the tectonic contraction and by the difference in structural style of the subsurface between the eastern Tuscan Nappe and the Umbria–Marche Apennine Ridge.

Description: We developed the WRF-GC model, an online coupling of the Weather Research and Forecasting (WRF) mesoscale meteorological model and the GEOS-Chem atmospheric chemistry model, for regional atmospheric chemistry and air quality modeling. WRF and GEOS-Chem are both open-source community models. WRF-GC offers regional modellers access to the latest GEOS-Chem chemical module, which is state of the science, well documented, traceable, benchmarked, actively developed by a large international user base, and centrally managed by a dedicated support team. At the same time, WRF-GC enables GEOS-Chem users to perform high-resolution forecasts and hindcasts for any region and time of interest. WRF-GC uses unmodified copies of WRF and GEOS-Chem from their respective sources; the coupling structure allows future versions of either one of the two parent models to be integrated into WRF-GC with relative ease. Within WRF-GC, the physical and chemical state variables are managed in distributed memory and translated between WRF and GEOS-Chem by the WRF-GC coupler at runtime. We used the WRF-GC model to simulate surface PM2.5 concentrations over China during 22 to 27 January 2015 and compared the results to surface observations and the outcomes from a GEOS-Chem Classic nested-China simulation. Both models were able to reproduce the observed spatiotemporal variations of regional PM2.5, but the WRF-GC model (r=0.68, bias =29 %) reproduced the observed daily PM2.5 concentrations over eastern China better than the GEOS-Chem Classic model did (r=0.72, bias =55 %). This was because the WRF-GC simulation, nudged with surface and upper-level meteorological observations, was able to better represent the pollution meteorology during the study period. The WRF-GC model is parallelized across computational cores and scales well on massively parallel architectures. In our tests where the two models were similarly configured, the WRF-GC simulation was 3 times more efficient than the GEOS-Chem Classic nested-grid simulation due to the efficient transport algorithm and the Message Passing Interface (MPI)-based parallelization provided by the WRF software framework. WRF-GC v1.0 supports one-way coupling only, using WRF-simulated meteorological fields to drive GEOS-Chem with no chemical feedbacks. The development of two-way coupling capabilities, i.e., the ability to simulate radiative and microphysical feedbacks of chemistry to meteorology, is under way. The WRF-GC model is open source and freely available from http://wrf.geos-chem.org (last access: 10 July 2020).

Description: Mixed-method bicultural research in Aotearoa New Zealand, including the weaving of Indigenous and other knowledge, is emerging within many academic disciplines. However, mātauranga Māori (the knowledge, culture, values, and world view of the Indigenous peoples of Aotearoa New Zealand) and Te Ao Māori (the Māori world) is poorly represented within geomorphological investigations. Here, we review international efforts to include Indigenous knowledge in geologic and geomorphic studies and provide an overview of the current state of mātauranga Māori within research endeavours in Aotearoa New Zealand. We review three theoretical frameworks (i.e. methodologies) for including mātauranga Māori in research projects and three models (i.e. methods) for including Māori values within research. We identify direct benefits to geomorphology and discuss how these frameworks and models can be adapted for use with Indigenous knowledge systems outside of Aotearoa New Zealand. The aim of this review is to encourage geomorphologists around the world to engage with local Indigenous peoples to develop new approaches to geomorphic research. In Aotearoa New Zealand, we hope to inspire geomorphologists to embark on research journeys in genuine partnership with Māori that promote toitū te mātauranga – the enduring protection, promotion and respect of mātauranga Māori.

Description: Suspended sediment load in rivers is highly uncertain because sediment production and transport at catchment scale are strongly variable in space and time, and they are affected by catchment hydrology, topography, and land cover. Among the main sources of this variability are the spatially distributed nature of overland flow as an erosion driver and of surface erodibility given by soil type and vegetation cover distribution. Temporal variability mainly results from the time sequence of rainfall intensity during storms and snowmelt leading to soil saturation and overland flow. We present a new spatially distributed soil erosion and suspended sediment transport module integrated into the computationally efficient physically based hydrological model TOPKAPI-ETH, with which we investigate the effects of the two erosion drivers – precipitation and surface erodibility – on catchment sediment fluxes in a typical pre-Alpine mesoscale catchment. By conducting a series of numerical experiments, we quantify the impact of spatial variability in the two key erosion drivers on erosion–deposition patterns, sediment delivery ratio, and catchment sediment yields. Main findings are that the spatial variability in erosion drivers affects sediment yield by (i) increasing sediment production due to a spatially variable precipitation, while decreasing it due to a spatially variable surface erodibility, (ii) favouring the clustering of sediment source areas in space by surface runoff generation, and (iii) decreasing their connectivity to the river network by magnifying sediment buffers. The results highlight the importance of resolving spatial gradients controlling hydrology and sediment processes when modelling sediment dynamics at the mesoscale, in order to capture the key effects of sediment sources, buffers, and hillslope hydrological pathways in determining the sediment signal.

Description: Atmospheric mineral dust influences Earth's radiative budget, cloud formation, and lifetime; has adverse health effects; and affects air quality through the increase of regulatory PM10 concentrations, making its real-time quantification in the atmosphere of strategic importance. Only few near-real-time techniques can discriminate dust aerosol in PM10 samples and they are based on the dust chemical composition. The online determination of mineral dust using aerosol absorption photometers offers an interesting and competitive alternative but remains a difficult task to achieve. This is particularly challenging when dust is mixed with black carbon, which features a much higher mass absorption cross section. We build on previous work using filter photometers and present here for the first time a highly time-resolved online technique for quantification of mineral dust concentration by coupling a high-flow virtual impactor (VI) sampler that concentrates coarse particles with an aerosol absorption photometer (Aethalometer, model AE33). The absorption of concentrated dust particles is obtained by subtracting the absorption of the submicron (PM1) aerosol fraction from the absorption of the virtual impactor sample (VI-PM1 method). This real-time method for detecting desert dust was tested in the field for a period of 2 months (April and May 2016) at a regional background site of Cyprus, in the Eastern Mediterranean. Several intense desert mineral dust events were observed during the field campaign with dust concentration in PM10 up to 45 µg m−3. Mineral dust was present most of the time during the campaign with an average PM10 of about 8 µg m−3. Mineral dust absorption was most prominent at short wavelengths, yielding an average mass absorption cross section (MAC) of 0.24±0.01 m2 g−1 at 370 nm and an absorption Ångström exponent of 1.41±0.29. This MAC value can be used as a site-specific parameter for online determination of mineral dust concentration. The uncertainty of the proposed method is discussed by comparing and validating it with different methods.

Description: Dry–rewetting perturbations are natural disturbances in the edaphic environment and particularly in dryland cultivation areas. The interaction of this disturbance with glyphosate-based herbicides (GBHs) deserves special attention in the soil environment due to the intensification of agricultural practices and the acceleration of climate change with an intensified water cycle. The objective of this study was to assess the response of microbial communities in a soil with a long history of GBHs to a secondary imposed perturbation (a single dry–rewetting event). A factorial microcosm study was conducted to evaluate the potential conditioning effect of an acute glyphosate exposure on the response to a following dry–rewetting event. A respiratory quotient (RQ) based on an ecologically relevant substrate (p-coumaric acid) and basal respiration was used as a physiological indicator. Similarly, DNA-based analyses were considered, including quantitative PCR (qPCR) of functional sensitive microbial groups linked to cycles of carbon (Actinobacteria) and nitrogen (ammonia-oxidizing microorganisms), qPCR of total bacteria and denaturing gradient gel electrophoresis (DGGE) of ammonia-oxidizing bacteria (AOB). Significant effects of herbicide and of dry–rewetting perturbations were observed in the RQ and in the copy number of the amoA gene of AOB, respectively. However, no significant interaction was observed between them when analyzing the physiological indicator and the copy number of the evaluated genes. PCR–DGGE results were not conclusive regarding a potential effect of dry–rewetting × herbicide interaction on AOB community structure, suggesting further analysis by deep sequencing of the amoA gene. The results of this study indicate that the perturbation of an acute glyphosate exposure in a soil with a long history of this herbicide does not have a conditioning effect on the response to a subsequent dry–rewetting disturbance according to a physiological indicator or the quantified bacterial/archaeal genes. This is particularly relevant for the sustainability of soils in rainfed agriculture, where frequent exposure to GBHs along with intensification of hydrological cycles are expected to occur. Further studies considering multiple dry–rewetting disturbances and in different soil types should be conducted to simulate those conditions and to validate our results.

Description: Storm tides are a major hazard for the German North Sea coasts. For coastal protection and economic activities, planning information on the probability and magnitude of extreme storm tides and their possible future changes is important. This study focuses on the most extreme events and examines whether they could have become more severe under slightly different conditions while still remaining within physical plausibility. In the face of a limited number of observational data on very severe events, an extensive set of model data is used to extract most extreme storm tide events for locations in the German Bight, in particular Borkum and the Ems estuary. The data set includes water levels and respective atmospheric conditions from a hindcast and future climate realizations without sea level rise describing today's and possible future conditions. A number of very severe events with water levels exceeding those measured near Borkum since 1906 are identified in the data set. A possible further amplification of the highest events is investigated by simulating these events for the North Sea with different phase lags between the astronomical tide given at the open model boundaries and the wind forcing. It is found that superposition of spring tide conditions, different timing of the astronomical high water and atmospheric conditions during the highest storm event would cause an enhancement of the highest water level up to about 50 cm. The water levels of the two highest events from the data set are used to analyse the effects in the Ems estuary using a high-resolution model of the German Bight. Additionally, the influences of an extreme river runoff and of sea level rise are studied. The extreme river runoff of 1200 m3 s−1 increases the highest water levels by several decimetres in the narrow upstream part of the Ems estuary. This effect diminishes downstream. The sea level rise increases the water level in the downstream part of the Ems estuary by the amount applied at the model boundary to the North Sea. In the upstream part, its influence on the water level decreases. This study may serve as a first step towards an impact assessment for severe storm tides and towards implications for coastal zone management in times of climate change.

Description: More than 300 non-dispersive infrared (NDIR) CO2 low-cost sensors labelled as LP8 were integrated into sensor units and evaluated for the purpose of long-term operation in the Carbosense CO2 sensor network in Switzerland. Prior to deployment, all sensors were calibrated in a pressure and climate chamber and in ambient conditions co-located with a reference instrument. To investigate their long-term performance and to test different data processing strategies, 18 sensors were deployed at five locations equipped with a reference instrument after calibration. Their accuracy during 19 to 25 months deployment was between 8 and 12 ppm. This level of accuracy requires careful sensor calibration prior to deployment, continuous monitoring of the sensors, efficient data filtering, and a procedure to correct drifts and jumps in the sensor signal during operation. High relative humidity (〉 ∼85 %) impairs the LP8 measurements, and corresponding data filtering results in a significant loss during humid conditions. The LP8 sensors are not suitable for the detection of small regional gradients and long-term trends. However, with careful data processing, the sensors are able to resolve CO2 changes and differences with a magnitude larger than about 30 ppm. Thereby, the sensor can resolve the site-specific CO2 signal at most locations in Switzerland. A low-power network (LPN) using LoRaWAN allowed for reliable data transmission with low energy consumption and proved to be a key element of the Carbosense low-cost sensor network.

Description: Peatlands store ∼ 20 %–30 % of the global soil organic carbon stock and are an important source of dissolved organic carbon (DOC) for inland waters. Recent improvements for in situ optical monitoring revealed that the DOC concentration in streams draining peatlands is highly variable, showing seasonal variation and short and intense DOC concentration peaks. This study aimed to statistically determine the variables driving stream DOC concentration variations at seasonal and event scales. Two mountainous peatlands (one fen and one bog) were monitored in the French Pyrenees to capture their outlet DOC concentration variability at a high-frequency rate (30 min). Abiotic variables including precipitation, stream temperature and water level, water table depth, and peat water temperature were also monitored at high frequency and used as potential predictors to explain DOC concentration variability. Results show that at both sites DOC concentration time series can be decomposed into a seasonal baseline interrupted by many short and intense peaks of higher concentrations. The DOC concentration baseline is driven, at the seasonal scale, by peat water temperature. At the event scale, DOC concentration increases are mostly driven by a rise in the water table within the peat at both sites. Univariate linear models between DOC concentration and peat water temperature or water table increases show greater efficiency at the fen site. Water recession times were derived from water level time series using master recession curve coefficients. They vary greatly between the two sites but also within one peatland site. They partly explain the differences between DOC dynamics in the studied peatlands, including peat porewater DOC concentrations and the links between stream DOC concentration and water table rise within the peatlands. This highlights that peatland complexes are composed of a mosaic of heterogeneous peat units distinctively producing or transferring DOC to streams.

Description: Submarine groundwater discharge (SGD) is an important gateway for nutrients and pollutants from land to sea. While understanding SGD is crucial for managing nearshore ecosystems and coastal freshwater reserves, studying this discharge is complicated by its occurrence at the limit between land and sea, a dynamic environment. This practical difficulty is exacerbated by the significant spatial and temporal variability. Therefore, to capture the magnitude of SGD, a variety of techniques and measurements, applied over multiple periods, is needed. Here, we combine several geophysical methods to detect zones of fresh submarine groundwater discharge (FSGD) in the intertidal zone, upper beach, dunes, and shallow coastal area. Both terrestrial electrical-resistivity tomography (ERT; roll-along) and marine continuous resistivity profiling (CRP) are used from the shallow continental shelf up to the dunes and combined with frequency domain electromagnetic (FDEM) mapping in the intertidal zone. In particular, we apply an estimation of robust apparent electrical conductivity (rECa) from FDEM data to provide reliable lateral and vertical discrimination of FSGD zones. The study area is a very dynamic environment along the North Sea, characterized by semi-diurnal tides between 3 and 5 m. CRP is usually applied in calmer conditions, but we prove that such surveys are possible and provide additional information to primarily land-bound ERT surveying. The 2D inversion models created from ERT and CRP data clearly indicate the presence of FSGD on the lower beach or below the low-water line. This discharge originates from a potable freshwater lens below the dunes and flows underneath a thick saltwater lens, present from the dunes to the lower sandy beach, which is fully observed with ERT. Freshwater outflow intensity has increased since 1980, due to a decrease of groundwater pumping in the dunes. FDEM mapping at two different times reveals discharge at the same locations, clearly displays the lateral variation of the zone of discharge, and suggests that FSGD is stronger at the end of winter compared to the beginning of autumn. ERT, CRP, and FDEM are complementary tools in the investigation of SGD. They provide a high-resolution 3D image of the saltwater and freshwater distribution in the phreatic coastal aquifer over a relatively large area, both off- and onshore.

Description: An optimized method is presented to determine dissolved free (DFCHO) and dissolved combined carbohydrates (DCCHO) in saline matrices, such as oceanic seawater, Arctic ice core samples or brine using a combination of a desalination with electro-dialysis (ED) and high-performance anion exchange chromatography coupled to pulsed amperometric detection (HPAEC-PAD). Free neutral sugars, such as glucose and galactose, were found with 95 %–98 % recovery rates. Free amino sugars and free uronic acids were strongly depleted during ED at pH=8, but an adjustment of the pH could result in higher recoveries (58 %–59 % for amino sugars at pH=11; 45 %–49 % for uronic acids at pH=1.5). The applicability of this method for the analysis of DCCHO was evaluated with standard solutions and seawater samples compared with another established desalination method using membrane dialysis. DFCHO in field samples from different regions on Earth ranged between 11 and 118 nM and DCCHO between 260 and 1410 nM. This novel method has the potential to contribute to a better understanding of biogeochemical processes in the oceans and sea–air transfer processes of organic matter into the atmosphere in future studies.

Description: A sensor comprised of an electronic circuit and a hybrid single and dual heat pulse probe was constructed and tested along with a novel signal processing procedure to determine changes in the effective dual-probe spacing radius over the time of measurement. The circuit utilized a proportional–integral–derivative (PID) controller to control heat inputs into the soil medium in lieu of a variable resistor. The system was designed for onboard signal processing and implemented USB, RS-232, and SDI-12 interfaces for machine-to-machine (M2M) exchange of data, thereby enabling heat inputs to be adjusted to soil conditions and data availability shortly after the time of experiment. Signal processing was introduced to provide a simplified single-probe model to determine thermal conductivity instead of reliance on late-time logarithmic curve fitting. Homomorphic and derivative filters were used with a dual-probe model to detect changes in the effective probe spacing radius over the time of experiment to compensate for physical changes in radius as well as model and experimental error. Theoretical constraints were developed for an efficient inverse of the exponential integral on an embedded system. Application of the signal processing to experiments on sand and peat improved the estimates of soil water content and bulk density compared to methods of curve fitting nominally used for heat pulse probe experiments. Applications of the technology may be especially useful for soil and environmental conditions under which effective changes in probe spacing radius need to be detected and compensated for over the time of experiment.

Description: The Copernicus Polar Ice and Snow Topography Altimeter (CRISTAL) mission is one of six high-priority candidate missions (HPCMs) under consideration by the European Commission to enlarge the Copernicus Space Component. Together, the high-priority candidate missions fill gaps in the measurement capability of the existing Copernicus Space Component to address emerging and urgent user requirements in relation to monitoring anthropogenic CO2 emissions, polar environments, and land surfaces. The ambition is to enlarge the Copernicus Space Component with the high-priority candidate missions in the mid-2020s to provide enhanced continuity of services in synergy with the next generation of the existing Copernicus Sentinel missions. CRISTAL will carry a dual-frequency synthetic-aperture radar altimeter as its primary payload for measuring surface height and a passive microwave radiometer to support atmospheric corrections and surface-type classification. The altimeter will have interferometric capabilities at Ku-band for improved ground resolution and a second (non-interferometric) Ka-band frequency to provide information on snow layer properties. This paper outlines the user consultations that have supported expansion of the Copernicus Space Component to include the high-priority candidate missions, describes the primary and secondary objectives of the CRISTAL mission, identifies the key contributions the CRISTAL mission will make, and presents a concept – as far as it is already defined – for the mission payload.

Description: It is widely accepted that the atmospheric boundary layer is drastically under-sampled in the vertical dimension. In recent years, the commercial availability of ground-based remote sensors combined with the widespread use of small, weather-sensing uncrewed aerial systems (WxUAS) has opened up many opportunities to fill this measurement gap. In July 2018, the University of Oklahoma (OU) deployed a state-of-the-art WxUAS, dubbed the CopterSonde, and the Collaborative Lower Atmospheric Mobile Profiling System (CLAMPS) in the San Luis Valley in south-central Colorado. Additionally, these systems were deployed to the Kessler Atmospheric and Ecological Field Station (KAEFS) in October 2018. The colocation of these various systems provided ample opportunity to compare and contrast kinematic and thermodynamic observations from different methodologies of boundary layer profiling, namely WxUAS, remote sensing, and the traditional in situ radiosonde. In this study, temperature, dew point temperature, wind speed, and wind direction from these platforms are compared statistically with data from the two campaigns. Moreover, we present select instances from the dataset to highlight differences between the measurement techniques. This analysis highlights strengths and weaknesses of planetary boundary layer profiling and helps lay the groundwork for developing highly adaptable systems that integrate remote and in situ profiling techniques.

Description: Für die Verortung der deutschen akademischen Humangeographie im Geflecht zwischen wissenschaftlicher Praxis, theoretischen Propositionen und Lehralltag war der Geographentag in Kiel im Jahre 1969 ein Meilenstein. Anlässlich der Erinnerung an den Ort und die Debatten vor 50 Jahren wurde wiederum Kiel im Jahre 2019 zu einem Ort der Reflektion. Der hier vorliegende Beitrag versucht in einer bewusst persönlich formulierten Art die Impulse, die von „Kiel 1969“ ausgingen, im universitären Alltag des Geographischen Institutes der TU München in den 1980er Jahren zu verorten und hierdurch gewissermaßen zu relativieren. Hierdurch entsteht ein differenziertes Bild von richtungsweisenden Veränderungen und verharrenden Strukturen, welche ineinander verwoben die damals überregional bekannte Münchener Sozialgeographie charakterisierten – und für die deutschsprachige Humangeographie über die speziell Münchner Zustände hinaus bezeichnend waren.

Description: An improved two-sphere integration (TSI) technique has been developed to quantify black carbon (BC) concentrations in the atmosphere and seasonal snow. The major advantage of this system is that it combines two distinct integrated spheres to reduce the scattering effect due to light-absorbing particles and thus provides accurate determinations of total light absorption from BC collected on Nuclepore filters. The TSI technique can be calibrated using a series of 15 filter samples of standard fullerene soot. This technique quantifies the mass of BC by separating the spectrally resolved total light absorption into BC and non-BC fractions. To assess the accuracy of the improved system, an empirical procedure for measuring BC concentrations with a two-step thermal–optical method is also applied. Laboratory results indicate that the BC concentrations determined using the TSI technique and theoretical calculations are well correlated (R2=0.99), whereas the thermal–optical method underestimates BC concentrations by 35 %–45 % compared to that measured by the TSI technique. Assessments of the two methods for atmospheric and snow samples revealed excellent agreement, with least-squares regression lines with slopes of 1.72 (r2=0.67) and 0.84 (r2=0.93), respectively. However, the TSI technique is more accurate in quantifications of BC concentrations in both the atmosphere and seasonal snow, with an overall lower uncertainty. Using the improved TSI technique, we find that light absorption at a wavelength of 550 nm due to BC plays a dominant role relative to non-BC light absorption in both the atmosphere (62.76 %–91.84 % of total light absorption) and seasonal snow (43.11 %–88.56 %) over northern China.

Description: Quantitative precipitation estimation with commercial microwave links (CMLs) is a technique developed to supplement weather radar and rain gauge observations. It is exploiting the relation between the attenuation of CML signal levels and the integrated rain rate along a CML path. The opportunistic nature of this method requires a sophisticated data processing using robust methods. In this study we focus on the processing step of rain event detection in the signal level time series of the CMLs, which we treat as a binary classification problem. This processing step is particularly challenging, because even when there is no rain, the signal level can show large fluctuations similar to that during rainy periods. False classifications can have a high impact on falsely estimated rainfall amounts. We analyze the performance of a convolutional neural network (CNN), which is trained to detect rainfall-specific attenuation patterns in CML signal levels, using data from 3904 CMLs in Germany. The CNN consists of a feature extraction and a classification part with, in total, 20 layers of neurons and 1.4×105 trainable parameters. With a structure inspired by the visual cortex of mammals, CNNs use local connections of neurons to recognize patterns independent of their location in the time series. We test the CNN's ability to recognize attenuation patterns from CMLs and time periods outside the training data. Our CNN is trained on 4 months of data from 800 randomly selected CMLs and validated on 2 different months of data, once for all CMLs and once for the 3104 CMLs not included in the training. No CMLs are excluded from the analysis. As a reference data set, we use the gauge-adjusted radar product RADOLAN-RW provided by the German meteorological service (DWD). The model predictions and the reference data are compared on an hourly basis. Model performance is compared to a state-of-the-art reference method, which uses the rolling standard deviation of the CML signal level time series as a detection criteria. Our results show that within the analyzed period of April to September 2018, the CNN generalizes well to the validation CMLs and time periods. A receiver operating characteristic (ROC) analysis shows that the CNN is outperforming the reference method, detecting on average 76 % of all rainy and 97 % of all nonrainy periods. From all periods with a reference rain rate larger than 0.6 mm h−1, more than 90 % was detected. We also show that the improved event detection leads to a significant reduction of falsely estimated rainfall by up to 51 %. At the same time, the quality of the correctly estimated rainfall is kept at the same level in regards to the Pearson correlation with the radar rainfall. In conclusion, we find that CNNs are a robust and promising tool to detect rainfall-induced attenuation patterns in CML signal levels from a large CML data set covering all of Germany.

Description: Per capita arable land is decreasing due to the rapidly increasing population, and fresh water is becoming scarce and more expensive. Therefore, farmers should continue to use technology and innovative solutions to improve efficiency, save input costs, and optimise environmental resources (such as water). In the case study presented in this paper, the Global Navigation Satellite System interferometric reflectometry (GNSS-IR) technique was used to monitor soil moisture during 66 d, from 3 December 2018 to 6 February 2019, in the installations of the Cajamar Centre of Experiences, Paiporta, Valencia, Spain. Two main objectives were pursued. The first was the extension of the technique to a multi-constellation solution using GPS, GLONASS, and GALILEO satellites, and the second was to test whether mass-market sensors could be used for this technique. Both objectives were achieved. At the same time that the GNSS observations were made, soil samples taken at 5 cm depth were used for soil moisture determination to establish a reference data set. Based on a comparison with that reference data set, all GNSS solutions, including the three constellations and the two sensors (geodetic and mass market), were highly correlated, with a correlation coefficient between 0.7 and 0.85.