Explore the vast range of titles available.

Climate models often exhibit spurious long-term changes independent of either internal variability or changes to external forcing. Such changes, referred to as model “drift,” may distort the estimate of forced change in transient climate simulations. The importance of drift is examined in comparison to historical trends over recent decades in the Coupled Model Intercomparison Project (CMIP). Comparison based on a selection of metrics suggests a significant overall reduction in the magnitude of drift from phase 3 of CMIP (CMIP3) to phase 5 of CMIP (CMIP5). The direction of both ocean and atmospheric drift is systematically biased in some models introducing statistically significant drift in globally averaged metrics. Nevertheless, for most models globally averaged drift remains weak compared to the associated forced trends and is often smaller than the difference between trends derived from different ensemble members or the error introduced by the aliasing of natural variability. An exception to this is metrics that include the deep ocean (e.g., steric sea level) where drift can dominate in forced simulations. In such circumstances drift must be corrected for using information from concurrent control experiments. Many CMIP5 models now include ocean biogeochemistry. Like physical models, biogeochemical models generally undergo long spinup integrations to minimize drift. Nevertheless, based on a limited subset of models, it is found that drift is an important consideration and must be accounted for. For properties or regions where drift is important, the drift correction method must be carefully considered. The use of a drift estimate based on the full control time series is recommended to minimize the contamination of the drift estimate by internal variability.

How much Antarctic ice shelf basal melt rates can increase in response to global warming remains an open question. Here we describe the response of the Southern Ocean and ice shelf cavities to an abrupt change to high-end atmospheric conditions plausible by the late 23rd century under the SSP5-8.5 scenario. To achieve this objective, we first present and evaluate a new 0.25∘ global configuration of the NEMO (Nucleus for European Modelling of the Ocean NEMO System Team, 2019) ocean and sea ice model. Our present-day simulations demonstrate good agreement with observational data for key variables such as temperature, salinity, and ice shelf melt rates, despite the remaining difficulties to simulate the interannual variability in the Amundsen Sea. The ocean response to the high-end atmospheric perturbation includes a strengthening and extension of the Ross and Weddell gyres and a quasi-disappearance of sea ice, with a subsequent decrease in production of High Salinity Shelf Water and increased intrusion of warmer water onto the continental shelves favoured by changes in baroclinic currents at the shelf break. We propose to classify the perturbed continental shelf as a “warm–fresh shelf”. This induces a substantial increase in ice shelf basal melt rates, particularly in the coldest seas, with a total basal mass loss rising from 1180 to 15 700 Gt yr−1 and an Antarctica averaged melt rate increasing from 0.8 to 10.6 m yr−1. In the perturbed simulation, most ice shelves around Antarctica experience conditions that are currently found in the Amundsen Sea, while the Amundsen Sea warms by 2 ∘C. These idealised projections can be used as a base to calibrate basal melt parameterisations used in long-term ice sheet projections.

Climate model projections have previously been used to compute ice shelf basal melt rates in ice sheet models, but the strategies employed – e.g., ocean input, parameterization, calibration technique, and corrections – have varied widely and are often ad hoc. Here, a methodology is proposed for the calculation of circum-Antarctic basal melt rates for floating ice, based on climate models, that is suitable for ISMIP6, the Ice Sheet Model Intercomparison Project for CMIP6 (6th Coupled Model Intercomparison Project). The past and future evolution of ocean temperature and salinity is derived from a climate model by estimating anomalies with respect to the modern day, which are added to a present-day climatology constructed from existing observational datasets. Temperature and salinity are extrapolated to any position potentially occupied by a simulated ice shelf. A simple formulation is proposed for a basal melt parameterization in ISMIP6, constrained by the observed temperature climatology, with a quadratic dependency on either the nonlocal or local thermal forcing. Two calibration methods are proposed: (1) based on the mean Antarctic melt rate (MeanAnt) and (2) based on melt rates near Pine Island's deep grounding line (PIGL). Future Antarctic mean melt rates are an order of magnitude greater in PIGL than in MeanAnt. The PIGL calibration and the local parameterization result in more realistic melt rates near grounding lines. PIGL is also more consistent with observations of interannual melt rate variability underneath Pine Island and Dotson ice shelves. This work stresses the need for more physics and less calibration in the parameterizations and for more observations of hydrographic properties and melt rates at interannual and decadal timescales.

Changes in precipitation over the European Alps are investigated with the regional climate model MAR (Modèle Atmosphérique Régional) applied with a 7 km resolution over the period 1903–2010 using the reanalysis ERA-20C as forcing. A comparison with several observational datasets demonstrates that the model is able to reproduce the climatology as well as both the interannual variability and the seasonal cycle of precipitation over the European Alps. The relatively high resolution allows us to estimate precipitation at high elevations. The vertical gradient of precipitation simulated by MAR over the European Alps reaches 33% km−1 (1.21 mm d−1 km−1) in summer and 38 % km−1 (1.15 mm d−1 km−1) in winter, on average, over 1971–2008 and shows a large spatial variability. A significant (p value < 0.05) increase in mean winter precipitation is simulated in the northwestern Alps over 1903–2010, with changes typically reaching 20 % to 40 % per century. This increase is mainly explained by a stronger simple daily intensity index (SDII) and is associated with less-frequent but longer wet spells. A general drying is found in summer over the same period, exceeding 20 % to 30 % per century in the western plains and 40 % to 50 % per century in the southern plains surrounding the Alps but remaining much smaller (<10 %) and not significant above 1500 m a.s.l. Below this level, the summer drying is explained by a reduction in the number of wet days, reaching 20 % per century over the northwestern part of the Alps and 30 % to 50 % per century in the southern part of the Alps. It is associated with shorter but more-frequent wet spells. The centennial trends are modulated over the last decades, with the drying occurring in the plains in winter also affecting high-altitude areas during this season and with a positive trend of autumn precipitation occurring only over the last decades all over the Alps. Maximum daily precipitation index (Rx1day) takes its highest values in autumn in both the western and the eastern parts of the southern Alps, locally reaching 50 to 70 mm d−1 on average over 1903–2010. Centennial maxima up to 250 to 300 mm d−1 are simulated in the southern Alps, in France and Italy, as well as in the Ticino valley in Switzerland. Over 1903–2010, seasonal Rx1day shows a general and significant increase at the annual timescale and also during the four seasons, reaching local values between 20 % and 40 % per century over large parts of the Alps and the Apennines. Trends of Rx1day are significant (p value < 0.05) only when considering long time series, typically 50 to 80 years depending on the area considered. Some of these trends are nonetheless significant when computed over 1970–2010, suggesting a recent acceleration of the increase in extreme precipitation, whereas earlier periods with strong precipitation also occurred, in particular during the 1950s and 1960s.

Oceanic melting beneath ice shelves is the main driver of the current mass loss of the Antarctic ice sheet and is mostly parameterised in stand-alone ice-sheet modelling. Parameterisations are crude representations of reality, and their response to ocean warming has not been compared to 3-D ocean–ice-sheet coupled models. Here, we assess various melting parameterisations ranging from simple scalings with far-field thermal driving to emulators of box and plume models, using a new coupling framework combining the ocean model NEMO and the ice-sheet model Elmer/Ice. We define six idealised one-century scenarios for the far-field ocean ranging from cold to warm, and representative of potential futures for typical Antarctic ice shelves. The scenarios are used to constrain an idealised geometry of the Pine Island glacier representative of a relatively small cavity. Melt rates and sea-level contributions obtained with the parameterised stand-alone ice-sheet model are compared to the coupled model results. The plume parameterisations give good results for cold scenarios but fail and underestimate sea level contribution by tens of percent for warm(ing) scenarios, which may be improved by adapting its empirical scaling. The box parameterisation with five boxes compares fairly well to the coupled results for almost all scenarios, but further work is needed to grasp the correct number of boxes. For simple scalings, the comparison to the coupled framework shows that a quadratic as opposed to linear dependency on thermal forcing is required. In addition, the quadratic dependency is improved when melting depends on both local and non-local, i.e. averaged over the ice shelf, thermal forcing. The results of both the box and the two quadratic parameterisations fall within or close to the coupled model uncertainty. All parameterisations overestimate melting for thin ice shelves while underestimating melting in deep water near the grounding line. Further work is therefore needed to assess the validity of these melting parameteriations in more realistic set-ups.

The representation of the El Niño–Southern Oscillation (ENSO) under historical forcing and future projections is analyzed in 34 models from the Coupled Model Intercomparison Project phase 5 (CMIP5). Most models realistically simulate the observed intensity and location of maximum sea surface temperature (SST) anomalies during ENSO events. However, there exist systematic biases in the westward extent of ENSO-related SST anomalies, driven by unrealistic westward displacement and enhancement of the equatorial wind stress in the western Pacific. Almost all CMIP5 models capture the observed asymmetry in magnitude between the warm and cold events (i.e., El Niños are stronger than La Niñas) and between the two types of El Niños: that is, cold tongue (CT) El Niños are stronger than warm pool (WP) El Niños. However, most models fail to reproduce the asymmetry between the two types of La Niñas, with CT stronger than WP events, which is opposite to observations. Most models capture the observed peak in ENSO amplitude around December; however, the seasonal evolution of ENSO has a large range of behavior across the models. The CMIP5 models generally reproduce the duration of CT El Niños but have biases in the evolution of the other types of events. The evolution of WP El Niños suggests that the decay of this event occurs through heat content discharge in the models rather than the advection of SST via anomalous zonal currents, as seems to occur in observations. No consistent changes are seen across the models in the location and magnitude of maximum SST anomalies, frequency, or temporal evolution of these events in a warmer world.

The ice sheet model intercomparison project for CMIP6 (ISMIP6) effort brings together the ice sheet and climate modeling communities to gain understanding of the ice sheet contribution to sea level rise. ISMIP6 conducts stand-alone ice sheet experiments that use space- and time-varying forcing derived from atmosphere–ocean coupled global climate models (AOGCMs) to reflect plausible trajectories for climate projections. The goal of this study is to recommend a subset of CMIP5 AOGCMs (three core and three targeted) to produce forcing for ISMIP6 stand-alone ice sheet simulations, based on (i) their representation of current climate near Antarctica and Greenland relative to observations and (ii) their ability to sample a diversity of projected atmosphere and ocean changes over the 21st century. The selection is performed separately for Greenland and Antarctica. Model evaluation over the historical period focuses on variables used to generate ice sheet forcing. For stage (i), we combine metrics of atmosphere and surface ocean state (annual- and seasonal-mean variables over large spatial domains) with metrics of time-mean subsurface ocean temperature biases averaged over sectors of the continental shelf. For stage (ii), we maximize the diversity of climate projections among the best-performing models. Model selection is also constrained by technical limitations, such as availability of required data from RCP2.6 and RCP8.5 projections. The selected top three CMIP5 climate models are CCSM4, MIROC-ESM-CHEM, and NorESM1-M for Antarctica and HadGEM2-ES, MIROC5, and NorESM1-M for Greenland. This model selection was designed specifically for ISMIP6 but can be adapted for other applications.

At the interface between the ocean and the ice shelf base, in the framework of the shear-controlled melt parameterisation, the ice melts due to combined actions of temperature, salinity and friction velocity. In the NEMO ocean model, the friction velocity is usually computed based on a constant drag coefficient and an ocean velocity averaged vertically within a distance from the ice, which is often referred to as the Losch layer. Instead, in this study, we use a logarithmic approach, where a constant hydrographic roughness length detetermines the drag coefficient through the law of the wall and the horizontal current speed is sampled in the first wet cell. The aim is to reduce the vertical resolution dependency, to homogeneise the sampling of horizontal current speed between the thermodynamic and the dynamic drag equation and to enable the use of a variable drag coefficient based on the subgrid-scale (or unresolved) ice shelf basal topography. The motivation behind a variable drag based on the topography comes from observations showing that regions with rough topographic features such as crevasses or basal melt channels experience more melts than flat ones. We compare different experiments in a configuration of the Amundsen Sea at 1/12°. We find that our approach is less sensitive (6 % melt rates difference) to a coarser vertical resolution, such as the one used in global Earth System Models, than the Losch layer approach (22 % melt rates difference). We also find that it succeeds in reproducing higher melt rates in rougher regions while keeping total ice shelf melt rate within the observed range. Finally, to assess the effect of increasing ice shelf damage, we tested the sensitivity of a higher hydrographic roughness length. If the roughness of all the ice shelf grid points were to increase to the highest value currently observed, the overall ice shelf melting would increase by 16 %. This suggests the possibility of a positive feedback in which more melting leads to more ice damage and increased roughness, in turn increasing melt rates.

A large spread exists in both Indian and Australian average monsoon rainfall and in their interannual variations diagnosed from various observational and reanalysis products. While the multi model mean monsoon rainfall from 59 models taking part in the Coupled Model Intercomparison Project (CMIP3 and CMIP5) fall within the observational uncertainty, considerable model spread exists. Rainfall seasonality is consistent across observations and reanalyses, but most CMIP models produce either a too peaked or a too flat seasonal cycle, with CMIP5 models generally performing better than CMIP3. Considering all North-Australia rainfall, most models reproduce the observed Australian monsoon-El Niño Southern Oscillation (ENSO) teleconnection, with the strength of the relationship dependent on the strength of the simulated ENSO. However, over the Maritime Continent, the simulated monsoon-ENSO connection is generally weaker than observed, depending on the ability of each model to realistically reproduce the ENSO signature in the Warm Pool region. A large part of this bias comes from the contribution of Papua, where moisture convergence seems to be particularly affected by this SST bias. The Indian summer monsoon-ENSO relationship is affected by overly persistent ENSO events in many CMIP models. Despite significant wind anomalies in the Indian Ocean related to Indian Ocean Dipole (IOD) events, the monsoon-IOD relationship remains relatively weak both in the observations and in the CMIP models. Based on model fidelity in reproducing realistic monsoon characteristics and ENSO teleconnections, we objectively select 12 “best” models to analyze projections in the rcp8.5 scenario. Eleven of these models are from the CMIP5 ensemble. In India and Australia, most of these models produce 5–20 % more monsoon rainfall over the second half of the twentieth century than during the late nineteenth century. By contrast, there is no clear model consensus over the Maritime Continent.

The disintegration of the ice shelves along the Antarctic Peninsula have spurred much discussion on the various processes leading to their eventual dramatic collapse, but without a consensus on an atmospheric forcing that could connect these processes. Here, using an atmospheric river detection algorithm along with a regional climate model and satellite observations, we show that the most intense atmospheric rivers induce extremes in temperature, surface melt, sea-ice disintegration, or large swells that destabilize the ice shelves with 40% probability. This was observed during the collapses of the Larsen A and B ice shelves during the summers of 1995 and 2002 respectively. Overall, 60% of calving events from 2000–2020 were triggered by atmospheric rivers. The loss of the buttressing effect from these ice shelves leads to further continental ice loss and subsequent sea-level rise. Under future warming projections, the Larsen C ice shelf will be at-risk from the same processes.