SHELLS: A thin-shell program for modeling neotectonics of regional or global lithosphere with faults

摘要: The thin-plate method of modeling neotectonics uses isostasy and vertical integration lithospheric strength to reduce three-dimensional problems two dimensions. We introduce new thin-shell continuum elements fault which satisfy the completeness compatibility requirements, permit extension these methods lithospheres spherical planets. Even a coarse grid with elastic rheology can reproduce low-order toroidal free oscillations. In program SHELLS, realistic frictional/dislocation-creep is handled by iteration; this converges solutions be numerically tested confirm that they both local equilibrium balance torques on whole plates. Because code incorporate natural plate shapes (with internal faults) rheologies, it yields models are readily their predictions geodetic velocities, stresses, slip rates.