Preprints
https://doi.org/10.5194/egusphere-2023-2491
https://doi.org/10.5194/egusphere-2023-2491
02 Nov 2023
 | 02 Nov 2023

A Fast Surrogate Model for 3D-Earth Glacial Isostatic Adjustment using Tensorflow (v2.8.10) Artificial Neural Networks

Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev

Abstract. Models of glacial isostatic adjustment (GIA) play a central role in the interpretation of various geologic and geodetic data to understand and simulate past and future changes in ice sheets and sea level, and infer rheological properties of the deep Earth. A relatively recent advance has been the development of models that include 3D Earth structure, as opposed to 1D, spherically symmetric structure. However, a major limitation in employing 3D GIA models is their high computational expense. As such, we have developed a method using artificial neural networks (ANNs) and the Tensorflow library to emulate the influence of 3D Earth models with the goal of more affordably constraining the parameter space of these models: specifically the radial (1D) viscosity profile upon which the lateral variations are added.

This study provides an initial “proof of concept” assessment of using ANNs to emulate the influence of lateral Earth structure on GIA model output. Our goal is to test whether the fast surrogate model can accurately predict the difference in these outputs (i.e., RSL and uplift rates) for the 3D case relative to the SS case. If so, the surrogate model can be used with a computationally efficient SS (Earth) GIA model to generate output that reproduces output from a 3D (Earth) GIA model. Evaluation of the surrogate model performance for deglacial RSL indicates that it is able to provide useful estimates of this field throughout the parameter space when trained on only ≈ 15 % (≈ 50) of the parameter vectors considered (330 in total). Our results indicate that the ANN:model misfits, while not negligible, are of a scale such that useful predictions of deglacial RSL changes can be made.

We applied the surrogate model in a model:data comparison exercise using RSL data distributed along the North American coasts from the Canadian Arctic to the US Gulf coast. We find that the surrogate model is able to successfully reproduce the data:model misfit values such that the region of minimum misfit either overlaps the 3D GIA model results, or is within two increments in the parameter space. The surrogate model can, therefore, be used to accurately explore this aspect of the 3D Earth model parameter space. While the 3D Earth models can outperform the SS Earth models for some regional subsets of the RSL data set, the SS Earth models still produce better fits overall. In summary, this work demonstrates the utility of machine learning in 3D Earth GIA modelling and so future work to expand on this analysis is warranted.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this preprint. The responsibility to include appropriate place names lies with the authors.
Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev

Status: final response (author comments only)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on egusphere-2023-2491', Anonymous Referee #1, 30 Dec 2023
    • AC1: 'Reply on RC1', Ryan Love, 15 Feb 2024
  • RC2: 'Comment on egusphere-2023-2491', Wouter van der Wal, 08 Jan 2024
    • AC2: 'Reply on RC2', Ryan Love, 15 Feb 2024
  • EC1: 'Comment on egusphere-2023-2491', Andrew Wickert, 26 Jan 2024
Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev

Data sets

Input Data for A Fast Surrogate Model for 3D-Earth Glacial Isostatic Adjustment using Tensorflow (v2.8.10) Artificial Neural Networks Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev https://zenodo.org/records/10042047

Model code and software

Supplemental Materials for A Fast Surrogate Model for 3D-Earth Glacial Isostatic Adjustment using Tensorflow (v2.8.10) Artificial Neural Networks Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev https://zenodo.org/records/10045463

Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev

Viewed

Total article views: 434 (including HTML, PDF, and XML)
HTML PDF XML Total Supplement BibTeX EndNote
308 99 27 434 45 12 13
  • HTML: 308
  • PDF: 99
  • XML: 27
  • Total: 434
  • Supplement: 45
  • BibTeX: 12
  • EndNote: 13
Views and downloads (calculated since 02 Nov 2023)
Cumulative views and downloads (calculated since 02 Nov 2023)

Viewed (geographical distribution)

Total article views: 429 (including HTML, PDF, and XML) Thereof 429 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 
Latest update: 22 May 2024
Download
Short summary
A relatively recent advance in glacial isostatic adjustment modelling has been the development of models that include 3D Earth structure, as opposed to 1D structure. However, a major limitation is the computational expense. We have developed a method using artificial neural networks to emulate the influence of 3D Earth models to affordably constrain the viscosity parameter space. Our results indicate that the misfits are of a scale such that useful predictions of relative sea level can be made.