Abstract: Surface melting on the Antarctic Ice Sheet has been monitored by satellite microwave radiometry for over 40 years. Despite this long perspective, our understanding of the microwave emission from wet snow is still limited, preventing the full exploitation of these observations to study supraglacial hydrology. Using the Snow Microwave Radiative Transfer (SMRT) model, this study investigates the sensitivity of microwave brightness temperature to snow liquid water content at frequencies from 1.4 to 37 GHz. We first determine the snowpack properties for eight selected coastal sites by retrieving profiles of density, grain size and ice layers from microwave observations when the snowpack is dry during wintertime. Second, a series of brightness temperature simulations is run with added water. The results show that (i) a small quantity of liquid water (≈0.5 kg m−2) can be detected, but the actual quantity cannot be retrieved out of the full range of possible water quantities; (ii) the detection of a buried wet layer is possible up to a maximum depth of 1 to 6 m depending on the frequency (6–37 GHz) and on the snow properties (grain size, density) at each site; (iii) surface ponds and water-saturated areas may prevent melt detection, but the current coverage of these waterbodies in the large satellite field of view is presently too small in Antarctica to have noticeable effects; and (iv) at 1.4 GHz, while the simulations are less reliable, we found a weaker sensitivity to liquid water and the maximal depth of detection is relatively shallow (<10 m) compared to the typical radiation penetration depth in dry firn (≈1000 m) at this low frequency. These numerical results pave the way for the development of improved multi-frequency algorithms to detect melt intensity and the depth of liquid water below the surface in the Antarctic snowpack.

Abstract: The energy budget of Arctic sea ice is strongly affected by the snow cover. Intensive sampling of snow properties was conducted near Qikiqtarjuaq in Baffin Bay on typical landfast sea ice during two melt seasons in 2015 and 2016. The sampling included stratigraphy, vertical profiles of snow specific surface area (SSA), density and irradiance, and spectral albedo (300–1100 nm). Both years featured four main phases: (I) dry snow cover, (II) surface melting, (III) ripe snowpack, and (IV) melt pond formation. Each phase was characterized by distinctive physical and optical properties. A high SSA value of 49.3 m2 kg−1 was measured during phase I on surface wind slabs together with a corresponding broadband albedo (300–3000 nm) of 0.87. Phase II was marked by alternating episodes of surface melting, which dramatically decreased the SSA below 3 m2 kg−1, and episodes of snowfall re-establishing pre-melt conditions. Albedo was highly time-variable, with minimum broadband values around 0.70. In phase III, continued melting led to a fully ripe snowpack composed of clustered rounded grains. Albedo began to decrease in the visible as snow thickness decreased but remained steady at longer wavelengths. Moreover, significant spatial variability appeared for the first time following snow depth heterogeneity. Spectral albedo was simulated by radiative transfer using measured SSA and density vertical profiles and estimated impurity contents based on limited measurements. Simulations were most of the time within 1 % of measurements in the visible and within 2 % in the infrared. Simulations allowed the calculations of albedo and of the spectral flux at the snow–ice interface. These showed that photosynthetically active radiation fluxes at the bottom of the snowpack durably exceeded 5 W m−2 (∼9.2 µmol m−2 s−1) only when the snowpack thickness started to decrease at the end of phase II.

Abstract:

Abstract. Arctic landscapes are covered in snow for at least 6 months of the year. The energy balance of the snow cover plays a key role in these environments, influencing the surface albedo, the thermal regime of the permafrost, and other factors. Our goal is to quantify all major heat fluxes above, within, and below a low-Arctic snowpack at a shrub tundra site on the east coast of Hudson Bay in eastern Canada. The study is based on observations from a flux tower that uses the eddy covariance approach and from profiles of temperature and thermal conductivity in the snow and soil. Additionally, we compared the observations with simulations produced using the Crocus snow model. We found that radiative losses due to negative longwave radiation are mostly counterbalanced by the sensible heat flux, whereas the latent heat flux is minimal. At the snow surface, the heat flux into the snow is similar in magnitude to the sensible heat flux. Because the snow cover stores very little heat, the majority of the upward heat flux in the snow is used to cool the soil. Overall, the model was able to reproduce the observed energy balance, but due to the effects of atmospheric stratification, it showed some deficiencies when simulating turbulent heat fluxes at an hourly timescale.

Abstract: The Regional Ocean Modeling System (ROMS), including an ice shelf component, has been applied on a circum-Antarctic domain to derive estimates of ice shelf basal melting. Significant improvements made compared to previous models of this scale are the inclusion of tides and a horizontal spatial resolution of 2 km, which is sufficient to resolve on-shelf heat transport by bathymetric troughs and eddy-scale circulation. We run the model with ocean–atmosphere–sea ice conditions from the year 2007 to represent nominal present-day climate. We force the ocean surface with buoyancy fluxes derived from sea ice concentration observations and wind stress from ERA-Interim atmospheric reanalysis. Boundary conditions are derived from the ECCO2 ocean state estimate; tides are incorporated as sea surface height and barotropic currents at the open boundary. We evaluate model results using satellite-derived estimates of ice shelf melting and established compilations of ocean hydrography. The Whole Antarctic Ocean Model (WAOM v1.0) qualitatively captures the broad scale difference between warm and cold regimes as well as many of the known characteristics of regional ice–ocean interaction. We identify a cold bias for some warm-water ice shelves and a lack of high-salinity shelf water (HSSW) formation. We conclude that further calibration and development of our approach are justified. At its current state, the model is ideal for addressing specific, process-oriented questions, e.g. related to tide-driven ice shelf melting at large scales.

Abstract: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2021, EFOS increased by 5.1 % relative to 2020, with fossil emissions at 10.1 ± 0.5 GtC yr−1 (9.9 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.1 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 10.9 ± 0.8 GtC yr−1 (40.0 ± 2.9 GtCO2). Also, for 2021, GATM was 5.2 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 3.5 ± 0.9 GtC yr−1, with a BIM of −0.6 GtC yr−1 (i.e. the total estimated sources were too low or sinks were too high). The global atmospheric CO2 concentration averaged over 2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest an increase in EFOS relative to 2021 of +1.0 % (0.1 % to 1.9 %) globally and atmospheric CO2 concentration reaching 417.2 ppm, more than 50 % above pre-industrial levels (around 278 ppm). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2021, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extratropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set. The data presented in this work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).