Modelling the impact of ocean warming on melting and water masses of ice shelves in the Eastern Weddell Sea

The Eastern Weddell Ice Shelves (EWIS) are believed to modify the water masses of the coastal current and thus preconditions the water mass formation in the southern and western Weddell Sea. We apply various ocean warming scenarios to investigate the impact on the temperature–salinity distribution and the sub-ice shelf melting in the Eastern Weddell Sea. In our numerical experiments, the warming is imposed homogeneously along the open inflow boundaries of the model domain, leading to a warming of the warm deep water (WDW) further downstream. Our modelling results indicate a weak quadratic dependence of the melt rate at the ice shelf base on the imposed amount of warming, which is consistent with earlier studies. The total melt rate has a strong dependence on the applied ocean warming depth. If the warming is restricted to the upper ocean (above 1,000  m), the water column (aside from the mixed surface layer) in the vicinity of the ice shelves stabilises. Hence, reduced vertical mixing will reduce the potential of Antarctic Bottom Water formation further downstream with consequences on the global thermohaline circulation. If the warming extends to the abyss, the WDW core moves significantly closer to the continental shelf break. This sharpens the Antarctic Slope Front and leads to a reduced density stratification. In contrast to the narrow shelf bathymetry in the EWIS region, a wider continental shelf (like in the southern Weddell Sea) partly protects ice shelves from remote ocean warming. Hence, the freshwater production rate of, e.g., the Filchner–Ronne Ice Shelf increases much less compared with the EWIS for identical warming scenarios. Our study therefore indicates that the ice-ocean interaction has a significant impact on the temperature-salinity distribution and the water column stability in the vicinity of ice shelves located along a narrow continental shelf. The effects of ocean warming and the impact of increased freshwater fluxes on the circulation are of the same order of magnitude and superimposed. Therefore, a consideration of this interaction in large-scale climate studies is essential.     

Ocean/ice shelf interaction in the southern Weddell Sea: results of a regional numerical helium/neon simulation

Ocean/ice interaction at the base of deep-drafted Antarctic ice shelves modifies the physical properties of inflowing shelf waters to become Ice Shelf Water (ISW). In contrast to the conditions at the atmosphere/ocean interface, the increased hydrostatic pressure at the glacial base causes gases embedded in the ice to dissolve completely after being released by melting. Helium and neon, with an extremely low solubility, are saturated in glacial meltwater by more than 1000%. At the continental slope in front of the large Antarctic caverns, ISW mixes with ambient waters to form different precursors of Antarctic Bottom Water. A regional ocean circulation model, which uses an explicit formulation of the ocean/ice shelf interaction to describe for the first time the input of noble gases to the Southern Ocean, is presented. The results reveal a long-term variability of the basal mass loss solely controlled by the interaction between waters of the continental shelf and the ice shelf cavern. Modeled helium and neon supersaturations from the Filchner–Ronne Ice Shelf front show a “low-pass” filtering of the inflowing signal due to cavern processes. On circumpolar scales, the simulated helium and neon distributions allow us to quantify the ISW contribution to bottom water, which spreads with the coastal current connecting the major formation sites in Ross and Weddell Seas.

Southern Ocean warming and increased ice shelf basal melting in the twenty-first and twenty-second centuries based on coupled ice-ocean finite-element modelling

We utilise a global finite-element sea ice–ocean model (FESOM), focused on the Antarctic marginal seas, to analyse projections of ice shelf basal melting in a warmer climate. Ice shelf–ocean interaction is described using a three-equation system with a diagnostic computation of temperature and salinity at the ice–ocean interface. A tetrahedral mesh with a minimumhorizontal resolution of 4 km and hybrid vertical coordinates is used. Ice shelf draft, cavity geometry, and global ocean bathymetry have been derived from the RTopo-1 data set. The model is forced with the atmospheric output from two climate models: (1) the Hadley Centre Climate Model (HadCM3) and (2) Max Planck Institute’s ECHAM5/MPI-OM coupled climate model. Results from experiments forced with their twentieth century output are used to evaluate the modelled present-day ocean state. Sea ice coverage is largely realistic in both simulations; modelled ice shelf basal melt rates compare well with observations in both cases, but are consistently smaller for ECHAM5/MPI-OM. Projections for future ice shelf basal melting are computed using atmospheric output for the Intergovernmental Panel on Climate Change (IPCC) scenarios E1 and A1B. In simulations forced with ECHAM5 data, trends in ice shelf basal melting are small. In contrast, decreasing convection along the Antarctic coast in HadCM3 scenarios leads to a decreasing salinity on the continental shelf and to intrusions of warm deep water of open ocean origin. In the case of the Filchner–Ronne Ice Shelf (FRIS), this water reaches deep into the cavity, so that basal melting increases by a factor of 4 to 6 compared to the present value of about 90 Gt/year. By the middle of the twenty-second century, FRIS becomes the dominant contributor to total ice shelf basal mass loss in these simulations. Our results indicate that the surface freshwater fluxes on the continental shelves may be crucial for the future of especially the large cold water ice shelves in the Southern Ocean.     

Ross ice shelf cavity circulation,residence time,and melting: Results from a model of oceanic chlorofluorocarbons

Despite their harmful effects in the upper atmosphere, anthropogenic chlorofluorocarbons dissolved in seawater are extremely useful for studying ocean circulation and ventilation, particularly in remote locations. Because they behave as a passive tracer in seawater, and their atmospheric concentrations are well-mixed, well-known, and have changed over time, they are ideal for gaining insight into the oceanographic characteristics of the isolated cavities found under Antarctic ice shelves, where direct observations are difficult to obtain. Here we present results from a modeling study of air–sea chlorofluorocarbon exchange and ocean circulation in the Ross Sea, Antarctica. We compare our model estimates of oceanic CFC-12 concentrations along an ice shelf edge transect to field data collected during three cruises spanning 16 yr. Our model produces chlorofluorocarbon concentrations that are quite similar to those measured in the field, both in magnitude and distribution, showing high values near the surface, decreasing with depth, and increasing over time. After validating modeled circulation and air–sea gas exchange through comparison of modeled temperature, salinity, and chlorofluorocarbons with field data, we estimate that the residence time of water in the Ross Ice Shelf cavity is approximately 2.2 yr and that basal melt rates for the ice shelf average 10 cm yr⁻¹. The model predicts a seasonal signature to basal melting, with highest melt rates in the spring and also the fall.     

Greenland and Antarctica Ice Sheet Mass Changes and Effects on Global Sea Level

Thirteen years of GRACE data provide an excellent picture of the current mass changes of Greenland and Antarctica, with mass loss in the GRACE period 2002–2015 amounting to 265 ± 25 GT/year for Greenland (including peripheral ice caps), and 95 ± 50 GT/year for Antarctica, corresponding to 0.72 and 0.26 mm/year average global sea level change. A significant acceleration in mass loss rate is found, especially for Antarctica, while Greenland mass loss, after a corresponding acceleration period, and a record mass loss in the summer of 2012, has seen a slight decrease in short-term mass loss trend. The yearly mass balance estimates, based on point mass inversion methods, have relatively large errors, both due to uncertainties in the glacial isostatic adjustment processes, especially for Antarctica, leakage from unmodelled ocean mass changes, and (for Greenland) difficulties in separating mass signals from the Greenland ice sheet and the adjacent Canadian ice caps. The limited resolution of GRACE affects the uncertainty of total mass loss to a smaller degree; we illustrate the “real” sources of mass changes by including satellite altimetry elevation change results in a joint inversion with GRACE, showing that mass change occurs primarily associated with major outlet glaciers, as well as a narrow coastal band. For Antarctica, the primary changes are associated with the major outlet glaciers in West Antarctica (Pine Island and Thwaites Glacier systems), as well as on the Antarctic Peninsula, where major glacier accelerations have been observed after the 2002 collapse of the Larsen B Ice Shelf.     

冰川均衡调整对南极冰质量平衡监测的影响及其不确定性

本文研究了新的全球冰川均衡调整(GIA)模型对南极冰盖质量平衡监测的影响,考虑现有冰川负荷模型和地幔黏滞度模型的差异,完整评估了结果的不确定性,最后结合GRACE和卫星测高的结果进行了对比分析.结果表明,GIA对GRACE监测的等效水柱变化有重大影响,较大的GIA影响出现在西南极,沿罗斯冰架-卡姆布冰流-罗尼冰架-南极...     

Mass budgets of the Lambert, Mellor and Fisher Glaciers and basal fluxes beneath their flowbands on Amery Ice Shelf

We used in situ measurements and remote-sensing data sets to evaluate the mass budgets of the Lambert, Mellor and Fisher Glaciers and the basal melting and freezing rates beneath their flowbands on the Amery Ice Shelf. Our findings show the Lambert and Mellor Glaciers upstream of the ANARE Lambert Glacier Basin (LGB) traverse may have positive imbalances of 3.9±2.1 Gt a-1 and 2.1±2.4 Gt a-1, respectively, while the Fisher Glacier is approximately in balance. The upstream region as a whole has a positive imbalance of 5.9±4.9 Gt a-1. The three same glaciers downstream of the ANARE LGB traverse line are in negative imbalance, where the whole downstream region has a negative imbalance of -8.5±5.8 Gt a-1. Overall the mass budgets of the Lambert, Mellor, and Fisher Glaciers are close to bal-ance, and the collective three-glacier system is also nearly in balance with a mass budget of -2.6±6.5 Gt a-1. The significant positive imbalances for the interior basin upstream of the ice-movement stations established in the early 1970s (GL line) reported previously are possibly due to an overestimate of the total accumulation and an underestimate of the ice flux through the GL line. The mean melting rate is -23.0±3.5 m ice a-1 near the southern grounding line, which decreases rapidly downstream, and transitions to refreezing at around 300 km from the southern extremity of the Amery Ice Shelf. Freezing rates along the flowbands are around 0.5±0.1 to 1.5±0.2 m ice a-1. The per-centage of ice lost from the interior by basal melting beneath the flowbands is about 80%±5%. The total basal melting and refreezing beneath the three flowbands is 50.3±7.5 Gt ice a-1 and 7.0±1.1 Gt ice a-1, respectively. We find a much larger total basal melting and net melting than the results for the whole Amery Ice Shelf derived from previous modeling and oceanographic measurements.     

Numerical study of three-dimensional shelf circulation on the Scotian Shelf using a shelf circulation model

A numerical shelf circulation model was developed for the Scotian Shelf, using a nested-grid setup consisting of a three-dimensional baroclinic inner model embedded inside a two-dimensional barotropic outer model. The shelf circulation model is based on the Princeton Ocean Model and driven by three-hourly atmospheric forcing provided by a numerical weather forecast model and by tidal forcing specified at the inner model's open boundaries based on pre-calculated tidal harmonic constants. The outer model simulates the depth-mean circulation forced by wind and atmospheric pressure fields over the northwest Atlantic Ocean with a horizontal resolution of 1/12°. The inner model simulates the three-dimensional circulation over the Gulf of St. Lawrence, the Scotian Shelf, and the adjacent slope with a horizontal resolution of 1/16°. The performance of the shelf circulation model is assessed by comparing model results with oceanographic observations made along the Atlantic coast of Nova Scotia and in the vicinity of Sable Island (on the Scotian Shelf) during two periods: October 2000–March 2001 and April–June 2002. Analysis of model results on Sable Island Bank indicates that tidal currents account for as much as ∼80% of the total variance of near-bottom currents, and currents driven by local winds account for ∼30% of the variance of the non-tidal near-bottom currents. Shelf waves generated remotely by winds and propagating into the region also play an important role in the near-bottom circulation on the bank.     

Long-term ocean simulations in FESOM: evaluation and application in studying the impact of Greenland Ice Sheet melting

The Finite Element Sea-ice Ocean Model (FESOM) is formulated on unstructured meshes and offers geometrical flexibility which is difficult to achieve on traditional structured grids. In this work, the performance of FESOM in the North Atlantic and Arctic Ocean on large time scales is evaluated in a hindcast experiment. A water-hosing experiment is also conducted to study the model sensitivity to increased freshwater input from Greenland Ice Sheet (GrIS) melting in a 0.1-Sv discharge rate scenario. The variability of the Atlantic Meridional Overturning Circulation (AMOC) in the hindcast experiment can be explained by the variability of the thermohaline forcing over deep convection sites. The model also reproduces realistic freshwater content variability and sea ice extent in the Arctic Ocean. The anomalous freshwater in the water-hosing experiment leads to significant changes in the ocean circulation and local dynamical sea level (DSL). The most pronounced DSL rise is in the northwest North Atlantic as shown in previous studies, and also in the Arctic Ocean. The released GrIS freshwater mainly remains in the North Atlantic, Arctic Ocean and the west South Atlantic after 120 model years. The pattern of ocean freshening is similar to that of the GrIS water distribution, but changes in ocean circulation also contribute to the ocean salinity change. The changes in Arctic and sub-Arctic sea level modify exchanges between the Arctic Ocean and subpolar seas, and hence the role of the Arctic Ocean in the global climate. Not only the strength of the AMOC, but also the strength of its decadal variability is notably reduced by the anomalous freshwater input. A comparison of FESOM with results from previous studies shows that FESOM can simulate past ocean state and the impact of increased GrIS melting well.

     

Long-term ocean simulations in FESOM: evaluation and application in studying the impact of Greenland Ice Sheet melting

The Finite Element Sea-ice Ocean Model (FESOM) is formulated on unstructured meshes and offers geometrical flexibility which is difficult to achieve on traditional structured grids. In this work, the performance of FESOM in the North Atlantic and Arctic Ocean on large time scales is evaluated in a hindcast experiment. A water-hosing experiment is also conducted to study the model sensitivity to increased freshwater input from Greenland Ice Sheet (GrIS) melting in a 0.1-Sv discharge rate scenario. The variability of the Atlantic Meridional Overturning Circulation (AMOC) in the hindcast experiment can be explained by the variability of the thermohaline forcing over deep convection sites. The model also reproduces realistic freshwater content variability and sea ice extent in the Arctic Ocean. The anomalous freshwater in the water-hosing experiment leads to significant changes in the ocean circulation and local dynamical sea level (DSL). The most pronounced DSL rise is in the northwest North Atlantic as shown in previous studies, and also in the Arctic Ocean. The released GrIS freshwater mainly remains in the North Atlantic, Arctic Ocean and the west South Atlantic after 120 model years. The pattern of ocean freshening is similar to that of the GrIS water distribution, but changes in ocean circulation also contribute to the ocean salinity change. The changes in Arctic and sub-Arctic sea level modify exchanges between the Arctic Ocean and subpolar seas, and hence the role of the Arctic Ocean in the global climate. Not only the strength of the AMOC, but also the strength of its decadal variability is notably reduced by the anomalous freshwater input. A comparison of FESOM with results from previous studies shows that FESOM can simulate past ocean state and the impact of increased GrIS melting well.     

Detection of crevasses over polar ice shelves using Satellite Laser Altimeter

Ice shelf breakups account for most mass loss from the Antarctic Ice Sheet as the consequence of the propagation of crevasses(or rift)in response to stress.Thus there is a pressing need for detecting crevasses’location and depth,to understand the mechanism of calving processes.This paper presents a method of crevasse detection using the ICESat-1/GLAS data.A case study was taken at the Amery Ice Shelf of Antarctica to verify the accuracy of geo-location and depth of crevasses detected.Moreover,based on the limited crevasse points,we developed a method to detect the peak stress points which can be used to track the location of the crack tips and to identify the possible high-risk area where an ice shelf begins to break up.The spatial and temporal distribution of crevasse depth and the spatial distribution of peak stress points of the Amery Ice Shelf were analyzed through 132 tracks in 16 campaign periods of ICESat-1/GLAS between 2003 and 2008.The results showed that the depth of the detected crevasse points ranged from 2 to 31.7 m,which were above the sea level;the crevasse that advected downstream to the front edge of an ice shelf has little possibility to directly result in breakups because the crevasse depth did not show any increasing trend over time;the local stress concentration is distributed mainly in the suture zones on the ice shelves.     

Numerical simulation of ice-ocean variability in the Barents Sea region

A dynamic–thermodynamic sea ice model has been coupled to a three-dimensional ocean general circulation model for the purpose of conducting ocean climate dynamical downscaling experiments for the Barents Sea region. To assess model performance and suitability for such an application, the coupled model has been used to conduct a hindcast for the period 1990–2002. A comparison with available observations shows that the model successfully tracks seasonal and inter-annual variability in the ocean temperature field and that the simulated horizontal and vertical distribution of temperature are in good agreement with observations. The model results follow the seasonal and inter-annual variability in sea ice cover in the region, with the exception that the model results show too much ice melting in the northern Barents Sea during summer. The spatial distribution of the winter simulated sea ice cover is in close agreement with observations. Modelled temperatures and ice concentrations in the central Barents Sea are biased too high and too low, respectively. The probable cause is too high inflow of Atlantic Water into the Barents. The seasonal and inter-annual fluctuations in temperature and sea ice cover in the central Barents are, however, in excellent agreement with observations. Salt release during the freezing process in the numerical simulation exhibits considerable inter-annual variability and tends to vary in an opposite manner to the net inflow volume flux at the western entrance of the Barents Sea. Overall, the model produces realistic ice-ocean seasonal and inter-annual variability and should prove to be a useful tool for dynamical downscaling applications.     

Evaluation of the glacial isostatic adjustment (GIA) models for Antarctica based on GPS vertical velocities

Due to the scarcity of data, modeling the glacial isostatic adjustment(GIA) for Antarctica is more difficult than it is for the ancient ice sheet area in North America and Northern Europe. Large uncertainties are observed in existing GIA models for Antarctica. Modern space-based geodetic measurements provide checks and constraints for GIA models. The present-day uplift velocities of global positioning system(GPS) stations at 73 stations in Antarctica and adjacent regions from 1996 to 2014 have been estimated using GAMIT/GLOBK version 10.5 with a colored noise model. To easily analyze the effect of difference sources on the vertical velocities, and for easy comparison with both GIA model predictions and GPS results from Argus et al.(2014) and Thomas et al.(2011), seven sub-regions are divided. They are the northern Antarctic Peninsula, the Filchner-Ronne Ice Shelf, the Amundsen Sea coast, the Ross Ice Shelf, Mount Erebus, inland Southwest Antarctica and the East Antarctic coast,respectively. The results show that the fast uplift in the north Antarctic Peninsula and Pine Island Bay regions may be caused by the elastic response to snow and ice mass loss. The fast subsidence near Mount Erebus may be related to the activity of a magma body. The uplift or subsidence near the East Antarctic coast is very slow while the uplift for the rest regions is mainly caused by GIA. By analyzing the correlation and the associated weighted root mean square(WRMS) between the GIA predictions and the GPS velocities, we found that the ICE-6G_C(VM5a) model and the Geruo 13 model show the most consistency with our GPS results, while the W12a and IJ05_R2 series models show poor consistency with our GPS results. Although improved greatly in recent years, the GIA modeling in Antarctica still lags behind the modeling of the North American. Some GPS stations, for example the Bennett Nunatak station(BENN), have observed large discrepancies between GIA predictions and GPS velocities.Because of the large uncertainties in calculating elastic responses due to the significant variations of ice and snow loads, the GPS velocities still cannot be used as a precise constraint on GIA models.     

Study on subtidal circulation and variability in the Gulf of St. Lawrence,Scotian Shelf,and Gulf of Maine using a nested-grid shelf circulation model

This paper examines the subtidal circulation and associated variability in the Gulf of St. Lawrence, the Scotian Shelf, and the Gulf of Maine (GSL-SS-GOM) in 1988–2004 based on results produced by a nested-grid shelf circulation model. The model has a fine-resolution child model (～ (1/12)°) embedded inside a coarse-resolution parent model (～ (1/4)°) of the northwest Atlantic. The combination of the semi-prognostic method and the spectral nudging method is used to reduce the model seasonal bias and drift. The child model reproduces the general features of the observed circulation and hydrography over the study region during the study period. The child model results demonstrate that the time-mean circulation in the GSL is affected by the time-mean atmospheric forcing and inflow through the Strait of Belle Isle. The temporal variability in atmospheric forcing affects the outflow through western Cabot Strait, which in turn affects the transport of the Nova Scotian Current and the gulf-wide cyclonic circulation in the GOM. The simulated seasonal variability of salinity in the top 30 m of the GSL-SS-GOM is mainly affected by the equatorward advection of low-salinity waters from the lower St. Lawrence Estuary to the GOM through the Scotian Shelf. The simulated intraseasonal variability of circulation in the GSL is affected by the variability in the estuarine circulation in response to the temporal variability in atmospheric forcing. On the Scotian Shelf, the intraseasonal variability is mainly driven by the variability of wind forcing and mesoscale and nonlinear dynamics over the shelf break and slope region. The interannual variability in the simulated temperature and salinity are spatially coherent in the intermediate waters in the GSL, which is caused partially by the local response to atmospheric variability and partially by variabilities over the southern Newfoundland Shelf that enter the GSL through the eastern Cabot Strait. By comparison, on the Scotian Shelf, the interannual variability of simulated circulation is affected by anomalies produced by the nonlinear dynamics which are advected equatorwards by the shelf break jet.     

Long-range detection of hydroacoustic signals from large icebergs in the Ross Sea, Antarctica

Hydroacoustic signals detected in late 2000 by seismic stations in Polynesia are shown to originate from huge icebergs which at the time were drifting in the Ross Sea after calving off the Ross Ice Shelf. The signals present a broad variety of spectral characteristics, most of them featuring prominent eigenfrequencies in the 4-7 Hz range, often complemented by overtones. Most epicenters, obtained by combining observations of distant hydroacoustic and regional seismic records, follow the spatio-temporal evolution of the drift of iceberg B-15B. Most of the signals are generated during a 36-day time window when it is speculated that B-15B collided with smaller icebergs or was scraping the ocean floor on the shallow continental shelf. We speculate on the possible physical nature of the resonator generating the signals, which could correspond to an elastic mode of the iceberg, or to the oscillation of fluid-filled cracks in the ice.     

Response of the Greenland and Antarctic Ice Sheets to Multi-Millennial Greenhouse Warming in the Earth System Model of Intermediate Complexity LOVECLIM

Calculations were performed with the Earth system model of intermediate complexity LOVECLIM to study the response of the Greenland and Antarctic ice sheets to sustained multi-millennial greenhouse warming. Use was made of fully dynamic 3D thermomechanical ice-sheet models bidirectionally coupled to an atmosphere and an ocean model. Two 3,000-year experiments were evaluated following forcing scenarios with atmospheric CO₂ concentration increased to two and four times the pre-industrial value, and held constant thereafter. In the high concentration scenario the model shows a sustained mean annual warming of up to 10°C in both polar regions. This leads to an almost complete disintegration of the Greenland ice sheet after 3,000 years, almost entirely caused by increased surface melting. Significant volume loss of the Antarctic ice sheet takes many centuries to initiate due to the thermal inertia of the Southern Ocean but is equivalent to more than 4 m of global sea-level rise by the end of simulation period. By that time, surface conditions along the East Antarctic ice sheet margin take on characteristics of the present-day Greenland ice sheet. West Antarctic ice shelves have thinned considerably from subshelf melting and grounding lines have retreated over distances of several 100 km, especially for the Ross ice shelf. In the low concentration scenario, corresponding to a local warming of 3?C4°C, polar ice-sheet melting proceeds at a much lower rate. For the first 1,200 years, the Antarctic ice sheet is even slightly larger than today on account of increased accumulation rates but contributes positively to sea-level rise after that. The Greenland ice sheet loses mass at a rate equivalent to 35 cm of global sea level rise during the first 1,000 years increasing to 150 cm during the last 1,000 years. For both scenarios, ice loss from the Antarctic ice sheet is still accelerating after 3,000 years despite a constant greenhouse gas forcing after the first 70?C140 years of the simulation.     

Reconstruction of the annual variations of thermohaline characteristics and water circulation on the northeastern Sakhalin Shelf

Methods of statistical analysis and GIS Sakhalin Shelf have been used to reconstruct the annual variations of water temperature, salinity, and density from the surface to the bed at 60 standard oceanographic stations on northeastern Sakhalin shelf (with the total duration of observations of ～40 years). Numerical simulation using Bergen University oceanic model has been applied to calculate the spatial fields of water temperature and salinity at standard horizons in nodes of a regular grid. The analysis of calculation results revealed new features of the hydrological regime and gave previously unknown thermohaline fields with maximal detail for full annual cycle. The obtained data were used to evaluate the mean monthly three-dimensional water circulation.     

Modeled steric and mass-driven sea level change caused by Greenland Ice Sheet melting

Meltwater from the Greenland Ice Sheet (GIS) has been a major contributor to sea level change in the recent past. Global and regional sea level variations caused by melting of the GIS are investigated with the finite element sea-ice ocean model (FESOM). We consider changes of local density (steric effects), mass inflow into the ocean, redistribution of mass, and gravitational effects. Five melting scenarios are simulated, where mass losses of 100, 200, 500, and 1000 Gt/yr are converted to a continuous volume flux that is homogeneously distributed along the coast of Greenland south of 75°N. In addition, a scenario of regional melt rates is calculated from daily ice melt characteristics. The global mean sea level modeled with FESOM increases by about 0.3 mm/yr if 100 Gt/yr of ice melts, which includes eustatic and steric sea level change. In the global mean the steric contribution is one order of magnitude smaller than the eustatic contribution. Regionally, especially in the North Atlantic, the steric contribution leads to strong deviations from the global mean sea level change. The modeled pattern mainly reflects the structure of temperature and salinity change in the upper ocean. Additionally, small steric variations occur due to local variability in the heat exchange between the atmosphere and the ocean. The mass loss has also affects on the gravitational attraction by the ice sheet, causing spatially varying sea level change mainly near the GIS, but also at greater distances. This effect is accounted for by using Green's functions.     

Preserved landscapes underneath the Antarctic Ice Sheet reveal the geomorphological history of Jutulstraumen Basin

The landscape of Antarctica, hidden beneath kilometre-thick ice in most places, has been shaped by the interactions between tectonic and erosional processes. The flow dynamics of the thick ice cover deepened pre-formed topographic depressions by glacial erosion, but also preserved the subglacial landscapes in regions with moderate to slow ice flow. Mapping the spatial variability of these structures provides the basis for reconstruction of the evolution of subglacial morphology. This study focuses on the Jutulstraumen Glacier drainage system in Dronning Maud Land, East Antarctica. The Jutulstraumen Glacier reaches the ocean via the Jutulstraumen Graben, which is the only significant passage for draining the East Antarctic Ice Sheet through the western part of the Dronning Maud Land mountain chain. We acquired new bed topography data during an airborne radar campaign in the region upstream of the Jutulstraumen Graben to characterise the source area of the glacier. The new data show a deep relief to be generally under-represented in available bed topography compilations. Our analysis of the bed topography, valley characteristics and bed roughness leads to the conclusion that much more of the alpine landscape that would have formed prior to the Antarctic Ice Sheet is preserved than previously anticipated. We identify an active and deeply eroded U-shaped valley network next to largely preserved passive fluvial and glacial modified landscapes. Based on the landscape classification, we reconstruct the temporal sequence by which ice flow modified the topography since the beginning of the glaciation of Antarctica.     

Response of a Coupled Ocean�CAtmosphere Model to Greenland Ice Melting

We investigate the transient response of the global coupled ocean?Catmosphere system to enhanced freshwater forcing representative of melting of the Greenland ice sheets. A 50-year long simulation by a coupled atmosphere?Cocean general circulation model (CGCM) is compared with another of the same length in which Greenland melting is prescribed. To highlight the importance of coupled atmosphere?Cocean processes, the CGCM results are compared with those of two other experiments carried out with the oceanic general circulation model (OGCM). In one of these OGCM experiments, the prescribed surface fluxes of heat, momentum and freshwater correspond to the unperturbed simulation by the CGCM; in the other experiment, Greenland melting is added to the freshwater flux. The responses by the CGCM and OGCM to the Greenland melting have similar patterns in the Atlantic, albeit the former having five times larger amplitudes in sea surface height anomalies. The CGCM shows likewise stronger variability in all state variables in all ocean basins because the impact of Greenland melting is quickly communicated to all ocean basins via atmospheric bridges. We conclude that the response of the global climate to Greenland ice melting is highly dependent on coupled atmosphere?Cocean processes. These lead to reduced latent heat flux into the atmosphere and an associated increase in net freshwater flux into the ocean, especially in the subpolar North Atlantic. The combined result is a stronger response of the coupled system to Greenland ice sheet melting.