Publications

For imaging of seismic discontinuities at depth, reverse time migration (RTM) is a powerful method to apply to recordings of seismic events. It is especially powerful when an extensive receiver array, numerous seismic sources, or both, permit adequate reconstruction of incident and scattered wavefields at depth. Reconstructing either the incident or scattered wavefield at depth becomes...

Kīlauea volcano hosts numerous pit craters that are inferred to have formed in competent bedrock (lava flows with minor tephra and other sediments), including Wood Valley Pit Crater. The Wood Valley Pit Crater is a 50-meter-deep, nearly circular pit that includes access to a cave entrance, which provides an opportunity to monitor cave climate throughout a cave that is ordinarily...

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between...

The purpose of this project was to monitor the cave climate of the Sunset Crater National Monument “Bonito Flow” Ice Cave. The main purpose of the climate monitoring was to determine if “Ice Cave” was still an ice cave, i.e., a cave that contains perennial ice. The data acquired from 2 March 2021 to 13 Dec. 2022 consisted of temperature, humidity, and pressure throughout the cave...

On Earth, caves are unique environments at the intersection of geology, climate, and biology. Given that the same terrestrial speleogenetic processes exist throughout the solar system, it would be surprising if caves beyond Earth did not exist. Thousands of potential cave entrances (or subsurface access points) have been identified from Earth’s Moon to Pluto’s moon, Charon. To date, our...

Fitting parametric seismological models to earthquake catalogs often comes with numerical challenges, especially when catalogs are small. An alternative way to quantify parameter values for a seismic region is by eliciting expert opinions on the seismological characteristics that each parameter corresponds to. For instance, expert beliefs on aftershock patterns can be formulated into...

When sunspots are large and numerous, intense magnetic storms are likely to occur on the Earth. Magnetic storms can generate electric fields in the Earth, and these fields can, in turn, interfere with electric power transmission grids that are grounded at the Earth’s surface. Across the contiguous United States, geoelectric hazards are highest in the Upper Midwest and in the East. These...

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

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R(t)=K(c+t)−p⁠, with parameters K, c, and p. The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced...

On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical...

The square-root-impedance (SRI) method is commonly used to approximate the seismic site amplifications computed using the full-resonance (FR) method for gradient shear-wave velocity (VS) profiles that are smoothly varying with depth. The SRI site amplifications have been observed to systematically underpredict the FR site amplifications by a ratio of FR/SRI amplifications around 1.05 to...

This paper introduces cross-fade sampling, a computationally efficient Markov Chain Monte Carlo simulation method that uses a semi-analytical approach to quickly solve Bayesian inverse problems that do not themselves have an analytical solution. Cross-fading is efficient in two ways. First, it requires fewer samples to obtain the same quality simulation of the target probability density...