Search USGS⌕ Search

USGS · 70277088

The digital archivist: Automating legacy macroseismic data processing using large language models

Abstract

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Aarnav Agrawal, Susan E. Hough, Mostafa Mousavi, Margaret Hellweg, William Ellsworth, Clara Yoon, Salvador Blanco. 2026-07-01. The digital archivist: Automating legacy macroseismic data processing using large language models. https://doi.org/10.1785/0220250362

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

Discover connections

Connections use source metadata and explicit phrase matches, not verified experimental comparisons.

KEEP EXPLORING

Related USGS reports

Amploc: A complete Python workflow for amplitude-based locations

Amplitude‐based locations are well suited for characterizing long‐duration seismic sources that may migrate with time, such as debris flows, avalanches, rockfall, and volcanic tremor. Although amplitude‐based location methods have existed for decades, their adoption in real‐time monitoring systems remains limited. We present a complete Python‐based workflow that includes station corrections, amplitude‐based location algorithms, and interactive review tools that are appropriate for real‐time use. The algorithms themselves are not new; however, Amploc integrates several previously published methods into a unified operational framework. The workflow also introduces a browser‐based graphical interface that allows users to review real‐time location results using formats familiar to many observatories. By integrating calibration, location algorithms, and visualization tools, Amploc lowers the barrier to operational use of amplitude‐based locations and improves characterization of phenomena that inform hazard assessments such as long‐period earthquakes, volcanic tremor, and debris flows.

British Columbia, Oregon, Washington↗

Global seasonal microseism intensity and hemispheric asymmetry

Earth’s microseism wavefield dominates seismic background levels at periods between approximately 4 and 20 s, and reflects periodic and secular variations in ocean swell energy. Ocean wave energy couples to the seismic wavefield via distinct primary microseism and secondary microseism (PM and SM) source mechanisms, which are excited by basal swell tractions and seafloor pressure variations due to crossing seas, respectively. This study examines annual amplitude variations for the globally dominant PM (14–20 s) and SM (4–10 s) period bands. Annual harmonic variations are represented by four‐term Fourier series fits to vertical‐component seasonally smoothed acceleration time series from 73 stations in the Global Seismographic and GEOSCOPE networks with over 20 yr of recording and at least 75% data completeness. These annual periodic functions fit between 14%–95% (PM) and 22%–97% (SM) of signal variance. Station annual peak‐to‐peak variations range between 1.2–14.3 dB (PM) and 1.5–20.6 dB (SM). An asymmetry in microseism features exists between the Northern (NH) and Southern (SH) Hemispheres. High‐latitude NH stations show highly correlated PM and SM annual amplitude variations. This character dominates the wider extratropical NH but diminishes at tropical latitudes, and widespread relative PM–SM decorrelation is observed in the SH. These hemispheric characteristics reflect systematic differences in both extratropical storm activity and ocean wave state. Greater annual variation and seasonal predictability in the NH reflect the influence of the large continental landmasses that enhance both storm intensity seasonality and SM‐generating coastal wave reflection. Notable clusters of low PM–SM correlation stations are also observed in continental Antarctica due to seasonal sea ice influences, and in East Asia, reflecting unusual PM and SM source responses to South Asian monsoonal and tropical cyclone ocean wave influences.

Seismological Research Letters↗

Digitizer Suite: The Albuquerque Seismological Laboratory Digitizer Testing Suite

Laboratory testing of digitizers and seismometers helps ensure that prior to deployment the instrumentation can produce high quality data and is operating within specifications. In this work we detail the software package called: the Albuquerque Seismological Laboratory (ASL) Digitizer Test Suite. This Java software package provides several algorithms to verify various performance parameters of digitizers commonly used for recording analog seismic instruments. The goal of these tests is not to be exhaustive, but to identify common failures that could compromise the integrity of seismic data being recorded on the digitizer. For example, Sandia National Laboratories (e.g., Slad and Merchant, 2018) routinely do comprehensive testing of digitizers for various monitoring missions. While these tests reports are valuable for comprehensively characterizing a recording system, it would be resource intensive to conduct such tests on every seismic recorder used in a network. We focus on tests that include ways to estimate the sensitivity, timing, self-noise, and clip-level of the digitizer, as well as the fidelity of the signal being recorded. The software is publicly available and provides a way for the community to verify the integrity of a digitizer using a minimum amount of outside equipment.

Seismological Research Letters↗