ALBERT

Description: We present the result of the third Marine Ice Sheet Model Intercomparison Project, MISMIP+. MISMIP+ is intended to be a benchmark for ice-flow models which include fast sliding marine ice streams and floating ice shelves and in particular a treatment of viscous stress that is sufficient to model buttressing, where upstream ice flow is restrained by a downstream ice shelf. A set of idealized experiments first tests that models are able to maintain a steady state with the grounding line located on a retrograde slope due to buttressing and then explore scenarios where a reduction in that buttressing causes ice stream acceleration, thinning, and grounding line retreat. The majority of participating models passed the first test and then produced similar responses to the loss of buttressing. We find that the most important distinction between models in this particular type of simulation is in the treatment of sliding at the bed, with other distinctions – notably the difference between the simpler and more complete treatments of englacial stress but also the differences between numerical methods – taking a secondary role.

Description: The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

Description: Satellite retrievals of column-averaged dry-air mole fractions of carbon dioxide (CO2) and methane (CH4), denoted XCO2 and XCH4, respectively, have been used in recent years to obtain information on natural and anthropogenic sources and sinks and for other applications such as comparisons with climate models. Here we present new data sets based on merging several individual satellite data products in order to generate consistent long-term climate data records (CDRs) of these two Essential Climate Variables (ECVs). These ECV CDRs, which cover the time period 2003–2018, have been generated using an ensemble of data products from the satellite sensors SCIAMACHY/ENVISAT and TANSO-FTS/GOSAT and (for XCO2) for the first time also including data from the Orbiting Carbon Observatory 2 (OCO-2) satellite. Two types of products have been generated: (i) Level 2 (L2) products generated with the latest version of the ensemble median algorithm (EMMA) and (ii) Level 3 (L3) products obtained by gridding the corresponding L2 EMMA products to obtain a monthly 5∘×5∘ data product in Obs4MIPs (Observations for Model Intercomparisons Project) format. The L2 products consist of daily NetCDF (Network Common Data Form) files, which contain in addition to the main parameters, i.e., XCO2 or XCH4, corresponding uncertainty estimates for random and potential systematic uncertainties and the averaging kernel for each single (quality-filtered) satellite observation. We describe the algorithms used to generate these data products and present quality assessment results based on comparisons with Total Carbon Column Observing Network (TCCON) ground-based retrievals. We found that the XCO2 Level 2 data set at the TCCON validation sites can be characterized by the following figures of merit (the corresponding values for the Level 3 product are listed in brackets) – single-observation random error (1σ): 1.29 ppm (monthly: 1.18 ppm); global bias: 0.20 ppm (0.18 ppm); and spatiotemporal bias or relative accuracy (1σ): 0.66 ppm (0.70 ppm). The corresponding values for the XCH4 products are single-observation random error (1σ): 17.4 ppb (monthly: 8.7 ppb); global bias: −2.0 ppb (−2.9 ppb); and spatiotemporal bias (1σ): 5.0 ppb (4.9 ppb). It has also been found that the data products exhibit very good long-term stability as no significant long-term bias trend has been identified. The new data sets have also been used to derive annual XCO2 and XCH4 growth rates, which are in reasonable to good agreement with growth rates from the National Oceanic and Atmospheric Administration (NOAA) based on marine surface observations. The presented ECV data sets are available (from early 2020 onwards) via the Climate Data Store (CDS, https://cds.climate.copernicus.eu/, last access: 10 January 2020) of the Copernicus Climate Change Service (C3S, https://climate.copernicus.eu/, last access: 10 January 2020).

Description: It is widely accepted that orbital variations are responsible for the generation of glacial cycles during the late Pleistocene. However, the relative contributions of the orbital forcing compared to CO2 variations and other feedback mechanisms causing the waxing and waning of ice sheets have not been fully understood. Testing theories of ice ages beyond statistical inferences, requires numerical modeling experiments that capture key features of glacial transitions. Here, we focus on the glacial buildup from Marine Isotope Stage (MIS) 7 to 6 covering the period from 240 to 170 ka (ka: thousand years before present). This transition from interglacial to glacial conditions includes one of the fastest Pleistocene glaciation–deglaciation events, which occurred during MIS 7e–7d–7c (236–218 ka). Using a newly developed three-dimensional coupled atmosphere–ocean–vegetation–ice sheet model (LOVECLIP), we simulate the transient evolution of Northern Hemisphere and Southern Hemisphere ice sheets during the MIS 7–6 period in response to orbital and greenhouse gas forcing. For a range of model parameters, the simulations capture the evolution of global ice volume well within the range of reconstructions. Over the MIS 7–6 period, it is demonstrated that glacial inceptions are more sensitive to orbital variations, whereas terminations from deep glacial conditions need both orbital and greenhouse gas forcings to work in unison. For some parameter values, the coupled model also exhibits a critical North American ice sheet configuration, beyond which a stationary-wave–ice-sheet topography feedback can trigger an unabated and unrealistic ice sheet growth. The strong parameter sensitivity found in this study originates from the fact that delicate mass imbalances, as well as errors, are integrated during a transient simulation for thousands of years. This poses a general challenge for transient coupled climate–ice sheet modeling, with such coupled paleo-simulations providing opportunities to constrain such parameters.

Description: Solar trackers are often used by spectrometers to measure atmospheric trace gas concentrations using direct sun spectroscopy. The ideal solar tracker should be sufficiently accurate, highly reliable, and with a longevity that exceeds the lifetime of the spectrometer that it serves. It should also be affordable, easy to use, and not too complex should maintenance be required. In this paper we present a design that fulfils these requirements using some simple innovations. Our altitude–azimuth design features a custom coaxial power transformer, enabling continuous 360∘ azimuth rotation. This increases reliability and avoids the need to reverse the tracker each day. In polar regions, measurements can continue uninterrupted through the summer polar day. Tracking accuracy is enhanced using a simple optical feedback technique that adjusts error offset variables while monitoring the edges of a focused solar image with four photodiodes. Control electronics are modular, and our software is written in Python, running as a web server on a recycled laptop with a Linux operating system. Over a period of 11 years we have assembled four such trackers. These are in use at Lauder (45∘ S), New Zealand, and Arrival Heights (78∘ S), Antarctica, achieving a history of good reliability even in polar conditions. Tracker accuracy is analysed regularly and can routinely produce a pointing accuracy of 0.02∘.

Description: The use of a boundary-layer parameterization of buttressing and ice flux across grounding lines in a two-dimensional ice-sheet model is improved by allowing general orientations of the grounding line. This and another modification to the model's grounding-line parameterization are assessed in three settings: rectangular fjord-like domains – the third Marine Ice Sheet Model Intercomparison Project (MISMIP+) and Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d) – and future simulations of West Antarctic ice retreat under Representative Concentration Pathway (RCP)8.5-based climates. The new modifications are found to have significant effects on the fjord-like results, which are now within the envelopes of other models in the MISMIP+ and MISMIP3d intercomparisons. In contrast, the modifications have little effect on West Antarctic retreat, presumably because dynamics in the wider major Antarctic basins are adequately represented by the model's previous simpler one-dimensional formulation. As future grounding lines retreat across very deep bedrock topography in the West Antarctic simulations, buttressing is weak and deviatoric stress measures exceed the ice yield stress, implying that structural failure at these grounding lines would occur. We suggest that these grounding-line quantities should be examined in similar projections by other ice models to better assess the potential for future structural failure.

Description: This work presents the latest release (v9.0) of the University of Leicester GOSAT Proxy XCH4 dataset. Since the launch of the GOSAT satellite in 2009, these data have been produced by the UK National Centre for Earth Observation (NCEO) as part of the ESA Greenhouse Gas Climate Change Initiative (GHG-CCI) and Copernicus Climate Change Services (C3S) projects. With now over a decade of observations, we outline the many scientific studies achieved using past versions of these data in order to highlight how this latest version may be used in the future. We describe in detail how the data are generated, providing information and statistics for the entire processing chain from the L1B spectral data through to the final quality-filtered column-averaged dry-air mole fraction (XCH4) data. We show that out of the 19.5 million observations made between April 2009 and December 2019, we determine that 7.3 million of these are sufficiently cloud-free (37.6 %) to process further and ultimately obtain 4.6 million (23.5 %) high-quality XCH4 observations. We separate these totals by observation mode (land and ocean sun glint) and by month, to provide data users with the expected data coverage, including highlighting periods with reduced observations due to instrumental issues. We perform extensive validation of the data against the Total Carbon Column Observing Network (TCCON), comparing to ground-based observations at 22 locations worldwide. We find excellent agreement with TCCON, with an overall correlation coefficient of 0.92 for the 88 345 co-located measurements. The single-measurement precision is found to be 13.72 ppb, and an overall global bias of 9.06 ppb is determined and removed from the Proxy XCH4 data. Additionally, we validate the separate components of the Proxy (namely the modelled XCO2 and the XCH4∕XCO2 ratio) and find these to be in excellent agreement with TCCON. In order to show the utility of the data for future studies, we compare against simulated XCH4 from the TM5 model. We find a high degree of consistency between the model and observations throughout both space and time. When focusing on specific regions, we find average differences ranging from just 3.9 to 15.4 ppb. We find the phase and magnitude of the seasonal cycle to be in excellent agreement, with an average correlation coefficient of 0.93 and a mean seasonal cycle amplitude difference across all regions of −0.84 ppb. These data are available at https://doi.org/10.5285/18ef8247f52a4cb6a14013f8235cc1eb (Parker and Boesch, 2020).

Description: We describe the change in operational instrument for the routine measurement of column-averaged dry-air mole fraction of several greenhouse gases (denoted Xgas) at the Lauder Total Carbon Column Observing Network (TCCON) site and the steps taken to demonstrate comparability between the two observation systems following a systematic methodology. Further, we intercompare retrieved Xgas values during an intensive intercomparison period during October and November 2018, when both instruments were performing optimally, and on subsequent, less frequent occasions. The average difference between the two observing systems was found to be well below the expected level of uncertainty for TCCON retrievals for all compared species. In the case of XCO2 the average difference was 0.0264±0.0465 % (0.11±0.19 µmol mol−1).