Comment on gmd-2020-440

As I think it is urgently necessary to consider more realistic stand conditions as well as management impacts into estimates of carbon and wood storages, I welcome the incorporation of such features into LPG-GUESS. I think that the results obtained for managed forest area are indeed advancing modelling science. Despite being not the first modelling approach addressing biomass development in dependence on forest management, I guess it is worthwhile publishing the respective components specifically for LPJ, in particular since the now available options are quite comprehensive and the question addressed is still not well investigated. In my opinion, the presentation is well structured and the language is fluent (some spelling errors may occur that are not specifically pinpointed but this can easily be corrected). Some specific recommendations related to references, figures and issues that might be better addressed are given below.

regarding standing volume and harvested volume in monoculture beech and spruce (space allowing, also pine and oak) stands located in a few European countries. Input to the model will be stand age, tree species and harvest intensity. Harvesting alternatives will be: no harvest, detailed harvest intensity every 5-10 years and automated harvest as used in the European simulations.
With this respect, I am particular puzzled about the example simulation, where beech seems to perform the superior growth relative to spruce in the PNV (for the first 100 years) as well as (and particularly in) all managed forest options. This seems to be counter-intuitive to what I would expect for a site in southern Sweden (see e.g. Bolte et al. 2010), although I am aware of only one extensive biomass study which gives almost the same biomass/production values for both species in this area (Nihlgard 1972). However, although the particular conditions seem to be more representative to Poland than to Sweden (see e.g. Jagodizinki et al. 2020), I might misjudge the growth and will happily accept if I am proven wrong. In any case, I think that particularly (the) illustrating example(s) need some data or references that support their realism.
REPLY: While it is our aim to improve the realism of interspecies competition in the model, we don't expect a DGVM to to a perfect job for every species in every region. Current developments of the model, out of scope for this article, that should be able to achieve a more realistic interspecies competition include an improved soil water-holding capacity description and growth limitations of nutrients other than nitrogen. Another important factor that influences competition at the tree seedling stage is herbivory, which is not yet currently implemented in the model. Figures 2 and 4, it is true that beech out-competes spruce when growing in PNV at the beginning of a vegetation succession (Fig.2b) and in mixed beech-spruce stands (Fig.3a). However, when grown in monocultures (Fig.4b), beach and spruce grows equally well until 1985 (when measured in kgC/m 2 ), when beech growth increases, and when measured in wood volume, spruce is ahead of beech all the time (ending up at 377 m 3 /ha vs. 302 m 3 /ha in 2000). Figure 4 illustrates the two different ways to achieve a pre-defined species mix in the model, but we could add a plot showing the standing volume development of spruce and beech although it might detract from the main message. The fact that the plot inside the red box Fig.4b is the biomass C development at the landscape level and that it's not possible to compare the growth of the two species in that plot can be better stressed in the figure caption.

At our south Swedish site in
We will add a short section in the discussion about present problems in interspecies competition between the European species in the model.
Second, it is unclear why old forests (one third of the area) are considered as PNV despite it is acknowledged by the authors that this leads to a simulated species composition widely divergent from the actual forests. It also leads to total forest estimates of carbon in the different compartments that are heavily influenced from this (unrealistic?) assumption. I guess that considering even a species composition from coarse data as is available from literature would largely decrease uncertainties related to the inventory as it is now.
REPLY: The representation of >140y forests as PNV is meant as a conservative compromise and should result in an underestimation of the C sink of European forests. Alternative setups using the GFAD age database would not necessarily improve the situation satisfactorily. Spinup using PNV, attaining a state of equilibrium of vegetation and soil carbon/nitrogen, is the standard procedure in LPJ-GUESS, has been evaluated extensively and has an uncertainty due to lack of land-use history which is to some extent predictable, while a spinup using the four GFAD forest classes would introduce further uncertainties. Thus, since we don't have any information about when to initialise the oldgrowth forest classes in GFAD and we believe that unsupported assumptions about age structure is more problematic than uncertainties regarding species coverage, we find it far from ideal but less problematic at this stage to leave the old-growth forest as PNV with a disturbance interval of 400 years. This also allows us to compare more directly with the global simulations of Pugh et al. (2020). Also, the >140y forests in GFAD contain disturbing artefacts, most likely as a result of including shrubland as forest. There are large (up to 20 % in eastern Europe) areas supposedly containing broadleaf evergreen forests and far too much needleleaf deciduous forests. The regrowth part of GFAD is not affected by these peculiarities, so in the European simulations, we like to focus on the regrowth part and treat >140y forest as a source for improvement for future setup refinements. We will further stress this in the text.
Also, the errors in NAI, GS and harvested volume resulting from species inconsistencies in >140y forests are expected to be modest since these are calculated from conversion factors for each European country and not from species-specific wood volumes.

Specific comments L17/18: How is the underestimation of carbon sink explained if stock and net growth is fairly well met?
REPLY: Two obvious reasons for this are a too-high soil carbon content because of a lack of land-use history in the setup that is leaking into the atmosphere after the creation of secondary stands, as well as an unrealistic lack of any removal of biomass in the >140y forests (see reply about this and suggestions for improvement to referee #2). We can add a short sentence about this in the abstract.

L29: Morin et al. 2017 not in ref. Do you mean Morin et al. 2018?
REPLY: Yes, thank you for spotting this error.

L41: Are the Pan et al. data really providing below-and aboveground carbon stocks for
Europe? The numbers given here don't seem to match the same values indicated for the same source in Table 2. Estimates from other sources seem to be also far different as e.g. app. 5/5 (Liski et al.),9/14 (Goodale et al.), and 9/8 PgC (Pilli et al.). Please check. Consider to introduce the given references into Table 2-4 if appropriate.
REPLY: Table S3 in Pan et al. 2011 provides carbon in total living biomass, dead wood, litter and soil. We included dead wood in the soil/litter compartment.We will add this detail in the table caption. We didn't manage to find discrepancies between L41 and table 3: The suggested references seem to include less than complete coverages of the study area (e.g. western Europe, excluding baltic countries, Belarus and Ukraine). Pan et al. includes data for soil and litter for all european countries, even though they have to rely on estimates from similar countries in 14 cases. The vakues also agree well with the other source of comparison in the manuscript, Forest Europe, taken into the account that soil data is missing there for a number of countries. We believe the message comes through even if the references used claim higher soil pools than others, i.e that the simulation setup probably overestimates soil carbon. REPLY (Revised text): In this study, we describe the implementation of expanded forest management capabilities in LPJ-GUESS v.4.0, which considers, in addition to detailed carbon-and water-cycle processes, nitrogen-cycling and nitrogen-limitation (Smith et al., 2014). With this, forest management in LPJ-GUESS is for the first time fully integrated in a model version capable of simulating a landscape containing a mosaic of land cover types like PNV, cropland, pasture and peatland and with a sophisticated land-use and land-cover change functionality.

L76ff: What is (are) the timestep(s) of the model?
REPLY (Entered into the text): Photosynthesis, respiration, phenology, soil carbon and nitrogen cycling and hydrology occur at a daily timestep, while biomass growth allocation and turnover, establishment and mortality occur at a yearly timestep.
L132ff: For the within-patch mixture of Fig.4a, the 60/40 ratio refers to while for the among-stand-types mixture of Fig. 4b the 60/40 ratio refers to groundcover area. I think that this should be indicated when presenting the different options.

L192/Tab. 1: consider changing 'N fertilization' to 'N fertilization/deposition' since the continuous kind of application looks more similar to a deposition regime (and is more likely to happen anyway).
REPLY: Nitrogen fertilisation, even though not currently available at defined points in time in the model, can be applied per management type on top of N deposition, which is applied as standard in all simulations, as described in the Methods. Although N leaching will be unrealistic, extra N will increase growth in the model considerably e.g. in Scandinavia.

L205/Fig. 5: This is a nice presentation indicating the potential options that can be selected. However, I wonder what happens, if 'old' would be selected together with 'small' or vice versa in the thinning rule options? What would happen? Wouldn't be a preference 'from above' or 'from below' be more suitable? What effect does this have on the simulation if biomass is taken e.g. preferentially from larger trees? I guess that the effect would depend on underlying assumptions about the biomass in different size groupswhat are these assumptions?
REPLY: Size preference overrides age preference, so there are no ambiguities possible. Age or size selection are most often interchangable in the model, but not always. But the terms 'from above' and 'from below' could be helpful to mention in the text or in the table. We will try to incorporate these terms. Since the cohorts suffer bioth intra-and interspecies competition for light in the layered canopy, the choice of cutting from below or above will be important for the success of the different cohorts and especially for the amount of regeneration of seedlings over time. The biomass in different size groups is a dynamic outcome of competition beween cohorts and species.
L215ff: I am a bit unhappy with the unit 'densmax' being trees per ha since it implies that the stand is harvested after a certain number of trees has reached. However, the number is actually not defined and can be very different dependent on the diameter distribution.

Densmax thus is an arbitrary density value in 1/m2 or similar unit. Similarly, I guess that 'dens' is also not indicating real tree numbers.
REPLY: The terms in the equations are exactly as in Bellassen et al. We think we should use the same terms to avoid confusion. The different cohorts of a pft are represented by a mean individual with a density in trees per m 2 . Tab REPLY: Agreed, this sounds akward. We should be able to use the wording 'harvest fraction' instead (since we already mentioned vegetation carbon).

L311: What does 'reductions in wood products and residuals' mean? Do you mean 'due to' instead of 'in'? Consider that the argument is repeated a couple of lines later.
REPLY: A better wording is probably 'not taking into account the fate of wood products and residues following removal from the site'. (Revised text) L330-332: irritating punctuation (needs ":" after forest, and "," after pubescens) REPLY: Revised text as suggested.
L486: The simulations are initialized and driven by the given information but not 'constrained', correct? The word would be correct, if repeated inventory data would be used to parameterize or adjust the simulations which I think has not been done.
REPLY: If using groundcover area and land-cover change functionality, we can constrain the simulations-In other cases, we can only constrain relative species mixes and harvest fractions, so yes, the wording is not accurate. (Revised text): 'The initialisation and harvest alternatives used are informed by available forest inventory data and harvest information.'

L494: As I take it, there are no observations but only estimates of 'mean growing stock' that are based on different information, sometimes on rather inhomogeneous sources, correct?
REPLY: That is probably true. We can change 'observations' to 'reported values' in most places in the text.
L526: You mean that the NPP decline that ORCHIDEE simulates is not simulated by LPJ, which might have various reasons that are caused by the model structure and processes. You don't mean that the decline is a function as such that is 'included' or not 'included', correct? Could such a decline be related to nutrient depletion? How is the effect of nutrient export -that is a main concern of a sustainable management -considered anyway? I think it should at least be part of the discussion (see e.g. Parolari and Porporato; Sverdrup et al.) REPLY: In some earlier versions of our forest management code, a reduction factor was applied to the light extiction coefficient after thinning, following Näslund et al. (1971, Stud.For.Suec. 89:124) and Hale (2003, For.Ecol.Manage. 179:341-349). This functionality was not included in the current version. But the sentence omits a lot of information and we are uncertain about by which mechanism this decline is achieved in ORCHIDEE, so we will change the wording to 'is not simulated in this version of LPJ-GUESS'.
Nitrogen depletion of the soil in previous land-use history has an effect on forest growth in the model as seen in Figure 2c. As discussed in the reply to referee #2's comment on line 107, a change in the amount of removed leaf nitrogen causes little changes in productivity and carbon pools at the European scale using the current version of LPJ-GUESS, but in further studies, the effect of different management practices on productivity at both stand and regional scale will be fully investigated. Obviously, the effect of nutrient export in different management schemes is potentially a very important use of forest management in LPJ-GUESS. Upcoming updates to the model such as improved soil water holding capacity description, growth limitation and cycling of nutrient other than nitrogen, will hopefully improve the usefulness of the model in this respect.
We will add a part of this to the discussion after consulting the suggested references. Referee #2: Including a more realistic forest management in DGVMs is important, and this manuscript clearly addresses and demonstrates this. Although everything is presented and explained, it was difficult to follow all approaches and concepts in the first read, perhaps because I am not familiar with LPJ-GUESS. Including (many) different options to achieve the same result, and extensively listing alternative options or exceptions are confusing and distract from the main message. An example of the first case is the presentation of 3 alternative ways of creating an (initial) age class distribution (Section 2.2.2., Figure 2), which requires already quite a good understanding of the concepts on patch/stand/management type in LPJ-GUESS. An example of the second case is section 2.3.2.1 (Species selection). It is full of words like "may", "can", "or", "optionally" and "possibly". I think it would be worthwhile to introduce early on in the paper predefined categories, and explain how these are implemented in LPJ-GUESS. For Europe, I would say three types of management are important to consider: age-class (even-aged) based forestry, uneven-aged forestry/continuous cover forestry, and pristine forests. Pristine forests would be implemented as Natural (PNV) in Figure 1,  REPLY: The suggestions of the referee concerning the presentation of the different management methods are all useful. We will add more background why the different options are needed, structure the text within the frame of the three major types of management and revise Table 1 according to the suggestions.