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

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↗
USGS2026

Regional variations in sea ice and primary productivity in the Bering Sea during Marine Isotope Stage 11

Marine Isotope Stage (MIS) 11 (424-374 ka) has long been an analog for Holocene climate, because it is the most recent interglacial period with similar orbital conditions. However, there is significant global and regional variability in the climate response to MIS 11 warmth. Here, we review sediment core records from across the Bering Sea to investigate changes in paleoceanographic conditions during Marine Isotope Stages 12-10. Sea ice was present over much of the Bering Sea during MIS 11, but today, none of the sites investigated are ever ice-covered. This suggests that sea ice regimes in the Bering Sea during MIS 11 were different to those of the Holocene. There are also regional differences in the response of sea ice to MIS 11 warming. At the Umnak Plateau, Southeastern Bering Sea, sea ice concentrations decline during deglaciation, but they remain high at the slope sites until Peak MIS 11. Sea ice re-advances over the Umnak Plateau during peak interglacial warmth, at the same time that it declines over the slope sites. Late MIS 11 is characterized by high concentrations of seasonal sea ice at the Umnak Plateau, whilst sea ice at the slope sites fluctuates between consolidated and unconsolidated ice cover. This east-west dichotomy may be explained by changes in the behavior of the Aleutian Low. Productivity increases dramatically during deglaciation due to increased upwelling and sea level rise bringing fresh nutrients into the oceans. This is characterized by increased diatom productivity, increased terrestrial carbon deposition, and laminations at all sites.

Alaska↗
USGS2026

Constraining sources of mid-Pleistocene to present explosive volcanism in the Gulf of Alaska using machine learning and compositional data analysis

The long-term eruptive record of a region helps better elucidate magmatic processes at depth as well as volcanic hazards at the surface. Typically, reconstructing such records is done using proximal tephrostratigraphy and linking individual tephras to source volcanoes. These records, however, can only be accurately constructed if they are both well preserved and correctly linked to source volcanoes—not a trivial task, especially in areas such as Alaska that have experienced numerous glaciation events since the Pleistocene. This ultimately necessitates another way of assessing the long-term volcanic record such that these histories may be better discerned. Here we present data from 70 marine core tephras from the Gulf of Alaska, which have a virtually uninterrupted depositional record going back through the mid-Pleistocene. We utilize compositional data analysis techniques to quantify 37 eruptions over the span of eight cores, machine learning classification and conformal prediction algorithms to assign the most probable volcanic source(s) to each eruption, and multivariate distance-based metrics when machine learning classification algorithms are inappropriate. We find that the Mount Katmai magmatic system is the most probable volcanic source for analyzed tephras and that large volcanic centers such as Mount Katmai, Fisher Caldera, and Emmons Lake volcanic center have produced nearly invariant incompatible trace element ratio magmas, allowing for them to be confidently identified in long-term tephra records. This highlights the utility of trace elements for tephra studies, especially when paired with compositional data analysis, multivariate statistical tests, and petrologically informed discriminants.

Alaska↗