the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Coupling framework (1.0) for the Úa (2023b) ice sheet model and the FESOM-1.4 z-coordinate ocean model in an Antarctic domain
Abstract. The rate at which the Antarctic Ice Sheet loses mass is to a large degree controlled by ice-ocean interactions underneath small ice shelves, with the most sensitive regions concentrated in even smaller areas near grounding lines and local pinning points. Sufficient horizontal resolution is key to resolving critical ice-ocean processes in these regions, but difficult to afford in large-scale models used to predict the coupled response of the entire Antarctic Ice Sheet and the global ocean to climate change. In this study we describe the implementation of a framework that couples the ice sheet flow model Úa with the Finite Element Sea Ice Ocean Model (FESOM-1.4) in a configuration using depth-dependent vertical coordinates. The novelty of this approach is the use of horizontally unstructured grids in both model components, allowing us to resolve critical processes directly, while keeping computational demands within the range of feasibility. We use the Marine Ice Sheet–Ocean Model Intercomparison Project framework to verify that ice retreat and readvance is reliably simulated, and inaccuracies in mass, heat and salt conservation are small compared to the forcing signal. Further, we demonstrate the capabilities of our approach for a global ocean/Antarctic Ice Sheet domain. In a 39-year hindcast simulation (1979–2018) we resolve retreat behaviour of Pine Island Glacier, a known challenge for coarser resolution models. We conclude that Úa-FESOM is well suited to improve predictions of the Antarctic Ice Sheet evolution over centennial time scales.
- Preprint
(2595 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2024-648', Anonymous Referee #1, 06 May 2024
This paper describes the coupling framework for the Úa ice sheet and FESOM ice-ocean models. This coupling is particularly interesting because both models use unstructured meshes and mesh refinement to resolve better local processes in otherwise global or regional setups. The paper is well-structured and written, and I expect it to be of substantial interest to GMD readers. I find the paper quite mature, and I recommend publishing it once the authors have addressed my suggestions and minor concerns.
Sequential coupling: You use a sequential coupling, which essentially doubles the runtime of the whole system, at least at first guess. Yet, there is no discussion of why you don't use a parallel approach. You need to give at least the reason or justification for why you went with the sequential approach. An estimate of the performance benefits of using a parallel coupling would also benefit the paper.
I don't like your names for sections 2.4 and 2.5, as they could be more descriptive; please consider giving these more descriptive names.
L116: The description of the background mesh left me hanging. How do you decide what regions "could possibly unground during the simulation period"? You describe this later, but more details are in order here. Even saying that it's simulation-specific and will be described in more detail later would suffice (if this is the case).
In section names for sections 3 and 3.3, you use "verification" when "evaluation" would be more fitting. You can only verify the model results if you have an analytical solution.
In section 3.1, I would appreciate more justification for choosing the IceOcean1ra experiment. I expect that there are several idealised experiments in the MISOMIP framework, so you should say why you chose this particular one.
L160: "... with the later described pan-Arctic setup" should be "... with the pan-Arctic setup described later".
L161: "summaries" should be "summarised".
In section 3.3, I miss a reference to what results other modelling groups get for the IceOcean1ra setup. Are your results similar to those of others? Is there a large spread in the results in general?
In section 4.4, you say that you use three different machines to run the model. Such a setup is very unusual, and it would be nice to have more details. Is everything automated, or is there some manual work involved? Do all the machines have access to a shared storage area, or do you need to copy data between machines? Why did you choose this setup?
You only mention Greenland in the "summary and conclusions" sections. You should either remove this or also mention it in the main text. As it is, it comes completely out of nowhere.
Citation: https://doi.org/10.5194/egusphere-2024-648-RC1 -
RC2: 'Comment on egusphere-2024-648', Anonymous Referee #2, 27 Jun 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-648/egusphere-2024-648-RC2-supplement.pdf
-
AC1: 'Comment on egusphere-2024-648', Ole Richter, 10 Sep 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-648/egusphere-2024-648-AC1-supplement.pdf
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
263 | 60 | 30 | 353 | 22 | 24 |
- HTML: 263
- PDF: 60
- XML: 30
- Total: 353
- BibTeX: 22
- EndNote: 24
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1