USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth and Planetary Science Letters”

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 361 records · Page 20Linked to original sources
USGS2018

Tracing enhanced oil recovery signatures in casing gases from the Lost Hills oil field using noble gases

Enhanced oil recovery (EOR) and hydraulic fracturing practices are commonly used methods to improve hydrocarbon extraction efficiency; however the environmental impacts of such practices remain poorly understood. EOR is particularly prevalent in oil fields throughout California where water resources are in high demand and disposal of high volumes of produced water may affect groundwater quality. Consequently, it is essential to better understand the fate of injected (EOR) fluids in California and other subsurface petroleum systems, as well as any potential effect on nearby aquifer systems. Noble gases can be used as tracers to understand hydrocarbon generation, migration, and storage conditions, as well as the relative proportions of oil and water present in the subsurface. In addition, a noble gas signature diagnostic of injected (EOR) fluids can be readily identified. We report noble gas isotope and concentration data in casing gases from oil production wells in the Lost Hills oil field, northwest of Bakersfield, California, and injectate gas data from the Fruitvale oil field, located within the city of Bakersfield. Casing and injectate gas data are used to: 1) establish pristine hydrocarbon noble-gas signatures and the processes controlling noble gas distributions, 2) characterize the noble gas signature of injectate fluids, 3) trace injectate fluids in the subsurface, and 4) construct a model to estimate EOR efficiency. Noble gas results range from pristine to significantly modified by EOR, and can be best explained using a solubility exchange model between oil and connate/formation fluids, followed by gas exsolution upon production. This model is sensitive to oil-water interaction during hydrocarbon expulsion, migration, and storage at reservoir conditions, as well as any subsequent modification by EOR.

California↗
USGS2018

Rapid late Miocene surface uplift of the Central Anatolian Plateau margin

The Central Anatolian Plateau (CAP), Turkey, is bordered to its south by a steep mountain belt that emerged ∼8–7 Ma ago from the Mediterranean Sea. Knowledge of the onset, duration and rate of surface uplift and orographic barrier formation along the plateau margin is crucial for understanding the geodynamic drivers of plateau uplift. We present a new comprehensive data set that includes 12 40 Ar/ 39 Ar ages and lacustrine carbonate δ 18 "> δ18 O data ( n = 637 "> n=637 ) from 13 sections in upper Oligocene to Pliocene continental basins of the CAP interior. We aim at documenting the development of a rain shadow and therefore the surface uplift history of the CAP and its southern margin (Tauride Mts.). In the rain shadow of the Tauride Mts. we observe a gradual 3.9‰ decrease of δ 18 "> δ18 O values of lacustrine carbonate between ∼11 and 5 Ma that we interpret to originate from a similar change in δ 18 "> δ18 O values of precipitation owing to the late Miocene development of an orographic barrier. Our stable isotope paleoaltimetry data show that by 5 Ma the southern CAP margin had reached similar-to-present elevations of ∼2 km. Surface uplift was coeval with ignimbritic magmatism, forearc shortening and distributed compression. We suggest that the removal of lithospheric mantle below Anatolia led to surface uplift of the CAP interior, which was followed by surface uplift of the southern CAP margin due to crustal thickening as a result of northward subduction of the African plate below central Anatolia.

Earth and Planetary Science Letters↗
USGS2018

Pliocene erosional pulse and glacier-landscape feedbacks in the western Alaska Range

Pliocene–Pleistocene glaciation modified the topography and erosion of most middle- and high-latitude mountain belts, because the evolution of catchment topography controls long-term glacier mass balance and erosion. Hence, characterizing how erosion rates change during repeated glaciations can help test hypothesized glacier erosion-landscape feedbacks across a range of settings. To better understand how glaciations and landscapes coevolve on geologic timescales, I quantify erosion rates in the glaciated western Alaska Range with low-temperature thermochronometric data and modeling. Zircon (U–Th)/He and apatite fission track data suggest mountain-building was underway by the early Miocene. In contrast, lower-temperature apatite (U–Th)/He age-elevation and grain age-kinetic data indicate that erosion accelerated coincident with regional Pliocene glaciation ca. 4 Ma. Furthermore, erosion rates calculated within an eroding half-space indicate slow erosion at a rate ≤0.3 km/m.y. before 4.2 Ma, an initial pulse of rapid erosion at a rate of 1.0–1.6 km/m.y. during 4.2–2.9 Ma, and more moderate erosion at a rate of 0.4–0.7 km/m.y. since 2.9 Ma. The initial erosion pulse suggests a significant transient landscape adjustment to the introduction of efficient glacial erosion. The subsequent decrease in Pleistocene erosion rates is consistent with a negative feedback between continuing glaciation and glacier size/erosivity: If glacial erosion outpaces rock uplift, glacier erosion decreases over time as topography, mass balance, valley gradients, and ice flux are reduced. These findings imply that in areas of moderate rock uplift rates, the onset of local Plio–Pleistocene glaciation may have been punctuated by an initial pulse of rapid landscape change, after which change became more gradual.

Alaska↗
USGS2019

Seismic velocity structure across the 2013 Craig, Alaska rupture from aftershock tomography: Implications for seismogenic conditions

The 2013 Craig, Alaska M W 7.5 earthquake ruptured along ∼150 km of the Queen Charlotte Fault (QCF), a right-lateral strike-slip plate boundary fault separating the Pacific and North American plates. Regional shear wave analyses suggest that the Craig earthquake rupturepropagated in the northward direction faster than the S-wave (supershear). Theoretical studies suggest that a bimaterial interface, such as that along the QCF, which separates oceanic and continental crust with differing elastic properties, can promote supershear rupture propagation. We deployed short-period ocean-bottom seismometers (OBS) as a part of a rapid-response effort less than four months after the Craig earthquake mainshock. During a 21-day period, 1,133 aftershocks were recorded by 8 OBS instruments. Aftershock spatial distribution indicates that the base of the seismogenic zone along the QCF approaches ∼25 km depth, consistent with a thermally-controlled fault rheology expected for igneous rocks at oceanic transform faults. The spatial distribution also provides supporting evidence for a previously hypothesized active strand of the QCF system within the Pacific Plate. Tomographic traveltime inversion for velocity structure indicates a low-velocity (V P and V S ) zone on the Pacific side of the plate boundary at 5–20 km depths, where NeogenePacific crust and upper mantle seismic velocities average ∼3–11% slower than the North American side, where the Paleozoic North American crust is seismically faster. Our results suggest that elastic properties along the studied portion of the QCF are different than those of a simple oceanic–continental plate boundary fault. In our study region, velocity structure across the QCF, while bimaterial, does not support faster material on the west side of the fault, which has been proposed as one possible explanation for northward supershear propagation during the Craig earthquake. Instead, we image low-velocity material on the west side of the fault. Explanations could include that part of the rupture was subshear, or that fault damage zone properties or fault smoothness are more important controls on supershear rupture than a bimaterial contrast.

Alaska↗
USGS2018

Eruptions in sync: Improved constraints on Kīlauea Volcano's hydraulic connection

Kīlauea Volcano is an archetype for the complex interactions that can occur between a volcano’s summit and flanks. Decades of monitoring at Kīlauea have demonstrated that magma rises beneath the summit and flows laterally at shallow depths to erupt along the rift zones. Kīlauea’s recent eruptions at Halema‘uma‘u and Pu‘u ‘Ō‘ō mark the first time in the historic record that long-term (>1 year) eruptions have been concurrent at the summit and a rift zone, offering a new opportunity to improve our understanding of the relationship between these two segments of the magmatic system. While magma supply rate beneath the summit has been shown in previous studies to be a primary control on magmatic system pressure and eruptive activity, the role of the eruptive vent has been less clear. Our study shows that a dynamic equilibrium is maintained between Kīlauea’s summit and East Rift Zone (ERZ) eruptive vent—and lava lake level fluctuations are closely coupled at the two eruption sites—providing new constraints on the hydraulic connection and ERZ conduit. We show that localized changes at the ERZ eruption site during 2010-2011 regulated summit behavior in an uprift direction over distances of ~20 km. Changes in the elevation and efficiency of the ERZ vent affect pressure in Kīlauea’s magmatic system and impact summit behavior. Thus, the hydraulic connection between the summit and rift zone is a “two-way street” that transmits both downrift- and uprift-directed changes. Our results support recent work at other volcanoes that shows a complex interplay between a volcano’s summit reservoir and flank conduit during flank eruptions, and suggest that explosive summit activity may in some cases be triggered by changes far away on a volcano’s rift.

Hawaii↗
USGS2019

Detrital K-feldspar Pb isotopic evaluation of extraregional sediment transported through an Eocene tectonic breach of southern California's Cretaceous batholith

Sedimentary provenance studies have come to be overwhelmingly based upon U–Pb geochronologic measurements performed with detrital zircon while alternative and potentially complementary approaches such as conglomerate clast studies and heavy mineral analysis have faded in importance. Measurement of Pb isotopic compositions in detrital K-feldspar is among the under-utilized approaches available to ascertain sedimentary source regions. While it has been long recognized that common Pb isotope compositions recorded by K-feldspar vary widely and reflect the crustal provinces from which the host basement rocks crystallized, use of the approach has suffered due to a lack of appropriate statistical models and ground truth compositional data from source regions. In this paper, we: (1) present high-throughput LA-ICPMS analysis protocols needed to generate statistically meaningful detrital K-feldspar Pb isotope data sets; (2) develop an interpretative approach based upon 208 Pb/ 206 Pb vs. 207 Pb/ 206 Pb that incorporate information from the U- and Th-decay systems into one two-dimensional plot that is amenable to analysis using two-dimensional Kolmogorov–Smirnoff statistical tests; (3) generate new Pb isotopic data from basement rocks from southwestern North America to improve knowledge of the Pb isotopic properties of potential source regions; and (4) generate new Pb isotopic data from Lower Eocene to Lower Miocene sedimentary rocks to evaluate changes in drainage patterns that occurred in response to deformation that affected the southern California margin. Through this case study, we demonstrate how our new analytical and interpretative methods could be profitably applied to future geochemical and provenance studies and tectonically driven re-organization of drainage patterns.

Earth and Planetary Science Letters↗
USGS2019

Hydrogen isotopes in high 3He/4He submarine basalts: Primordial vs. recycled water and the veil of mantle enrichment

The hydrogen isotope value ( δ D) of water indigenous to the mantle is masked by the early degassing and recycling of surface water through Earth's history. High 3 He/ 4 He ratios in some ocean island basalts, however, provide a clear geochemical signature of deep, primordial mantle that has been isolated within the Earth's interior from melting, degassing, and convective mixing with the upper mantle. Hydrogen isotopes were measured in high 3 He/ 4 He submarine basalt glasses from the Southeast Indian Ridge (SEIR) at the Amsterdam–St. Paul (ASP) Plateau ( δ D = −51 to −90‰, 3 He/ 4 He = 7.6 to 14.1 R A ) and in submarine glasses from Loihi seamount south of the island of Hawaii ( δ D = −70 to −90‰, 3 He/ 4 He = 22.5 to 27.8 R A ). These results highlight two contrasting patterns of δ D for high 3 He/ 4 He lavas: one trend toward high δ D of approximately −50‰, and another converging at δ D = −75‰. These same patterns are evident in a global compilation of previously reported δ D and 3 He/ 4 He results. We suggest that the high δ D values result from water recycled during subduction that is carried into the source region of mantle plumes at the core–mantle boundary where it is mixed with primordial mantle, resulting in high δ D and moderately high 3 He/ 4 He. Conversely, lower δ D values of −75‰, in basalts from Loihi seamount and also trace element depleted mid-ocean ridge basalts, imply a primordial Earth hydrogen isotopic value of −75‰ or lower. δ D values down to −100‰ also occur in the most trace element-depleted mid-ocean ridge basalts, typically in association with 87 Sr/ 86 Sr ratios near 0.703. These lower δ D values may be a result of multi-stage melting history of the upper mantle where minor D/H fractionation could be associated with hydrogen retention in nominally anhydrous residual minerals. Collectively, the predominance of δ D around −75‰ in the majority of mid-ocean ridge basalts and in high 3 He/ 4 He Loihi basalts is consistent with an origin of water on Earth that was dominated by accretion of chondritic material.

Earth and Planetary Science Letters↗
USGS2019

A constant slip rate for the western Qilian Shan frontal thrust during the last 200 ka consistent with GPS-derived and geological shortening rates

Active thrust faulting at the front of the Qilian Shan accommodates the northeastward growth of the Tibetan Plateau, however, the lifespan of individual faults and their slip history on different timescales remain largely unknown. Here, we show that the main range-bounding thrust fault of the western Qilian Shan has accrued tectonic slip at an almost constant rate during the last ∼200 ka, and possibly since fault initiation in the mid-Miocene. Our finding is based on 10Be ages from a flight of five deformed fluvial terraces along the Hongshuiba river, which constrain the vertical slip rate of the Qilian Shan frontal thrust to be m/ka during the last 200 ka. With a fault dip of ° constrained by seismic reflection data, we obtain a horizontal shortening rate of m/ka. This value is consistent with both the short-term shortening rate of mm/a derived from GPS data and the long-term shortening rate of m/ka, which is based on a balanced geological cross-section. The latter provides a total shortening estimate of km since the thrust fault initiated Ma ago. The agreement between the shortening rates on the range of timescales between 100 and 107 years suggests that the western Qilian Shan frontal thrust has slipped at a steady rate since its initiation and implies that this fault is the main structure responsible for the growth of the western Qilian Shan.

Tibetan Plateau, Qilian Shan↗
USGS2019

Comment on “Particle fluxes in groundwater change subsurface rock chemistry over geologic time”

Over the last decade, studies at the Shale Hills Critical Zone Observatory (Shale Hills) have greatly expanded knowledge of weathering in previously understudied, shale-mantled terrains, as well as Earth's Critical Zone as a whole. Among the many discoveries made was the importance of redistribution and losses of micron-sized particles during development of shale-derived soils. A geochemical fingerprint of this process for Al and Fe was illustrated quantitatively by Jin et al. (2010). Subsequent papers, too numerous to list in a Comment, built upon this new recognition by evaluating the spatial and temporal aspects element mobilization. Recently, Kim et al. (2018) examined the composition of suspended, generally micron-sized particles in the Shale Hills stream, along with the dissolved load, across seasons and ranges of discharge. One prominent conclusion from Kim et al. (2018) is that Zr is essentially immobile at Shale Hills. Such a broad conclusion is in direct contradiction with one from Bern and Yesavage (2018) that Zr has been mobilized from soils at Shale Hills, and the losses relative to soil parent material are significant (median 41%). The point is important, because assuming Zr immobility is necessary to index gains and losses of other elements using the open-chemical-system transport function ( τ ). Both papers draw upon patterns and calculations using elemental concentration data from Shale Hills and attempt to construct conceptual frameworks to explain the results. Here, the argument is made that the understanding of substantial Zr mobility from soils at Shale Hills described by Bern and Yesavage (2018) is more accurate. Additionally, issues with adaptations of the standard τ equations used in Kim et al. (2018) and some previous papers are also addressed.

Pennsylvania↗
USGS2019

Ten-million years of activity within the Eastern California Shear Zone from U-Pb dating of fault-zone opal

Reconstructions of long-term fault activity are essential for understanding both the mechanisms controlling fault behavior and accurate earthquake hazard assessments. Increasing evidence for temporal variations in strain accumulation suggests non-uniform strain rates over a range of historic to geologic timescales. The paucity of long-term records of fault activity has limited our ability to resolve these variations. We present a method for constraining long-term fault activity based on U-Pb dating of fault-related opal from secondary fault segments within the Eastern California Shear Zone (ECSZ). The presence of sheared and breccia-cemented opaline silica within well-exposed faults at near-surface conditions suggest that opal formation is associated with high-magnitude earthquakes capable of surface rupture (>6 M). Temporal constraints from sheared syntectonic opal (n=74) on related secondary faults from this study provide new insights on the timing of fault initiation, reactivation, and longevity. The oldest dates obtained indicate that ECSZ activity commenced at or before 10 Ma. Multiple deformation events dated within a single structure on episodically deposited and sheared opal (up to six generations), demonstrate that fault reactivation occurred over 105 year timescales (0.7-0.1 Ma). Relative probabilities of dated deformation events can be used to evaluate changes in fault activity in the past 2.5 Ma (n=60). This analysis indicates enhanced fault activity starting at 2 Ma and peaking around 1 Ma, possibly due to fault-interactions and distribution of deformation between the ECSZ and the San Andreas Fault.

California↗
USGS2019

Physicochemical models of effusive rhyolitic eruptions constrained with InSAR and DEM data: A case study of the 2011-2012 Cordon Caulle eruption

The 9 month long 2011-2012 eruption of Cordon Caulle (Southern Andes, Chile) is the best instrumentally recorded rhyolitic eruption to date and the first time that the effusion of a rhyolitic flow has been observed in detail. We use Interferometric Synthetic Aperture Radar (InSAR), with time-lapse DEMs and numerical models to study the dynamics of coupled magma reservoir deflation and lava effusion. InSAR recorded ~2.2-2.5 m of deflation after the first three days of the eruption, which can be modeled using a spheroidal magma reservoir at a depth of ∼5 km, ∼20 km long, and with a pressure drop of ~20-30 MPa. The source is elongated in the NW-SE direction and its large dimensions imply a large plumbing system spanning neighboring volcanoes and active throughout the eruption, with a slight change halfway through the effusive phase. TanDEM-X DEMs record the extrusion of both the rhyolitic lava flow and the intrusion of a shallow laccolith around the eruptive vent, with a total volume of ~1.2 km3 DRE. The laccolith was emplaced during the first month of the eruption, during both the eruption explosive and effusive stages. Both the reservoir pressure drop and the extruded volume time series follow quasi-exponential trends, and can be explained by a model that couples the reservoir pres- sure decrease, time- and pressured ependent variations in the magma properties inside of the reservoir, and conduit flow. This model predicts both the temporal evolution and amplitude of both time series, and a magma compressibility of ∼10^−10 Pa−1, half the compressibility of the magma of the sub-Plinian explosive phase. Further, we estimate that the reservoir contained 1-3 wt.% dissolved H2O at the onset of lava effusion, with no exsolved CO2 and H2O in the reservoir throughout the eruption. This is in accord with a magma that was significantly degassed after the explosive phase. These remaining volatiles might have been responsible for magma fragmentation, consistent with the hybrid explosive and effusive style observed during the waning of the eruption.

Cordón Caulle↗
USGS2019

Clustered BSRs: Evidence for gas hydrate-bearing turbidite complexes in folded regions, example from the Perdido Fold Belt, northern Gulf of Mexico

We describe previously undocumented but extensive gas hydrate accumulations in the mouth of Perdido Canyon in the northern Gulf of Mexico. The accumulations are located within central parts of structural domes (four-way closures) and are characterized by stacked, high-amplitude bottom simulating reflections (BSRs) that we call clustered BSRs. Seismic data from Perdido Canyon show two clustered BSRs associated with turbidite sequences within two dome folds formed from tectonic folding and salt diapir rise. The northwestern (NW) and southeastern (SE) clustered BSRs have aerial extents of ~25 km2 and 50 km2, respectively. Well log data confirm gas hydrate occurs above the NW clustered BSR, within a 225 m-thick consistently high-resistivity interval that we interpret as gas hydrate in near-vertical fractures and turbidite sands. The SE dome is only drilled at the edge of the BSR; nevertheless, the well log data indicate that a 30 m-thick gas hydrate accumulation is present. Gas chromatographic logs in both domes suggest a gradual transition from predominantly microbial gas below the BSR (500–1000 meters below seafloor (mbsf)) to thermogenic gas at 1000–2000 mbsf. Based on the well log data and seismic stratigraphic analysis, we find gas hydrate is concentrated in fractures in marine mud, as well as in the pores of submarine fan turbidities, where saturations reach as high as 75%. An estimate of the total gas hydrate-bound gas volume at standard temperature and pressure is between 0.04 and 0.17 trillion cubic meters (TCM) assuming average hydrate saturation of 5-20% in a ~45 m thick turbidite sand unit above the Perdido Canyon BSR area. Measured BSR extent and gas volume estimates indicate that the NW and SE reservoirs are among the largest gas hydrate occurrences known in the Gulf of Mexico.

Perdido Fold Belt↗
USGS2020

Plate boundary localization, slip-rates and rupture segmentation of the Queen Charlotte Fault based on submarine tectonic geomorphology

Linking fault behavior over many earthquake cycles to individual earthquake behavior is a primary goal in tectonic geomorphology, particularly across an entire plate boundary. Here, we examine the 1150-km-long, right-lateral Queen Charlotte-Fairweather fault system using comprehensive multibeam bathymetry data acquired along the Queen Charlotte Fault (QCF) offshore southeastern Alaska and western British Columbia. Fine-scale analysis of tectonic geomorphology allowed us to identify and reconstruct 184 strike-slip piercing points over a 630 km stretch of the QCF. Age constraints from glacial recession and offshore sedimentation patterns yield a consistent slip-rate of ∼50–57 mm/yr since ∼17–12 ka, the fastest rate for a continent-ocean strike-slip fault on Earth. These slip-rates equal or exceed estimates of Pacific-North America (PA-NA) relative motion from global plate reconstructions, indicating that PA-NA motion is highly localized. The QCF cuts the seafloor along a narrow and unusually straight trace for its entire length and multiple fault traces are observed only at local step-overs. The geometry and behavior of the QCF over many earthquake cycles is simple and typical of mature faults with relatively homogeneous stress fields. Since the QCF is the primary PA-NA plate boundary, we used the trace of the QCF to define the small circle path for relative plate motion and computed the associated Euler pole. Predicted along-strike obliquity variations based on the new pole agree with observed tectonic geomorphology and suggest that previous global plate reconstructions overestimated the degree of oblique convergence along the QCF. We also find that subtle, long-wavelength (75–150 km) bends and discrete step-overs appear to define the endpoints of M>7 earthquakes, suggesting that obliquity and resultant fault geometry may control rupture segmentation and asperity development. Lastly, the agreement between predicted obliquity and tectonic geomorphology along the entire length of QCF compelled a reevaluation of regional tectonic models. In the north, the eastern Yakatat Terrane appears to be translating northwest with the Pacific plate, and slip transferred from the QCF to the Fairweather Fault results in ∼20 mm/yr of convergence along the southern St. Elias mountains. In the south, we predict a reduced rate of convergence along the QCF west of Haida Gwaii (∼5–6 mm/yr of shortening, on average) relative to previous studies. Our results support a model for transpression and strike-slip partitioning along the edge of a hot and weak Pacific Plate, leading to crustal thickening and growth of the Queen Charlotte Terrace to the west of Haida Gwaii.

Alaska, British Columbia↗
USGS2019

Late Quaternary slip rate of the Central Sierra Madre fault, southern California: Implications for slip partitioning and earthquake hazard

The Sierra Madre fault system accommodates contraction within a large restraining bend area of the San Andreas fault along the northern margin of the Los Angeles metropolitan area in Southern California. Reverse slip along this fault system during earthquakes controls growth of the San Gabriel Mountains and poses a significant seismic hazard to the region. Here, we measure the late Quaternary slip rate of the Central Sierra Madre fault (CSMF) using analysis of high-resolution topography combined with cosmogenic 10Be surface exposure dating and post-IR IRSL geochronology. We mapped terrace and fan surfaces from three drainages that cross the CSMF and correlate them based on soil development and geomorphic position. Cosmogenic nuclide and luminescence ages are consistent amongst the three prominent surfaces offset ~5 to 28 m across the fault zone. We devised a new strategy to estimate surface ages, incorporating data from two dating methods at three locations, refined by inset age relationships, that yields surface age estimates of 54 +21/-13 ka, 36 ± 8 ka, and 12 ± 4 ka. Estimated slip for these geomorphic markers is more uncertain than the measured vertical separation due to uncertainties in fault dip and ranges from 7.5 +5.4/-3.1 m to 58.5 +46.3/-14.4 m. Incremental dip-slip rate estimates from different age surfaces and locations overlap within uncertainty, with median values ranging from 0.7 to 1.1 mm/yr. The average slip rate for all three generations of markers is 1.1 +1.2/-0.4 mm/yr. This late Quaternary slip rate for the CSMF is slower than estimates based on interseismic geodetic data, and emphasizes the importance of contraction distributed across multiple structures south of the Sierra Madre fault when assessed against the geodetic shortening budget. Despite being the central portion of the broader Sierra Madre fault system, the CSMF has a slip rate similar to or lower than neighboring sections, suggesting that slip transfer onto other nearby faults control the along-strike pattern of deformation rate. Paleoseismic evidence indicates that the last earthquake on the CSMF was in the early Holocene, and the slip rate we estimate suggests that the accumulated elastic strain could produce many meters of slip in future earthquakes.

California↗
USGS2019

What determines the downstream evolution of turbidity currents?

Seabed sediment flows called turbidity currents form some of the largest sediment accumulations, deepest canyons and longest channel systems on Earth. Only rivers transport comparable sediment volumes over such large areas; but there are far fewer measurements from turbidity currents, ensuring they are much more poorly understood. Turbidity currents differ fundamentally from rivers, as turbidity currents are driven by the sediment that they suspend. Fast turbidity currents can pick up sediment, and self-accelerate (ignite); whilst slow flows deposit sediment and dissipate. Self-acceleration cannot continue indefinitely, and flows might reach a near-uniform state (autosuspension). Here we show how turbidity currents evolve using the first detailed measurements from many locations along their pathway, which come from Monterey Canyon offshore California. All flows initially ignite. Typically, initially-faster flows then achieve near-uniform velocities (autosuspension), whilst slower flows dissipate. Fractional increases in initial velocity favour much longer runout, and a new model explains this bifurcating behaviour. However, the only flow during less-stormy summer months is anomalous as it self-accelerated, which is perhaps due to erosion of surficial-mud layer with fine sands mid-canyon. Turbidity current evolution is therefore highly sensitive to both initial velocities and seabed character.

Earth and Planetary Science Letters↗
USGS2020

Seismic and geodetic progression of the 2018 summit caldera collapse of Kīlauea Volcano

The 2018 eruption of Kīlauea volcano, Hawaiʻi, resulted in a major collapse of the summit caldera along with an effusive eruption in the lower East Rift Zone. The caldera collapse comprised 62 highly similar collapse cycles of strong ground deformation and earthquake swarms that ended with a magnitude 5 collapse event and one partial cycle that did not end with a collapse event. We analyzed geodetic and seismic data to better understand how the caldera collapse progressed over 3 months of activity, focusing on the cyclical activity. We identified 3 main phases of collapse: initial ring-fault activation and small explosions (Phase 1), an eastward shift in activity and freeing of the central piston (Phase 2), and a recoupling of the piston to the reservoir followed by relatively steady behavior until the eruption’s end (Phase 3). Additionally, we observed geodetic evidence of tangential motion from the localization of the main ring fault (Phase 2) and the formation of a major peripheral ring fault on the eastern side of the collapse caldera during Phase 3. Both geodetic and seismic parameters suggest that the collapse may have had an eastward-component of motion after the ring fault system had formed. The cyclical seismic and geodetic parameters show no obvious signs that the collapse was coming to an end, with the only notable change being a significant increase in the ratio of cyclical displacement to co-collapse displacement observed during the last complete cycle on GNSS stations outside the caldera region.

Hawaii↗
USGS2020

Mysterious tsunami in the Caribbean Sea following the 2010 Haiti earthquake possibly generated by dynamically triggered early aftershocks

Dynamically triggered offshore aftershocks, caused by passing seismic waves from main shocks located on land, are currently not considered in tsunami warnings. The M7.0 2010 Haiti earthquake epicenter was located on land 27 km north of the Caribbean Sea and its focal mechanism was oblique strike-slip. Nevertheless, a tsunami recorded on a Caribbean Deep-Ocean Assessment and Reporting of Tsunami (DART) buoy and a tide gauge produced runup heights of 1–3 m along Haiti southeast coast. Earthquake finite-fault model inversions of the DART waveform suggest that a reverse fault doublet with magnitudes of M6.8 and M6.5 located 85 km southwest of the epicenter may have excited the tsunami. This doublet collocates with dynamically triggered aftershocks, derived from back-projection analysis, that occurred 20-60 s after the main shock of the Haiti earthquake. The aftershocks are within a region of maximum dynamic strain predicted by the main shock, on a possibly tectonically active submarine ridge southwest of Haiti's Southern Peninsula. The agreement between the tsunami finite-fault source models and the seismic and tectonic evidence suggests that earthquakes on land, even strike-slip faults, can generate tsunamis by dynamically triggering offshore aftershocks.

Earth and Planetary Science Letters↗
USGS2021

Material failure and caldera collapse: Insights from the 2018 Kilauea eruption

The Failure Forecast Method (FFM) was introduced as an empirical model for forecasting catastrophic material failures related to natural hazards, such as landslides and volcanic eruptions, with mixed success. During the 2018 eruption of Kilauea volcano, Hawaii, the draining of the summit magma reservoir into the Lower East Rift Zone resulted in the formation of a new caldera at the summit. I tested the applicability of the FFM to caldera collapse by analyzing the cyclical earthquake swarms and ground deformation that occurred between 62 sudden major caldera collapse events. The progression of both the cumulative moment release of the cyclical earthquakes and the GNSS displacement show a major change in mid-June. In late May through early June, the progression of the parameters is consistent with strain localization or creep progression related to the development or activation of the ring fault system. From late June until the end of the eruption, parameter progression is roughly steady with initial accelerating increases in cumulative moment and displacement that shift to approximately linear progression. Analysis of repeating earthquake families in the cyclical swarms showed that the behavior of the repeaters was consistent with that of the cyclical swarms as a whole and suggested that each family undergoes its own progression of activation to termination. While the FFM analysis identified the system change in mid-June, it did not demonstrate an ability to forecast collapse events or the end of the eruption.

Hawaii↗