Page 42 - gt1506
P. 42

Moving lithospheric modeling forward: Attributes of a
                       community computer code

GSA TODAY | JUNE 2015  C.M. Cooper, Washington State University, School of the                into numerical models that are reliable, portable, and computa-
                       Environment, P.O. Box 624812, Pullman, Washington 99164-2812,          tionally efficient.
                       USA; Eric Mittelstaedt, University of Idaho, Dept. of Geological
                       Sciences, 875 Perimeter Drive, MS 3022, Moscow, Idaho 83844-             Lithospheric modelers are confronted with a broad range of
                       3022, USA; Claire A. Currie, University of Alberta, Dept. of Physics,  challenges to address these drivers. Scientifically, crucial geolog-
                       Edmonton, Alberta, Canada T6G 2E1; Jolante van Wijk, New               ical processes lack theoretical or empirical descriptions (e.g., vari-
                       Mexico Institute of Mining and Technology, Dept. of Earth &            able fault dip at depth, spacing in shear bands, localization of
                       Environmental Science, 801 Leroy Place, Socorro, New Mexico            deformation, and coupled deformation with melting and melt
                       87801, USA; Louise K. Kellogg, Lorraine Hwang, University of           migration). Incorporating the vast quantity of new data available
                       California Davis, Earth and Planetary Sciences, Computational          through such initiatives as the National Science Foundation’s
                       Infrastructure for Geodynamics, 2215 Earth and Physical Sciences,      EarthScope and data compilations such as Gplates (Qin et al.,
                       One Shields Avenue, Davis, California 95616, USA; and Ramon            2012) and PetDB (Lehnert et al., 2000) requires both the develop-
                       Arrowsmith, Arizona State University, School of Earth & Space          ment of new data-handling methods and an understanding of
                       Exploration, P.O. Box 876004, Tempe, Arizona 85287-6004, USA           their interrelationships. Added to these challenges are the difficul-
                                                                                              ties in implementing the numerical methods required to run the
                         We live on a planet with an active surface that is modified and      desired simulations, including modeling systems with large-
                       deformed at multiple temporal and spatial scales owing to diverse      magnitude variations in material properties occurring over short
                       processes occurring at plate boundaries and plate interiors. The       spatial scales; maintaining discrete material boundaries as the
                       processes of mid-ocean-ridge spreading, mountain building,             model evolves; and incorporating realistic fault evolution and
                       subduction of tectonic plates, mantle drag, intra-continental          faulting behavior. Lastly, extending models to three dimensions
                       deformation, earthquakes, and volcanism cross traditional disci-       increases the numerical and model complexity, an area that has
                       plinary boundaries (Fig. 1A). Understanding these lithospheric         seen limited development.
                       processes is valuable not only for intellectual curiosity and to
                       refine our working knowledge of plate tectonics, but also for            Modeling complex systems requires validation and verification
                       understanding threats to life, property, and infrastructure.           of software. Establishing and running benchmarks and test suites
                       Computer modeling and simulation are increasingly powerful             not only “proves” a code, it also provides important insight to the
                       tools that researchers employ to better understand lithospheric        researcher. Limits in parameter space and trade-offs between
                       deformation and unravel the complex feedbacks that drive the           different model specifications become better known.
                       evolution of Earth’s surface. The field is poised for a significant    Benchmarking performance helps to inform the use of computa-
                       advance to take advantage of recent expansions in computing            tional resources and to understand numerical uncertainty.
                       power, improved representation of idealized processes, increased
                       data availability, and better communication between software             The heterogeneity of the lithosphere translates to a heteroge-
                       developers and geoscientists.                                          neous approach to modeling lithospheric processes. Computational
                                                                                              approaches employed to address the key scientific interests of the
                         To move forward as a community, we must address key scien-           community tend to be based either in continuum, analytical, or
                       tific drivers motivating present and future lithospheric deforma-      discontinuous methods (Fig. 1B). Usage of these different math-
                       tion research. The scientific processes to incorporate include         ematical methods, several of which may be deployed in any one
                       melting and melt transport, strain localization and de-localiza-       code, depends on the maturity of the research area and the
                       tion, surface processes (e.g., erosion and deposition), and mantle-    specifics of the research question. Individual researchers will
                       lithosphere interaction. Understanding these requires the              often develop numerical techniques and modeling software
                       integration of results from seismic imaging, the earthquake cycle,     capable of solving specific geologic problems. While these efforts
                       plate boundary evolution, and more realistic Earth-like rheologies     often result in numerical codes that are powerful and apt for the
                                                                                              problem at hand, they often do not translate into a more universal
                                                                                              modeling tool.

                           GSA Today, v. 25, no. 5, doi: 10.1130/GSATG230GW.1.

                           E-mails: Cooper: cmcooper@wsu.edu; Mittelstaedt: emittelstaedt@uidaho.edu; Currie: claire.currie@ualberta.ca; van Wijk: jvanwijk@ees.nmt.edu; Kellogg: kellogg@
                           ucdavis.edu; Hwang: lorraine@geodynamics.org; Arrowsmith: ramon.arrowsmith@asu.edu.

                     42
   37   38   39   40   41   42   43   44   45   46   47