Atmospheric River Representation in the Energy Exascale Earth System Model (E3SM) Version 1.0
- 1Department of Geography, University of California, Berkeley, CA, USA
- 2Pacific Northwest National Laboratory, Richland, WA, USA
- 3Joint Institute for Regional Earth System Science and Engineering, University of California, Los Angeles, CA, USA
- 4Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA
- 1Department of Geography, University of California, Berkeley, CA, USA
- 2Pacific Northwest National Laboratory, Richland, WA, USA
- 3Joint Institute for Regional Earth System Science and Engineering, University of California, Los Angeles, CA, USA
- 4Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA
Abstract. The Energy Exascale Earth System Model (E3SM) Project is an ongoing, state-of-the-science Earth system modeling, simulation, and prediction project developed by the U.S. Department of Energy (DOE). With an emphasis on supporting DOE's energy mission, understanding and quantifying how well the model simulates water cycle processes is of particular importance. Here, we evaluate E3SM version v1.0 for its ability to represent atmospheric rivers (ARs), which play significant roles in water vapor transport and precipitation. The characteristics and precipitation associated with global ARs in E3SM at standard resolution (1° × 1°) are compared to the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA2). Global pattern of AR frequencies in E3SM show high degrees of correlation (>= 0.97) with MERRA2 and low mean absolute errors (< 1 %) annually, seasonally, and across different ensemble members. However, some large-scale condition biases exist leading to AR biases - most significant of which are: the double-ITCZ, a stronger and/or equatorward shifted subtropical jet during boreal and austral winter, and enhanced northern hemisphere westerlies during summer. By comparing atmosphere-only and fully coupled simulations, we attribute the sources of the biases to the atmospheric component or to a coupling response. Using relationships revealed in Dong et al. (2021), we provide evidence showing the stronger north Pacific jet in winter and enhanced northern hemisphere westerlies during summer associated with E3SM's double-ITCZ and related weaker AMOC, respectively, are the sources of much of the AR biases found in the coupled simulations.
Sol Kim et al.
Status: final response (author comments only)
-
CEC1: 'Comment on gmd-2021-364', Juan Antonio Añel, 29 Dec 2021
Dear authors,
After checking your manuscript, it has come to our attention that it does not comply with our Code and Data Policy.
https://www.geoscientific-model-development.net/policies/code_and_data_policy.htmlUnfortunately, we can not accept that tARget v3 is available only on personal request. This software is a critical part of your work, and therefore must be included with the manuscript. You must deposit it in one of the acceptable repositories following our policy, under a FLOSS (free-libre open source) license, and you must obtain a DOI for it. Then you must cite it in the manuscript.
Also, we can not accept that you restrict access to ARTMIP data. In the same way, you must add the relevant data for the work to the repository.
Please, address these issues as soon as possible to continue with the Discussions stage and the review process for your manuscript.Regards,
Juan A. Añel
Geosc. Mod. Dev. Executive Editor
-
AC1: 'Reply on CEC1', Sol Kim, 03 Jan 2022
Dear Juan A. Añel,
We thank you for bringing this to our attention. Our updated Code and data availability section should read:
Code and data availability. E3SM and MERRA2 datasets used in this study are available at https://esgf-node.llnl.gov/projects/e3sm/ and https://disc.gsfc.nasa.gov/datasets?project=MERRA-2 respectively. The tARget v3 algorithm is available at https://doi.org/10.25346/S6/B89KXF and outputs using tARget are available via ARTMIP (Shields et al. (2018b); https://doi.org/10.5065/D6R78D1M).
-
CEC2: 'Reply on AC1', Juan Antonio Añel, 04 Jan 2022
Dear authors,
Many thanks for making your code and data available in the Dataverse repository. However, I have found an additional problem: There is no license listed. If you do not include a license, the code is not FLOSS (free-libre open-source); it continues to be your property. Therefore, you should choose a FLOSS license and add it to your repository so that it applies to all the code. We recommend the GPLv3. You only need to include the file 'https://www.gnu.org/licenses/gpl-3.0.txt' as LICENSE.txt with your code. Also, you can choose other options that Zenodo provides: GPLv2, Apache License, MIT License, etc.
Also, we can not accept the UCAR/NCAR repository. I understand that the ARTMIP outputs are too big to be copied to another repository. However, you should store on a different one at least a minimum set of final fields/variables that you have used for your work, not complete model outputs.
Please, address these additional issues as soon as possible.Regards,
Juan A. Añel
Geosc. Mod. Dev. Executive Editor
-
AC2: 'Reply on CEC2', Sol Kim, 05 Jan 2022
Dear Juan A. Añel,
Thank you for guiding us through this process. We have added the license to the linked code. On second thought, adding the ARTMIP repository does not seem necessary - we do not use any of the outputs directly from ARTMIP in this work; it was included for readers to have access to previous version outputs of the tARget algorithm and for comparing to other atmospheric river detection algorithms. We only used the publically available E3SM and MERRA2 datasets in combination with the linked tARget v3 algorithm to produce all of our results. Given that both the code and data necessary to run our analysis are available, does this satisfy the policy? We remove mention of ARTMIP and our statement will read as below:
"Code and data availability. E3SM and MERRA2 datasets used in this study are publically available at https://esgf-node.llnl.gov/projects/e3sm/ and https://disc.gsfc.nasa.gov/datasets?project=MERRA-2 respectively. The tARget v3 algorithm is available at https://doi.org/10.25346/S6/B89KXF."
- CEC3: 'Reply on AC2', Juan Antonio Añel, 06 Jan 2022
-
AC2: 'Reply on CEC2', Sol Kim, 05 Jan 2022
-
CEC2: 'Reply on AC1', Juan Antonio Añel, 04 Jan 2022
-
AC1: 'Reply on CEC1', Sol Kim, 03 Jan 2022
-
RC1: 'Comment on gmd-2021-364', Anonymous Referee #1, 10 Jan 2022
The authors evaluated the simulations of atmospheric rivers in the model E3SM and explored the relevant physical and dynamical processes. However, I have some major concerns as well as many other comments listed below. Before the author solve these concerns, this manuscript might not be ready to publish.
Major Concerns:
(1) Lines 112-120: The authors compared the AR frequency between MERRA2 and E3SM. Although the authors clarified that “All % differences mentioned below … are absolute differences, not relative difference”, the description is misleading. For example, “there is a slight positive bias (1-3%) in the E3SM frequency near the edge between the tropics and subtropics …” I do not think 1-3% (absolute difference in AR frequency) is a “slight positive bias”. For example, while the AR frequency is below 10% over the southern California coastal region (Figs. 1a-b), the 1-3% bias (Fig. 1c) is large. In other words, the relative difference/bias of AR frequency over that region is larger than 10-30%. The bias in Chile is even larger (Fig. 1c). Please rewrite this part and the potential reasons for the large bias should be discussed.
(2) Line 122 and line 130: “… exhibits a close match …” and “… in agreement with …” Similar to my last comment, please be careful with these vague descriptions. For example, as the authors mentioned that in NDJFM over the west coast of North America, the AR frequency difference is 3-4% (line 125). That means the relative difference is roughly > 30-40% over the US West Coast (one of the most important area that affected by ARs) during the AR season (NDJFM). These model bias should be carefully examined and described.
(3) Many statistical analyses (such as Figs 2 and 4) are conducted over the global domain. Those results are important. However, many useful information/signals might be smoothed out using the global domain. It would be helpful if the authors could conduct similar analyses for some regions with high frequency or high impacts of ARs, such as the west coast of North America in NDJFM. I believe many readers will be interested in the analysis for that kind of region (with large model bias as well as high social impacts), rather than a smoothed result for a global domain.
(4) I have some concerns and suggestions for Fig. 4. (a) In my understanding, in each panel the total probability of E3SM or MERRA2 should be 1. However, it seems like the total probability is much lower than 1. Please clarify. (b) In panel a, please put “x104” closer to the unit or the numbers in X-axis. (c) Please keep consistency in the number of decimals for the values of medians in each panel. (d) In caption, “sets of line” should be “sets of lines”. (e) Please extend the Y-axis in panel e since it seems like the maximum value is higher than the probability of 8x10-3 (is the probability values at Y-axis correct?).
(5) Lines 227-233: “A large source of general E3SM precipitation biases come from … suggesting certain large-scale circulation biases may have larger influence on AR frequency than the frequency of non-AR storms.” These sentences are important to interpret the potential mechanism responsible for the model bias in AR precipitation. However, it is difficult to follow the logic. For example, how did the authors conclude that “… certain large-scale circulation biases may have larger influence on AR frequency than the frequency of non-AR storms” without analyzing the frequency and precipitation rate of non-AR storms?
(6) Fig. 5 shows the comparison of AR precipitation rate. How about the total AR precipitation amount? For example, over the US West Coast E3SM has positive bias in both AR precipitation rate and AR frequency in NDJFM. I am curious how large the bias of total AR precipitation will be over there.
(7) I suggest the authors go through the manuscript to improve the writing. This manuscript would be easier to read if the authors could improve the writing. I listed some issues below, but there are more places that could be improved.
Other Comments:
(1) Line 58: What is “standard resolution”?
(2) Lines 58-63: These two sentences provide a general background for the E3SMv1 performance, but they are vague. Please re-write and provide more details.
(3) Line 79: is “daily data” daily mean or instantaneous?
(4) Line 80: “Five ensemble members …”. Please clarify the difference between the ensemble members, as well as the motivation to use the five ensemble members.
(5) Line 87: “AMIP”, please spell out the full name when it is used for the first time.
(6) Line 88: “CMIP6 DECK simulations”, please define DECK.
(7) Lines 105-106: “This means the threshold is calculated separately for MERRA2 and the E3SM simulation.” Is there any large difference in the IVT threshold between MERRA2 and E3SM? The difference in the IVT threshold (85th percentile) could be a part of the model bias.
(8) Line 110: “All % differences …” should be “All percentage differences …”
(9) Line 152: “SDs are consistent with MERRA2.” This sentence is too vague.
(10) Line 161: “… AR frequencies are well < 1.0 %.” Difficult to understand.
(11) Line 162: “The seasonal SDs reveal sources of the higher annual SDs.” Higher than which SDs? I saw the annual SD is obviously lower than the NDJFM and MJJAS SDs in Fig. 3.
(12) Lines 165-168: “MJJAS SDs (Fig. 3b peak for ~1.5 % over various regions of … ” Difficult to follow. Please rewrite.
(13) Line 168: “In general, the northern hemisphere shows more internal variability.” This is an interesting result, but do the authors have any idea about the potential reasons?
(14) Line 170: “… using a single historical simulation …” Why did the authors use a single historical simulation? How did the authors select the single simulation?
(15) Line 171: “The distribution of all the ARs …” Do the authors mean the “characteristics”?
(16) Line 172-173: “All characteristics show strong similarities in shape and peak at the same values, barring magnitude of mean IVT (4e).” I do not understand this sentence, what are the “same values”? I do not understand the logic to mention 4e (Fig. 4e?) here either?
(17) Lines 175-177: These two sentences are difficult to follow. Please rewrite.
(18) Lines 199-200: “For annual AR precipitation, E3SM reproduces the …” Please be careful, there might be large differences in the distributions and magnitudes of AR precipitation if using different precipitation data. “Reproduce” might be too vague.
(19) Line 202, “The western coasts … produce higher rates of AR precipitation in E3SM.” “Produce” might not be suitable here.
(20) Line 226-227: “… model bias in the subtropical jet would affect precipitation …, as both are influenced by the jet and storm tracks.” Please provide reference.
(21) Lines 250-252: “Although Dong et al. (2021) looked at future projections of large-scale circulation and precipitation changes … which can explain the sources of some of the AR biases.” I do not understand the logic to use the results of future projections from Dong et al. (2021) to explain “the sources of some of the AR biases” in the historical simulations in this study. Many factors and even mechanisms may change under climate change.
(22) Fig. 8: It would be helpful if the authors could add contours to show the distribution of AR frequency in the historical simulations.
-
AC3: 'Reply on RC1', Sol Kim, 19 Feb 2022
The comment was uploaded in the form of a supplement: https://gmd.copernicus.org/preprints/gmd-2021-364/gmd-2021-364-AC3-supplement.pdf
-
AC3: 'Reply on RC1', Sol Kim, 19 Feb 2022
-
RC2: 'Comment on gmd-2021-364', Anonymous Referee #2, 31 Jan 2022
In this manuscript, the authors compare historical simulations of global AR frequency in an Earth System Model (E3SM) to AR frequency in a reanalysis dataset (MERRA2). The authors highlight differences in model depictions of global AR frequency and provide some physical insights into these biases. The manuscript is detailed, and I appreciate the authors’ attempt to put the E3SM biases into context. However, I feel that the manuscript could benefit from several major changes prior to publication. I base my recommendation on the general comments listed below.
General Comments:
The methodology section is rather sparse on details, especially those related to E3SM, MERRA2, and the AR detection algorithm. For E3SM and reanalysis, which output fields are obtained? Is the daily data instantaneous (if so, at what time?) or daily-averaged? For the algorithm, a few additional details (e.g., regarding geometric criteria) would be helpful.
Regarding the choice of AR detection algorithm, recent intercomparison studies have shown that climate model simulations of AR activity vary based on choice of AR detection algorithm. While the authors are justified in choosing the Guan and Waliser algorithm, I think the authors should acknowledge possible uncertainty in results owing to the choice of algorithm.
I agree with Referee #1 that the characterization of differences in AR frequency should be in a relative sense vs. an absolute sense. Given that ARs only occur for at most 15-20% of the time steps, these absolute differences are considerable. I agree that this section should be overhauled to reflect the relative differences.
Regarding the choice of domain, while I agree with Referee #1 that a regional domain would be useful for more targeted studies, I think that the global domain is appropriate for the scope of this paper, which I interpret as a first step in understanding depictions of AR frequency in E3SM.
While the authors provide a good initial overview of the biases seen in E3SM, the authors do not provide much insight into potential pathways toward improving the model. I am specifically thinking of the parameterization of convection, PBL turbulence, microphysics, etc., which could have an influence on the depiction of water vapor transport. At the very least, the authors should note that the model physics could play a role in depicting ARs and precipitation, even if this is left to future research.
Specific Comments:
Lines 40-54: The authors may consider adding a recent article, O’Brien et al. (2021), which analyzes changes in AR counts and size in a future climate within the CMIP5/6 models:
O’Brien, Travis Allen and Wehner, Michael F and Payne, Ashley E. and Shields, Christine A and Rutz, Jonathan J. and Leung, L. Ruby and Ralph, F. Martin and Marquardt Collow, Allison B. and Guan, Bin and Lora, Juan Manuel and et al., (2021) Increases in Future AR Count and Size: Overview of the ARTMIP Tier 2 CMIP5/6 Experiment. JGR A. https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2021JD036013
Lines 78-95: Which fields are analyzed for the model and reanalysis datasets? At the bare minimum, the specific humidity and wind fields would be required to calculate IVT. Are there additional fields downloaded?
Lines 96-106: Given the importance of the choice of AR detection algorithm to this study, I suggest that the authors provide some additional details about the detection algorithm (even though the details have already been published elsewhere).
Lines 109-110: How are the ensemble mean AR frequencies calculated? Are ARs detected for each ensemble member and the frequencies averaged? Or are the ARs detected using the ensemble-average IVT?
Figure 1: While the absolute differences are important, relative differences are likely more interesting given the relative infrequency of ARs, even over the midlatitude storm track.
Figure 2: This is a very interesting plot, but visually, the data points are difficult to discern. There is a lot of “empty” space in these diagrams. Would it be possible to either a) make the data points larger or b) only show a fraction of each diagram? Zooming in somehow would be very helpful for the reader to help discern the slight differences in correlations between ensemble members.
Lines 165-168: Do features like the Indian Monsoon, large-scale tropical convection, or TCs get detected as ARs?
Lines 171-173: Are the feature-averaged values (namely, feature-averaged IVT) weighted by latitude?
Line 178: Is the feature’s centroid based on area alone, or is there a weighting based on IVT intensity?
Lines 184-185: I don’t agree with this statement. While this may be the case here, tropical ARs could have weaker moisture transport due to a lack of strong winds aloft. I suggest the authors clarify this statement.
Figure 7: Although the details are provided in the figure caption, I suggest adding labels to these panels for quicker/easier interpretation.
Lines 322-323: Are these individual AMIP and fully-coupled simulations or an ensemble average of each?
Figures 8-10: As with Figure 7, it would be helpful if the panels were labeled. Also, is AMIP subtracted from fully-coupled, or vice-versa?
Technical Corrections:
Line 51: “but only a few”
Lines 54, 60, 78: Parentheses around citation
Line 87: Spell out the AMIP abbreviation
Line 179: Use “fewer” rather than “less” here
Line 321: Is this “southern Africa”, as opposed to the country?
Line 344: Use “occur” rather than "occurs" here
-
AC3: 'Reply on RC1', Sol Kim, 19 Feb 2022
In this manuscript, the authors compare historical simulations of global AR frequency in an Earth System Model (E3SM) to AR frequency in a reanalysis dataset (MERRA2). The authors highlight differences in model depictions of global AR frequency and provide some physical insights into these biases. The manuscript is detailed, and I appreciate the authors’ attempt to put the E3SM biases into context. However, I feel that the manuscript could benefit from several major changes prior to publication. I base my recommendation on the general comments listed below.
General Comments:
The methodology section is rather sparse on details, especially those related to E3SM, MERRA2, and the AR detection algorithm. For E3SM and reanalysis, which output fields are obtained? Is the daily data instantaneous (if so, at what time?) or daily-averaged? For the algorithm, a few additional details (e.g., regarding geometric criteria) would be helpful.
Regarding the choice of AR detection algorithm, recent intercomparison studies have shown that climate model simulations of AR activity vary based on choice of AR detection algorithm. While the authors are justified in choosing the Guan and Waliser algorithm, I think the authors should acknowledge possible uncertainty in results owing to the choice of algorithm.
I agree with Referee #1 that the characterization of differences in AR frequency should be in a relative sense vs. an absolute sense. Given that ARs only occur for at most 15-20% of the time steps, these absolute differences are considerable. I agree that this section should be overhauled to reflect the relative differences.
Regarding the choice of domain, while I agree with Referee #1 that a regional domain would be useful for more targeted studies, I think that the global domain is appropriate for the scope of this paper, which I interpret as a first step in understanding depictions of AR frequency in E3SM.
While the authors provide a good initial overview of the biases seen in E3SM, the authors do not provide much insight into potential pathways toward improving the model. I am specifically thinking of the parameterization of convection, PBL turbulence, microphysics, etc., which could have an influence on the depiction of water vapor transport. At the very least, the authors should note that the model physics could play a role in depicting ARs and precipitation, even if this is left to future research.
Specific Comments:
Lines 40-54: The authors may consider adding a recent article, O’Brien et al. (2021), which analyzes changes in AR counts and size in a future climate within the CMIP5/6 models:
O’Brien, Travis Allen and Wehner, Michael F and Payne, Ashley E. and Shields, Christine A and Rutz, Jonathan J. and Leung, L. Ruby and Ralph, F. Martin and Marquardt Collow, Allison B. and Guan, Bin and Lora, Juan Manuel and et al., (2021) Increases in Future AR Count and Size: Overview of the ARTMIP Tier 2 CMIP5/6 Experiment. JGR A. https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2021JD036013
Lines 78-95: Which fields are analyzed for the model and reanalysis datasets? At the bare minimum, the specific humidity and wind fields would be required to calculate IVT. Are there additional fields downloaded?
Lines 96-106: Given the importance of the choice of AR detection algorithm to this study, I suggest that the authors provide some additional details about the detection algorithm (even though the details have already been published elsewhere).
Lines 109-110: How are the ensemble mean AR frequencies calculated? Are ARs detected for each ensemble member and the frequencies averaged? Or are the ARs detected using the ensemble-average IVT?
Figure 1: While the absolute differences are important, relative differences are likely more interesting given the relative infrequency of ARs, even over the midlatitude storm track.
Figure 2: This is a very interesting plot, but visually, the data points are difficult to discern. There is a lot of “empty” space in these diagrams. Would it be possible to either a) make the data points larger or b) only show a fraction of each diagram? Zooming in somehow would be very helpful for the reader to help discern the slight differences in correlations between ensemble members.
Lines 165-168: Do features like the Indian Monsoon, large-scale tropical convection, or TCs get detected as ARs?
Lines 171-173: Are the feature-averaged values (namely, feature-averaged IVT) weighted by latitude?
Line 178: Is the feature’s centroid based on area alone, or is there a weighting based on IVT intensity?
Lines 184-185: I don’t agree with this statement. While this may be the case here, tropical ARs could have weaker moisture transport due to a lack of strong winds aloft. I suggest the authors clarify this statement.
Figure 7: Although the details are provided in the figure caption, I suggest adding labels to these panels for quicker/easier interpretation.
Lines 322-323: Are these individual AMIP and fully-coupled simulations or an ensemble average of each?
Figures 8-10: As with Figure 7, it would be helpful if the panels were labeled. Also, is AMIP subtracted from fully-coupled, or vice-versa?
Technical Corrections:
Line 51: “but only a few”
Lines 54, 60, 78: Parentheses around citation
Line 87: Spell out the AMIP abbreviation
Line 179: Use “fewer” rather than “less” here
Line 321: Is this “southern Africa”, as opposed to the country?
Line 344: Use “occur” rather than "occurs" here
-
AC3: 'Reply on RC1', Sol Kim, 19 Feb 2022
Sol Kim et al.
Sol Kim et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
664 | 154 | 21 | 839 | 8 | 6 |
- HTML: 664
- PDF: 154
- XML: 21
- Total: 839
- BibTeX: 8
- EndNote: 6
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1