Abstract Laboratory testing of digitizers and seismometers helps ensure that, before 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 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. Although these test 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.

The use of fiber-optic sensing systems in seismology has exploded in the past decade. Despite an ever-growing library of ground-breaking studies, questions remain about the potential of fiber-optic sensing technologies as tools for advancing if not revolutionizing earthquake-hazards-related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber-optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber-optic sensing to impact monitoring efforts by documenting the data completeness in a number of long-term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber-optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

Earth's long period background seismic wavefield is dominated by two distinct processes that couple ocean wave energy to a global microseism wavefield. We assess global microseism intensity in the secondary (4-10 s) and primary (14-20 s) bands, and across eight 2 s-wide period bands between 4 and 20 s. Robustly estimated primary and secondary secular amplitude trends are estimated at 73 globally distributed seismic station sites with continuous recording spanning at least 20 years, from as early as the late 1980s through October 2025. These trends are positive at significance for 61 (84%) and 46 (63%) stations with global average rates for vertical-component acceleration of 0.17 0.04 and 0.11 0.05%/yr, for the primary and secondary bands, respectively, with corresponding rates of energy increase of 0.27 0.08 and 0.15 0.09%/yr. Secular intensification is also observed within all 2 s period bands between 4 and 20 s. Amplitude histories for the longest primary microseism periods (18-20 s) correlate to near-antipodal distances, reflecting long-range teleconnections attributed to large-fetch storm systems, long-range swell and Rayleigh wave propagation, and geographically correlated El Ni & ntilde;o Southern Oscillation and other geographically extensive atmospheric influences on storms and waves. The lower average rates of intensification for the secondary microseism suggest that crossing wave systems in remote regions are either under-observed or are intensifying more slowly than the primary microseism, possibly due to increasing swell unidirectionality. Secular intensification is greatest at the longest primary microseism periods. This is consistent with a broadening of the global ocean wave spectrum by approximately 0.01%/yr which may reflect an increasing occurrence of large storm systems.

The 29 July 2025 Mw 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near-source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real-time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite-fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The Mw 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

The U.S. Geological Survey’s Geomagnetism Program is collaborating with the Earthquake Hazards Program and Global Seismographic Network Program to densify magnetic field observations. This collaboration focuses on the installation of magnetometers, or magnetic variometers, at existing seismic stations. Along with improving the density of space weather observations for hazard monitoring, these data can be used to correct colocated magnetic field induced noise in seismic data. Such corrections are especially useful during time periods of large magnetic storms where the magnetic field-induced instrument noise can be of similar amplitude to earthquake ground-motion records.

Seismology has been used as a tool for understanding the current physical properties of the interior of the Earth and its dynamic evolution with remarkable success over the last century. Much of this progress is due to the ever-expanding set of high-quality quantitative observations of teleseismic waveforms recorded at seismographic stations worldwide. In this work, we revisit historical seismological studies that helped first identify a core distinct in physical properties from the overlying mantle, followed by the detection of an inner core that was eventually verified to be solid based on normal-mode eigenperiods. Along with a brief overview of past studies of the Earth's inner core, we examine the reproducibility of these results and discuss how historical data compare against modern observations. After accounting for past normal-mode misidentifications, we confirm that introducing a solid inner core is required to afford significant improvements in fits to both radial modes and core-sensitive spheroidal overtones. Strong shear dissipation in the inner core of the radial reference Earth model, REM1D (Q mu= 89:54), fits the reference datasets of both normal-mode eigenperiods and quality factors accounting for physical dispersion. Because a liquid region would only have bulk dissipation, a narrow range of low Q mu values that are preferred by the reference datasets affords additional evidence of a solid inner core. In addition, we find that there is little systematic bias in the timing accuracy of historical data, although large variances exist. Investigations into the temperature, composition, and evolution of the inner core, as well as the reproducibility of past studies, can benefit from the reconciliation of historical and modern seismological datasets.

Abstract 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.

The finding of Yang et al. (2021; hereafter, YSR21) casts fundamental doubts on all the previous studies on the temporal changes of the inner core (IC) surface that are based on absolute arrival times, including Wen (2006; hereafter, WEN06) and its follow-up studies. The reports of small time shifts of the IC-reflected phase (PKiKP) likely misinterpreted station-related issues (clock errors and instrument changes) as evidence for rapid localized deformation at the IC surface. The comment by Zhang and Wen (2024; hereafter, ZW24) aims at defending the results and conclusions of WEN06 by questioning a tiny portion of the results (2 out of over 400 doublets and 3 out of over 400 problematic stations) in YSR21. In response, we acknowledge an accidental error in table S1 of YSR21, which is corrected in an erratum (Yang et al., 2023), although this does not affect the results or interpretations of YSR21. However, we strongly disagree with the criticisms in ZW24 due to many problems therein, especially the issue of irreproducibility due to its inconsistent procedures with YSR21. While the problems are related to the use of absolute arrivals in the previous studies by Wen’s group, key problems in this new comment are related to questionable data selection and inconsistent processing. We further caution the metadata updates, particularly when involving small temporal change signals. We focus on the main points in ZW24 and demonstrate that its arguments are problematic in many aspects and do not diminish the credibility of YSR21 or other related works by Yang and Song in support of the IC differential rotation.

Microelectromechanical system (MEMS) accelerometers are useful for seismological and engineering applications because of their ability to record unsaturated large seismic signals. Recent advances in MEMS technologies enable the design of instruments with improved capabilities that also allow the recording of small signals. As a result, MEMS can be useful across a broad dynamic range and for both major earthquakes and smaller magnitude events. Here, we analyze the performance of a MEMS-based accelerometer with the capability of real-time data transmission. We discuss the self-noise level, dynamic range, and sensitivity, along with the timing precision and data transmission latency, and compare these parameters to other MEMS-based instruments. We also summarize observations from a regional network deployed in southern Mexico over a period of 3+ yr for the purpose of earthquake early warning. In addition to the sensor evaluation, we present a large, openly available data set of strong motion data from the Mexican network that comprises continuous ground-motion records from 24 accelerometers in the period between 2017 and 2022.

The amplitude and frequency content of background seismic noise is highly variable with geographic location. Understanding the characteristics and behavior of background seismic noise as a function of location can inform approaches to improve network performance and in turn increase earthquake detection capabilities. Here, we calculate power spectral density estimates in one-hour windows for over 15 yr of vertical-component data from the nine-station Caribbean network (CU) and look at background noise within the 0.05 -300 s period range. We describe the most visually apparent features observed at the CU stations. One of the most prominent features occurs in the 0.75 -3 s band for which power levels are systematically elevated and decay as a function of proximity to the coastline. Further examination of this band on 1679 contiguous USArray Transportable Array stations reveals the same relationship. Such a relationship with coastal distance is not observed in the 4 -8 s range more typical of globally observed secondary microseisms. A simple surface-wave amplitude decay model fits the observed decay well with geometric spreading as the most important factor for stations near the coast ( <- 50 km). The model indicates that power levels are strongly influenced by proximity to coastline at 0.75 -3 s. This may be because power from nearshore wave action at 0.75 -3 s overwhelms more distant and spatially distributed secondary microseism generation. Application of this basic model indicates that a power reduction of- 25 dB can be achieved by simply installing the seismometer 25 km away from the coastline. This finding may help to inform future site locations and array design thereby improving network performance and data quality, and subsequently earthquake detection capabilities.

ABSTRACT Intermediate sized earthquakes (≈M4–6.5) are often measured using the teleseismic body-wave magnitude (mb). mb measurements are especially critical at the lower end of this range when teleseismic waveform modeling techniques (i.e., moment tensor analysis) are difficult. The U.S. Geological Survey National Earthquake Information Center (NEIC) determines the location and magnitude of all M 5 and greater earthquakes worldwide within 20 min of the rupture time, and therefore accurate mb magnitude estimates are essential to fulfill its mission. To better understand how network geometry and noise levels affect the global response capabilities, we developed a method to spatially estimate the minimum measurable mb. To do this, we compare expected mb amplitudes at every station to the station’s background noise level. We find that using NEIC’s current network geometry and these idealized thresholds, NEIC can potentially estimate mb magnitudes down to M 4.5 globally. Low-latitude regions in the Southern Hemisphere present the biggest opportunity to improve monitoring capabilities. However, logistically they also present the biggest hurdles for network operators. Finally, to test the resiliency of the network we removed the 20 most important stations and found the mb threshold remains mb 4.5. However, the region where only mb 4.5 and greater can be estimated increases and is again restricted to the Southern Hemisphere.

The U.S. Geological Survey (USGS) Global Seismographic Network (GSN) Program operates two thirds of the GSN, a network of state-of-the-art, digital seismological and geophysical sensors with digital telecommunications. This network serves as a multiuse scientific facility and a valuable resource for research, education, and monitoring. The other one third of the GSN is funded by the National Science Foundation (NSF), and the operations of this component are overseen by EarthScope. This collaboration between the USGS, EarthScope, and NSF has allowed for the development and operations of the GSN to be a truly multiuse network that provides near real-time open access data, facilitating fundamental discoveries by the Earth science community, supporting the earthquake hazards mission of the USGS, benefitting tsunami monitoring by the National Oceanic and Atmospheric Administration, and contributing to nuclear test monitoring and treaty verification. In this article, we describe the installation and evolution of the seismic networks operated by the USGS that ultimately led to the USGS portion of the GSN (100 stations under network codes IU, IC, and CU) as they are today and envision technological advances and opportunities to further improve the utility of the network in the future. This article focuses on the USGS-operated component of the GSN; a companion article on the GSN stations funded by the NSF and operated by the Cecil and Ida Green Institute of Geophysics and Planetary Physics, Scripps Institution of Oceanography, University of California at San Diego by Davis et al. (2023) appears in this volume.

The body-wave magnitude (mb) is a long-standing network-averaged, amplitude-based magnitude used to estimate the magnitude of seismic sources from teleseismic observations. The U.S. Geological Survey National Earthquake Information Center (NEIC) relies on mb in its global real-time earthquake monitoring mission. Although waveform modeling-based moment magnitudes are the modern standard to characterize earthquake size, mb is important because (1) in many cases, waveform modeling is not possible (e.g., low signal-to-noise events), (2) mb is applicable over a broad range of magnitudes, ∼M 4–7, and (3) there is a many decades-long history of estimating mb magnitudes. We use the NEIC Preliminary Determination of Epicenters earthquake catalog to investigate the uncertainty in NEIC station mb measurements. We show that mb measurements are spatially correlated, which can bias event mb, and we describe an empirical relation between this spatial correlation and station-to-station distance. We further describe an approach to mitigate bias from the spatial correlation. Accounting for the spatial covariance of observations can change the event mb from −0.15 to 0.07 mb units (10th to 90th percentile) for smaller events (mb≤4.5). These smaller events have the largest mb standard deviations ranging from 0.05 to 0.15 mb units (10th to 90th percentile).

Compilations of earthquake moment tensors from global and regional catalogs find pervasive non-double-couple (NDC) componentswith a mean deviation from a double-couple (DC) source of around 20%. Their distributions vary only slightly with magnitude, faulting mechanism, or geologic environments. This consistency suggests thatfor most earthquakes, especially smaller ones whose rupture processes are expected to be simpler, the NDC components are largely artifacts of the moment tensor inversion procedure. This possibility is also supported by the fact that NDC components for individual earthquakes with Mw<6.5 are only weakly correlated betweencatalogs. We explore this possibility by generating synthetic seismograms for the double-couple components of earthquakes around theworld using one Earth model and inverting them with a different Earth model. To match the waveforms with a different Earth model, the inversion changes the mechanisms to include a substantial NDC component while largely preserving the fault geometry (DC component). The resulting NDC components have a size and distribution similar to those reported for the earthquakes in the Global Centroid Moment Tensor (GCMT) catalog. The fact that numerical experiments replicate general features of the pervasive NDC components reported in moment tensor catalogs implies that these components are largely artifacts of the inversions not adequately accounting for the effects of laterally varying Earth structure.

Modern seismic data are collected, distributed, and analyzed using digital formats, and this has become a standard for the field. Although most modern seismometers still make use of analog electronic circuits, their data are converted from an analog voltage output to time-tagged counts by way of digitization. Although much of the digitization process is not complicated to conceptualize, there is a fair bit of jargon in digitizer specifications, and a few pitfalls that can arise in the processes of recording and analyzing ground-motion data. In this article, we review some of the fundamental physical properties of data acquisition systems and the basic steps in digitizing data from an analog instrument (specifically a seismometer). We then briefly discuss the digitization process and some of the key properties needed to make these data useful for seismological applications. Finally, we discuss some of the filtering processes that naturally arise from digitization and how it can affect the processing workflow. The end goal is to provide a user guide that will enable seismologists to have a working knowledge of the digitization process. We focus on aspects central to seismological applications and have tried to avoid getting bogged down in signal processing formalism.

Abstract Ocean waves excite continuous globally observable seismic signals. We use data from 52 globally distributed seismographs to analyze the vertical component primary microseism wavefield at 14–20 s period between the late 1980s and August 2022. This signal is principally composed of Rayleigh waves generated by ocean wave seafloor tractions at less than several hundred meters depth, and is thus a proxy for near-coastal swell activity. Here we show that increasing seismic amplitudes at 3σ significance occur at 41 (79%) and negative trends occur at 3σ significance at eight (15%) sites. The greatest absolute increase occurs for the Antarctic Peninsula with respective acceleration amplitude and energy trends ( ± 3σ) of 0.037 ± 0.008 nm s−2y−1 (0.36 ± 0.08% y−1) and 4.16 ± 1.07 nm2 s−2y−1 (0.58 ± 0.15% y−1), where percentage trends are relative to historical medians. The inferred global mean near-coastal ocean wave energy increase rate is 0.27 ± 0.03% y−1 for all data and is 0.35 ± 0.04% y−1 since 1 January 2000. Strongly correlated seismic amplitude station histories occur to beyond 50∘ of separation and show regional-to-global associations with El Niño and La Niña events.

ABSTRACT An increase in seismic stations also having microbarographs has led to increased interest in the field of seismoacoustics. A review of the recent advances in this field can be found in Dannemann Dugick et al. (2023). The goal of this note is to draw the attention of the readers of Dannemann Dugick et al. (2023) to several additional interactions between the solid Earth and atmosphere that have not been classically considered in the field of seismoacoustics. The 15 January 2022 Hunga Tonga–Hunga Ha‘api eruption produced acoustic gravity waves that were recorded globally. For example, the Lamb wave from this eruption produced early-arriving and long-lasting tsunami waves. This eruption also provided globally recorded coupling of atmospheric modes with solid Earth modes, providing another example of the complex interactions that can occur at the boundary between the atmosphere and the solid Earth. Even in the absence of large atmospheric signals, collocated pressure sensors at seismic stations can be a useful tool for estimating the local substructure, such at VS30, the average shear velocity of the upper 30 m. Finally, at low frequencies, it is possible to use pressure records to correct out atmospheric disturbances recorded on seismometers. We briefly review the aforementioned, nontraditional seismoacoustic topics that we feel are important to consider as part of the full suite of interactions occurring between the solid Earth and atmosphere.