ALBERT

Description: Previous studies show accelerations of West Antarctic glaciers, implying that basal melt rates of these glaciers were previously small and increased in the middle of the 20th century. This enhanced melting is a likely source of the observed Ross Sea (RS) freshening, but its long-term impact on the Southern Ocean hydrography has not been well investigated. Here, we conduct coupled sea ice–ice shelf–ocean simulations with different levels of ice shelf melting from West Antarctic glaciers. Freshening of RS shelf and bottom water is simulated with enhanced West Antarctic ice shelf melting, while no significant changes in shelf water properties are simulated when West Antarctic ice shelf melting is small. We further show that the freshening caused by glacial meltwater from ice shelves in the Amundsen and Bellingshausen seas can propagate further downstream along the East Antarctic coast into the Weddell Sea. The freshening signal propagates onto the RS continental shelf within a year of model simulation, while it takes roughly 5–10 and 10–15 years to propagate into the region off Cape Darnley and into the Weddell Sea, respectively. This advection of freshening modulates the shelf water properties and possibly impacts the production of Antarctic Bottom Water if the enhanced melting of West Antarctic ice shelves continues for a longer period.

Description: The Regional Antarctic ice and Global Ocean (RAnGO) model has been developed to study the interaction between the world ocean and the Antarctic ice sheet. The coupled model is based on a global implementation of the Finite Element Sea-ice Ocean Model (FESOM) with a mesh refinement in the Southern Ocean, particularly in its marginal seas and in the sub-ice-shelf cavities. The cryosphere is represented by a regional setup of the ice flow model RIMBAY comprising the Filchner–Ronne Ice Shelf and the grounded ice in its catchment area up to the ice divides. At the base of the RIMBAY ice shelf, melt rates from FESOM's ice-shelf component are supplied. RIMBAY returns ice thickness and the position of the grounding line. The ocean model uses a pre-computed mesh to allow for an easy adjustment of the model domain to a varying cavity geometry. RAnGO simulations with a 20th-century climate forcing yield realistic basal melt rates and a quasi-stable grounding line position close to the presently observed state. In a centennial-scale warm-water-inflow scenario, the model suggests a substantial thinning of the ice shelf and a local retreat of the grounding line. The potentially negative feedback from ice-shelf thinning through a rising in situ freezing temperature is more than outweighed by the increasing water column thickness in the deepest parts of the cavity. Compared to a control simulation with fixed ice-shelf geometry, the coupled model thus yields a slightly stronger increase in ice-shelf basal melt rates.

Description: Recent observation-based results show that the seasonal amplitude of surface ocean partial pressure of CO2 (pCO2) has been increasing on average at a rate of 2–3 µatm per decade (Landschützer et al., 2018). Future increases in pCO2 seasonality are expected, as marine CO2 concentration ([CO2]) will increase in response to increasing anthropogenic carbon emissions (McNeil and Sasse, 2016). Here we use seven different global coupled atmosphere–ocean–carbon cycle–ecosystem model simulations conducted as part of the Coupled Model Intercomparison Project Phase 5 (CMIP5) to study future projections of the pCO2 annual cycle amplitude and to elucidate the causes of its amplification. We find that for the RCP8.5 emission scenario the seasonal amplitude (climatological maximum minus minimum) of upper ocean pCO2 will increase by a factor of 1.5 to 3 over the next 60–80 years. To understand the drivers and mechanisms that control the pCO2 seasonal amplification we develop a complete analytical Taylor expansion of pCO2 seasonality in terms of its four drivers: dissolved inorganic carbon (DIC), total alkalinity (TA), temperature (T), and salinity (S). Using this linear approximation we show that the DIC and T terms are the dominant contributors to the total change in pCO2 seasonality. To first order, their future intensification can be traced back to a doubling of the annual mean pCO2, which enhances DIC and alters the ocean carbonate chemistry. Regional differences in the projected seasonal cycle amplitude are generated by spatially varying sensitivity terms. The subtropical and equatorial regions (40∘ S–40∘ N) will experience a ≈30–80 µatm increase in seasonal cycle amplitude almost exclusively due to a larger background CO2 concentration that amplifies the T seasonal effect on solubility. This mechanism is further reinforced by an overall increase in the seasonal cycle of T as a result of stronger ocean stratification and a projected shoaling of mean mixed layer depths. The Southern Ocean will experience a seasonal cycle amplification of ≈90–120 µatm in response to the mean pCO2-driven change in the mean DIC contribution and to a lesser extent to the T contribution. However, a decrease in the DIC seasonal cycle amplitude somewhat counteracts this regional amplification mechanism.

Description: The responses of carbon dioxide (CO2) and other climate variables to an emission pulse of CO2 into the atmosphere are often used to compute the Global Warming Potential (GWP) and Global Temperature change Potential (GTP), to characterize the response timescales of Earth System models, and to build reduced-form models. In this carbon cycle-climate model intercomparison project, which spans the full model hierarchy, we quantify responses to emission pulses of different magnitudes injected under different conditions. The CO2 response shows the known rapid decline in the first few decades followed by a millennium-scale tail. For a 100 Gt-C emission pulse added to a constant CO2 concentration of 389 ppm, 25 ± 9% is still found in the atmosphere after 1000 yr; the ocean has absorbed 59 ± 12% and the land the remainder (16 ± 14%). The response in global mean surface air temperature is an increase by 0.20 ± 0.12 °C within the first twenty years; thereafter and until year 1000, temperature decreases only slightly, whereas ocean heat content and sea level continue to rise. Our best estimate for the Absolute Global Warming Potential, given by the time-integrated response in CO2 at year 100 multiplied by its radiative efficiency, is 92.5 × 10−15 yr W m−2 per kg-CO2. This value very likely (5 to 95% confidence) lies within the range of (68 to 117) × 10−15 yr W m−2 per kg-CO2. Estimates for time-integrated response in CO2 published in the IPCC First, Second, and Fourth Assessment and our multi-model best estimate all agree within 15% during the first 100 yr. The integrated CO2 response, normalized by the pulse size, is lower for pre-industrial conditions, compared to present day, and lower for smaller pulses than larger pulses. In contrast, the response in temperature, sea level and ocean heat content is less sensitive to these choices. Although, choices in pulse size, background concentration, and model lead to uncertainties, the most important and subjective choice to determine AGWP of CO2 and GWP is the time horizon.

Description: 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. Each of the simulations begins with a geometry and velocity close to present day observations, and evolves according to variation in meteoric ice accumulation, ice shelf melting, and mesh resolution. 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 rates 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, 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. Sensitivity to mesh resolution is spurious, and we find that sub-kilometer resolution is needed along most regions of the grounding line to avoid systematic under-estimates of the retreat rate, although resolution requirements are more stringent in some regions – for example the Amundsen Sea Embayment – than others – such as the Möller and Institute ice streams.

Description: The largest uncertainty in projections of future sea-level change still results from the potentially changing dynamical ice discharge from Antarctica. While ice discharge can alter through a number of processes, basal ice-shelf melting induced by a warming ocean has been identified as a major if not the major cause for possible additional ice flow across the grounding line. Here we derive dynamic ice-sheet response functions for basal ice-shelf melting using experiments carried out within the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. As used here these response functions provide separate contributions for four different Antarctic drainage regions. Under the assumptions of linear-response theory we project future ice-discharge for each model, each region and each of the four Representative Concentration Pathways (RCP) using oceanic temperatures from 19 comprehensive climate models of the Coupled Model Intercomparison Project, CMIP-5, and two ocean models from the EU-project Ice2Sea. Uncertainty in the climatic forcing, the oceanic response and the ice-model differences is combined into an uncertainty range of future Antarctic ice-discharge induced from basal ice-shelf melt. The additional ice-loss (Table 6) is clearly scenario-dependent and results in a median of 0.07 m (66%-range: 0.04–0.10 m; 90%-range: −0.01–0.26 m) of global sea-level equivalent for the low-emission RCP-2.6 scenario and yields 0.1 m (66%-range: 0.06–0.14 m; 90%-range: −0.01–0.45 m) for the strongest RCP-8.5. If only models with an explicit representation of ice-shelves are taken into account the scenario dependence remains and the values change to: 0.05 m (66%-range: 0.03–0.08 m) for RCP-2.6 and 0.07 m (66%-range: 0.04–0.11 m) for RCP-8.5. These results were obtained using a time delay between the surface warming signal and the subsurface oceanic warming as observed in the CMIP-5 models. Without this time delay the ranges for all ice-models changes to 0.10 m (66%-range: 0.07–0.12 m; 90%-range: 0.01–0.28 m) for RCP-2.6 and 0.15 m (66%-range: 0.10–0.21 m; 90%-range: 0.02–0.53 m) for RCP-8.5. All probability distributions as provided in Fig. 10 are highly skewed towards high values.

Description: 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.

Description: We analyzed hydrographic data from the northwestern Weddell Sea continental shelf of the three austral winters 1989, 1997, and 2006 and two summers following the last winter cruise. During summer a thermal front exists at ~64° S separating cold southern waters from warm northern waters that have similar characteristics as the deep waters of the central basin of the Bransfield Strait. In winter, the whole continental shelf exhibits southern characteristics with high Neon (Ne) concentrations, indicating a significant input of glacial melt water. The comparison of the winter data from the shallow shelf off the tip of the Antarctic Peninsula, spanning a period of 17 yr, shows a salinity decrease of 0.09 for the whole water column, which has a residence time of

Description: We develop a conceptual coupled atmosphere-ocean-ecosystem model for the tropical Pacific to investigate the interaction between marine biota and the El Niño-Southern Oscillation (ENSO). Ocean and atmosphere are represented by a two-box model for the equatorial Pacific cold tongue and the warm pool, including a simplified mixed layer scheme. Marine biota are represented by a three-component (nutrient, phytoplankton, and zooplankton) ecosystem model. The atmosphere-ocean model exhibits an oscillatory state which qualitatively captures the main physics of ENSO. During an ENSO cycle, the variation of nutrient upwelling, and, to a small extent, the variation of photosynthetically available radiation force an ecosystem oscillation. The simplified ecosystem in turn, due to the effect of phytoplankton on the absorption of shortwave radiation in the water column, leads to (1) a warming of the tropical Pacific, (2) a reduction of the ENSO amplitude, and (3) a prolongation of the ENSO period. We qualitatively investigate these bio-physical coupling mechanisms using continuation methods. It is demonstrated that bio-physical coupling may play a considerable role in modulating ENSO variability.

Description: We develop the theory of cyclic Markov chains and apply it to the El Niño-Southern Oscillation (ENSO) predictability problem. At the core of Markov chain modelling is a partition of the state space such that the transition rates between different state space cells can be computed and used most efficiently. We apply a partition technique, which divides the state space into multidimensional cells containing an equal number of data points. This partition leads to mathematical properties of the transition matrices which can be exploited further such as to establish connections with the dynamical theory of unstable periodic orbits. We introduce the concept of most and least predictable states. The data basis of our analysis consists of a multicentury-long data set obtained from an intermediate coupled atmosphere-ocean model of the tropical Pacific. This cyclostationary Markov chain approach captures the spring barrier in ENSO predictability and gives insight also into the dependence of ENSO predictability on the climatic state.