USGS Science⌕ Search
Home/Archive/USGS

USGS · 70211326

Distribution of earthquakes on a branching fault system using integer programming and greedy sequential methods

Eric L. Geist·Thomas E. Parsons
View original source ↗Export metadata ↓

Abstract

A new global optimization method is used to determine the distribution of earthquakes on a complex, connected fault system. The method, integer programming, has been advanced in the field of operations research, but has not been widely applied to geophysical problems until recently. In this application, we determine the optimal distribution of earthquakes on mapped faults to minimize the global misfit in slip rates for multi-fault ruptures. Integer programming solves for a decision vector composed of every possible location that a sample of earthquakes can occur on every fault, subject to slip-rate uncertainty constraints. Step over connections are straightforward to include, whereas branching fault connections are not. To include branching ruptures, we distinguish between individual multi-fault rupture paths, as opposed to formulating the integer-programming problem based on individual faults as in previous studies. The new method is applied to the complex fault system in the San Francisco Bay Area as a case study. Results from the integer-programming method are compared to those from a local optimization method, termed the greedy-sequential method. Several experiments using these two methods indicate that shape of the on-fault magnitude distributions and which branching faults are involved in multi-fault ruptures depend on how much emphasis is placed on fitting the target slip rate. In cases where the underlying data are not strong enough to warrant chasing the target slip rate, it is better to focus on the distribution of feasible results that better represents the uncertainty in the solutions imposed by the data.

Explore related subjects

Geochemistry, Geophysics, GeosystemsArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Eric L. Geist, Thomas E. Parsons. 2020-08-27. Distribution of earthquakes on a branching fault system using integer programming and greedy sequential methods. https://doi.org/10.1029/2020gc008964

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Erupted volume during the 2014–2023 Semisopochnoi unrest and eruption constrained by multi-platform InSAR and volcanic gas observations

Volcano deformation data provide important constraints on magma dynamics during volcanic unrest and eruptions. From these observations, we can estimate subsurface pressure changes and relate them to volume change in the reservoir. However, because magma and rock are compressible, this volume change cannot be directly equated with eruptive volume. When we combine geodetic analyses with gas observations, petrologic constraints, volatile solubility modeling, and magma compressibility equations, however, we can estimate realistic lower and upper bounds of eruptive volume. Leveraging multi‐platform InSAR time series and satellite SO 2 data, we apply these methods to the 2014–2023 unrest and eruption of Semisopochnoi, Alaska. We find that this phreatomagmatic eruption sequence was fed by a spherical magma source at 6.5–7.5 km depth bsl, centered beneath the caldera, ∼5 km offset to the east from the vent. Volatile solubility modeling indicates that ∼50% of the initial dissolved CO 2 has exsolved at depths consistent with the geodetic source depth, increasing the magma compressibility. Combining geodetic model constraints with SO 2 observations, petrologic constraints, solubility models, and magma compressibility findings yields 12–32 Mm 3 as the most‐likely juvenile, dense‐rock equivalent eruptive volume, with a best estimate of 20 Mm 3 . Stress change modeling suggests that the local stress field encouraged conduit formation west of the pre‐eruptive inflation center, and that the 2014 M7.9 Rat Island earthquake, ∼65 km to the west of Semisopochnoi, reduced the compressive stress underneath the volcano, supporting prolonged eruptive activity.

Alaska↗
USGS2026

Pre-eruptive volatiles and magma storage depth for the 2018 eruption at Mount Veniaminof, Alaska

Mount Veniaminof is a frequently active, open-system arc volcano on the Alaska Peninsula, part of the Alaska-Aleutian arc. Mount Veniaminof's September–December 2018 eruption began with small lava flows and culminated in a larger explosive event in November that produced a 450 km long continuous ash plume—one of the volcano's most significant ash-producing events in recently observed eruptions. A temporally constrained tephra record was sampled after the eruption. We present volatile element concentrations (H, C, S, Cl, F) from a total of 65 olivine-hosted melt inclusions (MIs) across six tephra layers to assess potential variations in magma chemistry throughout the course of the eruption. No variation is observed in volatile concentrations across the tephra layers, and there is no correlation between volatile concentrations and the explosivity of eruptive pulses leading to the November explosive event. H 2 O concentrations within MI suggest moderately hydrous magmas up to 2.4 wt% H 2 O and a large range in CO 2 concentrations in the MI (14–2,180 ppm) from Secondary Ion Mass Spectrometry and vapor bubble reconstructions. H 2 O and CO 2 measurements correspond to entrapment pressures up to 400 MPa (entrapment depths ∼15 km), but with modal storage around 100–200 MPa. H 2 O contents can be reconciled with previous microlite number density estimates using a polybaric storage history.

Alaska↗
USGS2026

Hydrothermal manganese-oxide mineralization of a carbonate ooze, Samoan hotspot region, South Pacific Ocean

Low-temperature hydrothermal manganese oxides occur throughout the global oceans. However, the hydrothermal replacement of a carbonate ooze by manganese oxides is described here for the first time. The 24 samples dredged from three locations around the Territory of American Samoa in the South Pacific Ocean include six samples with remnant carbonate and volcaniclastic sediments, which we refer to as “low Mn,” and 18 samples in which the mineralization is pervasive and has replaced most original sediment, termed “high Mn.” The 18 high-Mn samples exhibit a mean Mn content of 51 wt.%. Higher Li and Mn contents and lower Fe contents in the high-Mn samples indicate a hydrothermal origin and distinguish these samples from hydrogenetic ferromanganese crusts. Mn-oxide layers are up to 90 mm thick, some with columns to 44 mm long and 10 mm wide, the magnitude of which has not been described previously. Samples are composed of birnessite and 10 Å phyllomanganate minerals. Textures of the thickest Mn layers indicate mineralization below the seabed from ascending fluids during multiple hydrothermal pulses. Mineralization took place by complete to partial replacement and cementation of foraminiferal sediments intermixed with volcaniclastic sediments in varying amounts. Our results highlight the production of carbonate sediment-hosted hydrothermal Mn oxides from multiple hydrothermal sources within the Samoan volcanic chain. The potential extensive distribution on the regional scale of this newly described mineralization process and unique element enrichments raise questions about its broader distribution globally and potential importance to hydrothermal processes, element mass balance, and seabed mineral resources.

Geochemistry, Geophysics, Geosystems↗