the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 29 Oct 2024
CitcomSVE 3.0: A Three-dimensional Finite Element Software Package for Modeling Load-induced Deformation for an Earth with Viscoelastic and Compressible Mantle
Abstract. Earth and other terrestrial and icy planetary bodies deform visco-elastically under various forces. Numerical modeling plays a critical role in understanding the nature of various dynamic deformation processes. This article introduces a newly developed, open-source package, CitcomSVE-3.0, which efficiently solves the visco-elastic deformation of planetary bodies. Based on its predecessor, CitcomSVE-2.1, CitcomSVE-3.0 is updated to account for elastic compressibility and depth-dependent density, which are particularly important in modeling horizontal displacement for visco-elastic deformation. We benchmark CitcomSVE-3.0 against a semi-analytical code for two types of surface loading problems: 1) single harmonic loads on the surface and 2) the glacial isostatic adjustment (GIA) problem with a realistic ice sheet loading history (ICE-6G_D), in which an updated version of sea level equations is incorporated. The benchmark results presented here demonstrate the accuracy and efficiency of this package. CitcomSVE shows a second-order accuracy in terms of spatial resolution. For a typical GIA modeling with 122-ky glaciation-deglaciation history, surface horizontal resolution of ~50 km, and time increment of 125 yr, it takes ~ 3 hours on 384 CPU cores to complete with less than 5 % errors in displacement rates.
- Preprint
(4379 KB) - Metadata XML
-
Supplement
(2715 KB) - BibTeX
- EndNote
Status: open (until 24 Dec 2024)
-
RC1: 'Comment on egusphere-2024-3200', Volker Klemann, 25 Nov 2024
reply
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-3200/egusphere-2024-3200-RC1-supplement.pdf
-
RC2: 'Comment on egusphere-2024-3200', Anonymous Referee #2, 26 Nov 2024
reply
The paper by Yuan et al. presents the next step of the GIA modelling setup using CitcomS. While the method has been published previously (SLE solver in Zhong et al. and compressibility in A et al.), here they present a benchmark of the two combinations. The paper summarizes the applied methods very well. However, I felt sometimes a bit lost in the paper as you introduce a method/set-up and then something else is mentioned afterwards. You come back to the method/set-up at a later point. I would suggest splitting the paper in two parts: 1) introducing and benchmarking your model to simple harmonic loads (where no SLE is involved), and 2) the benchmark with the SLE and the ICE-6G_D ice model, which means that you introduce the SLE part after you have presented the first benchmark.
A few other things could be done to convince the reader that your GIA modelling setup works:
- Show a comparison of the geoid change in the benchmark with the SLE and ICE-6G_D.
- Compare your radial and horizontal displacement rates as well as geoid and RSL rates to the results by Peltier et al. (which you can find on his website).
Minor comments:
L.73 - Huang et al. is not a coupled spectral-finite element method as it is purely based on the finite-element method. In addition, the method by Wu is also based on the finite-element method and only the potential calculations are done in the spherical harmonic domain. Not sure, if this should be then called coupled spectral-FE method. If I understand you correctly, you are doing the same by "coupling" the spherical harmonic domain with the FE method.
L.77 - "center of mass".
L.144 - You mention that the Lamé parameters and viscosity can vary laterally and radially. Is this also true for the density? I couldn't find such a statement.
L.183 - "the applied potential, which is only relevant"
L.185 - "on the unknown incremental"
L.195-202 - Please move this paragraph somewhere else. It doesn't fit here.
L.212 - What "given threshold error tolerance" is implemented in CitcomSVE 2.1 & 3.0?
L.301 - Interesting approach. Please write it out in an equation how O(t) is calculated.
L.324 - Why do use the density and Lamé parameters from the mantle for the crust?
Section 3.2.1 - You refer only to figures in the supplement. Either move the text to the supplement and summarize the results in the main article or move figures from the supplement to the main manuscript. I would prefer the former.
Section 3.2.1 - Just out of curisosity: Could you please plot K2 in Fig. S1 as well? So, the results of the 2nd iteration based on the approach by Kendall et al. No need to do this for the other figures.
L.419 - In the abstract you say ICE-6G_D, but here it is ICE-6G now. Please use the same and correct abbreviations throughout the text.
L.546 - Please use more representative names for your model types, e.g. GIA_R1 could be GIA_135 to indicate that this model has a resolution of 135 km and the reader doesn't have to go back and check what R1 was referring to.
Table 3 - I find it difficult to compare the CPU hours when you use different amounts of cores. Please unify this to see the direct difference.
Fig.S2/S3 - Use the same abbreviations for both figures.
Figures – Copernicus recommends the authors to use colorblind-friendly colors/colormaps, which is stated on the GMD website (https://www.geoscientific-model-development.net/submission.html#figurestables). This includes avoiding the parallel usage of red and green as well as the avoiding jet and rainbow colormaps. The link above provides examples where you can find suitable colormaps and find colors for your line plots that are suitable.
Citation: https://doi.org/10.5194/egusphere-2024-3200-RC2 -
CEC1: 'Comment on egusphere-2024-3200 - No compliance with the policy of the journal', Juan Antonio Añel, 04 Dec 2024
reply
Dear authors,
Unfortunately, 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.html
You have archived your code on two servers that do not comply with our policy, namely GitHub (the CitcomSVE model) and a server hosted by the University of Toronto (the ICE-6G model). In the case of GitHub, itself instructs authors to use other long-term archival and publishing alternatives, such as Zenodo. Therefore, the current situation with your manuscript is irregular. Please, publish your code in one of the appropriate repositories and reply to this comment with the relevant information (links and a permanent identifier for it (e.g. DOIs)) as soon as possible, as we can not accept manuscripts in Discussions that do not comply with our policy.Also, you must include a modified 'Code and Data Availability' section in a potentially reviewed manuscript, containing the DOI of the new repositories.
Please, note that if you do not fix this problem, we will have to reject your manuscript for publication in our journal.
Juan A. Añel
Geosci. Model Dev. Executive EditorCitation: https://doi.org/10.5194/egusphere-2024-3200-CEC1 -
AC1: 'Reply on CEC1', Tao Yuan, 05 Dec 2024
reply
Dear Juan A. Añel,
Thanks for informing us of our mistakes. We uploaded the CitcomSVE version used in this manuscript to Zenodo (https://doi.org/10.5281/zenodo.14277856). The ICE-6G model has also been uploaded into the same Zenodo repository (see the description of this repository). Thus, all the models (code) and data used or generated for this study are deposited in the repository.
We will update the Code and Data Availability section with the following statement in this revision: “The current version of CitcomSVE3.0 is available from GitHub: https://github.com/shjzhong/CitcomSVE under the GPL licence. The exact version of the model used to produce the results used in this paper is archived on Zenodo (10.5281/zenodo.13932410), as are input data (including the ice model and Earth model used in this paper) and scripts to run the model and produce the plots for all the calculations presented in this paper.”
Best,
Tao Yuan
Citation: https://doi.org/10.5194/egusphere-2024-3200-AC1 -
CEC2: 'Reply on AC1', Juan Antonio Añel, 05 Dec 2024
reply
Dear authors,
Thanks for addressing this issue so quickly. We can consider now the current version of your manuscript compliant with our policy.
Regards,
Juan A. Añel
Geosci. Model Dev. Executive Editor
Citation: https://doi.org/10.5194/egusphere-2024-3200-CEC2
-
CEC2: 'Reply on AC1', Juan Antonio Añel, 05 Dec 2024
reply
-
AC1: 'Reply on CEC1', Tao Yuan, 05 Dec 2024
reply
Data sets
Dataset for CitcomSVE 3.0: A Three-dimensional Finite Element Software Package for Modeling Load-induced Deformation for an Earth with Viscoelastic and Compressible Mantle Tao Yuan https://doi.org/10.5281/zenodo.13932411
Model code and software
CitcomSVE Shijie Zhong https://github.com/shjzhong/CitcomSVE
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 177 | 46 | 11 | 234 | 23 | 3 | 3 |
- HTML: 177
- PDF: 46
- XML: 11
- Total: 234
- Supplement: 23
- BibTeX: 3
- EndNote: 3
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1