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.

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, Köhler, Peter, Knorr, Gregor, Stap, Lennert B., Ganopolski, Andrey, de Boer, Bas, van de Wal, Roderik S. W., Barker, Stephen, Rüpke, Lars H.

We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state-dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land ice sheet growth and sea level fall) the multimillennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction-based extrapolation of S for the future, we exclude intervals with strong ΔTg-CO2 divergence and find that S is less state-dependent, or even constant state-independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2.

2022, Feldmann, Johannes, Reese, Ronja, Winkelmann, Ricarda, Levermann, Anders

Basal ice-shelf melting is the key driver of Antarctica's increasing sea-level contribution. In diminishing the buttressing force of the ice shelves that fringe the ice sheet, the melting increases the ice discharge into the ocean. Here we contrast the influence of basal melting in two different ice-shelf regions on the time-dependent response of an isothermal, inherently buttressed ice-sheet-shelf system. In the idealized numerical simulations, the basal-melt perturbations are applied close to the grounding line in the ice-shelf's (1) ice-stream region, where the ice shelf is fed by the fastest ice masses that stream through the upstream bed trough and (2) shear margins, where the ice flow is slower. The results show that melting below one or both of the shear margins can cause a decadal to centennial increase in ice discharge that is more than twice as large compared to a similar perturbation in the ice-stream region. We attribute this to the fact that melt-induced ice-shelf thinning in the central grounding-line region is attenuated very effectively by the fast flow of the central ice stream. In contrast, the much slower ice dynamics in the lateral shear margins of the ice shelf facilitate sustained ice-shelf thinning and thereby foster buttressing reduction. Regardless of the melt location, a higher melt concentration toward the grounding line generally goes along with a stronger response. Our results highlight the vulnerability of outlet glaciers to basal melting in stagnant, buttressing-relevant ice-shelf regions, a mechanism that may gain importance under future global warming.

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