2022, Willeit, Matteo, Ganopolski, Andrey, Robinson, Alexander, Edwards, Neil R.

The newly developed fast Earth system model CLIMBER-X is presented. The climate component of CLIMBER-X consists of a 2.5-D semi-empirical statistical-dynamical atmosphere model, a 3-D frictional-geostrophic ocean model, a dynamic-thermodynamic sea ice model and a land surface model. All the model components are discretized on a regular lat-long grid with a horizontal resolution of 5 ° ×5 °. The model has a throughput of ° ∼ 10 000 simulation years per day on a single node with 16 CPUs on a high-performance computer and is designed to simulate the evolution of the Earth system on temporal scales ranging from decades to >100000 years. A comprehensive evaluation of the model performance for the present day and the historical period shows that CLIMBER-X is capable of realistically reproducing many observed climate characteristics, with results that generally lie within the range of state-of-the-art general circulation models. The analysis of model performance is complemented by a thorough assessment of climate feedbacks and model sensitivities to changes in external forcings and boundary conditions. Limitations and applicability of the model are critically discussed. CLIMBER-X also includes a detailed representation of the global carbon cycle and is coupled to an ice sheet model, which will be described in separate papers. CLIMBER-X is available as open-source code and is expected to be a useful tool for studying past climate changes and for the investigation of the long-term future evolution of the climate.

2022, Kikstra, Jarmo S., Nicholls, Zebedee R. J., Smith, Christopher J., Lewis, Jared, Lamboll, Robin D., Byers, Edward, Sandstad, Marit, Meinshausen, Malte, Gidden, Matthew J., Rogelj, Joeri, Kriegler, Elmar, Peters, Glen P., Fuglestvedt, Jan S., Skeie, Ragnhild B., Samset, Bjørn H., Wienpahl, Laura, van Vuuren, Detlef P., van der Wijst, Kaj-Ivar, Al Khourdajie, Alaa, Forster, Piers M., Reisinger, Andy, Schaeffer, Roberto, Riahi, Keywan

While the Intergovernmental Panel on Climate Change (IPCC) physical science reports usually assess a handful of future scenarios, the Working Group III contribution on climate mitigation to the IPCC's Sixth Assessment Report (AR6 WGIII) assesses hundreds to thousands of future emissions scenarios. A key task in WGIII is to assess the global mean temperature outcomes of these scenarios in a consistent manner, given the challenge that the emissions scenarios from different integrated assessment models (IAMs) come with different sectoral and gas-to-gas coverage and cannot all be assessed consistently by complex Earth system models. In this work, we describe the "climate-assessment"workflow and its methods, including infilling of missing emissions and emissions harmonisation as applied to 1202 mitigation scenarios in AR6 WGIII. We evaluate the global mean temperature projections and effective radiative forcing (ERF) characteristics of climate emulators FaIRv1.6.2 and MAGICCv7.5.3 and use the CICERO simple climate model (CICERO-SCM) for sensitivity analysis. We discuss the implied overshoot severity of the mitigation pathways using overshoot degree years and look at emissions and temperature characteristics of scenarios compatible with one possible interpretation of the Paris Agreement. We find that the lowest class of emissions scenarios that limit global warming to "1.5 ° C (with a probability of greater than 50 %) with no or limited overshoot"includes 97 scenarios for MAGICCv7.5.3 and 203 for FaIRv1.6.2. For the MAGICCv7.5.3 results, "limited overshoot"typically implies exceedance of median temperature projections of up to about 0.1 ° C for up to a few decades before returning to below 1.5 ° C by or before the year 2100. For more than half of the scenarios in this category that comply with three criteria for being "Paris-compatible", including net-zero or net-negative greenhouse gas (GHG) emissions, median temperatures decline by about 0.3-0.4 ° C after peaking at 1.5-1.6 ° C in 2035-2055. We compare the methods applied in AR6 with the methods used for SR1.5 and discuss their implications. This article also introduces a "climate-assessment"Python package which allows for fully reproducing the IPCC AR6 WGIII temperature assessment. This work provides a community tool for assessing the temperature outcomes of emissions pathways and provides a basis for further work such as extending the workflow to include downscaling of climate characteristics to a regional level and calculating impacts.

2023, Niebsch, Jenny, Bloh, Werner von, Thonicke, Kirsten, Ramlau, Ronny

The increasing impacts of climate change require strategies for climate adaptation. Dynamic global vegetation models (DGVMs) are one type of multi-sectorial impact model with which the effects of multiple interacting processes in the terrestrial biosphere under climate change can be studied. The complexity of DGVMs is increasing as more and more processes, especially for plant physiology, are implemented. Therefore, there is a growing demand for increasing the computational performance of the underlying algorithms as well as ensuring their numerical accuracy. One way to approach this issue is to analyse the routines which have the potential for improved computational efficiency and/or increased accuracy when applying sophisticated mathematical methods. In this paper, the Farquhar–Collatz photosynthesis model under water stress as implemented in the Lund–Potsdam–Jena managed Land DGVM (4.0.002) was examined. We additionally tested the uncertainty of most important parameter of photosynthesis as an additional approach to improve model quality. We found that the numerical solution of a nonlinear equation, so far solved with the bisection method, could be significantly improved by using Newton’s method instead. The latter requires the computation of the derivative of the underlying function which is presented. Model simulations show a significantly lower number of iterations to solve the equation numerically and an overall run time reduction of the model of about 16 % depending on the chosen accuracy. Increasing the parameters θ and αC3 by 10 %, respectively, while keeping all other parameters at their original value, increased global gross primary production (GPP) by 2.384 and 9.542 GtC yr−1, respectively. The Farquhar–Collatz photosynthesis model forms the core component in many DGVMs and land surface models. An update in the numerical solution of the nonlinear equation in connection with adjusting globally important parameters to best known values can therefore be applied to similar photosynthesis models. Furthermore, this exercise can serve as an example for improving computationally costly routines while improving their mathematical accuracy.

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.

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.

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.

2018, Beyer, Sebastian, Kleiner, Thomas, Aizinger, Vadym, Rückamp, Martin, Humbert, Angelika

Subglacial hydrology plays an important role in ice sheet dynamics as it determines the sliding velocity. It also drives freshwater into the ocean, leading to undercutting of calving fronts by plumes. Modeling subglacial water has been a challenge for decades. Only recently have new approaches been developed such as representing subglacial channels and thin water sheets by separate layers of variable hydraulic conductivity. We extend this concept by modeling a confined-unconfined aquifer system (CUAS) in a single layer of an equivalent porous medium (EPM). The advantage of this formulation is that it prevents unphysical values of pressure at reasonable computational cost. We performed sensitivity tests to investigate the effect of different model parameters. The strongest influence of model parameters was detected in terms of governing the opening and closure of the system. Furthermore, we applied the model to the Northeast Greenland Ice Stream, where an efficient system independent of seasonal input was identified about 500 km downstream from the ice divide. Using the effective pressure from the hydrology model, the Ice Sheet System Model (ISSM) showed considerable improvements in modeled velocities in the coastal region.

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.

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 %.

2018, Calov, Reinhard, Beyer, Sebastian, Greve, Ralf, Beckmann, Johanna, Willeit, Matteo, Kleiner, Thomas, Rückamp, Martin, Humbert, Angelika, Ganopolski, Andrey

We introduce the coupled model of the Greenland glacial system IGLOO 1.0, including the polythermal ice sheet model SICOPOLIS (version 3.3) with hybrid dynamics, the model of basal hydrology HYDRO and a parameterization of submarine melt for marine-terminated outlet glaciers. The aim of this glacial system model is to gain a better understanding of the processes important for the future contribution of the Greenland ice sheet to sea level rise under future climate change scenarios. The ice sheet is initialized via a relaxation towards observed surface elevation, imposing the palaeo-surface temperature over the last glacial cycle. As a present-day reference, we use the 1961-1990 standard climatology derived from simulations of the regional atmosphere model MAR with ERA reanalysis boundary conditions. For the palaeo-part of the spin-up, we add the temperature anomaly derived from the GRIP ice core to the years 1961-1990 average surface temperature field. For our projections, we apply surface temperature and surface mass balance anomalies derived from RCP 4.5 and RCP 8.5 scenarios created by MAR with boundary conditions from simulations with three CMIP5 models. The hybrid ice sheet model is fully coupled with the model of basal hydrology. With this model and the MAR scenarios, we perform simulations to estimate the contribution of the Greenland ice sheet to future sea level rise until the end of the 21st and 23rd centuries. Further on, the impact of elevation-surface mass balance feedback, introduced via the MAR data, on future sea level rise is inspected. In our projections, we found the Greenland ice sheet to contribute between 1.9 and 13.0 cm to global sea level rise until the year 2100 and between 3.5 and 76.4 cm until the year 2300, including our simulated additional sea level rise due to elevation-surface mass balance feedback. Translated into additional sea level rise, the strength of this feedback in the year 2100 varies from 0.4 to 1.7 cm, and in the year 2300 it ranges from 1.7 to 21.8 cm. Additionally, taking the Helheim and Store glaciers as examples, we investigate the role of ocean warming and surface runoff change for the melting of outlet glaciers. It shows that ocean temperature and subglacial discharge are about equally important for the melting of the examined outlet glaciers.