A PYTHON PIPELINE FOR RAPID INVERSION OF GRAVITY DATA AIMED AT EARLY CAREER RESEARCHERS: APPLICATION TO THE LOST RIVER FAULT, IDAHO
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/m3 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.