ALBERT

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 terrestrial subsurface is the largest source of freshwater globally. The organic carbon contained within it and processes controlling its concentration remain largely unknown. The global median concentration of dissolved organic carbon (DOC) in groundwater is low compared to surface waters, suggesting significant processing in the subsurface. Yet the processes that remove this DOC in groundwater are not fully understood. The purpose of this study was to investigate the different sources and processes influencing DOC in a shallow anoxic coastal aquifer. Uniquely, this study combines liquid chromatography organic carbon detection with organic (δ13CDOC) carbon isotope geochemical analyses to fingerprint the various DOC sources that influence the concentration, carbon isotopic composition, and character with respect to distance from surface water sources, depth below surface, and inferred groundwater residence time (using 3H activities) in groundwater. It was found that the average groundwater DOC concentration was 5 times higher (5 mg L−1) than the global median concentration and that the concentration doubled with depth at our site, but the chromatographic character did not change significantly. The anoxic saturated conditions of the aquifer limited the rate of organic matter processing, leading to enhanced preservation and storage of the DOC sources from peats and palaeosols contained within the aquifer. All groundwater samples were more aromatic for their molecular weight in comparison to other lakes, rivers and surface marine samples studied. The destabilization or changes in hydrology, whether by anthropogenic or natural processes, could lead to the flux of up to 10 times more unreacted organic carbon from this coastal aquifer compared to deeper inland aquifers.

Description: One major challenge in the study of late-Quaternary extinctions (LQEs) is providing better estimates of past megafauna abundance. To show how megaherbivore population size varied before and after the last extinctions in interior Alaska, we use both a database of radiocarbon-dated bone remains (spanning 25–0 ka) and spores of the obligate dung fungus, Sporormiella, recovered from radiocarbon-dated lake-sediment cores (spanning 17–0 ka). Bone fossils show that the last stage of LQEs in the region occurred at about 13 ka ago, but the number of megaherbivore bones remains high into the Holocene. Sporormiella abundance also remains high into the Holocene and does not decrease with major vegetation changes recorded by arboreal pollen percentages. At two sites, the interpretation of Sporormiella was enhanced by additional dung fungal spore types (e.g., Sordaria). In contrast to many sites where the last stage of LQEs is marked by a sharp decline in Sporormiella abundance, in interior Alaska our results indicate the continuance of megaherbivore abundance, albeit with a major taxonomic turnover (including Mammuthus and Equus extinction) from predominantly grazing to browsing dietary guilds. This new and robust evidence implies that regional LQEs were not systematically associated with crashes of overall megaherbivore abundance.

Description: Multiple observation data sets – Interagency Monitoring of Protected Visual Environments (IMPROVE) network data, the Automated Smoke Detection and Tracking Algorithm (ASDTA), Hazard Mapping System (HMS) smoke plume shapefiles and aircraft acetonitrile (CH3CN) measurements from the NOAA Southeast Nexus (SENEX) field campaign – are used to evaluate the HMS–BlueSky–SMOKE (Sparse Matrix Operator Kernel Emission)–CMAQ (Community Multi-scale Air Quality Model) fire emissions and smoke plume prediction system. A similar configuration is used in the US National Air Quality Forecasting Capability (NAQFC). The system was found to capture most of the observed fire signals. Usage of HMS-detected fire hotspots and smoke plume information was valuable for deriving both fire emissions and forecast evaluation. This study also identified that the operational NAQFC did not include fire contributions through lateral boundary conditions, resulting in significant simulation uncertainties. In this study we focused both on system evaluation and evaluation methods. We discussed how to use observational data correctly to retrieve fire signals and synergistically use multiple data sets. We also addressed the limitations of each of the observation data sets and evaluation methods.

Description: Snow densification stores water in alpine regions and transforms snow into ice on the surface of glaciers. Despite its importance in determining snow-water equivalent and glacier-induced sea level rise, we still lack a complete understanding of the physical mechanisms underlying snow compaction. In essence, compaction is a rheological process, where the rheology evolves with depth due to variation in temperature, pressure, humidity, and meltwater. The rheology of snow compaction can be determined in a few ways, for example, through empirical investigations (e.g., Herron and Langway, 1980), by microstructural considerations (e.g., Alley, 1987), or by measuring the rheology directly, which is the approach we take here. Using a French-press or cafetière-à-piston compression stage, Wang and Baker (2013) compressed numerous snow samples of different densities. Here we derive a mixture theory for compaction and airflow through the porous snow to compare against these experimental data. We find that a plastic compaction law explains experimental results. Taking standard forms for the permeability and effective pressure as functions of the porosity, we show that this compaction mode persists for a range of densities and overburden loads. These findings suggest that measuring compaction in the lab is a promising direction for determining the rheology of snow through its many stages of densification.

Description: We reanalyzed mass balance records at Taku and Lemon Creek Glaciers to better understand the relative roles of hypsometry, local climate and dynamics as mass balance drivers. Over the 1946–2018 period, the cumulative mass balances diverged. Tidewater Taku Glacier advanced and gained mass at an average rate of +0.25 ± 0.28 m w.e. a–1, contrasting with retreat and mass loss of −0.60 ± 0.15 m w.e. a−1 at land-terminating Lemon Creek Glacier. The uniform influence of regional climate is demonstrated by strong correlations among annual mass balance and climate data. Regional warming trends forced similar statistically significant decreases in surface mass balance after 1989: −0.83 m w.e. a–1 at Taku Glacier and −0.81 m w.e. a–1 at Lemon Creek Glacier. Divergence in cumulative mass balance arises from differences in glacier hypsometry and local climate. Since 2013 negative mass balance and glacier-wide thinning prevailed at Taku Glacier. These changes initiated terminus retreat, which could increase dramatically if calving begins. The future mass balance trajectory of Taku Glacier hinges on dynamics, likely ending the historic dichotomy between these glaciers.

Description: Recent successes in passive remote sensing of far-red solar-induced chlorophyll fluorescence (SIF) have spurred the development and integration of canopy-level fluorescence models in global terrestrial biosphere models (TBMs) for climate and carbon cycle research. The interaction of fluorescence with photochemistry at the leaf and canopy scales provides opportunities to diagnose and constrain model simulations of photosynthesis and related processes, through direct comparison to and assimilation of tower, airborne, and satellite data. TBMs describe key processes related to the absorption of sunlight, leaf-level fluorescence emission, scattering, and reabsorption throughout the canopy. Here, we analyze simulations from an ensemble of process-based TBM–SIF models (SiB3 – Simple Biosphere Model, SiB4, CLM4.5 – Community Land Model, CLM5.0, BETHY – Biosphere Energy Transfer Hydrology, ORCHIDEE – Organizing Carbon and Hydrology In Dynamic Ecosystems, and BEPS – Boreal Ecosystems Productivity Simulator) and the SCOPE (Soil Canopy Observation Photosynthesis Energy) canopy radiation and vegetation model at a subalpine evergreen needleleaf forest near Niwot Ridge, Colorado. These models are forced with local meteorology and analyzed against tower-based continuous far-red SIF and gross-primary-productivity-partitioned (GPP) eddy covariance data at diurnal and synoptic scales during the growing season (July–August 2017). Our primary objective is to summarize the site-level state of the art in TBM–SIF modeling over a relatively short time period (summer) when light, canopy structure, and pigments are similar, setting the stage for regional- to global-scale analyses. We find that these models are generally well constrained in simulating photosynthetic yield but show strongly divergent patterns in the simulation of absorbed photosynthetic active radiation (PAR), absolute GPP and fluorescence, quantum yields, and light response at the leaf and canopy scales. This study highlights the need for mechanistic modeling of nonphotochemical quenching in stressed and unstressed environments and improved the representation of light absorption (APAR), distribution of light across sunlit and shaded leaves, and radiative transfer from the leaf to the canopy scale.

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: Characterizing the temporal uncertainty in palaeoclimate records is crucial for analysing past climate change, correlating climate events between records, assessing climate periodicities, identifying potential triggers and evaluating climate model simulations. The first global compilation of speleothem isotope records by the SISAL (Speleothem Isotope Synthesis and Analysis) working group showed that age model uncertainties are not systematically reported in the published literature, and these are only available for a limited number of records (ca. 15 %, n=107/691). To improve the usefulness of the SISAL database, we have (i) improved the database's spatio-temporal coverage and (ii) created new chronologies using seven different approaches for age–depth modelling. We have applied these alternative chronologies to the records from the first version of the SISAL database (SISALv1) and to new records compiled since the release of SISALv1. This paper documents the necessary changes in the structure of the SISAL database to accommodate the inclusion of the new age models and their uncertainties as well as the expansion of the database to include new records and the quality-control measures applied. This paper also documents the age–depth model approaches used to calculate the new chronologies. The updated version of the SISAL database (SISALv2) contains isotopic data from 691 speleothem records from 294 cave sites and new age–depth models, including age–depth temporal uncertainties for 512 speleothems. SISALv2 is available at https://doi.org/10.17864/1947.256 (Comas-Bru et al., 2020a).

Description: With semiconductor technology gradually approaching its physical and thermal limits, recent supercomputers have adopted major architectural changes to continue increasing the performance through more power-efficient heterogeneous many-core systems. Examples include Sunway TaihuLight that has four management processing elements (MPEs) and 256 computing processing elements (CPEs) inside one processor and Summit that has two central processing units (CPUs) and six graphics processing units (GPUs) inside one node. Meanwhile, current high-resolution Earth system models that desperately require more computing power generally consist of millions of lines of legacy code developed for traditional homogeneous multicore processors and cannot automatically benefit from the advancement of supercomputer hardware. As a result, refactoring and optimizing the legacy models for new architectures become key challenges along the road of taking advantage of greener and faster supercomputers, providing better support for the global climate research community and contributing to the long-lasting societal task of addressing long-term climate change. This article reports the efforts of a large group in the International Laboratory for High-Resolution Earth System Prediction (iHESP) that was established by the cooperation of Qingdao Pilot National Laboratory for Marine Science and Technology (QNLM), Texas A&M University (TAMU), and the National Center for Atmospheric Research (NCAR), with the goal of enabling highly efficient simulations of the high-resolution (25 km atmosphere and 10 km ocean) Community Earth System Model (CESM-HR) on Sunway TaihuLight. The refactoring and optimizing efforts have improved the simulation speed of CESM-HR from 1 SYPD (simulation years per day) to 3.4 SYPD (with output disabled) and supported several hundred years of pre-industrial control simulations. With further strategies on deeper refactoring and optimizing for remaining computing hotspots, as well as redesigning architecture-oriented algorithms, we expect an equivalent or even better efficiency to be gained on the new platform than traditional homogeneous CPU platforms. The refactoring and optimizing processes detailed in this paper on the Sunway system should have implications for similar efforts on other heterogeneous many-core systems such as GPU-based high-performance computing (HPC) systems.