the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Global and Northern-High-Latitude Terrestrial carbon sinks in the 21st century from CMIP6 experiments
Abstract. Climate warming is accelerating the changes in the global terrestrial ecosystems and particularly those in the northern high latitudes (NHL), and rendering the land-atmosphere carbon exchange highly uncertain. The Coupled Model Intercomparison Project Phase 6 (CMIP6) employs the most updated climate models to estimate terrestrial ecosystem carbon dynamics driven by a new set of socioeconomic and climate change pathways. By analyzing the future (2015–2100) carbon fluxes estimated by ten CMIP6 models, we quantitatively evaluated the projected magnitudes, trends and uncertainties of global and NHL carbon fluxes under four scenarios plus the role of NHL in the global terrestrial ecosystem carbon dynamics. Overall, the models suggest that the global and NHL terrestrial ecosystems will be consistent carbon sinks in the future, and the extent of the carbon sinks is projected to be larger under scenarios with higher radiative forcing. By the end of this century, the models by average estimate the NHL net ecosystem productivity (NEP) as 0.54±0.77, 1.01±0 .98, 0.97±1.62, and 1.05±1.83 PgC/yr under SSP126, SSP245, SSP370 and SSP585, respectively. The uncertainties are not substantially reduced compared with earlier results, e.g., the Coupled Climate Carbon Cycle Model Intercomparison Project (C4MIP). Although NHL contributes a small fraction of the global carbon sink (~13 %), the relative uncertainties of NHL NEP are much larger than the global level. Our results provide insights into future carbon flux evolutions under future scenarios and highlight the urgent need to constrain the large uncertainties associated with model projections for making better climate mitigation strategies.
-
Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
-
Preprint
(1417 KB)
-
Supplement
(1534 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(1417 KB) - Metadata XML
-
Supplement
(1534 KB) - BibTeX
- EndNote
- Final revised paper
Journal article(s) based on this preprint
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2022-417', Xu Lian, 22 Jun 2022
The authors investigated the future trajectories, patterns and uncertainties of northern ecosystem carbon fluxes using an ensemble of CMIP6 climate models. They found that under future warmer climates, the CMIP6 models project an overall enhanced net carbon uptake by the ecosystems, resulted from a tradeoff between NPP and RH increases. The spread of estimates across individual models is larger than that of the global average. The paper is methodologically sound, well written, and the results are nicely presented. I have a number of line comments on the manuscript, as given below. I would like to recommend its publication with the following possible revisions or clarifications.
Line comments:
Lines 113-114: It’s better to use “ten models” directly, rather than “nine models with ten datasets”.
Line 123: I don’t understand why land surface temperature, rather than 2-m temperature, is used in this analysis. When people say ecosystem response to temperature, they often refer to air temperature.
Line 129: Do you mean the original annual outputs from models, or the annual values aggregated from original monthly outputs?
Line 130: If the model data are resampled to 1-degree global grids, this is supposed to be 360x180? Please clarify.
Lines 145-151: The model ensemble mean of global NEP is strongly higher than the estimate by the Global Carbon Project. What’s the implications for the NHL NEP and future projections of NEP changes?
Lines 167-172: It’s interesting to look at how the MHL mean NEP compares with the global mean, and whether this difference is contributed more by NPP or RH.
Figure 2 and 3: Please also include the multi-model ensemble mean result.
Citation: https://doi.org/10.5194/egusphere-2022-417-RC1 -
AC1: 'Reply on RC1', Han Qiu, 08 Sep 2022
Reviewer 1
The authors investigated the future trajectories, patterns and uncertainties of northern ecosystem carbon fluxes using an ensemble of CMIP6 climate models. They found that under future warmer climates, the CMIP6 models project an overall enhanced net carbon uptake by the ecosystems, resulted from a tradeoff between NPP and RH increases. The spread of estimates across individual models is larger than that of the global average. The paper is methodologically sound, well written, and the results are nicely presented. I have a number of line comments on the manuscript, as given below. I would like to recommend its publication with the following possible revisions or clarifications.
Line comments:
Lines 113-114: It’s better to use “ten models” directly, rather than “nine models with ten datasets”.
Response: we used ‘ten models’ instead in the revised manuscript.
Line 123: I don’t understand why land surface temperature, rather than 2-m temperature, is used in this analysis. When people say ecosystem response to temperature, they often refer to air temperature.
Response: Thanks for pointing out this. In the revision, we replaced the land surface temperature with 2-m air temperature in our analyses and reported the new results in the updated manuscript. We found no essential changes to the results and findings with the new variable. But it is indeed true that the use of 2-m air temperature will potentially make our report more useful.
Line 129: Do you mean the original annual outputs from models, or the annual values aggregated from original monthly outputs?
Response: We meant the annual values aggregated from original monthly outputs. We clarified this in the revised manuscript.
Line 130: If the model data are resampled to 1-degree global grids, this is supposed to be 360x180? Please clarify.
Response: The ‘1-degree grid’ was an approximation. In practice we resampled the model data based on the model grids of BCC-CSM2-MR whose resolution is in the middle (globally 320x160). We clarified it in the revised manuscript.
Lines 145-151: The model ensemble mean of global NEP is strongly higher than the estimate by the Global Carbon Project. What’s the implications for the NHL NEP and future projections of NEP changes?
Response: Thanks for pointing out that. We focused on NEP in our analyses while Global Carbon Project (GCP) reported Net Biosphere Productivity (NBP) which took disturbance-induced carbon fluxes from NEP and thus was lower. To make a meaningful comparison between CMIP6 and GCP data, we added historical CMIP6 NBP into analysis in the revision. According to the updated Fig.1a, the CMIP6 estimated lower global NBP than GCP data, and the CMIP6 NBP is closer to GCP data than NEP,
Lines 167-172: It’s interesting to look at how the NHL mean NEP compares with the global mean, and whether this difference is contributed more by NPP or RH.
Response: Thanks for the suggestion. We looked into the mean carbon fluxes and added the following sentences:
“Except SSP126, similarly positive but generally smaller trends were found for RH at global scales (Figure S4, Table 2) with the rates as 87.15, 173.39, 254.43 and 318.31 Tg C/year2 under the four scenarios. The NHL RH trends are 18.64, 36.27, 55.39 and 72.56 Tg C/year2, normalized by the area, the growth rates are 0.44, 1.33, 1.99 and 2.62 g C/year2 for global NPP over the four scenarios respectively. The area-normalized growth rates in the NHL NPP are 0.54, 1.37, 2.03 and 2.63 g C/year2, respectively. Area-normalized global RH growth rates are 0.59, 1.17, 1.72 and 2.15 g C/year2 while the area-normalized NHL RH growth rates are 0.62, 1.20, 1.84 and 2.41 g C/year2 under the four scenarios, respectively. These results indicate that a faster average growing NPP and RH in the NHL than the global average. The fast-growing RH cancelled a large part of the NPP growth and resulted in small growing NEPs.”
Figure 2 and 3: Please also include the multi-model ensemble mean result.
Response: we added multi-model ensemble mean results in the updated Figure 2 and 3.
-
AC1: 'Reply on RC1', Han Qiu, 08 Sep 2022
-
RC2: 'Comment on egusphere-2022-417', Anonymous Referee #2, 05 Jul 2022
General Comments
The manuscript by Qiu and colleagues analyzes future projections from models included in the Coupled Model Intercomparison Project Phase 6 (CMIP6), with a focus on carbon fluxes in the northern high latitudes (NHL) compared to rest of the globe. The authors find that the CMIP6 models project terrestrial ecosystems to remain carbon sinks, with net sequestration rates increasing with global temperature changes, and that NHL ecosystems are a relatively minor fraction of the net sink but considerably more uncertain than the global mean.
Overall, this is a straightforward and relatively simple paper documenting CMIP6 output in the NHL domain. It reads somewhat like a report and contains little novelty, using familiar tools and approaches with largely unsurprising results, especially for readers following the evolution of these models and how they represent the carbon cycle. Nonetheless, the data being reported is important to the scientific community. The manuscript would benefit from additional context and discussion, grammar checks, and potentially additional analyses. It’s also unclear to me if the authors are actually analyzing the net carbon sink (Net Ecosystem Productivity vs. Net Biome Productivity).
Specific comments
My main technical question is why the authors chose to analyze Net Ecosystem Productivity, which is the balance between photosynthesis and ecosystem respiration, versus Net Biome Productivity, which includes other fluxes such as fire and harvest disturbance and is a land model’s attempt at the complete carbon cycle (it would also be helpful to include lateral fluxes if possible). I realize the participating models represent different aspects of the carbon cycle, which may make comparisons more challenging, but the manuscript essentially equates NEP to net carbon balance, which is incorrect (although the difference may be minor). It may be worth adding Supplementary analyses to address this: does the inclusion of these other fluxes change the results shown, even for a subset of the models? Furthermore, is direct comparison of CMIP6 NEP to values from the Global Carbon Project valid given these differences? Does that partially account for the large differences shown between the two in Figure 1? The authors do not provide details on how they made this comparison, so it’s not possible for the reader to tell.
Personally, I also believe a manuscript such as this that discusses future carbon budgets would benefit from a historical comparison to data. This may be beyond the scope, but how should readers interpret the large differences between CMIP6 models and the GCP, and what does this ultimately mean for the future projections? Do models that tend to represent historical properties better (e.g., upscaled carbon fluxes from FLUXCOM, global biomass, LAI, etc.) tend to project higher or lower carbon sinks? These are important questions, some of which could be addressed in the Discussion and brining in past literature. Along similar lines, the Discussion could benefit from more context and interpretation. Why haven’t uncertainties been reduced throughout CMIP versions, how would including permafrost carbon and disturbances change the assessment, what are the major recommendations on ways forward considering this?
Why did the authors analyze land surface temperature as opposed to a property such as 2 m air temperature, which is much more commonly used as a metric and benchmark? Land surface temperature accounts for not only the climate changes but also land surface responses, and in that way seems to add unnecessary complexity.
Finally, the manuscript would benefit from a thorough grammar check throughout. I addressed a minor fraction of grammatical errors in my comments below, but many more remain.
Technical corrections
L 16: What’s the meaning of the word ‘extent’ here - spatial distribution, magnitude, or other?
L 20: do you mean to say CMIP5 here?
L20: NHL was defined and used previously as a plural noun, but here singular. Please remain consistent.
-Maybe define what domain NHL is referring to in the abstract?
L 32: The land carbon sink changes by a large amount interannually owing to annual climate oscillations, disturbances, etc.
L 45: this reference is almost a decade old now, and the literature it cites is over a decade; consider adding newer references for warming rates
L 76: What does the word ‘devoted’ mean here?
L 77: missing parenthesis
L 86: Can the authors expand on the ‘newly updated data’ they’re referring to?
L 114: Incorrect use of “i.e.”. Could say “…in this study, including…”
L 126: Did the authors account for non-land fractions of grid cells in their area-weighting?
L 138. The readers would benefit from more details on how these sensitivity analyses were conducted
Figure 1: 2095-2100 is a short period to use to calculate standard deviations. More typical would be something like a 20 year time period
L 189-190: Change ‘huge’ to something like ‘large’
Fig 1: Difficult to see the GCP values with all the light blue bars. Possible to change that into a shaded time series as well?
L 273: ‘Minimization’? Do the authors mean mineralization?
L 281: What does the word ‘special’ mean here?
L 287: The carbon balance will also significantly be impacted by disturbances, mentioned in the introduction but mostly not included in the CMIP6 models. This point should be emphasized.
L 291: Particularly poor grammar
L 302: “plant functional types”
L 305: What is meant by ‘compensation’ here?
L 332: Hyphen after Northern not necessary
Citation: https://doi.org/10.5194/egusphere-2022-417-RC2 -
AC2: 'Reply on RC2', Han Qiu, 08 Sep 2022
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2022/egusphere-2022-417/egusphere-2022-417-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Han Qiu, 08 Sep 2022
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2022-417', Xu Lian, 22 Jun 2022
The authors investigated the future trajectories, patterns and uncertainties of northern ecosystem carbon fluxes using an ensemble of CMIP6 climate models. They found that under future warmer climates, the CMIP6 models project an overall enhanced net carbon uptake by the ecosystems, resulted from a tradeoff between NPP and RH increases. The spread of estimates across individual models is larger than that of the global average. The paper is methodologically sound, well written, and the results are nicely presented. I have a number of line comments on the manuscript, as given below. I would like to recommend its publication with the following possible revisions or clarifications.
Line comments:
Lines 113-114: It’s better to use “ten models” directly, rather than “nine models with ten datasets”.
Line 123: I don’t understand why land surface temperature, rather than 2-m temperature, is used in this analysis. When people say ecosystem response to temperature, they often refer to air temperature.
Line 129: Do you mean the original annual outputs from models, or the annual values aggregated from original monthly outputs?
Line 130: If the model data are resampled to 1-degree global grids, this is supposed to be 360x180? Please clarify.
Lines 145-151: The model ensemble mean of global NEP is strongly higher than the estimate by the Global Carbon Project. What’s the implications for the NHL NEP and future projections of NEP changes?
Lines 167-172: It’s interesting to look at how the MHL mean NEP compares with the global mean, and whether this difference is contributed more by NPP or RH.
Figure 2 and 3: Please also include the multi-model ensemble mean result.
Citation: https://doi.org/10.5194/egusphere-2022-417-RC1 -
AC1: 'Reply on RC1', Han Qiu, 08 Sep 2022
Reviewer 1
The authors investigated the future trajectories, patterns and uncertainties of northern ecosystem carbon fluxes using an ensemble of CMIP6 climate models. They found that under future warmer climates, the CMIP6 models project an overall enhanced net carbon uptake by the ecosystems, resulted from a tradeoff between NPP and RH increases. The spread of estimates across individual models is larger than that of the global average. The paper is methodologically sound, well written, and the results are nicely presented. I have a number of line comments on the manuscript, as given below. I would like to recommend its publication with the following possible revisions or clarifications.
Line comments:
Lines 113-114: It’s better to use “ten models” directly, rather than “nine models with ten datasets”.
Response: we used ‘ten models’ instead in the revised manuscript.
Line 123: I don’t understand why land surface temperature, rather than 2-m temperature, is used in this analysis. When people say ecosystem response to temperature, they often refer to air temperature.
Response: Thanks for pointing out this. In the revision, we replaced the land surface temperature with 2-m air temperature in our analyses and reported the new results in the updated manuscript. We found no essential changes to the results and findings with the new variable. But it is indeed true that the use of 2-m air temperature will potentially make our report more useful.
Line 129: Do you mean the original annual outputs from models, or the annual values aggregated from original monthly outputs?
Response: We meant the annual values aggregated from original monthly outputs. We clarified this in the revised manuscript.
Line 130: If the model data are resampled to 1-degree global grids, this is supposed to be 360x180? Please clarify.
Response: The ‘1-degree grid’ was an approximation. In practice we resampled the model data based on the model grids of BCC-CSM2-MR whose resolution is in the middle (globally 320x160). We clarified it in the revised manuscript.
Lines 145-151: The model ensemble mean of global NEP is strongly higher than the estimate by the Global Carbon Project. What’s the implications for the NHL NEP and future projections of NEP changes?
Response: Thanks for pointing out that. We focused on NEP in our analyses while Global Carbon Project (GCP) reported Net Biosphere Productivity (NBP) which took disturbance-induced carbon fluxes from NEP and thus was lower. To make a meaningful comparison between CMIP6 and GCP data, we added historical CMIP6 NBP into analysis in the revision. According to the updated Fig.1a, the CMIP6 estimated lower global NBP than GCP data, and the CMIP6 NBP is closer to GCP data than NEP,
Lines 167-172: It’s interesting to look at how the NHL mean NEP compares with the global mean, and whether this difference is contributed more by NPP or RH.
Response: Thanks for the suggestion. We looked into the mean carbon fluxes and added the following sentences:
“Except SSP126, similarly positive but generally smaller trends were found for RH at global scales (Figure S4, Table 2) with the rates as 87.15, 173.39, 254.43 and 318.31 Tg C/year2 under the four scenarios. The NHL RH trends are 18.64, 36.27, 55.39 and 72.56 Tg C/year2, normalized by the area, the growth rates are 0.44, 1.33, 1.99 and 2.62 g C/year2 for global NPP over the four scenarios respectively. The area-normalized growth rates in the NHL NPP are 0.54, 1.37, 2.03 and 2.63 g C/year2, respectively. Area-normalized global RH growth rates are 0.59, 1.17, 1.72 and 2.15 g C/year2 while the area-normalized NHL RH growth rates are 0.62, 1.20, 1.84 and 2.41 g C/year2 under the four scenarios, respectively. These results indicate that a faster average growing NPP and RH in the NHL than the global average. The fast-growing RH cancelled a large part of the NPP growth and resulted in small growing NEPs.”
Figure 2 and 3: Please also include the multi-model ensemble mean result.
Response: we added multi-model ensemble mean results in the updated Figure 2 and 3.
-
AC1: 'Reply on RC1', Han Qiu, 08 Sep 2022
-
RC2: 'Comment on egusphere-2022-417', Anonymous Referee #2, 05 Jul 2022
General Comments
The manuscript by Qiu and colleagues analyzes future projections from models included in the Coupled Model Intercomparison Project Phase 6 (CMIP6), with a focus on carbon fluxes in the northern high latitudes (NHL) compared to rest of the globe. The authors find that the CMIP6 models project terrestrial ecosystems to remain carbon sinks, with net sequestration rates increasing with global temperature changes, and that NHL ecosystems are a relatively minor fraction of the net sink but considerably more uncertain than the global mean.
Overall, this is a straightforward and relatively simple paper documenting CMIP6 output in the NHL domain. It reads somewhat like a report and contains little novelty, using familiar tools and approaches with largely unsurprising results, especially for readers following the evolution of these models and how they represent the carbon cycle. Nonetheless, the data being reported is important to the scientific community. The manuscript would benefit from additional context and discussion, grammar checks, and potentially additional analyses. It’s also unclear to me if the authors are actually analyzing the net carbon sink (Net Ecosystem Productivity vs. Net Biome Productivity).
Specific comments
My main technical question is why the authors chose to analyze Net Ecosystem Productivity, which is the balance between photosynthesis and ecosystem respiration, versus Net Biome Productivity, which includes other fluxes such as fire and harvest disturbance and is a land model’s attempt at the complete carbon cycle (it would also be helpful to include lateral fluxes if possible). I realize the participating models represent different aspects of the carbon cycle, which may make comparisons more challenging, but the manuscript essentially equates NEP to net carbon balance, which is incorrect (although the difference may be minor). It may be worth adding Supplementary analyses to address this: does the inclusion of these other fluxes change the results shown, even for a subset of the models? Furthermore, is direct comparison of CMIP6 NEP to values from the Global Carbon Project valid given these differences? Does that partially account for the large differences shown between the two in Figure 1? The authors do not provide details on how they made this comparison, so it’s not possible for the reader to tell.
Personally, I also believe a manuscript such as this that discusses future carbon budgets would benefit from a historical comparison to data. This may be beyond the scope, but how should readers interpret the large differences between CMIP6 models and the GCP, and what does this ultimately mean for the future projections? Do models that tend to represent historical properties better (e.g., upscaled carbon fluxes from FLUXCOM, global biomass, LAI, etc.) tend to project higher or lower carbon sinks? These are important questions, some of which could be addressed in the Discussion and brining in past literature. Along similar lines, the Discussion could benefit from more context and interpretation. Why haven’t uncertainties been reduced throughout CMIP versions, how would including permafrost carbon and disturbances change the assessment, what are the major recommendations on ways forward considering this?
Why did the authors analyze land surface temperature as opposed to a property such as 2 m air temperature, which is much more commonly used as a metric and benchmark? Land surface temperature accounts for not only the climate changes but also land surface responses, and in that way seems to add unnecessary complexity.
Finally, the manuscript would benefit from a thorough grammar check throughout. I addressed a minor fraction of grammatical errors in my comments below, but many more remain.
Technical corrections
L 16: What’s the meaning of the word ‘extent’ here - spatial distribution, magnitude, or other?
L 20: do you mean to say CMIP5 here?
L20: NHL was defined and used previously as a plural noun, but here singular. Please remain consistent.
-Maybe define what domain NHL is referring to in the abstract?
L 32: The land carbon sink changes by a large amount interannually owing to annual climate oscillations, disturbances, etc.
L 45: this reference is almost a decade old now, and the literature it cites is over a decade; consider adding newer references for warming rates
L 76: What does the word ‘devoted’ mean here?
L 77: missing parenthesis
L 86: Can the authors expand on the ‘newly updated data’ they’re referring to?
L 114: Incorrect use of “i.e.”. Could say “…in this study, including…”
L 126: Did the authors account for non-land fractions of grid cells in their area-weighting?
L 138. The readers would benefit from more details on how these sensitivity analyses were conducted
Figure 1: 2095-2100 is a short period to use to calculate standard deviations. More typical would be something like a 20 year time period
L 189-190: Change ‘huge’ to something like ‘large’
Fig 1: Difficult to see the GCP values with all the light blue bars. Possible to change that into a shaded time series as well?
L 273: ‘Minimization’? Do the authors mean mineralization?
L 281: What does the word ‘special’ mean here?
L 287: The carbon balance will also significantly be impacted by disturbances, mentioned in the introduction but mostly not included in the CMIP6 models. This point should be emphasized.
L 291: Particularly poor grammar
L 302: “plant functional types”
L 305: What is meant by ‘compensation’ here?
L 332: Hyphen after Northern not necessary
Citation: https://doi.org/10.5194/egusphere-2022-417-RC2 -
AC2: 'Reply on RC2', Han Qiu, 08 Sep 2022
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2022/egusphere-2022-417/egusphere-2022-417-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Han Qiu, 08 Sep 2022
Peer review completion
Journal article(s) based on this preprint
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
380 | 160 | 12 | 552 | 29 | 5 | 6 |
- HTML: 380
- PDF: 160
- XML: 12
- Total: 552
- Supplement: 29
- BibTeX: 5
- EndNote: 6
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
Dalei Hao
Yelu Zeng
Xuesong Zhang
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(1417 KB) - Metadata XML
-
Supplement
(1534 KB) - BibTeX
- EndNote
- Final revised paper