USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geochemistry Geophysics Geosystems”

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 19 recordsLinked to original sources
USGS2022

A late Cenozoic kinematic model for fault motion within greater Cascadia

Widely accepted tectonic reconstructions indicate at least 100 km of coast-parallel northwestward translation of the Sierra Nevada block of California and 15–20° clockwise rotation of most of Oregon since the current phase of Basin and Range extension began ∼17 Ma. These reconstructions require at least 100 km of convergence between the central Coast Range of Oregon and rigid North America in mainland British Columbia, yet there is little discussion of how such convergence might be distributed. This study offers a kinematic model of the distribution of such deformation, constrained by geodesy, paleomagnetism, and fault offsets in Nevada, California and Oregon. The model includes differential rotation across the thrust faults of the Yakima fold and thrust belt (YFTB), compressive right-lateral faulting in the Washington Cascade Range, substantial thrust faulting within the Puget Lowland, and oroclinal bending and doming in the Olympic Mountains. Shortening across YTFB along 120°W longitude is modeled as 47 km, across Puget Lowland at 123°W (Olympia-Bellingham) is 94 km, and total shortening between the central Oregon Coast Range and northern Washington (Corvallis-Bellingham) is 125 km. Current motion of the coastal regions above the Cascadia subduction zone results from both permanent deformation of the continent and elastic coupling to the subducting plate. Permanent deformation in the model is based on extrapolating geodesy from east of 120°W or south of 40°N, indicating a very uniform convergence velocity with the Juan de Fuca plate for northernmost California and Oregon near 31 mm/yr at N61°E.

California, Idaho, Montana, Nevada, Oregon, Washin↗
USGS2014

Relationship between the Cascadia fore-arc mantle wedge, nonvolcanic tremor, and the downdip limit of seismogenic rupture

Great earthquakes anticipated on the Cascadia subduction fault can potentially rupture beyond the geodetically and thermally inferred locked zone to the depths of episodic tremor and slip (ETS) or to the even deeper fore-arc mantle corner (FMC). To evaluate these extreme rupture limits, we map the FMC from southern Vancouver Island to central Oregon by combining published seismic velocity structures with a model of the Juan de Fuca plate. These data indicate that the FMC is somewhat shallower beneath Vancouver Island (36–38 km) and Oregon (35–40 km) and deeper beneath Washington (41–43 km). The updip edge of tremor follows the same general pattern, overlying a slightly shallower Juan de Fuca plate beneath Vancouver Island and Oregon (∼30 km) and a deeper plate beneath Washington (∼35 km). Similar to the Nankai subduction zone, the best constrained FMC depths correlate with the center of the tremor band suggesting that ETS is controlled by conditions near the FMC rather than directly by temperature or pressure. Unlike Nankai, a gap as wide as 70 km exists between the downdip limit of the inferred locked zone and the FMC. This gap also encompasses a ∼50 km wide gap between the inferred locked zones and the updip limit of tremor. The separation of these features offers a natural laboratory for determining the key controls on downdip rupture limits.

Cascadia fore-arc mantle wedge↗
USGS2014

Modified expression for bulb-tracer depletion—Effect on argon dating standards

40 Ar/ 39 Ar geochronology depends critically on well-calibrated standards, often traceable to first-principles K-Ar age calibrations using bulb-tracer systems. Tracer systems also provide precise standards for noble-gas studies and interlaboratory calibration. The exponential expression long used for calculating isotope tracer concentrations in K-Ar age dating and calibration of 40 Ar/ 39 Ar age standards may provide a close approximation of those values, but is not correct. Appropriate equations are derived that accurately describe the depletion of tracer reservoirs and concentrations of sequential tracers. In the modified expression the depletion constant is not in the exponent, which only varies as integers by tracer-number. Evaluation of the expressions demonstrates that systematic error introduced through use of the original expression may be substantial where reservoir volumes are small and resulting depletion constants are large. Traditional use of large reservoir to tracer volumes and the resulting small depletion constants have kept errors well less than experimental uncertainties in most previous K-Ar and calibration studies. Use of the proper expression, however, permits use of volumes appropriate to the problems addressed.

Geochemistry, Geophysics, Geosystems↗
USGS2014

Soil CO2 emissions as a proxy for heat and mass flow assessment, Taupō Volcanic Zone, New Zealand

The quantification of heat and mass flow between deep reservoirs and the surface is important for understanding magmatic and hydrothermal systems. Here, we use high-resolution measurement of carbon dioxide flux (φCO 2 ) and heat flow at the surface to characterize the mass (CO 2 and steam) and heat released to the atmosphere from two magma-hydrothermal systems. Our soil gas and heat flow surveys at Rotokawa and White Island in the Taupō Volcanic Zone, New Zealand, include over 3000 direct measurements of φCO 2 and soil temperature and 60 carbon isotopic values on soil gases. Carbon dioxide flux was separated into background and magmatic/hydrothermal populations based on the measured values and isotopic characterization. Total CO 2 emission rates (ΣCO 2 ) of 441 ± 84 t d −1 and 124 ± 18 t d −1 were calculated for Rotokawa (2.9 km 2 ) and for the crater floor at White Island (0.3 km 2 ), respectively. The total CO 2 emissions differ from previously published values by +386 t d −1 at Rotokawa and +25 t d −1 at White Island, demonstrating that earlier research underestimated emissions by 700% (Rotokawa) and 25% (White Island). These differences suggest that soil CO 2 emissions facilitate more robust estimates of the thermal energy and mass flux in geothermal systems than traditional approaches. Combining the magmatic/hydrothermal-sourced CO 2 emission (constrained using stable isotopes) with reservoir H 2 O:CO 2 mass ratios and the enthalpy of evaporation, the surface expression of thermal energy release for the Rotokawa hydrothermal system (226 MW t ) is 10 times greater than the White Island crater floor (22.5 MW t ).

Taupō Volcanic Zone↗
USGS2014

Correlated patterns in hydrothermal plume distribution and apparent magmatic budget along 2500 km of the Southeast Indian Ridge

Multiple geological processes affect the distribution of hydrothermal venting along a mid-ocean ridge. Deciphering the role of a specific process is often frustrated by simultaneous changes in other influences. Here we take advantage of the almost constant spreading rate (65–71 mm/yr) along 2500 km of the Southeast Indian Ridge (SEIR) between 77°E and 99°E to examine the spatial density of hydrothermal venting relative to regional and segment-scale changes in the apparent magmatic budget. We use 227 vertical profiles of light backscatter and (on 41 profiles) oxidation-reduction potential along 27 first and second-order ridge segments on and adjacent to the Amsterdam-St. Paul (ASP) Plateau to map p h , the fraction of casts detecting a plume. At the regional scale, venting on the five segments crossing the magma-thickened hot spot plateau is almost entirely suppressed ( p h  = 0.02). Conversely, the combined p h (0.34) from all other segments follows the global trend of p h versus spreading rate. Off the ASP Plateau, multisegment trends in p h track trends in the regional axial depth, high where regional depth increases and low where it decreases. At the individual segment scale, a robust correlation between p h and cross-axis inflation for first-order segments shows that different magmatic budgets among first-order segments are expressed as different levels of hydrothermal spatial density. This correlation is absent among second-order segments. Eighty-five percent of the plumes occur in eight clusters totaling ∼350 km. We hypothesize that these clusters are a minimum estimate of the length of axial melt lenses underlying this section of the SEIR.

Geochemistry, Geophysics, Geosystems↗
USGS2014

Eruptive and tectonic history of the Endeavour Segment, Juan de Fuca Ridge, based on AUV mapping data and lava flow ages

High-resolution bathymetric surveys from autonomous underwater vehicles ABE and D. Allan B. were merged to create a coregistered map of 71.7 km2 of the Endeavour Segment of the Juan de Fuca Ridge. Radiocarbon dating of foraminifera in cores from three dives of remotely operated vehicle Doc Ricketts provide minimum eruption ages for 40 lava flows that are combined with the bathymetric data to outline the eruptive and tectonic history. The ages range from Modern to 10,700 marine-calibrated years before present (yr BP). During a robust magmatic phase from >10,700 yr BP to ~4300 yr BP, flows erupted from an axial high and many flowed >5 km down the flanks; some partly buried adjacent valleys. Axial magma chambers (AMCs) may have been wider than today to supply dike intrusions over a 2 km wide axial zone. Summit Seamount formed by ~4770 yr BP and was subsequently dismembered during a period of extension with little volcanism starting ~4300 yr BP. This tectonic phase with only rare volcanic eruptions lasted until ~2300 yr BP and may have resulted in near-solidification of the AMCs. The axial graben formed by crustal extension during this period of low magmatic activity. Infrequent eruptions occurred on the flanks between 2620–1760 yr BP and within the axial graben since ~1750 yr BP. This most recent phase of limited volcanic and intense hydrothermal activity that began ~2300 yr BP defines a hydrothermal phase of ridge development that coincides with the present-day 1 km wide AMCs and overlying hydrothermal vent fields.

Geochemistry, Geophysics, Geosystems↗
USGS2014

Low-frequency earthquakes reveal punctuated slow slip on the deep extent of the Alpine Fault, New Zealand

We present the first evidence of low-frequency earthquakes (LFEs) associated with the deep extension of the transpressional Alpine Fault beneath the central Southern Alps of New Zealand. Our database comprises a temporally continuous 36 month-long catalog of 8760 LFEs within 14 families. To generate this catalog, we first identify 14 primary template LFEs within known periods of seismic tremor and use these templates to detect similar events in an iterative stacking and cross-correlation routine. The hypocentres of 12 of the 14 LFE families lie within 10 km of the inferred location of the Alpine Fault at depths of approximately 20–30 km, in a zone of high P-wave attenuation, low P-wave speeds, and high seismic reflectivity. The LFE catalog consists of persistent, discrete events punctuated by swarm-like bursts of activity associated with previously and newly identified tremor periods. The magnitudes of the LFEs range between M L – 0.8 and M L 1.8, with an average of M L 0.5. We find that the frequency-magnitude distribution of the LFE catalog both as a whole and within individual families is not consistent with a power law, but that individual families' frequency-amplitude distributions approximate an exponential relationship, suggestive of a characteristic length-scale of failure. We interpret this LFE activity to represent quasi-continuous slip on the deep extent of the Alpine Fault, with LFEs highlighting asperities within an otherwise steadily creeping region of the fault.

Alpine Fault↗
USGS2014

Modeling ash fall distribution from a Yellowstone supereruption

We used the volcanic ash transport and dispersion model Ash3d to estimate the distribution of ashfall that would result from a modern-day Plinian supereruption at Yellowstone volcano. The simulations required modifying Ash3d to consider growth of a continent-scale umbrella cloud and its interaction with ambient wind fields. We simulated eruptions lasting 3 days, 1 week, and 1 month, each producing 330 km 3 of volcanic ash, dense-rock equivalent (DRE). Results demonstrate that radial expansion of the umbrella cloud is capable of driving ash upwind (westward) and crosswind (N-S) in excess of 1500 km, producing more-or-less radially symmetric isopachs that are only secondarily modified by ambient wind. Deposit thicknesses are decimeters to meters in the northern Rocky Mountains, centimeters to decimeters in the northern Midwest, and millimeters to centimeters on the East, West, and Gulf Coasts. Umbrella cloud growth may explain the extremely widespread dispersal of the ∼640 ka and 2.1 Ma Yellowstone tephra deposits in the eastern Pacific, northeastern California, southern California, and South Texas.

Geochemistry, Geophysics, Geosystems↗
USGS2014

Segmentation of plate coupling, fate of subduction fluids, and modes of arc magmatism in Cascadia, inferred from magnetotelluric resistivity

Five magnetotelluric (MT) profiles have been acquired across the Cascadia subduction system and transformed using 2-D and 3-D nonlinear inversion to yield electrical resistivity cross sections to depths of ∼200 km. Distinct changes in plate coupling, subduction fluid evolution, and modes of arc magmatism along the length of Cascadia are clearly expressed in the resistivity structure. Relatively high resistivities under the coasts of northern and southern Cascadia correlate with elevated degrees of inferred plate locking, and suggest fluid- and sediment-deficient conditions. In contrast, the north-central Oregon coastal structure is quite conductive from the plate interface to shallow depths offshore, correlating with poor plate locking and the possible presence of subducted sediments. Low-resistivity fluidized zones develop at slab depths of 35–40 km starting ∼100 km west of the arc on all profiles, and are interpreted to represent prograde metamorphic fluid release from the subducting slab. The fluids rise to forearc Moho levels, and sometimes shallower, as the arc is approached. The zones begin close to clusters of low-frequency earthquakes, suggesting fluid controls on the transition to steady sliding. Under the northern and southern Cascadia arc segments, low upper mantle resistivities are consistent with flux melting above the slab plus possible deep convective backarc upwelling toward the arc. In central Cascadia, extensional deformation is interpreted to segregate upper mantle melts leading to underplating and low resistivities at Moho to lower crustal levels below the arc and nearby backarc. The low- to high-temperature mantle wedge transition lies slightly trenchward of the arc.

Geochemistry, Geophysics, Geosystems↗
USGS2014

Cogenetic late Pleistocene rhyolite and cumulate diorites from Augustine Volcano revealed by SIMS 238U-230Th dating of zircon, and implications for silicic magma generation by extraction from mush

Augustine Volcano, a frequently active andesitic island stratocone, erupted a late Pleistocene rhyolite pumice fall that is temporally linked through zircon geochronology to cumulate dioritic blocks brought to the surface in Augustine's 2006 eruption. Zircon from the rhyolite yield a 238 U- 230 Th age of ∼25 ka for their unpolished rims, and their interiors yield a bimodal age populations at ∼26 ka and a minority at ∼41 ka. Zircon from dioritic blocks, ripped from Augustine's shallow magmatic plumbing system and ejected during the 2006 eruption, have interiors defining a ∼26 ka age population that is indistinguishable from that for the rhyolite; unpolished rims on the dioritic zircon are dominantly younger (≤12 ka) indicating subsequent crystallization. Zircon from rhyolite and diorite overlap in U, Hf, Ti, and REE concentrations although diorites also contain a second population of high-U, high temperature grains. Andesites that brought dioritic blocks to the surface in 2006 contain zircon with young (≤9 ka) rims and a scattering of older ages, but few zircon that crystallized during the 26 ka interval. Both the Pleistocene-age rhyolite and the 2006 dioritic inclusions plot along a whole-rock compositional trend distinct from mid-Holocene–present andesites and dacites, and the diorites, rhyolite, and two early Holocene dacites define linear unmixing trends often oblique to the main andesite array and consistent with melt (rhyolite) extraction from a mush (dacites), leaving behind a cumulate amphibole-bearing residue (diorites). Rare zircon antecrysts up to ∼300 ka from all rock types indicate that a Quaternary center has been present longer than preserved surficial deposits.

Alaska↗
USGS2015

The emergence of volcanic oceanic islands on a slow-moving plate: The example of Madeira Island, NE Atlantic

The transition from seamount to oceanic island typically involves surtseyan volcanism. However, the geological record at many islands in the NE Atlantic—all located within the slow-moving Nubian plate—does not exhibit evidence for an emergent surtseyan phase but rather an erosive unconformity between the submarine basement and the overlying subaerial shield sequences. This suggests that the transition between seamount and island may frequently occur by a relative fall of sea level through uplift, eustatic changes, or a combination of both, and may not involve summit volcanism. In this study, we explore the consequences for island evolutionary models using Madeira Island (Portugal) as a case study. We have examined the geologic record at Madeira using a combination of detailed fieldwork, biostratigraphy, and 40 Ar/ 39 Ar geochronology in order to document the mode, timing, and duration of edifice emergence above sea level. Our study confirms that Madeira's subaerial shield volcano was built upon the eroded remains of an uplifted seamount, with shallow marine sediments found between the two eruptive sequences and presently located at 320–430 m above sea level. This study reveals that Madeira emerged around 7.0–5.6 Ma essentially through an uplift process and before volcanic activity resumed to form the subaerial shield volcano. Basal intrusions are a likely uplift mechanism, and their emplacement is possibly enhanced by the slow motion of the Nubian plate relative to the source of partial melting. Alternating uplift and subsidence episodes suggest that island edifice growth may be governed by competing dominantly volcanic and dominantly intrusive processes.

Geochemistry, Geophysics, Geosystems↗
USGS2015

Episodic Holocene eruption of the Salton Buttes rhyolites, California, from paleomagnetic, U-Th, and Ar/Ar dating

In the Salton Trough, CA, five rhyolite domes form the Salton Buttes: Mullet Island, Obsidian Butte, Rock Hill, North and South Red Hill, from oldest to youngest. Results presented here include 40 Ar/ 39 Ar anorthoclase ages, 238 U- 230 Th zircon crystallization ages, and comparison of remanent paleomagnetic directions with the secular variation curve, which indicate that all domes are Holocene. 238 U- 230 Th zircon crystallization ages are more precise than but within uncertainty of 40 Ar/ 39 Ar anorthoclase ages, suggesting that zircon crystallization proceeded until shortly before eruption in all cases except one. Remanent paleomagnetic directions require three eruption periods: (1) Mullet Island, (2) Obsidian Butte, and (3) Rock Hill, North Red Hill, and South Red Hill. Borehole cuttings logs document up to two shallow tephra layers. North and South Red Hills likely erupted within 100 years of each other, with a combined 238 U- 230 Th zircon isochron age of: 2.83 ± 0.60 ka (2 sigma); paleomagnetic evidence suggests this age predates eruption by hundreds of years (1800 cal BP). Rock Hill erupted closely in time to these eruptions. The Obsidian Butte 238 U- 230 Th isochron age (2.86 ± 0.96 ka) is nearly identical to the combined Red Hill age, but its Virtual Geomagnetic Pole position suggests a slightly older age. The age of aphyric Mullet Island dome is the least well constrained: zircon crystals are resorbed and the paleomagnetic direction is most distinct; possible Mullet Island ages include ca. 2300, 5900, 6900, and 7700 cal BP. Our results constrain the duration of Salton Buttes volcanism to between ca. 5900 and 500 years.

Geochemistry, Geophysics, Geosystems↗
USGS2015

Effects of simple acid leaching of crushed and powdered geological materials on high-precision Pb isotope analyses

We present new results of simple acid leaching experiments on the Pb isotope composition of USGS standard reference material powders and on ocean island basalt whole rock splits and powders. Rock samples were leached with cold 6 N HCl in an ultrasonic bath, then on a hot plate, and washed with ultrapure H 2 O before sample digestion in HF‐HNO 3 and chromatographic purification of Pb. Lead isotope analyses were measured by Tl‐doped MC‐ICPMS. Intrasession and intersession analytical reproducibilities of repeated analyses of both synthetic Pb solutions and Pb from single digests of chemically processed natural samples were generally better than 100 ppm (2 SD). The comparison of leached and unleached samples shows that leaching consistently removes variable amounts of contaminants that differ in Pb isotopic composition for different starting materials. For repeated digests of a single sample, analyses of leached samples reproduce better than those of unleached ones, confirming that leaching effectively removes most of the heterogeneously distributed extraneous Pb. Nevertheless, the external reproducibility of leached samples is still up to an order of magnitude worse than that of Pb solution standards (∼100 ppm). More complex leaching methods employed by earlier studies yield Pb isotope ratios within error of those produced by our method and at similar levels of reproducibility, demonstrating that our simple leaching method is as effective as more complex leaching techniques. Therefore, any Pb isotope heterogeneity among multiple leached digests of samples in excess of the external reproducibility is attributed to inherent isotopic heterogeneity of the sample. The external precision of ∼100 ppm (2 SD) achieved for Pb isotope ratio determination by Tl‐doped MC‐ICPMS is thus sufficient for most rocks. The full advantage of the most precise Pb isotope analytical methods is only realized in cases where the natural isotopic heterogeneity among samples in a studied suite is substantially below 100 ppm.

Geochemistry, Geophysics, Geosystems↗
USGS2015

Pressure disequilibria induced by rapid valve closure in noble gas extraction lines

Pressure disequilibria during rapid valve closures can affect calculated molar quantities for a range of gas abundance measurements (e.g., K-Ar geochronology, (U-Th)/He geochronology, noble gas cosmogenic chronology). Modeling indicates this effect in a system with a 10 L reservoir reaches a bias of 1% before 1000 pipette aliquants have been removed from the system, and a bias of 10% before 10,000 aliquants. Herein we explore the causes and effects of this problem, which is the result of volume changes during valve closure. We also present a solution in the form of an electropneumatic pressure regulator that can precisely control valve motion. This solution reduces the effect to ∼0.3% even after 10,000 aliquants have been removed from a 10 L reservoir.

Geochemistry, Geophysics, Geosystems↗
USGS2015

Age of the Lava Creek supereruption and magma chamber assembly at Yellowstone based on 40 Ar/ 39 Ar and U-Pb dating of sanidine and zircon crystals

The last supereruption from the Yellowstone Plateau formed Yellowstone caldera and ejected the >1000 km 3 of rhyolite that composes the Lava Creek Tuff. Tephra from the Lava Creek eruption is a key Quaternary chronostratigraphic marker, in particular for dating the deposition of mid Pleistocene glacial and pluvial deposits in western North America. To resolve the timing of eruption and crystallization history for the Lava Creek magma, we performed (1) 40 Ar/ 39 Ar dating of single sanidine crystals to delimit eruption age and (2) ion microprobe U-Pb and trace-element analyses of the crystal faces and interiors of single zircons to date the interval of zircon crystallization and characterize magmatic evolution. Sanidines from the two informal members composing Lava Creek Tuff yield a preferred 40 Ar/ 39 Ar isochron date of 631.3 ± 4.3 ka. Crystal faces on zircons from both members yield a weighted mean 206 Pb/ 238 U date of 626.5 ± 5.8 ka, and have trace element concentrations that vary with the eruptive stratigraphy. Zircon interiors yield a mean 206 Pb/ 238 U date of 659.8 ± 5.5 ka, and reveal reverse and/or oscillatory zoning of trace element concentrations, with many crystals containing high U concentration cores that likely grew from highly evolved melt. The occurrence of distal Lava Creek tephra in stratigraphic sequences marking the Marine Isotope Stage 16–15 transition supports the apparent eruption age of ∼631 ka. The combined results reveal that Lava Creek zircons record episodic heating, renewed crystallization, and an overall up-temperature evolution for Yellowstone's subvolcanic reservoir in the 10 3 −10 4 year interval before eruption.

Wyoming↗
USGS2015

Late Pleistocene ages for the most recent volcanism and glacial-pluvial deposits at Big Pine volcanic field, California, USA, from cosmogenic 36 Cl dating

The Big Pine volcanic field is one of several Quaternary volcanic fields that poses a potential volcanic hazard along the tectonically active Owens Valley of east-central California, and whose lavas are interbedded with deposits from Pleistocene glaciations in the Sierra Nevada Range. Previous geochronology indicates an ∼1.2 Ma history of volcanism, but the eruption ages and distribution of volcanic products associated with the most-recent eruptions have been poorly resolved. To delimit the timing and products of the youngest volcanism, we combine field mapping and cosmogenic 36 Cl dating of basaltic lava flows in the area where lavas with youthful morphology and well-preserved flow structures are concentrated. Field mapping and petrology reveal approximately 15 vents and 6 principal flow units with variable geochemical composition and mineralogy. Cosmogenic 36 Cl exposure ages for lava flow units from the top, middle, and bottom of the volcanic stratigraphy indicate eruptions at ∼17, 27, and 40 ka, revealing several different and previously unrecognized episodes of late Pleistocene volcanism. Olivine to plagioclase-pyroxene phyric basalt erupted from several vents during the most recent episode of volcanism at ∼17 ka, and produced a lava flow field covering ∼35 km 2 . The late Pleistocene 36 Cl exposure ages indicate that moraine and pluvial shoreline deposits that overlie or modify the youngest Big Pine lavas reflect Tioga stage glaciation in the Sierra Nevada and the shore of paleo-Owens Lake during the last glacial cycle.

Geochemistry, Geophysics, Geosystems↗
USGS2015

Tremor-genic slow slip regions may be deeper and warmer and may slip slower than non-tremor-genic regions

Slow slip events (SSEs) are observed worldwide and often coincide with tectonic tremor. Notable examples of SSEs lacking observed tectonic tremor, however, occur beneath Kīlauea Volcano, Hawaii, the Boso Peninsula, Japan, near San Juan Bautista on the San Andreas Fault, California, and recently in Central Ecuador. These SSEs are similar to other worldwide SSEs in many ways (e.g., size or duration), but lack the concurrent tectonic tremor observed elsewhere; instead, they trigger swarms of regular earthquakes. We investigate the physical conditions that may distinguish these non-tremor-genic SSEs from those associated with tectonic tremor, including slip velocity, pressure, temperature, fluids, and fault asperities, although we cannot eliminate the possibility that tectonic tremor may be obscured in highly attenuating regions. Slip velocities of SSEs at Kīlauea Volcano (∼10 −6 m/s) and Boso Peninsula (∼10 −7 m/s) are among the fastest SSEs worldwide. Kīlauea Volcano, the Boso Peninsula, and Central Ecuador are also among the shallowest SSEs worldwide, and thus have lower confining pressures and cooler temperatures in their respective slow slip zones. Fluids also likely contribute to tremor generation, and no corresponding zone of high v p /v s has been noted at Kīlauea or Boso. We suggest that the relatively faster slip velocities at Kīlauea Volcano and the Boso Peninsula result from specific physical conditions that may also be responsible for triggering swarms of regular earthquakes adjacent to the slow slip, while different conditions produce slower SSE velocities elsewhere and trigger tectonic tremor.

Alaska, California, Hawaii↗
USGS2015

A random-walk algorithm for modeling lithospheric density and the role of body forces in the evolution of the Midcontinent Rift

This paper develops a Monte Carlo algorithm for extracting three-dimensional lithospheric density models from geophysical data. Empirical scaling relationships between velocity and density create a 3D starting density model, which is then iteratively refined until it reproduces observed gravity and topography. This approach permits deviations from uniform crustal velocity-density scaling, which provide insight into crustal lithology and prevent spurious mapping of crustal anomalies into the mantle. We test this algorithm on the Proterozoic Midcontinent Rift (MCR), north-central U.S. The MCR provides a challenge because it hosts a gravity high overlying low shear-wave velocity crust in a generally flat region. Our initial density estimates are derived from a seismic velocity/crustal thickness model based on joint inversion of surface-wave dispersion and receiver functions. By adjusting these estimates to reproduce gravity and topography, we generate a lithospheric-scale model that reveals dense middle crust and eclogitized lowermost crust within the rift. Mantle lithospheric density beneath the MCR is not anomalous, consistent with geochemical evidence that lithospheric mantle was not the primary source of rift-related magmas and suggesting that extension occurred in response to far-field stress rather than a hot mantle plume. Similarly, the subsequent inversion of normal faults resulted from changing far-field stress that exploited not only warm, recently faulted crust but also a gravitational potential energy low in the MCR. The success of this density modeling algorithm in the face of such apparently contradictory geophysical properties suggests that it may be applicable to a variety of tectonic and geodynamic problems.

Iowa, Kansas, Michigan, Minnesota, Nebraska, Wisco↗