@article{WinkelmannMartinHaseloffetal.2011, author = {Winkelmann, Ricarda and Martin, Maria A. and Haseloff, Monika and Albrecht, Torsten and Bueler, Ed and Khroulev, C. and Levermann, Anders}, title = {The Potsdam parallel ice sheet model (PISM-PIK) - Part 1: Model description}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {5}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {3}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-5-715-2011}, pages = {715 -- 726}, year = {2011}, abstract = {We present the Potsdam Parallel Ice Sheet Model (PISM-PIK), developed at the Potsdam Institute for Climate Impact Research to be used for simulations of large-scale ice sheet-shelf systems. It is derived from the Parallel Ice Sheet Model (Bueler and Brown, 2009). Velocities are calculated by superposition of two shallow stress balance approximations within the entire ice covered region: the shallow ice approximation (SIA) is dominant in grounded regions and accounts for shear deformation parallel to the geoid. The plug-flow type shallow shelf approximation (SSA) dominates the velocity field in ice shelf regions and serves as a basal sliding velocity in grounded regions. Ice streams can be identified diagnostically as regions with a significant contribution of membrane stresses to the local momentum balance. All lateral boundaries in PISM-PIK are free to evolve, including the grounding line and ice fronts. Ice shelf margins in particular are modeled using Neumann boundary conditions for the SSA equations, reflecting a hydrostatic stress imbalance along the vertical calving face. The ice front position is modeled using a subgrid-scale representation of calving front motion (Albrecht et al., 2011) and a physically-motivated calving law based on horizontal spreading rates. The model is tested in experiments from the Marine Ice Sheet Model Intercomparison Project (MISMIP). A dynamic equilibrium simulation of Antarctica under present-day conditions is presented in Martin et al. (2011).}, language = {en} } @article{MartinWinkelmannHaseloffetal.2011, author = {Martin, Maria A. and Winkelmann, Ricarda and Haseloff, M. and Albrecht, Tanja and Bueler, Ed and Khroulev, C. and Levermann, Anders}, title = {The Potsdam parallel ice sheet model (PISM-PIK) - Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {5}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {3}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-5-727-2011}, pages = {727 -- 740}, year = {2011}, abstract = {We present a dynamic equilibrium simulation of the ice sheet-shelf system on Antarctica with the Potsdam Parallel Ice Sheet Model (PISM-PIK). The simulation is initialized with present-day conditions for bed topography and ice thickness and then run to steady state with constant present-day surface mass balance. Surface temperature and sub-shelf basal melt distribution are parameterized. Grounding lines and calving fronts are free to evolve, and their modeled equilibrium state is compared to observational data. A physically-motivated calving law based on horizontal spreading rates allows for realistic calving fronts for various types of shelves. Steady-state dynamics including surface velocity and ice flux are analyzed for whole Antarctica and the Ronne-Filchner and Ross ice shelf areas in particular. The results show that the different flow regimes in sheet and shelves, and the transition zone between them, are captured reasonably well, supporting the approach of superposition of SIA and SSA for the representation of fast motion of grounded ice. This approach also leads to a natural emergence of sliding-dominated flow in stream-like features in this new 3-D marine ice sheet model.}, language = {en} } @article{AlbrechtMartinHaseloffetal.2011, author = {Albrecht, Tanja and Martin, M. and Haseloff, M. and Winkelmann, Ricarda and Levermann, Anders}, title = {Parameterization for subgrid-scale motion of ice-shelf calving fronts}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {5}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-5-35-2011}, pages = {35 -- 44}, year = {2011}, abstract = {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.}, language = {en} } @article{ScheweLevermannMeinshausen2011, author = {Schewe, Jacob and Levermann, Anders and Meinshausen, Malte}, title = {Climate change under a scenario near 1.5 degrees C of global warming: monsoon intensification, ocean warming and steric sea level rise}, series = {Earth system dynamics}, volume = {2}, journal = {Earth system dynamics}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2190-4979}, doi = {10.5194/esd-2-25-2011}, pages = {25 -- 35}, year = {2011}, abstract = {We present climatic consequences of the Representative Concentration Pathways (RCPs) using the coupled climate model CLIMBER-3 alpha, which contains a statistical-dynamical atmosphere and a three-dimensional ocean model. We compare those with emulations of 19 state-of-the-art atmosphere-ocean general circulation models (AOGCM) using MAGICC6. The RCPs are designed as standard scenarios for the forthcoming IPCC Fifth Assessment Report to span the full range of future greenhouse gas (GHG) concentrations pathways currently discussed. The lowest of the RCP scenarios, RCP3-PD, is projected in CLIMBER-3 alpha to imply a maximal warming by the middle of the 21st century slightly above 1.5 degrees C and a slow decline of temperatures thereafter, approaching today's level by 2500. We identify two mechanisms that slow down global cooling after GHG concentrations peak: The known inertia induced by mixing-related oceanic heat uptake; and a change in oceanic convection that enhances ocean heat loss in high latitudes, reducing the surface cooling rate by almost 50\%. Steric sea level rise under the RCP3-PD scenario continues for 200 years after the peak in surface air temperatures, stabilizing around 2250 at 30 cm. This contrasts with around 1.3 m of steric sea level rise by 2250, and 2 m by 2500, under the highest scenario, RCP8.5. Maximum oceanic warming at intermediate depth (300-800 m) is found to exceed that of the sea surface by the second half of the 21st century under RCP3-PD. This intermediate-depth warming persists for centuries even after surface temperatures have returned to present-day values, with potential consequences for marine ecosystems, oceanic methane hydrates, and ice-shelf stability. Due to an enhanced land-ocean temperature contrast, all scenarios yield an intensification of monsoon rainfall under global warming.}, language = {en} } @article{SchleussnerFrielerMeinshausenetal.2011, author = {Schleussner, Carl-Friedrich and Frieler, Katja and Meinshausen, Malte and Yin, J. and Levermann, Anders}, title = {Emulating Atlantic overturning strength for low emission scenarios consequences for sea-level rise along the North American east coast}, series = {Earth system dynamics}, volume = {2}, journal = {Earth system dynamics}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2190-4979}, doi = {10.5194/esd-2-191-2011}, pages = {191 -- 200}, year = {2011}, abstract = {In order to provide probabilistic projections of the future evolution of the Atlantic Meridional Overturning Circulation (AMOC), we calibrated a simple Stommel-type box model to emulate the output of fully coupled three-dimensional atmosphere-ocean general circulation models (AOGCMs) of the Coupled Model Intercomparison Project (CMIP). Based on this calibration to idealised global warming scenarios with and without interactive atmosphere-ocean fluxes and freshwater perturbation simulations, we project the future evolution of the AMOC mean strength within the covered calibration range for the lower two Representative Concentration Pathways (RCPs) until 2100 obtained from the reduced complexity carbon cycle-climate model MAGICC 6. For RCP3-PD with a global mean temperature median below 1.0 degrees C warming relative to the year 2000, we project an ensemble median weakening of up to 11\% compared to 22\% under RCP4.5 with a warming median up to 1.9 degrees C over the 21st century. Additional Greenland meltwater of 10 and 20 cm of global sea-level rise equivalent further weakens the AMOC by about 4.5 and 10 \%, respectively. By combining our outcome with a multi-model sea-level rise study we project a dynamic sea-level rise along the New York City coastline of 4 cm for the RCP3-PD and of 8 cm for the RCP4.5 scenario over the 21st century. We estimate the total steric and dynamic sea-level rise for New York City to be about 24 cm until 2100 for the RCP3-PD scenario, which can hold as a lower bound for sea-level rise projections in this region, as it does not include ice sheet and mountain glacier contributions.}, language = {en} } @article{LevermannAlbrechtWinkelmannetal.2012, author = {Levermann, Anders and Albrecht, Tanja and Winkelmann, Ricarda and Martin, Maria A. and Haseloff, Monika and Joughin, I.}, title = {Kinematic first-order calving law implies potential for abrupt ice-shelf retreat}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {6}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-6-273-2012}, pages = {273 -- 286}, year = {2012}, abstract = {Recently observed large-scale disintegration of Antarctic ice shelves has moved their fronts closer towards grounded ice. In response, ice-sheet discharge into the ocean has accelerated, contributing to global sea-level rise and emphasizing the importance of calving-front dynamics. The position of the ice front strongly influences the stress field within the entire sheet-shelf-system and thereby the mass flow across the grounding line. While theories for an advance of the ice-front are readily available, no general rule exists for its retreat, making it difficult to incorporate the retreat in predictive models. Here we extract the first-order large-scale kinematic contribution to calving which is consistent with large-scale observation. We emphasize that the proposed equation does not constitute a comprehensive calving law but represents the first-order kinematic contribution which can and should be complemented by higher order contributions as well as the influence of potentially heterogeneous material properties of the ice. When applied as a calving law, the equation naturally incorporates the stabilizing effect of pinning points and inhibits ice shelf growth outside of embayments. It depends only on local ice properties which are, however, determined by the full topography of the ice shelf. In numerical simulations the parameterization reproduces multiple stable fronts as observed for the Larsen A and B Ice Shelves including abrupt transitions between them which may be caused by localized ice weaknesses. We also find multiple stable states of the Ross Ice Shelf at the gateway of the West Antarctic Ice Sheet with back stresses onto the sheet reduced by up to 90 \% compared to the present state.}, language = {en} } @article{SchneidervonDeimlingMeinshausenLevermannetal.2012, author = {Schneider von Deimling, Thomas and Meinshausen, Malte and Levermann, Anders and Huber, Veronika and Frieler, Katja and Lawrence, D. M. and Brovkin, Victor}, title = {Estimating the near-surface permafrost-carbon feedback on global warming}, series = {Biogeosciences}, volume = {9}, journal = {Biogeosciences}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1726-4170}, doi = {10.5194/bg-9-649-2012}, pages = {649 -- 665}, year = {2012}, abstract = {Thawing of permafrost and the associated release of carbon constitutes a positive feedback in the climate system, elevating the effect of anthropogenic GHG emissions on global-mean temperatures. Multiple factors have hindered the quantification of this feedback, which was not included in climate carbon-cycle models which participated in recent model intercomparisons (such as the Coupled Carbon Cycle Climate Model Intercomparison Project - (CMIP)-M-4). There are considerable uncertainties in the rate and extent of permafrost thaw, the hydrological and vegetation response to permafrost thaw, the decomposition timescales of freshly thawed organic material, the proportion of soil carbon that might be emitted as carbon dioxide via aerobic decomposition or as methane via anaerobic decomposition, and in the magnitude of the high latitude amplification of global warming that will drive permafrost degradation. Additionally, there are extensive and poorly characterized regional heterogeneities in soil properties, carbon content, and hydrology. Here, we couple a new permafrost module to a reduced complexity carbon-cycle climate model, which allows us to perform a large ensemble of simulations. The ensemble is designed to span the uncertainties listed above and thereby the results provide an estimate of the potential strength of the feedback from newly thawed permafrost carbon. For the high CO2 concentration scenario (RCP8.5), 33-114 GtC (giga tons of Carbon) are released by 2100 (68\% uncertainty range). This leads to an additional warming of 0.04-0.23 degrees C. Though projected 21st century permafrost carbon emissions are relatively modest, ongoing permafrost thaw and slow but steady soil carbon decomposition means that, by 2300, about half of the potentially vulnerable permafrost carbon stock in the upper 3 m of soil layer (600-1000 GtC) could be released as CO2, with an extra 1-4\% being released as methane. Our results also suggest that mitigation action in line with the lower scenario RCP3-PD could contain Arctic temperature increase sufficiently that thawing of the permafrost area is limited to 9-23\% and the permafrost-carbon induced temperature increase does not exceed 0.04-0.16 degrees C by 2300.}, language = {en} } @article{ScheweLevermann2012, author = {Schewe, Jacob and Levermann, Anders}, title = {A statistically predictive model for future monsoon failure in India}, series = {Environmental research letters}, volume = {7}, journal = {Environmental research letters}, number = {4}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1748-9326}, doi = {10.1088/1748-9326/7/4/044023}, pages = {9}, year = {2012}, abstract = {Indian monsoon rainfall is vital for a large share of the world's population. Both reliably projecting India's future precipitation and unraveling abrupt cessations of monsoon rainfall found in paleorecords require improved understanding of its stability properties. While details of monsoon circulations and the associated rainfall are complex, full-season failure is dominated by large-scale positive feedbacks within the region. Here we find that in a comprehensive climate model, monsoon failure is possible but very rare under pre-industrial conditions, while under future warming it becomes much more frequent. We identify the fundamental intraseasonal feedbacks that are responsible for monsoon failure in the climate model, relate these to observational data, and build a statistically predictive model for such failure. This model provides a simple dynamical explanation for future changes in the frequency distribution of seasonal mean all-Indian rainfall. Forced only by global mean temperature and the strength of the Pacific Walker circulation in spring, it reproduces the trend as well as the multidecadal variability in the mean and skewness of the distribution, as found in the climate model. The approach offers an alternative perspective on large-scale monsoon variability as the result of internal instabilities modulated by pre-seasonal ambient climate conditions.}, language = {en} } @article{ScheweLevermannCheng2012, author = {Schewe, Jacob and Levermann, Anders and Cheng, Hai}, title = {A critical humidity threshold for monsoon transitions}, series = {Climate of the past : an interactive open access journal of the European Geosciences Union}, volume = {8}, journal = {Climate of the past : an interactive open access journal of the European Geosciences Union}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1814-9324}, doi = {10.5194/cp-8-535-2012}, pages = {535 -- 544}, year = {2012}, abstract = {Monsoon systems around the world are governed by the so-called moisture-advection feedback. Here we show that, in a minimal conceptual model, this feedback implies a critical threshold with respect to the atmospheric specific humidity q(o) over the ocean adjacent to the monsoon region. If q(o) falls short of this critical value q(o)(c), monsoon rainfall over land cannot be sustained. Such a case could occur if evaporation from the ocean was reduced, e.g. due to low sea surface temperatures. Within the restrictions of the conceptual model, we estimate q(o)(c) from present-day reanalysis data for four major monsoon systems, and demonstrate how this concept can help understand abrupt variations in monsoon strength on orbital timescales as found in proxy records.}, language = {en} } @article{MengelLevermannSchleussneretal.2012, author = {Mengel, Matthias and Levermann, Anders and Schleussner, Carl-Friedrich and Born, Andreas}, title = {Enhanced Atlantic subpolar gyre variability through baroclinic threshold in a coarse resolution model}, series = {Earth system dynamics}, volume = {3}, journal = {Earth system dynamics}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2190-4979}, doi = {10.5194/esd-3-189-2012}, pages = {189 -- 197}, year = {2012}, abstract = {Direct observations, satellite measurements and paleo records reveal strong variability in the Atlantic subpolar gyre on various time scales. Here we show that variations of comparable amplitude can only be simulated in a coupled climate model in the proximity of a dynamical threshold. The threshold and the associated dynamic response is due to a positive feedback involving increased salt transport in the subpolar gyre and enhanced deep convection in its centre. A series of sensitivity experiments is performed with a coarse resolution ocean general circulation model coupled to a statistical-dynamical atmosphere model which in itself does not produce atmospheric variability. To simulate the impact of atmospheric variability, the model system is perturbed with freshwater forcing of varying, but small amplitude and multi-decadal to centennial periodicities and observational variations in wind stress. While both freshwater and wind-stress-forcing have a small direct effect on the strength of the subpolar gyre, the magnitude of the gyre's response is strongly increased in the vicinity of the threshold. Our results indicate that baroclinic self-amplification in the North Atlantic ocean can play an important role in presently observed SPG variability and thereby North Atlantic climate variability on multi-decadal scales.}, language = {en} }