Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research- Solid Earth”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3Linked to original sources
USGS2018

Combining multiphase groundwater flow and slope stability models to assess stratovolcano flank collapse in the Cascade Range

Hydrothermal alteration can create low‐permeability zones, potentially resulting in elevated pore‐fluid pressures, within a volcanic edifice. Strength reduction by rock alteration and high pore‐fluid pressures have been suggested as a mechanism for edifice flank instability. Here we combine numerical models of multiphase heat transport and groundwater flow with a slope‐stability code that incorporates three‐dimensional distributions of strength and pore‐water pressure to address the following questions: (1) What permeability distributions and contrasts produce elevated pore‐fluid pressures in a stratovolcano? (2) What are the effects of these elevated pressures on flank stability? (3) Finally, what are the effects of magma intrusion on potential flank failure in an edifice? Simulation results show that under a range of plausible parameters, water tables in a stratovolcano can be elevated or perched. These elevated water tables result in universally lower stability (lower factor of safety) compared with equivalent dry edifices, indicating a higher likelihood of flank collapse. Low‐permeability (<1 × 10 −17 m 2 ) layers such as altered pyroclastic deposits or breccias can result in locally saturated regions (perched water) and lower factors of safety near the ground surface but may actually reduce liquid water saturation and pore pressures in the core of the edifice and thus may favor small, shallow collapses over larger, deeper collapses. Magma intrusion into the base of the edifice increases pore‐fluid pressures and decreases the factor of safety. However, the shear strength of edifice rocks also exerts a significant control on stability, so both mechanical properties and pore‐fluid pressures are important for stability assessments.

Washington↗
USGS2000

Seismovolcanic signals at Deception Island volcano, Antarctica: Wave field analysis and source modeling

The seismovolcanic signals associated with the volcanic activity of Deception Island (Antarctica), recorded during three Antarctic summers (1994–1995, 1995–1996 and 1996–1997), are analyzed using a dense small-aperture (500 m) seismic array. The visual and spectral classification of the seismic events shows the existence of long-period and hybrid isolated seismic events, and of low-frequency, quasi-monochromatic and spasmodic continuous tremors. All spectra have the highest amplitudes in the frequency band between 1 and 4 Hz, while hybrids and spasmodic tremors have also significant amplitudes in the high-frequency band (4–10 Hz). The array analysis indicates that almost all the well-correlated low-frequency signals share similar array parameters (slowness and back azimuth) and have the same source area, close to the array site. The polarization analysis shows that phases at high-frequency are mostly composed of P waves, and those phases dominated by low frequencies can be interpreted as surface waves. No clear shear waves are evidenced. From the energy evaluation, we have found that the reduced displacement values for surface and body waves are confined in a narrow interval. Volcano-tectonic seismicity is located close to the array, at a depth shallower than 1 km. The wave-field properties of the seismovolcanic signals allow us to assume a unique source model, a shallow resonating fluid-filled crack system at a depth of some hundreds of meters. All of the seismic activity is interpreted as the response of a reasonably stable stationary geothermal process. The differences observed in the back azimuth between low and high frequencies are a near-field effect. A few episodes of the degassification process in an open conduit were observed and modeled with a simple organ pipe.

Journal of Geophysical Research: Solid Earth↗
USGS2005

The U-tube: A novel system for acquiring borehole fluid samples from a deep geologic CO2 sequestration experiment

A novel system has been deployed to obtain geochemical samples of water and gas, at in situ pressure, during a geologic CO 2 sequestration experiment conducted in the Frio brine aquifer in Liberty County, Texas. Project goals required high-frequency recovery of representative and uncontaminated aliquots of a rapidly changing two-phase fluid (supercritical CO 2 and brine) fluid from 1.5 km depth. The data sets collected, using both the liquid and gas portions of the downhole samples, provide insights into the coupled hydrogeochemical issues affecting CO 2 sequestration in brine-filled formations. While the basic premise underlying the U-tube sampler is not new, the system is unique because careful consideration was given to the processing of the recovered two-phase fluids. In particular, strain gauges mounted beneath the high-pressure surface sample cylinders measured the ratio of recovered brine to supercritical CO 2 . A quadrupole mass spectrometer provided real-time gas analysis for perfluorocarbon and noble gas tracers that were injected along with the CO 2 . The U-tube successfully acquired frequent samples, facilitating accurate delineation of the arrival of the CO 2 plume, and on-site analysis revealed rapid changes in geochemical conditions.

Journal of Geophysical Research Solid Earth↗
USGS2008

Effects of supershear rupture speed on the high-frequency content of S waves investigated using spontaneous dynamic rupture models and isochrone theory

In this paper we achieve three goals: (1) We demonstrate that crack tips governed by friction laws, including slip weakening, rate- and state-dependent laws, and thermal pressurization of pore fluids, propagating at supershear speed have slip velocity functions with reduced high-frequency content compared to crack tips traveling at subshear speeds. This is demonstrated using a fully dynamic, spontaneous, three-dimensional earthquake model, in which we calculate fault slip velocity at nine points (locations) distributed along a quarter circle on the fault where the rupture is traveling at supershear speed in the in-plane direction and subshear speed in the antiplane direction. This holds for a fault governed by the linear slip-weakening constitutive equation, by slip weakening with thermal pressurization of pore fluid, and by rate- and state-dependent laws with thermal pressurization. The same is also true even assuming a highly heterogeneous initial shear stress field on the fault. (2) Using isochrone theory, we derive a general expression for the spectral characteristics and geometric spreading of two pulses arising from supershear rupture, the well-known Mach wave, and a second lesser known pulse caused by rupture acceleration. (3) We demonstrate that the Mach cone amplification of high frequencies overwhelms the de-amplification of high-frequency content in the slip velocity functions in supershear ruptures. Consequently, when earthquake ruptures travel at supershear speed, a net enhancement of high-frequency radiation is expected, and the alleged “low” peak accelerations observed for the 2002 Denali and other large earthquakes are probably not caused by diminished high-frequency content in the slip velocity function, as has been speculated.

Journal of Geophysical Research Solid Earth↗
USGS2008

Noise in GPS displacement measurements from southern California and southern Nevada

Time series of position changes estimated from data from 236 continuously recording GPS receivers operating in Southern California and Southern Nevada are evaluated for noise models that characterize their temporal correlations. The lengths of the time series range between 3.5 and 10 years. After adjusting these data for postseismic deformation, offsets, and annual periodicities, I find that about one-half of the time series have temporal correlations that are categorized as either flicker or random-walk noise. The remaining time series can be best categorized as either a combination of flicker and random-walk; power law noise; first-order Gauss-Markov plus random-walk noise; or power law plus broadband, seasonal noise. A variety of geodetic monuments are used in Southern California and Nevada, including deeply braced designs, cement piers, pins drilled in outcrop, and buildings. When I evaluate the noise for each time series in terms of an estimate of the standard error in velocity, I find that the sites with the smallest errors are those located in Nevada using deeply braced monuments. Sites that are installed within regions of active pumping, both for groundwater and oil, had the largest standard errors in velocity. Comparison of monument stability, as measured by standard error in rate, with average, annual rainfall nearby indicates a marginally significant correlation. In addition, even though regional filtering removed much of the common-mode signals in these time series, there still remains a common-mode seasonal signal which can and should be removed.

California, Nevada↗
USGS2008

A global search for stress shadows

[1] Debate continues regarding the relative proportion of earthquakes triggered by passing seismic waves versus static stress changes from a main shock. Static stress changes are expected to have long‐term effects on earthquake probabilities, whereas dynamic stress changes due to the passing of seismic waves should not. Both mechanisms are expected to raise seismicity rates in some areas, but only static stress change calculations predict rate decrease shadows. Thus, identification of post‐main‐shock earthquake suppression is diagnostic of a static stress change process. We note that in principle, static stress change theory predicts suppression of particular earthquake mechanisms in a shadow zone rather than an overall rate reduction. A stress shadow can therefore be characterized by a change in the average earthquake focal mechanism before and after a main shock that results from suppression of a given mechanism type. We examined average mechanisms from ±2° radii and 5‐year periods before and after 119 M s ≥ 7 main shock earthquakes drawn from the Harvard Centroid Moment Tensor (CMT) catalog. Significant average mechanism changes caused by earthquake suppression were found in only two cases. However, by stacking the data, we were able to resolve statistically significant suppression of particular post‐main‐shock focal mechanisms. This indicates that, while static stress shadows are subtle, they are indeed present in the global catalog.

Journal of Geophysical Research Solid Earth↗
USGS2008

Exponential decline of aftershocks of the M7.9 1868 great Kau earthquake, Hawaii, through the 20th century

The remarkable catalog of Hawaiian earthquakes going back to the 1820s is based on missionary diaries, newspaper accounts, and instrumental records and spans the great M 7.9 Kau earthquake of April 1868 and its aftershock sequence. The earthquake record since 1868 defines a smooth curve complete to M 5.2 of the declining rate into the 21st century, after five short volcanic swarms are removed. A single aftershock curve fits the earthquake record, even with numerous M 6 and 7 main shocks and eruptions. The timing of some moderate earthquakes may be controlled by magmatic stresses, but their overall long‐term rate reflects one of aftershocks of the Kau earthquake. The 1868 earthquake is, therefore, the largest and most controlling stress event in the 19th and 20th centuries. We fit both the modified Omori (power law) and stretched exponential (SE) functions to the earthquakes. We found that the modified Omori law is a good fit to the M ≥ 5.2 earthquake rate for the first 10 years or so and the more rapidly declining SE function fits better thereafter, as supported by three statistical tests. The switch to exponential decay suggests that a possible change in aftershock physics may occur from rate and state fault friction, with no change in the stress rate, to viscoelastic stress relaxation. The 61‐year exponential decay constant is at the upper end of the range of geodetic relaxation times seen after other global earthquakes. Modeling deformation in Hawaii is beyond the scope of this paper, but a simple interpretation of the decay suggests an effective viscosity of 10 19 to 10 20 Pa s pertains in the volcanic spreading of Hawaii's flanks. The rapid decline in earthquake rate poses questions for seismic hazard estimates in an area that is cited as one of the most hazardous in the United States.

Hawaii↗
USGS2009

Three‐dimensional model of Hellenic Arc deformation and origin of the Cretan uplift

[1] The Hellenic Arc of Greece is the most seismically active part of Europe, but little is know about its mechanics. We modeled deformation along the arc using a finite element model. The model was intended to capture large‐scale 3‐D structure of Nubian plate subduction beneath the Aegean block and its deformational consequences. The shape of the interface was developed using mapped traces at the surface and earthquake hypocenters at depth. Model block motions were constrained by recent compilations of GPS velocity vectors. We simulated a 10 ka period of convergence between Nubia and the Aegean and calculated the strain field in the overriding plate as well as the spatial distribution and orientation of differential stress (∣ σ 1 − σ 3 ∣). From these calculations we derived testable quantities such as the expected seismic moment rate on the interplate contact, uplift pattern, and distribution of strain modes. Our relatively simple model broadly reproduced observed uplift patterns, earthquake activity, and loci of extension and contraction. The model showed a localization of uplift near the island of Crete, where the fastest Aegean uplift rates are well documented. Comparison of calculated expected seismic moment and observed earthquake catalogs implies a nearly fully coupled interplate contact. On the basis of our modeling results, we suggest that south Aegean deformation is driven primarily by the fast moving (∼33 mm a −1 ) Aegean upper plate overriding a nearly stalled (∼5 mm a −1 ) Nubian lower plate. This tectonic setting thus more closely resembles a continental thrust than it does a typical oceanic subduction zone.

Aegean Sea, Hellenic Arc↗
USGS2009

Constraints on the stress state of the San Andreas fault with analysis based on core and cuttings from SAFOD drilling phases 1 and 2

Analysis of field data has led different investigators to conclude that the San Andreas Fault (SAF) has either anomalously low frictional sliding strength ( μ < 0.2) or strength consistent with standard laboratory tests ( μ > 0.6). Arguments for the apparent weakness of the SAF generally hinge on conceptual models involving intrinsically weak gouge or elevated pore pressure within the fault zone. Some models assert that weak gouge and/or high pore pressure exist under static conditions while others consider strength loss or fluid pressure increase due to rapid coseismic fault slip. The present paper is composed of three parts. First, we develop generalized equations, based on and consistent with the Rice (1992) fault zone model to relate stress orientation and magnitude to depth‐dependent coefficient of friction and pore pressure. Second, we present temperature‐ and pressure‐dependent friction measurements from wet illite‐rich fault gouge extracted from San Andreas Fault Observatory at Depth (SAFOD) phase 1 core samples and from weak minerals associated with the San Andreas Fault. Third, we reevaluate the state of stress on the San Andreas Fault in light of new constraints imposed by SAFOD borehole data. Pure talc ( μ ≈0.1) had the lowest strength considered and was sufficiently weak to satisfy weak fault heat flow and stress orientation constraints with hydrostatic pore pressure. Other fault gouges showed a systematic increase in strength with increasing temperature and pressure. In this case, heat flow and stress orientation constraints would require elevated pore pressure and, in some cases, fault zone pore pressure in excess of vertical stress.

California↗
USGS2010

Geodetic evidence for en echelon dike emplacement and concurrent slow slip during the June 2007 intrusion and eruption at Kilauea volcano, Hawaii

A series of complex events at Kīlauea Volcano, Hawaii, 17 June to 19 June 2007, began with an intrusion in the upper east rift zone (ERZ) and culminated with a small eruption (1500 m 3 ). Surface deformation due to the intrusion was recorded in unprecedented detail by Global Positioning System (GPS) and tilt networks as well as interferometric synthetic aperture radar (InSAR) data acquired by the ENVISAT and ALOS satellites. A joint nonlinear inversion of GPS, tilt, and InSAR data yields a deflationary source beneath the summit caldera and an ENE-striking uniform-opening dislocation with ~2 m opening, a dip of ∼80° to the south, and extending from the surface to ~2 km depth. This simple model reasonably fits the overall pattern of deformation but significantly misfits data near the western end of an inferred dike-like source. Three more complex dike models are tested that allow for distributed opening including (1) a dike that follows the surface trace of the active rift zone, (2) a dike that follows the symmetry axis of InSAR deformation, and (3) two en echelon dike segments beneath mapped surface cracks and newly formed steaming areas. The en echelon dike model best fits near-field GPS and tilt data. Maximum opening of 2.4 m occurred on the eastern segment beneath the eruptive vent. Although this model represents the best fit to the ERZ data, it still fails to explain data from a coastal tiltmeter and GPS sites on Kīlauea's southwestern flank. The southwest flank GPS sites and the coastal tiltmeter exhibit deformation consistent with observations of previous slow slip events beneath Kīlauea's south flank, but inconsistent with observations of previous intrusions. Slow slip events at Kīlauea and elsewhere are thought to occur in a transition zone between locked and stably sliding zones of a fault. An inversion including slip on a basal decollement improves fit to these data and suggests a maximum of ~15 cm of seaward fault motion, comparable to previous slow-slip events.

Hawai'i↗
USGS2010

The 2007 M5.4 Alum Rock, California, earthquake: Implications for future earthquakes on the central and southern Calaveras Fault

[1] The similarity of seismograms recorded by two seismic stations demonstrate that the 31 October 2007 moment magnitude M 5.4 Alum Rock earthquake is a repeat of a 1955 M L 5.5 earthquake. Both occurred on Oppenheimer et al.'s (1990) Zone V “stuck patch” on the central Calaveras fault, providing new support for their model of Calaveras fault earthquake activity. We suggest that Zone V fails only in a family of recurring M ∼ 5.4–5.5 earthquakes. The 1955 and 2007 earthquakes are the penultimate and ultimate Zone V events. Earthquakes in 1891 and 1864 are possible earlier Zone V events. The next Zone V event is not expected in the next few decades, assuming a time‐dependent recurrence model: the mean forecast date is 2064 (2035–2104, 95% confidence range). We further suggest that Zones I, II, III, and IV fail in recurring M ∼ 5.1–5.3, M ∼ 5.6–5.8, M ∼ 6.1–6.3, and M ∼ 4.9–5.0 earthquakes, respectively. If our earthquake recurrence model is correct, the next Zone I event is overdue and could occur anytime, and M 5–6 earthquakes should not occur on Zones II, III, and IV before 2014, 2012, and 2026, respectively. We cannot rule out the possibility that Zone VI, which lies at the southern end of the Mission Seismic Trend, where the southern Hayward and central Calaveras faults appear to connect at depth, fails aseismically or in large events on the southern Hayward fault, such as last occurred in 1868, or in large events on the adjoining northern Calaveras fault segment.

California↗
USGS2010

Postseismic relaxation following the 1989 MS7.1 Loma Prieta earthquake, central California

The postseismic relaxation (postseismic displacement less displacement that would have occurred at the preseismic rate) measured by GPS and leveling following the 1989 M S 7.1 Loma Prieta earthquake is reexamined. The temporal dependence of the relaxation over the first 1200 days postseismic is well described by 1 − e − t / τ , where τ = 414 ± 92 days. (Quoted uncertainties are standard deviations.) That temporal dependence appears to be a linear function of the cumulative number of M > 2.5 aftershocks that have occurred. The relaxation is attributed to afterslip (1.56 ± 0.20 m dextral strike slip and 0.60 ± 0.04 m reverse slip) on the downdip extension (depth, 16–21 km) of the coseismic rupture plus a collapse (0.11 ± 0.02 m fault-normal displacement) of the rupture zone (depth, 5–16 km). Because the postseismic uplift was determined by leveling over a route with relief in excess of 1000 m, an allowance (18 ± 2 ppm of height above the base elevation) for excess (beyond corrections already applied) height-dependent error in the measured uplift was estimated simultaneously with the afterslip and collapse parameters. This new solution for afterslip and collapse on the plane of the rupture provides an alternative explanation to the suggestion by Bürgmann et al. (1997) that the postseismic deformation was due to afterslip on the coseismic rupture and on a nearby, shallow, 30° dipping thrust.

California↗
USGS2019

The 12 November 2017 Mw 7.3 Ezgeleh–Sarpolzahab (Iran) earthquake and active tectonics of the Lurestan arc

The 12 November 2017 M w 7.3 Ezgeleh‐Sarpolzahab earthquake is the largest instrumentally recorded earthquake in the Zagros Simply Folded Belt by a factor of ∼10 in seismic moment. Exploiting local, regional, and teleseismic data and synthetic aperture radar interferometry imagery, we characterize the rupture, its aftershock sequence, background seismicity, and regional tectonics. The mainshock ruptured slowly (∼2 km/s), unilaterally southward, for ∼40 km along an oblique (dextral‐thrust) fault that dips ∼14°E beneath the northwestern Lurestan arc. Slip is confined to basement depths of ∼12–18 km, resolvably beneath the sedimentary cover which is ∼8 km thick in this area. The gentle dip angle and basement location allow for a broad slip area, explaining the large magnitude relative to earthquakes in the main Fars arc of the Zagros, where shallower, steeper faults are limited in rupture extent by weak sedimentary layers. Early aftershocks concentrate around the southern and western edges of the mainshock slip area and therefore cluster in the direction of rupture propagation, implying a contribution from dynamic triggering. A cluster of events ∼100 km to the south near Mandali (Iraq) reactivated the ∼50° dipping Zagros Foredeep Fault. The basement fault responsible for the Ezgeleh‐Sarpolzahab earthquake probably accounts for the ∼1 km elevation contrast between the Lurestan arc and the Kirkuk embayment but is distinct from sections of the Mountain Front Fault that define frontal escarpments elsewhere in the Zagros. It may be related to a seismic interface underlying the central and southern Lurestan arc, and a key concern is whether or not the more extensive regional structure is also seismogenic.

Journal of Geophysical Research: Solid Earth↗
USGS2019

Surface imaging functions for elastic reverse time migration

Reverse time migration is often used to interpret acoustic or three‐component seismic recordings by creating an image of subsurface seismic reflectors. Here I describe elastic reverse time migration imaging functions that are cast as waveform misfit sensitivity kernels of contrasts in material parameters across hypothetical seismic discontinuities, that is, specular reflectors. The proposed “surface” imaging functions are theoretically applicable to either reflected or converted waves in order to estimate the location and reflectivity of these discontinuities. The surface imaging functions, as well as volumetric sensitivity kernels that target point diffractors, are tested on sets of synthetic surface array recordings that sample the 3‐D seismic wavefield on simple 2‐D structures generated using a 2.5‐D spectral element method. These tests illustrate that in contrast with the volumetric sensitivity kernels, the reflectivity is generally dominated by a high‐amplitude peak that coincides with input locations of discontinuities. Passive recordings of microseismicity, shot gathers, or a combination thereof can be potentially interpreted with the new surface imaging functions to yield useful reflectivity images.

Journal of Geophysical Research Solid Earth↗
USGS2019

Insights into pāhoehoe lava emplacement using visible and thermal structure-from-motion photogrammetry

We present the evolution over 3 months of a 2016–2017 pāhoehoe flow at Kīlauea as it changed from a narrow sheet flow into a compound lava field fed by a stable system of tubes. The portion of the flow located on Kīlauea's coastal plain was characterized using helicopter‐based visible and thermal structure‐from‐motion photogrammetry to construct a series of georeferenced digital surface models and thermal maps on eight different days. Results reveal key influences on the emplacement and evolution of such long‐lived pāhoehoe flows. This region of the flow grew by ~12 × 10 6 m 3 with a near‐constant time‐average discharge rate of 1.2–2.7 m 3 /s. The development of two tube systems is captured and shows an initial nascent tube enhanced by a narrow topographic confinement, which later inflated and created a topographic inversion that modulated the emplacement of a second flow lobe with its own tube system. The analysis of breakouts at various stages of the field's life suggests that the evolution of the thermal and morphological properties of the flow surface reflect its maturity. Thermal properties of breakouts were used to expand the empirical relationship of breakout cooling to longer timescales. This study contributes to the long‐term development and validation of more accurate predictive models for pāhoehoe, required during the management of long‐lasting lava flow crises in Hawai'i and elsewhere.

Hawaii↗
USGS2020

Submarine landslide kinematics derived from high-resolution imaging in Port Valdez, Alaska

Submarine landslides caused by strong ground shaking during the M9.2 1964 Great Alaska earthquake generated a tsunami that destroyed much of the old town of Valdez, Alaska, and was responsible for 32 deaths at that location. We explore structural details of the 1964 landslide deposit, as well as landslide deposits from earlier events, in order to characterize kinematics of the landslide process. We present a new high‐resolution seismic reflection data set that images the 1964 landslide deposit and six pre‐1964 deposits with great detail. These deposits are represented by thick packages (~7–23 m) of debris within >500 m of fjord sedimentation above basement. Internal slide structures are associated with distinctive landslide failure mechanisms, including detailed erosional and depositional features and structures resolved within both landslide blocks and distal debris flow layers. Based on comparisons of deposit volume from subbottom structure and differenced bathymetry, we refine prior interpretations of the source of failed material. New data show evidence for basal erosion and reworking of fjord‐floor sedimentation. Additionally, material comprising the 1964 landslide appears to have been translated and deformed by lateral thrusting, rather than having been sourced entirely from upslope evacuation zones. Taking into account these complexities in depositional patterns, we show variations in slide size through Holocene time and relate the history of landslides to the paleoseismic record. Collectively, these new observations demonstrate that Port Valdez has a repeated history of large submarine landslides, which are likely associated with large megathrust earthquakes.

Alaska↗
USGS2020

Magma intrusion and volatile ascent beneath Norris Geyser Basin, Yellowstone National Park

Recent activity has provided new insights into the causes of surface deformation in and around the Yellowstone Caldera, a topic that has been debated since the discovery of caldera floor uplift more than four decades ago. An episode of unusually rapid uplift (>15 cm/yr) centered near Norris Geyser Basin along the north caldera rim began in late 2013 and continued until a M w 4.9 earthquake on 30 March 2014; thereafter, uplift abruptly switched to subsidence. Uplift at rates of several centimeters per year resumed in 2016 and continued at least through the end of 2018. Modeling of Global Positioning System and interferometric synthetic aperture radar data suggests an evolving process of deep magma intrusion during 1996–2001 followed by volatile ascent and accumulation at shallow levels, perhaps as shallow as a few hundred meters depth. The depth of shallow volatile accumulation appears to have shallowed from the 2014 to the 2016 deformation episode, and frequent eruptions of Steamboat Geyser since March 2018 are likely a surface manifestation of this ongoing process. Hydrothermal explosion features are prominent in the Norris Geyser Basin area, and the apparent shallow nature of the volatile accumulation implies an increased risk of hydrothermal explosions.

Wyoming↗
USGS2019

Separating sea and slow slip signals on the seafloor

Seafloor pressure measurements hold promise for estimating vertical displacements from transient slow slip events on submarine faults. We assess the accuracy of pressure offset estimates that evolve over days to weeks and the confidence with which they may be attributed to tectonic deformation or to the ocean water column. One common approach to resolve this ambiguity assumes water column pressures vary insignificantly over the study region and are represented by stable reference site pressures. Assessing the validity of this assumption requires independent evidence. Correlations between pressures and colocated temperatures collected during the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment suggest temperatures might provide a useful independent proxy for water column pressures. We compared offsets estimated using several methods, with temperature and other proxies. The use of a temperature proxy was unsuccessful, because seafloor temperatures did not track the seasonal signal that contributes significantly to seafloor pressure changes over the slow slip event period. Regardless of the estimation method, offsets varied within a few cm around some uncertain reference level. Commonly used statistical measures are shown not to be reliable indicators of offset accuracy since offsets contribute minimally to the total variance. Offsets estimated using identical methods but with seafloor pressures simulated using a regional ocean model were larger than those derived from the data but had a similar pattern. Since the model simulates only water column processes, this suggests a significant fraction of the estimated pressure offsets are due to seasonal water column signal and are not of tectonic origin.

Journal of Geophysical Research- Solid Earth↗