General comment

The authors implemented eco-evolutionary optimality principles into the LPJ-GUESS vegetation model, by merging the model with an adapted version of the P-model for photosynthesis. They then evaluated the performance of the model against in-situ eddy covariance fluxtower measurements and global data products of GPP, as well as existing global biome distribution maps. While the metrics show only a small (but significant) improvement compared to the standard version of LPJ-GUESS, these new simulations are based on highly reasonable optimality constraints, which allow the model to use less (often poorly constrained) parameters and assumptions. I think this is a very interesting approach and a very useful and meaningful contribution to the development of LPJ-GUESS and dynamic vegetation models in general. I do have a few (minor) questions and concerns, especially regarding the different values for leaf-to-canopy scaling factor and comparison with a nitrogen-limited version of standard LPJ-GUESS. However, I could see the article being published close to its current form.

Specific comments

• L194-216: These lines are sligthly confusing at first read, because it seems Eqs. (3) and (4) should be connected like Eqs. (1) and (2) are. However, when puzzling everything together further in the section it becomes clear that the relationship between GPP_indiv and I_abs,indiv is more complex in LPJ-GUESS. Maybe consider rephrasing a bit to make this clear earlier on? Maybe starting from Eq. (4) could help, and then introduce Eq. (3) only later when it's needed. And maybe don't mention Eqs (1) and (2) at all, because they are not really used? Just a suggestion.
• L342: Why is this extremely rare? Aren't such drought conditions rather common, especially when you do global simulations? Is this also extremely rare in Standard LPJ-GUESS? To which degree does the P-model capture the diversity in stomatal regulation (e.g. isohydry)?
• L349-350: Regarding the SMART scheme producing somewhat better results, is this in general, or in your simulations? What do you have to back up this statement? E.g. if in general, do you have a reference for this?
• L359: Correct me if I'm wrong, but as far as I know the leaf-to-canopy scaling factor (α_a) can have quite a large impact on the simulated GPP and carbon stocks in LPJ-GUESS. I understand that you set this parameter a-priori (ie. you did not tune it afterward to get a good result) even though the new value is still "an initial guess that appears more plausible" (L668-669). Still, is it fair to compare two versions of LPJ-GUESS that have a different value for this? How much of the model difference is due to this parameter, compared to the difference introduced by the EEO approach alone? See also my next comment for a similar concern.
• L371-373: The standard simulations have N limitation activated, but your EEO-model does not. Similarly to the previous comment, I wonder how this affects the model comparison, e.g. how much of the difference in model output can be explained by this? Doesn't this difference obscure the effect of the EEO approach on its own? I understand that you want to compare the "default versions" of both models (ie. the standard model has N-limitation by default, and your model does not) and that is totally valid (same for alpha_a). However, it also seems a bit "unfair" to compare a model with N-limitation to a model without it? See also next comment. I also see you discuss this in Sect 4.5 where you give good reasons to doubt the causal relationship between soil N and V_cmax, and you provide future pragmatic approach to include N limitation. But this also makes me wonder whether it would make sense to (also) compare your model against a version of LPJ-GUESS where N-limitation on V_c,max is turned off?
• L444: Similar to the previous comment, the standard model version produces a higher annual total, even though it is N limited. I wonder what N limitatation would have on the P-model? Maybe it will underestimate GPP? See also Table 2 (and text from L534) which shows that C stocks, leaf area and NPP were overall a bit lower in EEO version, despite not having N limitation.
• L690: Is this really a fundamental part of the P-model on its own?

Technical corrections

• L124-125: Maybe rephrase? A bit hard to read in one breath
• L37,L135: The term "mitochondrial respiration" is only used twice in this paper, while it's called dark respiration otherwise. Consider harmonising this terminology throughout the paper.
• L145: Typo in "LPJ-GUES"
• L165: It could be nice to refer to the three EEO hypotheses in the later part of this section, so we know which part of the theory (which formula) corresponds to which hypothesis, ie. from line 196
• L218: What is ξ? It has not been defined before. I do see that you define ξopt in eq. (6) however. Or did you maybe mean the ratio χ on this line, as in L174-175?
• L324: Word and maybe a comma missing?
• L473-474: ...but for many of these biomes still better than the standard model, right? (from Table B2)
• L476: The reference to Fig. 6 comes before the first reference to Fig. 5 as far as I can tell. Maybe reorder the figures?
• L804: I guess the bold values in Table B2 mean the "better" value. Maybe add that to the caption.
• L551: "lower left and upper right quadrants"?

The authors applied an eco-evolutionary optimality (EEO) approach for calculation of assimilation rates in the LPJ-GUESS DGVM. Additionally, an acclimation process was implemented. This alternative approach has two major advantages: It saves computation time and gives better agreement between observed and modelled GPP. They used for LPJ-GUESS mode in individual mode and see an improvement in model performance and better agreement with observational data. In general, the paper is well written and structured and fits into the scop of the journal.

The model is run with modelling individual trees. it might be interesting to see the performance in run time and quality of the average individual model mode. This mode is significantly faster and is best suited for coupling to Earth system models. `The authors state, however, that acclimation must be calculated on an individual level. A short discussion how this approach can be applied to average individual models should be added.

The authors compare the LPJ-GUESS model with EEO photosynthesis and acclimation to the unchanged model. It would be interesting to see the impact of different acclimation times on the model performance to see how important the addition of acclimation is.

More specific remarks:

Line 292: Soil moisture is denoted as Theta, identical to the co-limitation parameter Theta in Table 1.

Line 279: One question. Does this correspond to the original photosynthesis model setting the co-limitation factor Theta=1? Are there any justifcati9on to assume this?

Please add unit of GPP and R_d in the figure captions.

The authors use 500 yrs of spin up of the model. The carbon and nitrogen pool should be in equilibrium after the spin up period It might be in particular for the boreal and polar regions that the pools are in disequilibrium because of the low temperatures leading to very low decomposition rates This might not be important for evaluation of GPP values bit for soil pools it is important. Also nitrogen limitation might change.

The author claims that an error in the fAPAR values propagate directly to the GPP value due to its linear relationship. The fAPAR dynamics is also affected by the phenology. LPJ-GUESS has an annual update of the carbon pools distributing the assimilated carbon to all vegetation carbon pools. The phenology function determines the intra-annual the intra-annual variability of leaf carbon changing the fAPAR value. A discussion should be added.

The authors claim to have a faster model using the EEO approach. Pleas gives some numbers for a typical configuration on an HPC cluster.