A PYTHON PIPELINE FOR RAPID INVERSION OF GRAVITY DATA AIMED AT EARLY CAREER RESEARCHERS: APPLICATION TO THE LOST RIVER FAULT, IDAHO

Gravity data and modeling can be used to create a 3D view of basin and fault geometries, a valuable tool for seismic hazard assessment. Additionally, access to a high-speed, easy-to-use modeling workflow will enable more researchers to supplement their work with a new method and allow students to be exposed to modeling earlier in their academic careers. The primary goals of this work are to create a widely usable 3D gravity and basin model, to create an introduction to modeling teaching module for undergraduate geology students, and to apply the model to the northernmost two segments of the Lost River Fault in Idaho in order to constrain its dip and the sediment thickness of the overlying basin.

The model is written in Python, which is popular among geoscientists for being versatile, beginner-friendly, and open-source. It uses several powerful libraries, including pandas for importing and organizing data, pygmt for producing maps, and multiprocessing to decrease runtime by a factor of ~2-8. The model reads a CSV file containing the UTM coordinates and Bouguer gravity anomaly values, and the user configures parameters within the code, such as grid spacing and density contrast. Without any further input from the user, the code will generate several gravity maps, a sediment thickness map, and a basin floor slope map. It will also create a directory structure for saving the modeled data and figures to minimize the work required by the user before running the model. The ease of use and high speed make this a valuable tool for researchers and for introducing students in courses such as structural geology or geophysics to modeling with Python. The user can test different parameters and quickly compare the results to learn how the model is affected by each.

We apply this model to the northern Lost River Fault with an assumed density contrast of -450 kg/m³ between the Quaternary alluvium basin fill and Paleozoic basement rock. The modeled basin sediment thickness varies rapidly over short distances with steep gradients along the basin-bounding faults and in the middle of the valley. This pattern suggests that the basin floor retains paleotopographic features and that there are buried Eocene Challis volcanics with densities similar to those of the sediment fill, as evidenced by exposures at the surface.