Search USGS⌕ Search

Geology topics

K.M. Larson

Publications and source records attributed to K.M. Larson.

7 recordsLinked to original sources

Spatiotemporal evolution of dike opening and décollement slip at Kīlauea Volcano, Hawai'i

Rapid changes in ground tilt and GPS positions on Kīlauea Volcano, Hawai'i, are interpreted as resulting from a shallow, two‐segment dike intrusion into the east rift zone that began at 1217 UTC (0217 HST) on 17 June 2007 and lasted almost 3 days. As a result of the intrusion, a very small volume of basalt (about 1500 m 3 ) erupted on 19 June. Northward tilt at a coastal tiltmeter, subsidence of south flank GPS sites, southeastward displacements at southwestern flank GPS sites, and a swarm of flank earthquakes suggest that a slow slip event occurred on the décollement beneath Kīlauea's south flank concurrent with the rift intrusion. We use 4 min GPS positions that include estimates of time‐dependent tropospheric gradients and ground tilt data to study the spatial and temporal relationships between the two inferred shallow, steeply dipping dike segments extending from the surface to about 2 km depth and décollement slip at 8 km depth. We invert for the temporal evolution of distributed dike opening and décollement slip in independent inversions at each time step using a nonnegative least squares algorithm. On the basis of these inversions, the intrusion occurred in two stages that correspond spatially and temporally with concentrated rift zone seismicity. The dike opening began on the western of the two segments before jumping to the eastern segment, where the majority of opening accumulated. Dike opening preceded the start of décollement slip at an 84% confidence level; the latter is indicated by the onset of northward tilt of a coastal tiltmeter. Displacements at southwest flank GPS sites began about 18 h later and are interpreted as resulting from slow slip on the southwestern flank. Additional constraints on the evolution of the intrusion and décollement slip come from inversion of an Envisat interferogram that spans the intrusion until 0822 UTC on 18 June 2007, combined with GPS and tilt data. This inversion shows that up to 0822 UTC on 18 June, décollement slip is only required in a limited region offshore of Ka'ena Point. A similar inversion of the complete event, which includes GPS and tilt data up to 21 June and a second Envisat interferogram spanning the complete intrusion until 21 June, shows décollement slip spread westward across the south flank. This may suggest westward migration of the décollement slip as the event progressed.

Hawaii↗

Stability and uncertainty of finite-fault slip inversions: Application to the 2004 Parkfield, California, earthquake

The 2004 Parkfield, California, earthquake is used to investigate stability and uncertainty aspects of the finite-fault slip inversion problem with different a priori model assumptions. We utilize records from 54 strong ground motion stations and 13 continuous, 1-Hz sampled, geodetic instruments. Two inversion procedures are compared: a linear least-squares subfault-based methodology and a nonlinear global search algorithm. These two methods encompass a wide range of the different approaches that have been used to solve the finite-fault slip inversion problem. For the Parkfield earthquake and the inversion of velocity or displacement waveforms, near-surface related site response (top 100 m, frequencies above 1 Hz) is shown to not significantly affect the solution. Results are also insensitive to selection of slip rate functions with similar duration and to subfault size if proper stabilizing constraints are used. The linear and nonlinear formulations yield consistent results when the same limitations in model parameters are in place and the same inversion norm is used. However, the solution is sensitive to the choice of inversion norm, the bounds on model parameters, such as rake and rupture velocity, and the size of the model fault plane. The geodetic data set for Parkfield gives a slip distribution different from that of the strong-motion data, which may be due to the spatial limitation of the geodetic stations and the bandlimited nature of the strong-motion data. Cross validation and the bootstrap method are used to set limits on the upper bound for rupture velocity and to derive mean slip models and standard deviations in model parameters. This analysis shows that slip on the northwestern half of the Parkfield rupture plane from the inversion of strong-motion data is model dependent and has a greater uncertainty than slip near the hypocenter.

California↗

Slip history of the 2003 San Simeon earthquake constrained by combining 1-Hz GPS, strong motion, and teleseismic data

The slip history of the 2003 San Simeon earthquake is constrained by combining strong motion and teleseismic data, along with GPS static offsets and 1-Hz GPS observations. Comparisons of a 1-Hz GPS time series and a co-located strong motion data are in very good agreement, demonstrating a new application of GPS. The inversion results for this event indicate that the rupture initiated at a depth of 8.5 km and propagated southeastwards with a speed ???3.0 km/sec, with rake vectors forming a fan structure around the hypocenter. We obtained a peak slip of 2.8 m and total seismic moment of 6.2 ?? 1018 Nm. We interpret the slip distribution as indicating that the hanging wall rotates relative to the footwall around the hypocenter, in a sense that appears consistent with the shape of the mapped fault trace. Copyright 2004 by the American Geophysical Union.

Geophysical Research Letters↗

Crustal deformation measurements in Guerrero, Mexico

GPS measurements of crustal deformation in Guerrero, southern Mexico, include surveys collected between 1992 and 2001 as well as continuous GPS measurements at a few sites. These geodetic observations are used to calculate interseismic deformation rates and assess the presence and possible location of transient deformation during the period encompassing 1992.25 to 2001.75. The data are used to examine transient deformation in 1998 previously described from data at a single site by Lowry et al. [2001]. Survey measurements and continuous data from a site near Popocate??petl volcano confirm the 1998 transient, and survey data also suggest another transient occurred following the 14 September 1995 (Mw = 7.3) Copala earthquake. All of the available GPS position estimates have been inverted for a combined model of slip during each event plus the steady state slip on the plate interface. Modeling of the steady state deformation rates confirms that the Guerrero seismic gap is partially frictionally locked at depths shallower than about 25 km and accumulating strain that may eventually be released in a great earthquake. The data also suggest that there is frictional coupling to much greater (>40 km) depths, which releases more frequently in aseismic slip events. The locations and sizes of the transient events are only partially constrained by the available data. However, the transient models which best fit the GPS coordinate time series suggest that aseismic slip was centered downdip of the seismogenic portion of the plate-bounding thrust in both events, and the moment release had equivalent magnitudes Mw = 7.1 + 1.3/-1.0 in 1995-1996 and 7.1 + 0.4/-0.1 in 1998. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Using 1-Hz GPS data to measure deformations caused by the denali fault earthquake

The 3 November 2002 moment magnitude 7.9 Denali fault earthquake generated large, permanent surface displacements in Alaska and large-amplitude surface waves throughout western North America. We find good agreement between strong ground-motion records integrated to displacement and 1-hertz Global Positioning System (GPS) position estimates collected ??? 140 kilometers from the earthquake epicenter. One-hertz GPS receivers also detected seismic surface waves 750 to 3800 kilometers from the epicenter, whereas these waves saturated many of the seismic instruments in the same region. High-frequency GPS increases the dynamic range and frequency bandwidth of ground-motion observations, providing another tool for studying earthquake processes.

Science↗

Volcano monitoring using the Global Positioning System: Filtering strategies

Permanent Global Positioning System (GPS) networks are routinely used for producing improved orbits and monitoring secular tectonic deformation. For these applications, data are transferred to an analysis center each day and routinely processed in 24-hour segments. To use GPS for monitoring volcanic events, which may last only a few hours, real-time or near real-time data processing and subdaily position estimates are valuable. Strategies have been researched for obtaining station coordinates every 15 min using a Kalman filter; these strategies have been tested on data collected by a GPS network on Kilauea Volcano. Data from this network are tracked continuously, recorded every 30 s, and telemetered hourly to the Hawaiian Volcano Observatory. A white noise model is heavily impacted by data outages and poor satellite geometry, but a properly constrained random walk model fits the data well. Using a borehole tiltmeter at Kilauea's summit as ground-truth, solutions using different random walk constraints were compared. This study indicates that signals on the order of 5 mm/h are resolvable using a random walk standard deviation of 0.45 cm/√h. Values lower than this suppress small signals, and values greater than this have significantly higher noise at periods of 1–6 hours.

Hawai'i↗

Crustal displacements due to continental water loading

The effects of long-wavelength (> 100 km), seasonal variability in continental water storage on vertical crustal motions are assessed. The modeled vertical displacements (??rM) have root-mean-square (RMS) values for 1994-1998 as large as 8 mm, with ranges up to 30 mm, and are predominantly annual in character. Regional strains are on the order of 20 nanostrain for tilt and 5 nanostrain for horizontal deformation. We compare ??rM with observed Global Positioning System (GPS) heights (??rO) (which include adjustments to remove estimated effects of atmospheric pressure and annual tidal and non-tidal ocean loading) for 147 globally distributed sites. When the ??rO time series are adjusted by ??rM, their variances are reduced, on average, by an amount equal to the variance of the ??rM. Of the ??rO time series exhibiting a strong annual signal, more than half are found to have an annual harmonic that is in phase and of comparable amplitude with the annual harmonic in the ??rM. The ??rM time series exhibit long-period variations that could be mistaken for secular tectonic trends or post-glacial rebound when observed over a time span of a few years.

Geophysical Research Letters↗