This year's ISC Student Cluster Competition (ISC20-SCC) is special - it is for the first time (to our knolwedge) a virtual one. The other exciting fact about this high-profile competition, where computer science students from all over the planet compete in optimizing performance of real-world HPC applications, is that Elmer/Ice is included in this list of applications. The challenge involves a steady-state high-resolution model full-Stokes solution of the Greenland ice sheet
ISC20-SCC will end on June 17 with the award ceremony taking place on June 24. We wish all participants good luck in the competition!
The flow of large ice sheets and glaciers can be simulated by solving the full Stokes equations using a finite element method. The simulation is particularly sensitive to the discretization of the grounding line, which separates ice resting on bedrock and ice floating on water, and is moving with time. The boundary conditions at the ice base are enforced by Nitsche's method and a subgrid treatment of the grounding line element. Simulations with the method in two dimensions for an advancing and a retreating grounding line illustrate the performance of the method. The computed grounding line position is compared to previously published data with a fine mesh, showing that similar accuracy is obtained using subgrid modeling with more than 20-times-coarser meshes. This subgrid scheme is implemented in the two-dimensional version of the open-source code Elmer/Ice.
Read more: Cheng, G., P. Lötstedt and L. von Sydow, 2020. A full Stokes subgrid scheme in two dimensions for simulation of grounding line migration in ice sheets using Elmer/Ice (v8.3), Geosci. Model Dev., 13, 2245–2258, doi:10.5194/gmd-13-2245-2020
Interpretation of greenhouse gas records in polar ice cores requires a good understanding of the mechanisms controlling gas trapping in polar ice, and therefore of the processes of densification and pore closure in firn (compacted snow). Current firn densification models are based on a macroscopic description of the firn and rely on empirical laws and/or idealized geometries to obtain the equations governing the densification and pore closure.
Here, we propose a physically-based methodology explicitly representing the porous structure and its evolution over time. In order to handle the complex geometry and topological changes that occur during firn densification, we rely on a Level-Set representation of the interface between the ice and the pores. Two mechanisms are considered for the displacement of the interface: (i) mass surface diffusion driven by local pore curvature and (ii) ice dislocation creep. For the latter, ice is modeled as a viscous material and the flow velocities are solutions of the Stokes equations.
First applications show that the model is able to densify firn and split pores. Using the model in cold and arid conditions of the Antarctic plateau, we show that gas trapping models do not have to consider the reduced compressibility of closed pores compared to open pores in the deepest part of firns. Our results also suggest that the mechanism of curvature-driven surface diffusion does not result in pore splitting, and that ice creep has to be taken into account for pores to close. Future applications of this type of model could help quantify the evolution and closure of firn porous networks for various accumulation and temperature conditions.
Read More: Fourteau, K., F. Gillet-Chaulet, P. Martinerie and X. Faïn, 2020. A Micro-Mechanical Model for the Transformation of Dry Polar Firn Into Ice Using the Level-Set Method, Frontiers in Earth Science, 8, doi:10.3389/feart.2020.00101