Search Results

Now showing 1 - 10 of 51
  • Item
    Description and evaluation of the process-based forest model 4C v2.2 at four European forest sites
    (Katlenburg-Lindau : Copernicus, 2020) Lasch-Born, Petra; Suckow, Felicitas; Reyer, Christopher P. O.; Gutsch, Martin; Kollas, Chris; Badeck, Franz-Werner; Bugmann, Harald K. M.; Grote, Rüdiger; Fürstenau, Cornelia; Lindner, Marcus; Schaber, Jörg
    The process-based model 4C (FORESEE) has been developed over the past 20 years to study climate impacts on forests and is now freely available as an open-source tool. The objective of this paper is to provide a comprehensive description of this 4C version (v2.2) for scientific users of the model and to present an evaluation of 4C at four different forest sites across Europe. The evaluation focuses on forest growth as well as carbon (net ecosystem exchange, gross primary production), water (actual evapotranspiration, soil water content), and heat fluxes (soil temperature) using data from the PROFOUND database. We applied different evaluation metrics and compared the daily, monthly, and annual variability of observed and simulated values. The ability to reproduce forest growth (stem diameter and biomass) differs from site to site and is best for a pine stand in Germany (Peitz, model efficiency ME=0.98). 4C is able to reproduce soil temperature at different depths in Sorø and Hyytiälä with good accuracy (for all soil depths ME > 0.8). The dynamics in simulating carbon and water fluxes are well captured on daily and monthly timescales (0.51 < ME < 0.983) but less so on an annual timescale (ME < 0). This model–data mismatch is possibly due to the accumulation of errors because of processes that are missing or represented in a very general way in 4C but not with enough specific detail to cover strong, site-specific dependencies such as ground vegetation growth. These processes need to be further elaborated to improve the projections of climate change on forests. We conclude that, despite shortcomings, 4C is widely applicable, reliable, and therefore ready to be released to the scientific community to use and further develop the model.
  • Item
    The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) contribution to C4MIP: Quantifying committed climate changes following zero carbon emissions
    (Katlenburg-Lindau : Copernicus, 2019) Jones, Chris D.; Frölicher, Thomas L.; Koven, Charles; MacDougall, Andrew H.; Matthews, H. Damon; Zickfeld, Kirsten; Rogelj, Joeri; Tokarska, Katarzyna B.; Gillett, Nathan P.; Ilyina, Tatiana; Meinshausen, Malte; Mengis, Nadine; Séférian, Roland; Eby, Michael; Burger, Friedrich A.
    The amount of additional future temperature change following a complete cessation of CO2 emissions is a measure of the unrealized warming to which we are committed due to CO2 already emitted to the atmosphere. This “zero emissions commitment” (ZEC) is also an important quantity when estimating the remaining carbon budget – a limit on the total amount of CO2 emissions consistent with limiting global mean temperature at a particular level. In the recent IPCC Special Report on Global Warming of 1.5 ∘C, the carbon budget framework used to calculate the remaining carbon budget for 1.5 ∘C included the assumption that the ZEC due to CO2 emissions is negligible and close to zero. Previous research has shown significant uncertainty even in the sign of the ZEC. To close this knowledge gap, we propose the Zero Emissions Commitment Model Intercomparison Project (ZECMIP), which will quantify the amount of unrealized temperature change that occurs after CO2 emissions cease and investigate the geophysical drivers behind this climate response. Quantitative information on ZEC is a key gap in our knowledge, and one that will not be addressed by currently planned CMIP6 simulations, yet it is crucial for verifying whether carbon budgets need to be adjusted to account for any unrealized temperature change resulting from past CO2 emissions. We request only one top-priority simulation from comprehensive general circulation Earth system models (ESMs) and Earth system models of intermediate complexity (EMICs) – a branch from the 1 % CO2 run with CO2 emissions set to zero at the point of 1000 PgC of total CO2 emissions in the simulation – with the possibility for additional simulations, if resources allow. ZECMIP is part of CMIP6, under joint sponsorship by C4MIP and CDRMIP, with associated experiment names to enable data submissions to the Earth System Grid Federation. All data will be published and made freely available.
  • Item
    Drivers of Pine Island Glacier speed-up between 1996 and 2016
    (Katlenburg-Lindau : Copernicus, 2021-1-7) De Rydt, Jan; Reese, Ronja; Paolo, Fernando S.; Gudmundsson, G. Hilmar
    Pine Island Glacier in West Antarctica is among the fastest changing glaciers worldwide. Over the last 2 decades, the glacier has lost in excess of a trillion tons of ice, or the equivalent of 3 mm of sea level rise. The ongoing changes are thought to have been triggered by ocean-induced thinning of its floating ice shelf, grounding line retreat, and the associated reduction in buttressing forces. However, other drivers of change, such as large-scale calving and changes in ice rheology and basal slipperiness, could play a vital, yet unquantified, role in controlling the ongoing and future evolution of the glacier. In addition, recent studies have shown that mechanical properties of the bed are key to explaining the observed speed-up. Here we used a combination of the latest remote sensing datasets between 1996 and 2016, data assimilation tools, and numerical perturbation experiments to quantify the relative importance of all processes in driving the recent changes in Pine Island Glacier dynamics. We show that (1) calving and ice shelf thinning have caused a comparable reduction in ice shelf buttressing over the past 2 decades; that (2) simulated changes in ice flow over a viscously deforming bed are only compatible with observations if large and widespread changes in ice viscosity and/or basal slipperiness are taken into account; and that (3) a spatially varying, predominantly plastic bed rheology can closely reproduce observed changes in flow without marked variations in ice-internal and basal properties. Our results demonstrate that, in addition to its evolving ice thickness, calving processes and a heterogeneous bed rheology play a key role in the contemporary evolution of Pine Island Glacier.
  • Item
    Modeling the response of Greenland outlet glaciers to global warming using a coupled flow line-plume model
    (Katlenburg-Lindau : Copernicus, 2019) Beckmann, Johanna; Perrette, Mahé; Beyer, Sebastian; Calov, Reinhard; Willeit, Matteo; Ganopolski, Andrey
    In recent decades, the Greenland Ice Sheet has experienced an accelerated mass loss, contributing to approximately 25 % of contemporary sea level rise (SLR). This mass loss is caused by increased surface melt over a large area of the ice sheet and by the thinning, retreat and acceleration of numerous Greenland outlet glaciers. The latter is likely connected to enhanced submarine melting that, in turn, can be explained by ocean warming and enhanced subglacial discharge. The mechanisms involved in submarine melting are not yet fully understood and are only simplistically incorporated in some models of the Greenland Ice Sheet. Here, we investigate the response of 12 representative Greenland outlet glaciers to atmospheric and oceanic warming using a coupled line–plume glacier–flow line model resolving one horizontal dimension. The model parameters have been tuned for individual outlet glaciers using present-day observational constraints. We then run the model from present to the year 2100, forcing the model with changes in surface mass balance and surface runoff from simulations with a regional climate model for the RCP8.5 scenario, and applying a linear ocean temperature warming with different rates of changes representing uncertainties in the CMIP5 model experiments for the same climate change scenario. We also use different initial temperature–salinity profiles obtained from direct measurements and from ocean reanalysis data. Using different combinations of submarine melting and calving parameters that reproduce the present-day state of the glaciers, we estimate uncertainties in the contribution to global SLR for individual glaciers. We also perform a sensitivity analysis of the three forcing factors (changes in surface mass balance, ocean temperature and subglacial discharge), which shows that the roles of the different forcing factors are diverse for individual glaciers. We find that changes in ocean temperature and subglacial discharge are of comparable importance for the cumulative contribution of all 12 glaciers to global SLR in the 21st century. The median range of the cumulative contribution to the global SLR for all 12 glaciers is about 18 mm (the glaciers' dynamic response to changes of all three forcing factors). Neglecting changes in ocean temperature and subglacial discharge (which control submarine melt) and investigating the response to changes in surface mass balance only leads to a cumulative contribution of 5 mm SLR. Thus, from the 18 mm we associate roughly 70 % with the glaciers' dynamic response to increased subglacial discharge and ocean temperature and the remaining 30 % (5 mm) to the response to increased surface mass loss. We also find a strong correlation (correlation coefficient 0.74) between present-day grounding line discharge and their future contribution to SLR in 2100. If the contribution of the 12 glaciers is scaled up to the total present-day discharge of Greenland, we estimate the midrange contribution of all Greenland glaciers to 21st-century SLR to be approximately 50 mm. This number adds to SLR derived from a stand-alone ice sheet model (880 mm) that does not resolve outlet glaciers and thus increases SLR by over 50 %. This result confirms earlier studies showing that the response of the outlet glaciers to global warming has to be taken into account to correctly assess the total contribution of Greenland to sea level change.
  • Item
    Timescales of outlet-glacier flow with negligible basal friction: Theory, observations and modeling
    (Katlenburg-Lindau : Copernicus, 2023) Feldmann, Johannes; Levermann, Anders
    The timescales of the flow and retreat of Greenland's and Antarctica's outlet glaciers and their potential instabilities are arguably the largest uncertainty in future sea-level projections. Here we derive a scaling relation that allows the comparison of the timescales of observed complex ice flow fields with geometric similarity. The scaling relation is derived under the assumption of fast, laterally confined, geometrically similar outlet-glacier flow over a slippery bed, i.e., with negligible basal friction. According to the relation, the time scaling of the outlet flow is determined by the product of the inverse of (1) the fourth power of the width-To-length ratio of its confinement, (2) the third power of the confinement depth and (3) the temperature-dependent ice softness. For the outflow at the grounding line of streams with negligible basal friction, this means that the volume flux is proportional to the ice softness and the bed depth, but goes with the fourth power of the gradient of the bed and with the fifth power of the width of the stream. We show that the theoretically derived scaling relation is supported by the observed velocity scaling of outlet glaciers across Greenland as well as by idealized numerical simulations of marine ice-sheet instabilities (MISIs) as found in Antarctica. Assuming that changes in the ice-flow velocity due to ice-dynamic imbalance are proportional to the equilibrium velocity, we combine the scaling relation with a statistical analysis of the topography of 13 MISI-prone Antarctic outlets. Under these assumptions, the timescales in response to a potential destabilization are fastest for Thwaites Glacier in West Antarctica and Mellor, Ninnis and Cook Glaciers in East Antarctica; between 16 and 67 times faster than for Pine Island Glacier. While the applicability of our results is limited by several strong assumptions, the utilization and potential further development of the presented scaling approach may help to constrain timescale estimates of outlet-glacier flow, augmenting the commonly exploited and comparatively computationally expensive approach of numerical modeling.
  • Item
    Potential yield simulated by global gridded crop models: using a process-based emulator to explain their differences
    (Katlenburg-Lindau : Copernicus, 2021-3-23) Ringeval, Bruno; Müller, Christoph; Pugh, Thomas A. M.; Mueller, Nathaniel D.; Ciais, Philippe; Folberth, Christian; Liu, Wenfeng; Debaeke, Philippe; Pellerin, Sylvain
    How global gridded crop models (GGCMs) differ in their simulation of potential yield and reasons for those differences have never been assessed. The GGCM Intercomparison (GGCMI) offers a good framework for this assessment. Here, we built an emulator (called SMM for simple mechanistic model) of GGCMs based on generic and simplified formalism. The SMM equations describe crop phenology by a sum of growing degree days, canopy radiation absorption by the Beer–Lambert law, and its conversion into aboveground biomass by a radiation use efficiency (RUE). We fitted the parameters of this emulator against gridded aboveground maize biomass at the end of the growing season simulated by eight different GGCMs in a given year (2000). Our assumption is that the simple set of equations of SMM, after calibration, could reproduce the response of most GGCMs so that differences between GGCMs can be attributed to the parameters related to processes captured by the emulator. Despite huge differences between GGCMs, we show that if we fit both a parameter describing the thermal requirement for leaf emergence by adjusting its value to each grid-point in space, as done by GGCM modellers following the GGCMI protocol, and a GGCM-dependent globally uniform RUE, then the simple set of equations of the SMM emulator is sufficient to reproduce the spatial distribution of the original aboveground biomass simulated by most GGCMs. The grain filling is simulated in SMM by considering a fixed-in-time fraction of net primary productivity allocated to the grains (frac) once a threshold in leaves number (nthresh) is reached. Once calibrated, these two parameters allow for the capture of the relationship between potential yield and final aboveground biomass of each GGCM. It is particularly important as the divergence among GGCMs is larger for yield than for aboveground biomass. Thus, we showed that the divergence between GGCMs can be summarized by the differences in a few parameters. Our simple but mechanistic model could also be an interesting tool to test new developments in order to improve the simulation of potential yield at the global scale.
  • Item
    Patterns in Mongolian nomadic household movement derived from GPS trajectories
    (New York, NY [u.a.] : Elsevier, 2020) Teickner, Henning; Knoth, Christian; Bartoschek, Thomas; Kraehnert, Kati; Vigh, Melinda; Purevtseren, Myagmartseren; Sugar, Munkhnaran; Pebesma, Edzer
    This paper presents an approach for a quantitative analysis of movement patterns of nomadic households based on GPS trajectories. We distributed GPS loggers to 400 Mongolian herder households who carried them over a 9-month period, continuously recording position data every 30min. A total of 142of the resulting trajectories fulfilled our data quality criteria and were considered during the analysis. Based on this data, we derive summary indicators describing key parameters of the households’ mobility including measures of distance and number of movements as well as shape characteristics of the trajectories. We conduct an explorative statistical analysis of these summary indicators to investigate patterns in the nomadic mobility. We identify three movement strategies based on the number of different campsite locations and the distances traveled between campsites. We also compare the results to the existing literature on the mobility of Mongolian herders. Our findings show that GPS-based studies present a suitable framework to quantitatively analyze different movement strategies of nomadic herders.
  • Item
    Simulating the effect of tillage practices with the global ecosystem model LPJmL (version 5.0-tillage)
    (Katlenburg-Lindau : Copernicus, 2019) Lutz, Femke; Herzfeld, Tobias; Heinke, Jens; Rolinski, Susanne; Schaphoff, Sibyll; von Bloh, Werner; Stoorvogel, Jetse J.; Müller, Christoph
    The effects of tillage on soil properties, crop productivity, and global greenhouse gas emissions have been discussed in the last decades. Global ecosystem models have limited capacity to simulate the various effects of tillage. With respect to the decomposition of soil organic matter, they either assume a constant increase due to tillage or they ignore the effects of tillage. Hence, they do not allow for analysing the effects of tillage and cannot evaluate, for example, reduced tillage or no tillage (referred to here as “no-till”) practises as mitigation practices for climate change. In this paper, we describe the implementation of tillage-related practices in the global ecosystem model LPJmL. The extended model is evaluated against reported differences between tillage and no-till management on several soil properties. To this end, simulation results are compared with published meta-analyses on tillage effects. In general, the model is able to reproduce observed tillage effects on global, as well as regional, patterns of carbon and water fluxes. However, modelled N fluxes deviate from the literature values and need further study. The addition of the tillage module to LPJmL5 opens up opportunities to assess the impact of agricultural soil management practices under different scenarios with implications for agricultural productivity, carbon sequestration, greenhouse gas emissions, and other environmental indicators.
  • Item
    Stabilizing effect of mélange buttressing on the marine ice-cliff instability of the West Antarctic Ice Sheet
    (Katlenburg-Lindau : Copernicus, 2022) Schlemm, Tanja; Feldmann, Johannes; Winkelmann, Ricarda; Levermann, Anders
    Owing to global warming and particularly high regional ocean warming, both Thwaites and Pine Island Glaciers in the Amundsen region of the Antarctic Ice Sheet could lose their buttressing ice shelves over time. We analyse the possible consequences using the parallel ice sheet model (PISM), applying a simple cliff-calving parameterization and an ice mélange-buttressing model. We find that the instantaneous loss of ice-shelf buttressing, due to enforced ice-shelf melting, initiates grounding-line retreat and triggers marine ice sheet instability (MISI). As a consequence, the grounding line progresses into the interior of the West Antarctic Ice Sheet and leads to a sea level contribution of 0.6m within 100a. By subjecting the exposed ice cliffs to cliff calving using our simplified parameterization, we also analyse marine ice cliff instability (MICI). In our simulations it can double or even triple the sea level contribution depending on the only loosely constrained parameter that determines the maximum cliff-calving rate. The speed of MICI depends on this upper bound of the calving rate, which is given by the ice mélange buttressing the glacier. However, stabilization of MICI may occur for geometric reasons. Because the embayment geometry changes as MICI advances into the interior of the ice sheet, the upper bound on calving rates is reduced and the progress of MICI is slowed down. Although we cannot claim that our simulations bear relevant quantitative estimates of the effect of ice-mélange buttressing on MICI, the mechanism has the potential to stop the instability. Further research is needed to evaluate its role for the past and future evolution of the Antarctic Ice Sheet.
  • Item
    The role of history and strength of the oceanic forcing in sea level projections from Antarctica with the Parallel Ice Sheet Model
    (Katlenburg-Lindau : Copernicus, 2020) Reese, Ronja; Levermann, Anders; Albrecht, Torsten; Seroussi, Hélène; Winkelmann, Ricarda
    Mass loss from the Antarctic Ice Sheet constitutes the largest uncertainty in projections of future sea level rise. Ocean-driven melting underneath the floating ice shelves and subsequent acceleration of the inland ice streams are the major reasons for currently observed mass loss from Antarctica and are expected to become more important in the future. Here we show that for projections of future mass loss from the Antarctic Ice Sheet, it is essential (1) to better constrain the sensitivity of sub-shelf melt rates to ocean warming and (2) to include the historic trajectory of the ice sheet. In particular, we find that while the ice sheet response in simulations using the Parallel Ice Sheet Model is comparable to the median response of models in three Antarctic Ice Sheet Intercomparison projects – initMIP, LARMIP-2 and ISMIP6 – conducted with a range of ice sheet models, the projected 21st century sea level contribution differs significantly depending on these two factors. For the highest emission scenario RCP8.5, this leads to projected ice loss ranging from 1.4 to 4.0 cm of sea level equivalent in simulations in which ISMIP6 ocean forcing drives the PICO ocean box model where parameter tuning leads to a comparably low sub-shelf melt sensitivity and in which no surface forcing is applied. This is opposed to a likely range of 9.1 to 35.8 cm using the exact same initial setup, but emulated from the LARMIP-2 experiments with a higher melt sensitivity, even though both projects use forcing from climate models and melt rates are calibrated with previous oceanographic studies. Furthermore, using two initial states, one with a previous historic simulation from 1850 to 2014 and one starting from a steady state, we show that while differences between the ice sheet configurations in 2015 seem marginal at first sight, the historic simulation increases the susceptibility of the ice sheet to ocean warming, thereby increasing mass loss from 2015 to 2100 by 5 % to 50 %. Hindcasting past ice sheet changes with numerical models would thus provide valuable tools to better constrain projections. Our results emphasize that the uncertainty that arises from the forcing is of the same order of magnitude as the ice dynamic response for future sea level projections.