Search

The article presents a study which focuses on the use of Challenging Minisatellite Payload (CHAMP) Global Positioning System (GPS) radio occultation (RO) data in the investigation of temperature trends in Antarctica. The seasonal temperature trends are examined at 100 hectopascals and 500 hectopascals pressure levels. Variations in the temperature at spatial and temporal domains are highlighted.

Regional climate models (RCMs) and reanalysis datasets provide valuable information for assessing the vulnerability of ice shelves to collapse over Antarctica, which is important for future global sea level rise estimates. Within this context, this paper examines variability in snowfall, near-surface air temperature and melt across products from the Met Office Unified Model (MetUM), Regional Atmospheric Climate Model (RACMO) and Modèle Atmosphérique Régional (MAR) RCMs, as well as the ERA-Interim and ERA5 reanalysis datasets. Seasonal and trend decomposition using LOESS (STL) is applied to split the monthly time series at each model grid cell into trend, seasonal and residual components. Significant systematic differences between outputs are shown for all variables in the mean and in the seasonal and residual standard deviations, occurring at both large and fine spatial scales across Antarctica. Results imply that differences in the atmospheric dynamics, parametrisation, tuning and surface schemes between models together contribute more significantly to large-scale variability than differences in the driving data, resolution, domain specification, ice sheet mask, digital elevation model and boundary conditions. Despite significant systematic differences, high temporal correlations are found for snowfall and near-surface air temperature across all products at fine spatial scales. For melt, only moderate correlation exists at fine spatial scales between different RCMs and low correlation between RCM and reanalysis outputs. Root mean square deviations (RMSDs) between all outputs in the monthly time series for each variable are shown to be significant at fine spatial scales relative to the magnitude of annual deviations. Correcting for systematic differences results in significant reductions in RMSDs, suggesting the importance of observations and further development of bias-correction techniques. [ABSTRACT FROM AUTHOR]

HydroGFD3 (Hydrological Global Forcing Data) is a data set of bias-adjusted reanalysis data for daily precipitation and minimum, mean, and maximum temperature. It is mainly intended for large-scale hydrological modelling but is also suitable for other impact modelling. The data set has an almost global land area coverage, excluding the Antarctic continent and small islands, at a horizontal resolution of 0.25 ∘ , i.e. about 25 km. It is available for the complete ERA5 reanalysis time period, currently 1979 until 5 d ago. This period will be extended back to 1950 once the back catalogue of ERA5 is available. The historical period is adjusted using global gridded observational data sets, and to acquire real-time data, a collection of several reference data sets is used. Consistency in time is attempted by relying on a background climatology and only making use of anomalies from the different data sets. Precipitation is adjusted for mean bias as well as the number of wet days in a month. The latter is relying on a calibrated statistical method with input only of the monthly precipitation anomaly such that no additional input data about the number of wet days are necessary. The daily mean temperature is adjusted toward the monthly mean of the observations and applied to 1 h time steps of the ERA5 reanalysis. Daily mean, minimum, and maximum temperature are then calculated. The performance of the HydroGFD3 data set is on par with other similar products, although there are significant differences in different parts of the globe, especially where observations are uncertain. Further, HydroGFD3 tends to have higher precipitation extremes, partly due to its higher spatial resolution. In this paper, we present the methodology, evaluation results, and how to access the data set at 10.5281/zenodo.3871707. [ABSTRACT FROM AUTHOR]

Atmospheric rivers (ARs) are the primary mechanism for transporting water vapor from low latitudes to polar regions, playing a significant role as drivers of extreme weather, such as heavy precipitation and heat waves in both the Arctic and Antarctica. With the rapidly growing interest in polar ARs during the past decade, it is imperative to establish an objective framework to quantify the strength and impact of these ARs for both scientific research and practical application. The AR scale introduced by Ralph et al. (2019) ranks ARs based on the duration of AR conditions and the intensity. However, the thresholds of integrated water vapor transport (IVT) used to rank ARs are selected based on the IVT climatology at middle latitudes. These thresholds are insufficient for polar regions due to the substantially lower temperature and moisture content. In this study, we analyze the IVT climatology in polar regions, focusing on the coasts of Antarctica and Greenland. Then we introduce an extended version of the AR scale tuned to polar regions by adding lower IVT thresholds of 100, 150, and 200 kg m-1 s-1 to the standard AR scale, which starts at 250 kg m-1 s-1. The polar AR scale is utilized to examine AR frequency, seasonality, trends, and associated precipitation and surface melt over the Antarctic and Greenland coasts. The polar AR scale better characterizes the strength and impacts of ARs in the Antarctic and Arctic regions, and has the potential to enhance communications across observation, research, and forecasts for polar regions. [ABSTRACT FROM AUTHOR]

The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 timescale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al. (2013) focuses on the interval 120-800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception, respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on d15N data employed as the Datice background scenario, the AICC2012 presents a slightly improved timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the seesaw mechanism, with maximum differences of about 600 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. It is expected that the future contribution of both other long ice core records and other types of chronological constraints to the Datice tool will lead to further refinements in the ice core chronologies beyond the AICC2012 chronology. For the time being however, we recommend that AICC2012 be used as the preferred chronology for the Vostok, EDC, EDML and TALDICE ice core records, both over the last glacial cycle (this study), and beyond (following Bazin et al., 2013). The ages for NGRIP in AICC2012 are virtually identical to those of GICC05 for the last 60.2 ka, whereas the ages beyond are independent of those in GICC05modelext (as in the construction of AICC2012, the GICC05modelext was included only via the background scenarios and not as age markers). As such, where issues of phasing between Antarctic records included in AICC2012 and NGRIP are involved, the NGRIP ages in AICC2012 should therefore be taken to avoid introducing false offsets. However for issues involving only Greenland ice cores, there is not yet a strong basis to recommend superseding GICC05modelext as the recommended age scale for Greenland ice cores. [ABSTRACT FROM AUTHOR]

Eighteen General Circulation Models (GCMs) are compared to reference data for the present, the Mid- Holocene (MH) and the Last Glacial Maximum (LGM) for the Antarctic region. The climatology produced by a regional climate model is taken as a reference climate for the present. GCM results for the past are compared to ice-core data. The goal of this study is to find the best GCM that can be used to drive an ice sheet model that simulates the evolution of the Antarctic Ice Sheet. Because temperature and precipitation are the most important climate variables when modelling the evolution of an ice sheet, these two variables are considered in this paper. This is done by ranking the models according to how well their output corresponds with the references. In general, present-day temperature is simulated well, but precipitation is overestimated compared to the reference data. Another finding is that model biases play an important role in simulating the past, as they are often larger than the change in temperature or precipitation between the past and the present. Considering the results for the present-day as well as for the MH and the LGM, the best performing models are HadCM3 and MIROC 3.2.2. [ABSTRACT FROM AUTHOR]

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves. [ABSTRACT FROM AUTHOR]

Water stable isotope records from ice cores (δ18 O and δ D) are a critical tool for constraining long-term temperature variability at high latitudes. However, precipitation in Antarctica consists of semi-continuous small events and intermittent extreme events. In regions of high accumulation, this can bias ice core records towards recording the synoptic climate conditions present during extreme precipitation events. In this study we utilise a combination of ice core data, reanalysis products, and models to understand how precipitation intermittency impacts the temperature records preserved in an ice core from Mount Brown South in East Antarctica. Extreme precipitation events represent only the largest 10 % of all precipitation events, but they account for 52 % of the total annual snowfall at this site, leading to an overrepresentation of these events in the ice core record. Extreme precipitation events are associated with high-pressure systems in the mid-latitudes that cause increased transport of warm and moist air from the southern Indian Ocean to the ice core site. Warm temperatures associated with these events result in a +4.8 ∘ C warm bias in the mean annual temperature when weighted by daily precipitation, and water isotopes in the Mount Brown South ice core are shown to be significantly correlated with local temperature when this precipitation-induced temperature bias is included. The Mount Brown South water isotope record spans more than 1000 years and will provide a valuable regional reconstruction of long-term temperature and hydroclimate variability in the data-sparse southern Indian Ocean region. [ABSTRACT FROM AUTHOR]

Atmospheric rivers (ARs) transport large amounts of moisture from the mid- to high-latitudes and they are a primary driver of the most extreme snowfall events, along with surface melting, in Antarctica. In this study, we characterize the climatology and surface impacts of ARs on West Antarctica, focusing on the Amundsen Sea Embayment and Marie Byrd Land. First, we develop a climatology of ARs in this region, using an Antarctic-specific AR detection tool combined with the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) and the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5) atmospheric reanalyses. We find that while ARs are infrequent (occurring 3 % of the time), they cause intense precipitation in short periods of time and account for 11 % of the annual surface accumulation. They are driven by the coupling of a blocking high over the Antarctic Peninsula with a low-pressure system known as the Amundsen Sea Low. Next, we use observations from automatic weather stations on Thwaites Eastern Ice Shelf with the firn model SNOWPACK and interferometric reflectometry (IR) to examine a case study of three ARs that made landfall in rapid succession from 2 to 8 February 2020, known as an AR family event. While accumulation dominates the surface impacts of the event on Thwaites Eastern Ice Shelf (> 100 kgm-2 or millimeters water equivalent), we find small amounts of surface melt as well (< 5 kgm-2). The results presented here enable us to quantify the past impacts of ARs on West Antarctica's surface mass balance (SMB) and characterize their interannual variability and trends, enabling a better assessment of future AR-driven changes in the SMB. [ABSTRACT FROM AUTHOR]

This study presents airborne in situ and satellite remote sensing climatologies of cirrus clouds and humidity. The climatologies serve as a guide to the properties of cirrus clouds, with the new in situ database providing detailed insights into boreal midlatitudes and the tropics, while the satellite-borne data set offers a global overview. To this end, an extensive, quality-checked data archive, the Cirrus Guide II in situ database, is created from airborne in situ measurements during 150 flights in 24 campaigns. The archive contains meteorological parameters, ice water content (IWC), ice crystal number concentration (Nice), ice crystal mean mass radius (Rice), relative humidity with respect to ice (RH ice), and water vapor mixing ratio (H2O) for each of the flights. Depending on the parameter, the database has been extended by about a factor of 5–10 compared to earlier studies. As one result of our investigation, we show that the medians of Nice , Rice , and RH ice have distinct patterns in the IWC– T parameter space. Lookup tables of these variables as functions of IWC and T can be used to improve global model cirrus representation and remote sensing retrieval methods. Another outcome of our investigation is that across all latitudes, the thicker liquid-origin cirrus predominate at lower altitudes, while at higher altitudes the thinner in situ-origin cirrus prevail. Further, examination of the radiative characteristics of in situ-origin and liquid-origin cirrus shows that the in situ-origin cirrus only slightly warm the atmosphere, while liquid-origin cirrus have a strong cooling effect. An important step in completing the Cirrus Guide II is the provision of the global cirrus Nice climatology, derived by means of the retrieval algorithm DARDAR-Nice from 10 years of cirrus remote sensing observations from satellite. The in situ measurement database has been used to evaluate and improve the satellite observations. We found that the global median Nice from satellite observations is almost 2 times higher than the in situ median and increases slightly with decreasing temperature. Nice medians of the most frequently occurring cirrus sorted by geographical regions are highest in the tropics, followed by austral and boreal midlatitudes, Antarctica, and the Arctic. Since the satellite climatologies enclose the entire spatial and temporal Nicepage12570 occurrence, we could deduce that half of the cirrus are located in the lowest, warmest (224–242 K) cirrus layer and contain a significant amount of liquid-origin cirrus. A specific highlight of the study is the in situ observations of cirrus and humidity in the Asian monsoon anticyclone and the comparison to the surrounding tropics. In the convectively very active Asian monsoon, peak values of Nice and IWC of 30 cm-3 and 1000 ppmv are detected around the cold point tropopause (CPT). Above the CPT, ice particles that are convectively injected can locally add a significant amount of water available for exchange with the stratosphere. We found IWCs of up to 8 ppmv in the Asian monsoon in comparison to only 2 ppmv in the surrounding tropics. Also, the highest RH ice values (120 %–150 %) inside of clouds and in clear sky are observed around and above the CPT. We attribute this to the high H2O mixing ratios (typically 3–5 ppmv) observed in the Asian monsoon compared to 1.5 to 3 ppmv found in the tropics. Above the CPT, supersaturations of 10 %–20 % are observed in regions of weak convective activity and up to about 50 % in the Asian monsoon. This implies that the water available for transport into the stratosphere might be higher than the expected saturation value. [ABSTRACT FROM AUTHOR]

An increasing number of Southern Ocean models now include Antarctic ice-shelf cavities, and simulate thermodynamics at the ice-shelf/ocean interface. This adds another level of complexity to Southern Ocean simulations, as ice shelves interact directly with the ocean and indirectly with sea ice. Here, we present the first model intercomparison and evaluation of present-day ocean/sea-ice/iceshelf interactions, as simulated by two models: a circumpolar Antarctic configuration of MetROMS (ROMS: Regional Ocean Modelling System coupled to CICE: Community Ice CodE) and the global model FESOM (Finite Element Sea-ice Ocean Model), where the latter is run at two different levels of horizontal resolution. From a circumpolar Antarctic perspective, we compare and evaluate simulated ice-shelf basal melting and sub-ice-shelf circulation, as well as sea-ice properties and Southern Ocean water mass characteristics as they influence the sub-ice-shelf processes. Despite their differing numerical methods, the two models produce broadly similar results and share similar biases in many cases. Both models reproduce many key features of observations but struggle to reproduce others, such as the high melt rates observed in the small warm-cavity ice shelves of the Amundsen and Bellingshausen seas. Several differences in model design show a particular influence on the simulations. For example, FESOM's greater topographic smoothing can alter the geometry of some ice-shelf cavities enough to affect their melt rates; this improves at higher resolution, since less smoothing is required. In the interior Southern Ocean, the vertical coordinate system affects the degree of water mass erosion due to spurious diapycnal mixing, with MetROMS' terrain-following coordinate leading to more erosion than FESOM's z coordinate. Finally, increased horizontal resolution in FESOM leads to higher basal melt rates for small ice shelves, through a combination of stronger circulation and small-scale intrusions of warm water from offshore. [ABSTRACT FROM AUTHOR]

We use the BISICLES adaptive mesh ice sheet model to carry out one, two, and three century simulations of the fast-flowing ice streams of the West Antarctic Ice Sheet, deploying sub-kilometer resolution around the grounding line since coarser resolution results in substantial underestimation of the response. Each of the simulations begins with a geometry and velocity close to present-day observations, and evolves according to variation in meteoric ice accumulation rates and oceanic ice shelf melt rates. Future changes in accumulation and melt rates range from no change, through anomalies computed by atmosphere and ocean models driven by the E1 and A1B emissions scenarios, to spatially uniform melt rate anomalies that remove most of the ice shelves over a few centuries. We find that variation in the resulting ice dynamics is dominated by the choice of initial conditions and ice shelf melt rate and mesh resolution, although ice accumulation affects the net change in volume above flotation to a similar degree. Given sufficient melt rates, we compute grounding line retreat over hundreds of kilometers in every major ice stream, but the ocean models do not predict such melt rates outside of the Amundsen Sea Embayment until after 2100. Within the Amundsen Sea Embayment the largest single source of variability is the onset of sustained retreat in Thwaites Glacier, which can triple the rate of eustatic sea level rise. [ABSTRACT FROM AUTHOR]

Using the original setup described in Gallée et al. (2013), the MAR regional climate model including a coupled snowpack/aeolian snow transport parameterization, was run at a fine spatial (5 km horizontal and 2m vertical) resolution over 1 summer month in coastal Adélie Land. Different types of feedback were taken into account in MAR including drag partitioning caused by surface roughness elements. Model outputs are compared with observations made at two coastal locations, D17 and D47, situated respectively 10 and 100 km inland. Wind speed was correctly simulated with positive values of the Nash test (0.60 for D17 and 0.37 for D47) but wind velocities above 10 m s-1 were underestimated at both D17 and D47; at D47, the model consistently underestimated wind velocity by 2 m s-1. Aeolian snow transport events were correctly reproduced with the right timing and a good temporal resolution at both locations except when the maximum particle height was less than 1 m. The threshold friction velocity, evaluated only at D17 for a 7-day period without snowfall, was overestimated. The simulated aeolian snow mass fluxes between 0 and 2m at D47 displayed the same variations but were underestimated compared to the second-generation FlowCapt™ values, as was the simulated relative humidity at 2m above the surface. As a result, MAR underestimated the total aeolian horizontal snow transport for the first 2 m above the ground by a factor of 10 compared to estimations by the second-generation FlowCapt™. The simulation was significantly improved at D47 if a 1-order decrease in the magnitude of z0 was accounted for, but agreement with observations was reduced at D17. Our results suggest that z0 may vary regionally depending on snowpack properties, which are involved in different types of feedback between aeolian transport of snow and z0. [ABSTRACT FROM AUTHOR]

The historical simulations of sea ice during 1979 to 2005 by the Coupled Model Intercomparison Project Phase 5 (CMIP5) are compared with satellite observations, Global Ice-Ocean Modeling and Assimilation System (GIOMAS) output data and Pan-Arctic Ice Ocean Modeling and Assimilation System (PIOMAS) output data in this study. Forty-nine models, almost all of the CMIP5 climate models and earth system models with historical simulation, are used. For the Antarctic, multi-model ensemble mean (MME) results can give good climatology of sea ice extent (SIE), but the linear trend is incorrect. The linear trend of satellite-observed Antarctic SIE is 1.29 (±0.57)×105 km2 decade–1; only about 1/7 CMIP5 models show increasing trends, and the linear trend of CMIP5 MME is negative with the value of –3.36 (±0.15)×105 km2 decade–1. For the Arctic, both climatology and linear trend are better reproduced. Sea ice volume (SIV) is also evaluated in this study, and this is a first attempt to evaluate the SIV in all CMIP5 models. Compared with the GIOMAS and PIOMAS data, the SIV values in both the Antarctic and the Arctic are too small, especially for the Antarctic in spring and winter. The GIOMAS Antarctic SIV in September is 19.1×103 km3, while the corresponding Antarctic SIV of CMIP5 MME is 13.0×103 km3 (almost 32% less). The Arctic SIV of CMIP5 in April is 27.1×103 km3, which is also less than that from PIOMAS SIV (29.5×103 km3). This means that the sea ice thickness simulated in CMIP5 is too thin, although the SIE is fairly well simulated. [ABSTRACT FROM AUTHOR]

Concentrations of oxygen (O2) and other dissolved gases in the oceanic mixed layer are often used to calculate air-sea gas exchange fluxes. The mixed layer depth (Zmix) may be defined using criteria based on temperature or density differences to a reference depth near the ocean surface. However, temperature criteria fail in regions with strong haloclines such as the Southern Ocean where heat, freshwater and momentum fluxes interact to establish mixed layers. Moreover, the time scales of air-sea exchange differ for gases and heat, so that Zmix defined using oxygen may be different than Zmix defined using temperature or density. Here, we propose to define an O2-based mixed layer depth, Zmix (O2), as the depth where the relative difference between the O2 concentration and a reference value at a depth equivalent to 10 dbar equals 0.5%. This definition was established by analysis of O2 profiles from the Bellingshausen Sea (west of the Antarctic Peninsula) and corroborated by visual inspection. Comparisons of Zmix (O2) with Zmix based on potential temperature differences, i.e., Zmix (0.2 &deg:C) and Zmix (0.5 &deg:C), and potential density differences, i.e., Zmix (0.03 kgm-3) and Zmix (0.125 kgm-3), showed that Zmix (O2) closely follows Zmix (0.03 kgm-3). Further comparisons with published Zmix climatologies and Zmixderived from World Ocean Atlas 2005 data were also performed. To establish Zmixfor use with biological production estimates in the absence of O2 profiles, we suggest using Zmix (0.03 kgm-3), which is also the basis for the climatology by de Boyer Montégut et al. (2004). [ABSTRACT FROM AUTHOR]

The article presents a study on a firn densification model that simulates Antarctic firn density profiles as well as its temporal evolution and surface height. It mentions the model's utilization of an ameliorated firn densification expression that is adjusted to fit depth-density measurements. It also notes its steady-state and time-dependent model versions which respectively simulate strong variability in firn layer thickness and operate in Antarctic places with various climate conditions.

The Sub-Millimetre Radiometer (SMR) on board the Odin satellite, launched in February 2001, observes thermal emissions of stratospheric nitric acid (HNO3) originating from the Earth limb in a band centred at 544.6 GHz. Height-resolved measurements of the global distribution of nitric acid in the stratosphere were performed approximately on two observation days per week. An HNO3 climatology based on more than 7 years of observations from August 2001 to April 2009 covering the vertical range between typically ~19 and 45 km (~1.5-60 hPa or ~500-1800K in terms of potential temperature) was created. The study highlights the spatial and seasonal variation of nitric acid in the stratosphere, characterised by a pronounced seasonal cycle at middle and high latitudes with maxima during late fall and minima during spring, strong denitrification in the lower stratosphere of the Antarctic polar vortex during winter (the irreversible removal of NOy by the sedimentation of cloud particles containing HNO3), as well as large quantities of HNO3 formed every winter at high-latitudes in the middle and upper stratosphere. A strong inter-annual variability is observed in particular at high latitudes. A comparison with a stratospheric HNO3 climatology, based on over 7 years of UARS/MLS (Upper Atmosphere Research Satellite/Microwave Limb Sounder) measurements from the 1990s, shows good consistency and agreement of the main morphological features in the potential temperature range ~465 to ~960 K, if the different characteristics of the data sets such as the better altitude resolution of Odin/SMR as well as the slightly different altitude ranges are considered. Odin/SMR reaches higher up and UARS/MLS lower down in the stratosphere. An overview from 1991 to 2009 of stratospheric nitric acid is provided (with a short gap between 1998 and 2001), if the global measurements of both experiments are taken together. [ABSTRACT FROM AUTHOR]