Search Results

Now showing 1 - 10 of 14
  • Item
    A generic pixel-to-point comparison for simulated large-scale ecosystem properties and ground-based observations: An example from the Amazon region
    (Katlenburg-Lindau : Copernicus, 2018) Rammig, Anja; Heinke, Jens; Hofhansl, Florian; Verbeeck, Hans; Baker, Timothy R.; Christoffersen, Bradley; Ciais, Philippe; De Deurwaerder, Hannes; Fleischer, Katrin; Galbraith, David; Guimberteau, Matthieu; Huth, Andreas; Johnson, Michelle; Krujit, Bart; Langerwisch, Fanny; Meir, Patrick; Papastefanou, Phillip; Sampaio, Gilvan; Thonicke, Kirsten; von Randow, Celso; Zang, Christian; Rödig, Edna
    Comparing model output and observed data is an important step for assessing model performance and quality of simulation results. However, such comparisons are often hampered by differences in spatial scales between local point observations and large-scale simulations of grid cells or pixels. In this study, we propose a generic approach for a pixel-to-point comparison and provide statistical measures accounting for the uncertainty resulting from landscape variability and measurement errors in ecosystem variables. The basic concept of our approach is to determine the statistical properties of small-scale (within-pixel) variability and observational errors, and to use this information to correct for their effect when large-scale area averages (pixel) are compared to small-scale point estimates. We demonstrate our approach by comparing simulated values of aboveground biomass, woody productivity (woody net primary productivity, NPP) and residence time of woody biomass from four dynamic global vegetation models (DGVMs) with measured inventory data from permanent plots in the Amazon rainforest, a region with the typical problem of low data availability, potential scale mismatch and thus high model uncertainty. We find that the DGVMs under- and overestimate aboveground biomass by 25% and up to 60%, respectively. Our comparison metrics provide a quantitative measure for model-data agreement and show moderate to good agreement with the region-wide spatial biomass pattern detected by plot observations. However, all four DGVMs overestimate woody productivity and underestimate residence time of woody biomass even when accounting for the large uncertainty range of the observational data. This is because DGVMs do not represent the relation between productivity and residence time of woody biomass correctly. Thus, the DGVMs may simulate the correct large-scale patterns of biomass but for the wrong reasons. We conclude that more information about the underlying processes driving biomass distribution are necessary to improve DGVMs. Our approach provides robust statistical measures for any pixel-to-point comparison, which is applicable for evaluation of models and remote-sensing products.
  • Item
    A protocol for an intercomparison of biodiversity and ecosystem services models using harmonized land-use and climate scenarios
    (Katlenburg-Lindau : Copernicus, 2018) Kim, HyeJin; Rosa, Isabel M. D.; Alkemade, Rob; Leadley, Paul; Hurtt, George; Popp, Alexander; van Vuuren, Detlef P.; Anthoni, Peter; Arneth, Almut; Baisero, Daniele; Caton, Emma; Chaplin-Kramer, Rebecca; Chini, Louise; De Palma, Adriana; Di Fulvio, Fulvio; Di Marco, Moreno; Espinoza, Felipe; Ferrier, Simon; Fujimori, Shinichiro; Gonzalez, Ricardo E.; Gueguen, Maya; Guerra, Carlos; Harfoot, Mike; Harwood, Thomas D.; Hasegawa, Tomoko; Haverd, Vanessa; Havlík, Petr; Hellweg, Stefanie; Hill, Samantha L. L.; Hirata, Akiko; Hoskins, Andrew J.; Janse, Jan H.; Jetz, Walter; Johnson, Justin A.; Krause, Andreas; Leclère, David; Martins, Ines S.; Matsui, Tetsuya; Merow, Cory; Obersteiner, Michael; Ohashi, Haruka; Poulter, Benjamin; Purvis, Andy; Quesada, Benjamin; Rondinini, Carlo; Schipper, Aafke M.; Sharp, Richard; Takahashi, Kiyoshi; Thuiller, Wilfried; Titeux, Nicolas; Visconti, Piero; Ware, Christopher; Wolf, Florian; Pereira, Henrique M.
    To support the assessments of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES), the IPBES Expert Group on Scenarios and Models is carrying out an intercomparison of biodiversity and ecosystem services models using harmonized scenarios (BES-SIM). The goals of BES-SIM are (1) to project the global impacts of land-use and climate change on biodiversity and ecosystem services (i.e., nature's contributions to people) over the coming decades, compared to the 20th century, using a set of common metrics at multiple scales, and (2) to identify model uncertainties and research gaps through the comparisons of projected biodiversity and ecosystem services across models. BES-SIM uses three scenarios combining specific Shared Socio-economic Pathways (SSPs) and Representative Concentration Pathways (RCPs)-SSP1xRCP2.6, SSP3xRCP6.0, SSP5xRCP8.6-to explore a wide range of land-use change and climate change futures. This paper describes the rationale for scenario selection, the process of harmonizing input data for land use, based on the second phase of the Land Use Harmonization Project (LUH2), and climate, the biodiversity and ecosystem services models used, the core simulations carried out, the harmonization of the model output metrics, and the treatment of uncertainty. The results of this collaborative modeling project will support the ongoing global assessment of IPBES, strengthen ties between IPBES and the Intergovernmental Panel on Climate Change (IPCC) scenarios and modeling processes, advise the Convention on Biological Diversity (CBD) on its development of a post-2020 strategic plans and conservation goals, and inform the development of a new generation of nature-centred scenarios.
  • Item
    The GGCMI Phase 2 emulators: Global gridded crop model responses to changes in CO2, temperature, water, and nitrogen (version 1.0)
    (Katlenburg-Lindau : Copernicus, 2020) Franke, James A.; Müller, Christoph; Elliott, Joshua; Ruane, Alex C.; Jägermeyr, Jonas; Snyder, Abigail; Dury, Marie; Falloon, Pete D.; Folberth, Christian; François, Louis; Hank, Tobias; Izaurralde, R. Cesar; Jacquemin, Ingrid; Jones, Curtis; Li, Michelle; Liu, Wenfeng; Olin, Stefan; Phillips, Meridel; Pugh, Thomas A. M.; Reddy, Ashwan; Williams, Karina; Wang, Ziwei; Zabel, Florian; Moyer, Elisabeth J.
    Statistical emulation allows combining advantageous features of statistical and process-based crop models for understanding the effects of future climate changes on crop yields. We describe here the development of emulators for nine process-based crop models and five crops using output from the Global Gridded Model Intercomparison Project (GGCMI) Phase 2. The GGCMI Phase 2 experiment is designed with the explicit goal of producing a structured training dataset for emulator development that samples across four dimensions relevant to crop yields: Atmospheric carbon dioxide (CO2) concentrations, temperature, water supply, and nitrogen inputs (CTWN). Simulations are run under two different adaptation assumptions: That growing seasons shorten in warmer climates, and that cultivar choice allows growing seasons to remain fixed. The dataset allows emulating the climatological-mean yield response of all models with a simple polynomial in mean growing-season values. Climatological-mean yields are a central metric in climate change impact analysis; we show here that they can be captured without relying on interannual variations. In general, emulation errors are negligible relative to differences across crop models or even across climate model scenarios; errors become significant only in some marginal lands where crops are not currently grown. We demonstrate that the resulting GGCMI emulators can reproduce yields under realistic future climate simulations, even though the GGCMI Phase 2 dataset is constructed with uniform CTWN offsets, suggesting that the effects of changes in temperature and precipitation distributions are small relative to those of changing means. The resulting emulators therefore capture relevant crop model responses in a lightweight, computationally tractable form, providing a tool that can facilitate model comparison, diagnosis of interacting factors affecting yields, and integrated assessment of climate impacts. © 2020 EDP Sciences. All rights reserved.
  • Item
    Modelling the role of livestock grazing in C and N cycling in grasslands with LPJmL5.0-grazing
    (Katlenburg-Lindau : Copernicus, 2023) Heinke, Jens; Rolinski, Susanne; Müller, Christoph
    To represent the impact of grazing livestock on carbon (C) and nitrogen (N) dynamics in grasslands, we implement a livestock module into LPJmL5.0-tillage, a global vegetation and crop model with explicit representation of managed grasslands and pastures, forming LPJmL5.0-grazing. The livestock module uses lactating dairy cows as a generic representation of grazing livestock. The new module explicitly accounts for forage quality in terms of dry-matter intake and digestibility using relationships derived from compositional analyses for different forages. Partitioning of N into milk, feces, and urine as simulated by the new livestock module shows very good agreement with observation-based relationships reported in the literature. Modelled C and N dynamics depend on forage quality (C:N ratios in grazed biomass), forage quantity, livestock densities, manure or fertilizer inputs, soil, atmospheric CO2 concentrations, and climate conditions. Due to the many interacting relationships, C sequestration, GHG emissions, N losses, and livestock productivity show substantial variation in space and across livestock densities. The improved LPJmL5.0-grazing model can now assess the effects of livestock grazing on C and N stocks and fluxes in grasslands. It can also provide insights about the spatio-temporal variability of grassland productivity and about the trade-offs between livestock production and environmental impacts.
  • Item
    SURFER v2.0: a flexible and simple model linking anthropogenic CO2 emissions and solar radiation modification to ocean acidification and sea level rise
    (Katlenburg-Lindau : Copernicus, 2022) Martínez Montero, Marina; Crucifix, Michel; Couplet, Victor; Brede, Nuria; Botta, Nicola
    We present SURFER, a novel reduced model for estimating the impact of CO2 emissions and solar radiation modification options on sea level rise and ocean acidification over timescales of several thousands of years. SURFER has been designed for the analysis of CO2 emission and solar radiation modification policies, for supporting the computation of optimal (CO2 emission and solar radiation modification) policies and for the study of commitment and responsibility under uncertainty. The model is based on a combination of conservation laws for the masses of atmospheric and oceanic carbon and for the oceanic temperature anomalies, and of ad-hoc parameterisations for the different sea level rise contributors: ice sheets, glaciers and ocean thermal expansion. It consists of 9 loosely coupled ordinary differential equations, is understandable, fast and easy to modify and calibrate. It reproduces the results of more sophisticated, high-dimensional earth system models on timescales up to millennia.
  • Item
    LandInG 1.0: a toolbox to derive input datasets for terrestrial ecosystem modelling at variable resolutions from heterogeneous sources
    (Katlenburg-Lindau : Copernicus, 2023) Ostberg, Sebastian; Müller, Christoph; Heinke, Jens; Schaphoff, Sibyll
    We present the Land Input Generator (LandInG) version 1.0, a new toolbox for generating input datasets for terrestrial ecosystem models (TEMs) from diverse and partially conflicting data sources. While LandInG 1.0 is applicable to process data for any TEM, it is developed specifically for the open-source dynamic global vegetation, hydrology, and crop growth model LPJmL (Lund-Potsdam-Jena with managed Land). The toolbox documents the sources and processing of data to model inputs and allows for easy changes to the spatial resolution. It is designed to make inconsistencies between different sources of data transparent so that users can make their own decisions on how to resolve these should they not be content with the default assumptions made here. As an example, we use the toolbox to create input datasets at 5 and 30 arcmin spatial resolution covering land, country, and region masks, soil, river networks, freshwater reservoirs, irrigation water distribution networks, crop-specific annual land use, fertilizer, and manure application. We focus on the toolbox describing the data processing rather than only publishing the datasets as users may want to make different choices for reconciling inconsistencies, aggregation, spatial extent, or similar. Also, new data sources or new versions of existing data become available continuously, and the toolbox approach allows for incorporating new data to stay up to date.
  • Item
    Implementing the nitrogen cycle into the dynamic global vegetation, hydrology, and crop growth model LPJmL (version 5.0)
    (Katlenburg-Lindau : Copernicus, 2018) von Bloh, Werner; Schaphoff, Sibyll; Müller, Christoph; Rolinski, Susanne; Waha, Katharina; Zaehle, Sönke
    The well-established dynamical global vegetation, hydrology, and crop growth model LPJmL is extended with a terrestrial nitrogen cycle to account for nutrient limitations. In particular, processes of soil nitrogen dynamics, plant uptake, nitrogen allocation, response of photosynthesis and maintenance respiration to varying nitrogen concentrations in plant organs, and agricultural nitrogen management are included in the model. All new model features are described in full detail and the results of a global simulation of the historic past (1901-2009) are presented for evaluation of the model performance. We find that the implementation of nitrogen limitation significantly improves the simulation of global patterns of crop productivity. Regional differences in crop productivity, which had to be calibrated via a scaling of the maximum leaf area index, can now largely be reproduced by the model, except for regions where fertilizer inputs and climate conditions are not the yield-limiting factors. Furthermore, it can be shown that land use has a strong influence on nitrogen losses, increasing leaching by 93 %.
  • Item
    The Earth system model CLIMBER-X v1.0 - Part 2: The global carbon cycle
    (Katlenburg-Lindau : Copernicus, 2023) Willeit, Matteo; Ilyina, Tatiana; Liu, Bo; Heinze, Christoph; Perrette, Mahé; Heinemann, Malte; Dalmonech, Daniela; Brovkin, Victor; Munhoven, Guy; Börker, Janine; Hartmann, Jens; Romero-Mujalli, Gibran; Ganopolski, Andrey
    The carbon cycle component of the newly developed Earth system model of intermediate complexity CLIMBER-X is presented. The model represents the cycling of carbon through the atmosphere, vegetation, soils, seawater and marine sediments. Exchanges of carbon with geological reservoirs occur through sediment burial, rock weathering and volcanic degassing. The state-of-the-art HAMOCC6 model is employed to simulate ocean biogeochemistry and marine sediment processes. The land model PALADYN simulates the processes related to vegetation and soil carbon dynamics, including permafrost and peatlands. The dust cycle in the model allows for an interactive determination of the input of the micro-nutrient iron into the ocean. A rock weathering scheme is implemented in the model, with the weathering rate depending on lithology, runoff and soil temperature. CLIMBER-X includes a simple representation of the methane cycle, with explicitly modelled natural emissions from land and the assumption of a constant residence time of CH4 in the atmosphere. Carbon isotopes 13C and 14C are tracked through all model compartments and provide a useful diagnostic for model-data comparison. A comprehensive evaluation of the model performance for the present day and the historical period shows that CLIMBER-X is capable of realistically reproducing the historical evolution of atmospheric CO2 and CH4 but also the spatial distribution of carbon on land and the 3D structure of biogeochemical ocean tracers. The analysis of model performance is complemented by an assessment of carbon cycle feedbacks and model sensitivities compared to state-of-the-art Coupled Model Intercomparison Project Phase 6 (CMIP6) models. Enabling an interactive carbon cycle in CLIMBER-X results in a relatively minor slow-down of model computational performance by ∼ 20 % compared to a throughput of ∼ 10 000 simulation years per day on a single node with 16 CPUs on a high-performance computer in a climate-only model set-up. CLIMBER-X is therefore well suited to investigating the feedbacks between climate and the carbon cycle on temporal scales ranging from decades to >100000 years.
  • Item
    Grounding-line flux formula applied as a flux condition in numerical simulations fails for buttressed Antarctic ice streams
    (Katlenburg-Lindau : Copernicus, 2018) Reese, Ronja; Winkelmann, Ricarda; Gudmundsson, G. Hilmar
    Currently, several large-scale ice-flow models impose a condition on ice flux across grounding lines using an analytically motivated parameterisation of grounding-line flux. It has been suggested that employing this analytical expression alleviates the need for highly resolved computational domains around grounding lines of marine ice sheets. While the analytical flux formula is expected to be accurate in an unbuttressed flow-line setting, its validity has hitherto not been assessed for complex and realistic geometries such as those of the Antarctic Ice Sheet. Here the accuracy of this analytical flux formula is tested against an optimised ice flow model that uses a highly resolved computational mesh around the Antarctic grounding lines. We find that when applied to the Antarctic Ice Sheet the analytical expression provides inaccurate estimates of ice fluxes for almost all grounding lines. Furthermore, in many instances direct application of the analytical formula gives rise to unphysical complex-valued ice fluxes. We conclude that grounding lines of the Antarctic Ice Sheet are, in general, too highly buttressed for the analytical parameterisation to be of practical value for the calculation of grounding-line fluxes.
  • Item
    The future sea-level contribution of the Greenland ice sheet: A multi-model ensemble study of ISMIP6
    (Katlenburg-Lindau : Copernicus, 2020) Goelzer, Heiko; Nowicki, Sophie; Payne, Anthony; Larour, Eric; Seroussi, Helene; Lipscomb, William H.; Gregory, Jonathan; Abe-Ouchi, Ayako; Shepherd, Andrew; Simon, Erika; Agosta, Cécile; Alexander, Patrick; Aschwanden, Andy; Barthel, Alice; Calov, Reinhard; Chambers, Christopher; Choi, Youngmin; Cuzzone, Joshua; Dumas, Christophe; Edwards, Tamsin; Felikson, Denis; Fettweis, Xavier; Golledge, Nicholas R.; Greve, Ralf; Humbert, Angelika; Huybrechts, Philippe; Le clec'h, Sebastien; Lee, Victoria; Leguy, Gunter; Little, Chris; Lowry, Daniel P.; Morlighem, Mathieu; Nias, Isabel; Quiquet, Aurelien; Rückamp, Martin; Schlegel, Nicole-Jeanne; Slater, Donald A.; Smith, Robin S.; Straneo, Fiammetta; Tarasov, Lev; van de Wal, Roderik; van den Broeke, Michiel
    The Greenland ice sheet is one of the largest contributors to global mean sea-level rise today and is expected to continue to lose mass as the Arctic continues to warm. The two predominant mass loss mechanisms are increased surface meltwater run-off and mass loss associated with the retreat of marine-terminating outlet glaciers. In this paper we use a large ensemble of Greenland ice sheet models forced by output from a representative subset of the Coupled Model Intercomparison Project (CMIP5) global climate models to project ice sheet changes and sea-level rise contributions over the 21st century. The simulations are part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6).We estimate the sea-level contribution together with uncertainties due to future climate forcing, ice sheet model formulations and ocean forcing for the two greenhouse gas concentration scenarios RCP8.5 and RCP2.6. The results indicate that the Greenland ice sheet will continue to lose mass in both scenarios until 2100, with contributions of 90-50 and 32-17mm to sea-level rise for RCP8.5 and RCP2.6, respectively. The largest mass loss is expected from the south-west of Greenland, which is governed by surface mass balance changes, continuing what is already observed today. Because the contributions are calculated against an unforced control experiment, these numbers do not include any committed mass loss, i.e. mass loss that would occur over the coming century if the climate forcing remained constant. Under RCP8.5 forcing, ice sheet model uncertainty explains an ensemble spread of 40 mm, while climate model uncertainty and ocean forcing uncertainty account for a spread of 36 and 19 mm, respectively. Apart from those formally derived uncertainty ranges, the largest gap in our knowledge is about the physical understanding and implementation of the calving process, i.e. the interaction of the ice sheet with the ocean. © Author(s) 2020.