Open-source modular solutions for flexural isostasy: gFlex v1.0

Abstract

Abstract. Isostasy is one of the oldest and most widely applied concepts in the geosciences, but the geoscientific community lacks a coherent, easy-to-use tool to simulate flexure of a realistic (i.e., laterally heterogeneous) lithosphere under an arbitrary set of surface loads. Such a model is needed for studies of mountain building, sedimentary basin formation, glaciation, sea-level change, and other tectonic, geodynamic, and surface processes. Here I present gFlex (for GNU flexure), an open-source model that can produce analytical and finite difference solutions for lithospheric flexure in one (profile) and two (map view) dimensions. To simulate the flexural isostatic response to an imposed load, it can be used by itself or within GRASS GIS for better integration with field data. gFlex is also a component with the Community Surface Dynamics Modeling System (CSDMS) and Landlab modeling frameworks for coupling with a wide range of Earth-surfacerelated models, and can be coupled to additional models within Python scripts. As an example of this in-script coupling, I simulate the effects of spatially variable lithospheric thickness on a modeled Iceland ice cap. Finite difference solutions in gFlex can use any of five types of boundary conditions: 0-displacement, 0-slope (i.e., clamped); 0-slope, 0-shear; 0-moment, 0-shear (i.e., broken plate); mirror symmetry; and periodic. Typical calculations with gFlex require 1 s to ∼ 1 min on a personal laptop computer. These characteristics -multiple ways to run the model, multiple solution methods, multiple boundary conditions, and short compute time -make gFlex an effective tool for flexural isostatic modeling across the geosciences.

Tectono‐sedimentary evolution of the northern Iranian Plateau: insights from middle–late Miocene foreland‐basin deposits

“…To further test this hypothesis we performed flexural basin modelling (e.g. Wickert, 2015) using our stratigraphic information and assuming a 60-km-wide mountain range with~2.0 km of elevation above the surrounding basins (e.g. the proto Takab Range Complex; Appendix S5).…”

Widespread glacial erosion on the Scandinavian passive margin

“…We removed the offshore late Pliocene-Quaternary deposits and filled in onshore valleys and fjords using our erosion estimates from geophysical relief. We calculated flexural isostatic adjustments from these load changes using the open-source gFlex version 1.0 model (Wickert, 2016). We adopted a Young's modulus of 70 GPa, Pois-son ratio of 0.25, and densities of 1029 kg/m 3 , 2300 kg/m 3 , 2670 kg/m 3 , and 3300 kg/m 3 for water, sediment, eroded bedrock, and mantle, respectively.…”

Short communication: Landlab v2.0: a software package for Earth surface dynamics

“…Writing effective unit tests that ensure Landlab components reliably solve their equations under a variety of initial and boundary conditions is not a trivial task. When a set of equations that a component solves have an analytical solution then the numerics of a component can be verified based on the ability to reproduce such a relationship (e.g., stream power erosion produces a known slope-area relationship, Willgoose et al, 1991). When such analytical predictions do not exist-as is often the case-a more detailed analysis of the equations must be performed in order to create a full verification test.…”