Search Results

Now showing 1 - 2 of 2
  • 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
    Accelerated photosynthesis routine in LPJmL4
    (Katlenburg-Lindau : Copernicus, 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.