the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
IPSL-Perm-LandN: improving the IPSL Earth System Model to represent permafrost carbon-nitrogen interactions
Abstract. Permafrost soils have the potential to release large amounts of soil carbon to the atmosphere under climate change. However, in the Sixth Coupled Model Intercomparison Project (CMIP6), only two Earth System Models (ESM) represented permafrost carbon, both sharing the same land surface model. This makes future permafrost carbon dynamics highly uncertain and underscores the urgent need to include permafrost carbon in ESMs to enable more reliable future projections of climate change and remaining carbon budget estimates. Here, we present IPSL-Perm-LandN, an improved version of the Institut Pierre-Simon Laplace (IPSL) ESM (used for CMIP6) aiming at better representing high-latitude land ecosystems. The main developments are the inclusion of an explicit nitrogen cycle and of key permafrost physical and biogeochemical processes. The latent heat associated with soil water freeze/thaw is taken into account in the energy budget, as well as soil thermal insulation by soil organic matter and a surface organic layer (e.g. litter or moss). Soil organic carbon and nitrogen are vertically resolved with a depth-dependent decomposition dynamics, a key feature for representing the effect of gradual permafrost thaw on soil biogeochemistry. Cryoturbation is represented as a diffusion process that buries organic matter in the deeper soil layers. Compared to the previous version of the model used for CMIP6, we show that the extent of the permafrost region has improved significantly and that the simulated active layer thickness in the Arctic is in better agreement with observations. Permafrost soil carbon stocks have increased 20-fold to reach 1006 PgC in the top 3 meters of soil, which is consistent with observation-based estimates. We simulate that the permafrost region has been a net carbon sink over the past 150 years (+0.32±0.04 PgC.yr-1 on average between 2005 and 2014), primarily due to carbon uptake from boreal forests. This is comparable with recent pan-Arctic carbon balance estimates, when accounting for unrepresented processes in our model (fire and riverine carbon losses). Overall, the inclusion of permafrost processes has improved the response of the model to anthropogenic perturbations in high latitudes over the past century, marking a step forward in the representation of Arctic ecosystems.
- Preprint
(7155 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 18 Nov 2025)
- RC1: 'Comment on egusphere-2025-3656', Anonymous Referee #1, 10 Nov 2025 reply
-
RC2: 'Comment on egusphere-2025-3656', Anonymous Referee #2, 13 Nov 2025
reply
This manuscript describes the Earth System Model IPSL-Perm-LandN, a development based on IPSL’s CMIP6 model IPSL-CM6A-LR. Model improvements discussed here are based on a better physical representation of permafrost as well as an incorporation of vertically resolved soil carbon and an inclusion of soil nitrogen. A better representation of the terrestrial carbon cycle in general and permafrost carbon in particular in ESMs is becoming increasingly important especially since the CMIP framework is moving towards emission driven simulations instead of concentration driven simulations, and the permafrost carbon representation has been identified as a major source of uncertainty in climate projections. This study provides very relevant and important results on improvements mad in the IPSL’s contribution to CMIP and addresses important gaps in the CMIP6 contribution.
The manuscript is very well written, and results are presented in a clear and understandable manner. It addresses aspects of overall model performance, including atmosphere and ocean variables, as the model is discussed as a whole, even though improvements were restricted to the land model. The introduction would profit from a slightly broader view at the importance of representing the terrestrial carbon cycle under emission driven simulations, and the discussion would profit from a broader look at limitations of the chosen modelling set up.
I would like to echo one of the questions of reviewer 1 here: Will this be the model set up for IPSL’s contribution to CMIP7? If so, this should be emphasized in the conclusion section! Quite a number of LSMs that were used within the CMIP6 ESMs have had capabilities to represent land carbon processes better than it was done in CMIP6, but they were not included in the model set ups used for the simulations. A lot of LSM groups have made considerable improvements to their models since CMIP6, but to what degree these improvements will be part of the CMIP7 efforts remains to be seen, so if this model is IPSL’s contribution to CMIP7, that would be a nice signal to the ESM community.
I have one general comment concerning the results and discussion section with regard to permafrost dynamics and permafrost carbon. You describe a number of choices in the moel description that probably have technical reasons, but don’t seem straight forward. Their implications should be discussed in the results and discussion section.
- I assume there are technical reasons why there is only a 2m soil column for soil water, and that you wanted to do carbon and nitrogen pools for greater depth, but what are the consequences on latent heat and heat flux in general that soil moisture is only computed above 2m depths? ALT can be considerably deeper than 2m.
- Why is there not representation of moss and lichens as PFTs? And what processes are missing because you only consider their physical properties, but no biology?
- The soil column has 18 layers extending to 90m depth, and soil moisture is only done for 11 layers in the upper 2m. You assume a soil moisture below that, how does that effect decomposition, and are these deep layers even relevant? Soil carbon is often assumed to be limited to around the upper 3m of the soil column, what do your carbon and nitrogen pools in the deepest layer look like? And I assume that the soil moisture assumption for layers below 2m is only used for calculating decomposition, so mass is conserved, but does that make any sense in 90m depth? What are the consequences of these assumptions for your top layer carbon and nitrogen pools?
I have some more minor specific comments and questions to the authors listed below.
Specific Comments:
Lines 29-32: The permafrost-carbon feedback overall remains a major source of uncertainty, which is especially important in the light of the move towards emission driven simulations, where the actual CO2 concentrations in the atmosphere will to a large degree be determined by ESM's own carbon cycles, which should be included here. See eg Steinert, N. J. and Sanderson, B. M.: Normalizing the permafrost carbon feedback contribution to the Transient Climate Response to Cumulative Carbon Emissions and the Zero Emissions Commitment, Earth Syst. Dynam., 16, 1711–1721, https://doi.org/10.5194/esd-16-1711-2025, 2025 and Sanderson, B. M., Booth, B. B. B., Dunne, J., Eyring, V., Fisher, R. A., Friedlingstein, P., Gidden, M. J., Hajima, T., Jones, C. D., Jones, C. G., King, A., Koven, C. D., Lawrence, D. M., Lowe, J., Mengis, N., Peters, G. P., Rogelj, J., Smith, C., Snyder, A. C., Simpson, I. R., Swann, A. L. S., Tebaldi, C., Ilyina, T., Schleussner, C.-F., Séférian, R., Samset, B. H., van Vuuren, D., and Zaehle, S.: The need for carbon-emissions-driven climate projections in CMIP7, Geosci. Model Dev., 17, 8141–8172, https://doi.org/10.5194/gmd-17-8141-2024, 2024.
Line 44-46: There is an overview paper on the permafrost representation of the CMIP6 ESMs on a more general level than in Arora et al that should be mentioned here, Matthes, H., Damseaux, A., Westermann, S., Beer, C., Boone, A., Burke, E., ... & Wieder, W. R. (2025). Advances in Permafrost Representation: Biophysical Processes in Earth System Models and the Role of Offline Models. Permafrost and Periglacial Processes, 36(2), 302-318.
Line 56-59: The lack of resolving vertical soil carbon also prevents representation of abrupt permafrost thawing through eg fire or thaw slumps, which is estimated to be a major source of permafrost carbon loss (Turetsky, M.R., Abbott, B.W., Jones, M.C. et al. Carbon release through abrupt permafrost thaw. Nat. Geosci. 13, 138–143 (2020). https://doi.org/10.1038/s41561-019-0526-0), even though most LSMs in ESMs miss representation of the underlying processes as well. It would still be good to mention this here, since including vertically resolved soil carbon can be seen as creating a basis for inclusion of at least some of these fast processes.
Line 179: Here, and everywhere else you show units, you separate them by a period. While that increases readability, I don’t think it’s standard notation. It should probably be changed.
Line 305ff: Maybe I misunderstand what you are describing here. There is one mineral nitrogen pool, which is associated with the surface layer. If that pool is not sufficient, nitrogen is taken from the atmosphere. Nitrogen pools however are vertically resolved. So is mineralized nitrogen available to all soil layers out of the pool in the surface layer? And can nitrogen pools in all layers draw from the atmosphere through that mechanism as well? And if so, what consequences does that have?
Line 355: Does that definition of the permafrost area used here mean that the area for which cryoturbation is calculated is dynamic?
Line 360ff: If the SoilGrids data has vertical layers until 2m depth, how do you compute values for your layers below 2m? There is no information to interpolate, you could only extrapolate, and how does that make sense when there is hardly any carbon stored below 3m depth in most of the permafrost area?
Line 369: This is more of a comment: The equilibrium approach is the one everybody uses, but even in 1850, the carbon pools were not in equilibrium.
Line 500: While the realistic warming rates are good, when looking at permafrost dynamics and thawing related carbon release, what matters at the end are absolute temperatures, since freeze/thaw processes are coupled to an absolute temperature. From your latitudinal plot, it does seem that the underestimation is mainly originating from the tropics, but I still wonder what implications are, and if there is a seasonality to the bias that could impact permafrost?
Line 520: Does the overestimation of snow fall actually lead to an overestimation of snow cover? And what about the seasonality? Could that impact your permafrost dynamics?
Line 577: The overestimation of the southern boundary of permafrost might be associated with the underestimated surface air temperatures, depending on the seasonality of that bias.
Line 578ff: Is the underestimation of permafrost at the southern edge of the Canadian permafrost associated with maybe a spatial pattern in air temperature?
Line 603: I acknowledge that attribution of effects like this are difficult, but do you have any ideas why improving the terrestrial carbon cycle representation lead to a decrease in performance when it comes to GPP?
Citation: https://doi.org/10.5194/egusphere-2025-3656-RC2
Data sets
IPSL-Perm-LandN: improving the IPSL Earth System Model to represent permafrost carbon-nitrogen interactions Rémi Gaillard, Patricia Cadule, Philippe Peylin, Nicolas Vuichard, and Bertrand Guenet https://doi.org/10.5281/zenodo.16739216
Model code and software
IPSL-Perm-LandN: improving the IPSL Earth System Model to represent permafrost carbon-nitrogen interactions Rémi Gaillard, Patricia Cadule, Philippe Peylin, Nicolas Vuichard, and Bertrand Guenet https://doi.org/10.5281/zenodo.16739216
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 978 | 62 | 10 | 1,050 | 13 | 10 |
- HTML: 978
- PDF: 62
- XML: 10
- Total: 1,050
- BibTeX: 13
- EndNote: 10
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
General comments
This manuscript presents IPSL-Perm-LandN, an improved version of IPSL-CM6A-LR (used in CMIP6), incorporating a representation of key permafrost physical and biogeochemical processes and an explicit nitrogen cycle, for a better representation of high-latitude land ecosystems. The authors present the methodological details on the implementation of the model modifications for the relevant processes. In addition to relevant metrics and evaluations regarding permafrost dynamics and land carbon-cycle dynamics, the manuscript also provides an assessment and evaluation of atmospheric physics (here limited to temperature and precipitation) and ocean physics.
The manuscript provides relevant and novel results with regard to permafrost physics and carbon dynamics in historical simulations based on model modifications that address an important research gap. As the majority of CMIP climate models lack an adequate representation of permafrost physics, an assessment of the permafrost carbon-climate feedback was limited in CMIP6 and therefore AR6. The next iteration of climate models is in need of the inclusion of the relevant processes. This manuscript provides a description of including these processes in one of the participating CMIP models, which is not only a step towards a better representation of permafrost in CMIP, but could also serve as a blueprint for other climate models to consider the relevant processes.
I do have a few questions and comments, which I consider overall minor, that I would like the authors to address.
The manuscript is well written, structured and generally reads fluently. While the manuscript is quite long, the structuring and flow of the text help its understanding. However, several of the sentences in the Results Section (often the first sentence of the paragraph) do not explicitly state that they refer to the IPSL-Perm-LandN version (which I assume from the context and interpreting the figures), and therefore leave it up to interpretation for the reader. I suggest to carefully go through the Results Section to add these clarifications. I have indicated a few in my ‘Specific Comments’, but there may be more occurrences.
I commend the effort to assemble and present all model results in comparison to observations, and also evaluate the atmosphere and ocean to contextualize the results on permafrost modeling. Beyond that, the manuscript provides a valuable and concisely presented reference for other studies to compare their results to, with respect to both modeling results and observational products.
The methodological modifications presented here and the assumptions taken generally look good to me. Minor comments towards the methods are also specified in 'Specific Comments’.
I appreciate that the simulations presented here are run over the historical period, also because this lends itself to the evaluation with respect to observational products. Nonetheless, the efforts to improve the wider suite of CMIP models with regard to permafrost representation are often made with the motivation to address future changes of the climate. A comment in that direction was made in the, somewhat short, Conclusions Section, but I would appreciate some more discussion on how the improvements made in IPSL-Perm-LandN could/will improve the future assessment of permafrost changes and the associated permafrost carbon-climate feedback. Further, have the authors done any future experiments of note that could be presented here? Or will those be addressed in a separate follow-up paper?
Will the modifications of IPSL-Perm-LandN be incorporated in the next IPSL model version (i.e., for CMIP7)? And how does this operationally affect the spin-up/initialization strategy presented here, given that soil carbon initialization was done offline here?
Going into the manuscript, I expected an evaluation of the effect of the model changes on the permafrost carbon-climate feedback (even in the historical period). However, this assessment is not (directly) possible in the CO2-concentration driven simulations conducted here. What was the reasoning behind choosing to run CO2-concentration driven simulations? And can the results presented in Section 4.7 help to infer any of these results? Equally, with regards to this, the reader may appreciate some more discussion in the Conclusions section.
Specific comments
Technical corrections
References: