Refine
Year of publication
Document Type
- Article (87)
- Postprint (8)
- Other (3)
- Moving Images (1)
Is part of the Bibliography
- yes (99)
Keywords
- climate change (9)
- model (5)
- sea-level rise (4)
- Greenland (3)
- West Antarctica (3)
- climate impacts (3)
- impact (3)
- tropical cyclones (3)
- Disaster impact analysis (2)
- chains (2)
Institute
- Institut für Physik und Astronomie (75)
- Institut für Geowissenschaften (13)
- Mathematisch-Naturwissenschaftliche Fakultät (7)
- Institut für Informatik und Computational Science (2)
- Extern (1)
- Fachgruppe Volkswirtschaftslehre (1)
- Institut für Umweltwissenschaften und Geographie (1)
- Potsdam Institute for Climate Impact Research (PIK) e. V. (1)
- Wirtschaftswissenschaften (1)
- Zentrum für Umweltwissenschaften (1)
A parameterization for the motion of ice-shelf fronts on a Cartesian grid in finite-difference land-ice models is presented. The scheme prevents artificial thinning of the ice shelf at its edge, which occurs due to the finite resolution of the model. The intuitive numerical implementation diminishes numerical dispersion at the ice front and enables the application of physical boundary conditions to improve the calculation of stress and velocity fields throughout the ice-sheet-shelf system. Numerical properties of this subgrid modification are assessed in the Potsdam Parallel Ice Sheet Model (PISM-PIK) for different geometries in one and two horizontal dimensions and are verified against an analytical solution in a flow-line setup.
Recent observations and modeling studies emphasize the crucial role of fracture mechanics for the stability of ice shelves and thereby the evolution of ice sheets. Here we introduce a macroscopic fracture-density field into a prognostic continuum ice-flow model and compute its evolution incorporating the initiation and growth of fractures as well as their advection with two-dimensional ice flow. To a first approximation, fracture growth is assumed to depend on the spreading rate only, while fracture initiation is defined in terms of principal stresses. The inferred fracture-density fields compare well with observed elongate surface structures. Since crevasses and other deep-reaching fracture structures have been shown to influence the overall ice-shelf dynamics, we propose the fracture-density field introduced here be used as a measure for ice softening and decoupling of the ice flow in fracture-weakened zones. This may yield more accurate and realistic velocity patterns in prognostic simulations. Additionally, the memory of past fracture events links the calving front to the upstream dynamics. Thus the fracture-density field proposed here may be employed in fracture-based calving parameterizations. The aim of this study is to introduce the field and investigate which of the observed surface structures can be reproduced by the simplest physically motivated fracture source terms.
Antarctic ice-discharge constitutes the largest uncertainty in future sea-level projections. Floating ice shelves, fringing most of Antarctica, exert retentive forces onto the ice flow. While abrupt ice-shelf retreat has been observed, it is generally considered a localized phenomenon. Here we show that the disintegration of an ice shelf may induce the spontaneous retreat of its neighbor. As an example, we reproduce the spontaneous but gradual retreat of the Larsen B ice front as observed after the disintegration of the adjacent Larsen A ice shelf. We show that the Larsen A collapse yields a change in spreading rate in Larsen B via their connecting ice channels and thereby causes a retreat of the ice front to its observed position of the year 2000, prior to its collapse. This mechanism might be particularly relevant for the role of East Antarctica and the Antarctic Peninsula in future sea level.
Floating ice shelves can exert a retentive and hence stabilizing force onto the inland ice sheet of Antarctica. However, this effect has been observed to diminish by the dynamic effects of fracture processes within the protective ice shelves, leading to accelerated ice flow and hence to a sea-level contribution. In order to account for the macroscopic effect of fracture processes on large-scale viscous ice dynamics (i.e., ice-shelf scale) we apply a continuum representation of fractures and related fracture growth into the prognostic Parallel Ice Sheet Model (PISM) and compare the results to observations. To this end we introduce a higher order accuracy advection scheme for the transport of the two-dimensional fracture density across the regular computational grid. Dynamic coupling of fractures and ice flow is attained by a reduction of effective ice viscosity proportional to the inferred fracture density. This formulation implies the possibility of non-linear threshold behavior due to self-amplified fracturing in shear regions triggered by small variations in the fracture-initiation threshold. As a result of prognostic flow simulations, sharp across-flow velocity gradients appear in fracture-weakened regions. These modeled gradients compare well in magnitude and location with those in observed flow patterns. This model framework is in principle expandable to grounded ice streams and provides simple means of investigating climate-induced effects on fracturing (e. g., hydro fracturing) and hence on the ice flow. It further constitutes a physically sound basis for an enhanced fracture-based calving parameterization.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
Simulations of the glacial-interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. <br /> We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a "cookbook" for the growing community of PISM users and paleo-ice sheet modelers in general. <br /> For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (approximate to 210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation-age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model-data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016).
This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4 +/- 4.1m (or 6.5 +/- 2.0 x 10(6) km(3)), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.
Ten ice-sheet models are used to study sensitivity of the Greenland and Antarctic ice sheets to prescribed changes of surface mass balance, sub-ice-shelf melting and basal sliding. Results exhibit a large range in projected contributions to sea-level change. In most cases, the ice volume above flotation lost is linearly dependent on the strength of the forcing. Combinations of forcings can be closely approximated by linearly summing the contributions from single forcing experiments, suggesting that nonlinear feedbacks are modest. Our models indicate that Greenland is more sensitive than Antarctica to likely atmospheric changes in temperature and precipitation, while Antarctica is more sensitive to increased ice-shelf basal melting. An experiment approximating the Intergovernmental Panel on Climate Change's RCP8.5 scenario produces additional first-century contributions to sea level of 22.3 and 8.1 cm from Greenland and Antarctica, respectively, with a range among models of 62 and 14 cm, respectively. By 200 years, projections increase to 53.2 and 26.7 cm, respectively, with ranges of 79 and 43 cm. Linear interpolation of the sensitivity results closely approximates these projections, revealing the relative contributions of the individual forcings on the combined volume change and suggesting that total ice-sheet response to complicated forcings over 200 years can be linearized.
The 8.2 ka event : abrupt transition of the subpolar gyre toward a modern North Atlantic circulation
(2010)
Climate model simulations of the 8.2 ka event show an abrupt strengthening of the Atlantic subpolar gyre that allows us to connect two major but apparently contradictory climate events of the early Holocene: the freshwater outburst from proglacial lakes and the onset of Labrador Sea water formation. The 8.2 ka event is the largest climatic signal of our present interglacial with a widespread cooling in the North Atlantic region about 8200 years before present. It coincides with a meltwater outburst from North American proglacial lakes that is believed to have weakened the Atlantic meridional overturning circulation and northward heat transport, followed by a recovery of the deep ocean circulation and rising temperatures after a few centuries. Marine proxy data, however, date the onset of deep water formation in Labrador Sea to the same time. The subsequent strengthening of the slope current system created a regional signal recorded as an abrupt and persistent surface temperature decrease. Although similarities in timing are compelling, a mechanism to reconcile these apparently contradictory events was missing. Our simulations show that an abrupt and persistent strengthening of the Atlantic subpolar gyre provides a plausible explanation. The intense freshwater pulse triggered a transition of the gyre circulation into a different mode of operation, stabilized by internal feedbacks and persistent after the cessation of the perturbation. As a direct consequence, deep water formation around its center intensifies. This corresponds to the modern flow regime and stabilizes the meridional overturning circulation, possibly contributing to the Holocene's climatic stability.
We investigate the sensitivity of a coarse resolution coupled climate model to the representation of the overflows over the Greenland-Scotland ridge. This class of models suffers from a poor representation of the water mass exchange between the Nordic Seas and the North Atlantic, a crucial part of the large-scale oceanic circulation. We revisit the explicit representation of the overflows using a parameterisation by hydraulic constraints and compare it with the enhancement of the overflow transport by artificially deepened passages over the Greenland-Scotland ridge, a common practice in coarse resolution models. Both configurations increase deep water formation in the Nordic Seas and represent the large-scale dynamics of the Atlantic realistically in contrast to a third model version with realistic sill depths but without the explicit overflow transport. The comparison of the hydrography suggests that for the unperturbed equilibrium the Nordic Seas are better represented with the parameterised overflows. As in previous studies, we do not find a stabilising effect of the overflow parameterisation on the Atlantic meridional overturning circulation but merely on the overflow transport. As a consequence the surface air temperature in the Nordic Seas is less sensitive to anomalous surface fresh water forcing. Special attention is paid to changes in the subpolar gyre circulation. We find it sensitive to the overflow transport and the density of these water masses through baroclinic adjustments. The analysis of the governing equations confirms the presence of positive feedbacks inherent to the subpolar gyre and allows us to isolate the influence of the overflows on its dynamics.
The Atlantic subpolar gyre (SPG) is one of the main drivers of decadal climate variability in the North Atlantic. Here we analyze its dynamics in pre-industrial control simulations of 19 different comprehensive coupled climate models. The analysis is based on a recently proposed description of the SPG dynamics that found the circulation to be potentially bistable due to a positive feedback mechanism including salt transport and enhanced deep convection in the SPG center. We employ a statistical method to identify multiple equilibria in time series that are subject to strong noise and analyze composite fields to assess whether the bistability results from the hypothesized feedback mechanism. Because noise dominates the time series in most models, multiple circulation modes can unambiguously be detected in only six models. Four of these six models confirm that the intensification is caused by the positive feedback mechanism.