SEARCH · USGS Science

Results for “The Earth Observer”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,225 records · Page 68Linked to original sources

USGS2002

Aftershock triggering by complete Coulomb stress changes

We examine the correlation between seismicity rate change following the 1992, M 7.3, Landers, California, earthquake and characteristics of the complete Coulomb failure stress (CFS) changes (ΔCFS( t )) that this earthquake generated. At close distances the time-varying “dynamic” portion of the stress change depends on how the rupture develops temporally and spatially and arises from radiated seismic waves and from permanent coseismic fault displacement. The permanent “static” portion (ΔCFS) depends only on the final coseismic displacement. ΔCFS diminishes much more rapidly with distance than the transient, dynamic stress changes. A common interpretation of the strong correlation between ΔCFS and aftershocks is that load changes can advance or delay failure. Stress changes may also promote failure by physically altering properties of the fault or its environs. Because it is transient, ΔCFS( t ) can alter the failure rate only by the latter means. We calculate both ΔCFS and the maximum positive value of ΔCFS( t ) (peak ΔCFS( t )) using a reflectivity program. Input parameters are constrained by modeling Landers displacement seismograms. We quantify the correlation between maps of seismicity rate changes and maps of modeled ΔCFS and peak ΔCFS( t ) and find agreement for both models. However, rupture directivity, which does not affect ΔCFS, creates larger peak ΔCFS( t ) values northwest of the main shock. This asymmetry is also observed in seismicity rate changes but not in ΔCFS. This result implies that dynamic stress changes are as effective as static stress changes in triggering aftershocks and may trigger earthquakes long after the waves have passed.

Journal of Geophysical Research B: Solid Earth↗

USGS2004

Stress interaction between subduction earthquakes and forearc strike-slip faults: Modeling and application to the northern Caribbean plate boundary

Strike-slip faults in the forearc region of a subduction zone often present significant seismic hazard because of their proximity to population centers. We explore the interaction between thrust events on the subduction interface and strike-slip faults within the forearc region using three-dimensional models of static Coulomb stress change. Model results reveal that subduction earthquakes with slip vectors subparallel to the trench axis enhance the Coulomb stress on strike-slip faults adjacent to the trench but reduce the stress on faults farther back in the forearc region. In contrast, subduction events with slip vectors perpendicular to the trench axis enhance the Coulomb stress on strike-slip faults farther back in the forearc, while reducing the stress adjacent to the trench. A significant contribution to Coulomb stress increase on strike-slip faults in the back region of the forearc comes from "unclamping" of the fault, i.e., reduction in normal stress due to thrust motion on the subduction interface. We argue that although Coulomb stress changes from individual subduction earthquakes are ephemeral, their cumulative effects on the pattern of lithosphere deformation in the forearc region are significant. We use the Coulomb stress models to explain the contrasting deformation pattern between two adjacent segments of the Caribbean subduction zone. Subduction earthquakes with slip vectors nearly perpendicular to the Caribbean trench axis is dominant in the Hispaniola segment, where the strike-slip faults are more than 60 km inland from the trench. In contrast, subduction slip motion is nearly parallel to the Caribbean trench axis along the Puerto Rico segment, where the strike-slip fault is less than 15 km from the trench. This observed jump from a strike-slip fault close to the trench axis in the Puerto Rico segment to the inland faults in Hispaniola is explained by different distributions of Coulomb stress in the forearc region of the two segments, as a result of the change from the nearly trench parallel slip on the Puerto Rico subduction interface to the more perpendicular subduction slip beneath Hispaniola. The observations and modeling suggest that subduction-induced strike-slip seismic hazard to Puerto Rico may be smaller than previously assumed but the hazard to Hispaniola remains high. Copyright 2004 by the American Geophysical Union.

Puerto Rico↗

USGS2007

Sorted bed forms as self-organized patterns: 2. complex forcing scenarios

We employ a numerical model to study the development of sorted bed forms under a variety of hydrodynamic and sedimentary conditions. Results indicate that increased variability in wave height decreases the growth rate of the features and can potentially give rise to complicated, a priori unpredictable, behavior. This happens because the system responds to a change in wave characteristics by attempting to self-organize into a patterned seabed of different geometry and spacing. The new wavelength might not have enough time to emerge before a new change in wave characteristics occurs, leading to less regular seabed configurations. The new seabed configuration is also highly dependent on the preexisting morphology, which further limits the possibility of predicting future behavior. For the same reasons, variability in the mean current magnitude and direction slows down the growth of features and causes patterns to develop that differ from classical sorted bed forms. Spatial variability in grain size distribution and different types of net sediment aggradation/degradation can also result in the development of sorted bed forms characterized by a less regular shape. Numerical simulations qualitatively agree with observed geometry (spacing and height) of sorted bed forms. Also in agreement with observations is that at shallower depths, sorted bed forms are more likely to be affected by changes in the forcing conditions, which might also explain why, in shallow waters, sorted bed forms are described as ephemeral features. Finally, simulations indicate that the different sorted bed form shapes and patterns observed in the field might not necessarily be related to diverse physical mechanisms. Instead, variations in sorted bed form characteristics may result from variations in local hydrodynamic and/or sedimentary conditions.

Journal of Geophysical Research F: Earth Surface↗

USGS1997

Permeability and porosity of the Illinois UPH 3 drillhole granite and a comparison with other deep drillhole rocks

Permeability, porosity, and volumetric strain measurements were conducted on granite cores obtained at depths of 0.7 to 1.6 km from the Illinois UPH 3 drillhole at effective confining pressures from 5 to 100 MPa. Initial permeabilities were in the range of 10-17 to 10-19 m2 and dropped rapidly with applied pressure to values between 10-20 and 10-24 m2 at 100 MPa, typical of other deep granite core samples. These values are several decades lower than equivalent weathered surface granites at comparable effective confining pressures, where weathering products in cracks and pores inhibit crack closure with applied pressure. Permeabilities of the Illinois cores were inversely related to sample depth, suggesting that stress relief and thermal microfractures induced during core retrieval dominated the fluid flow. Thus these samples provide an upper bound on in situ matrix permeability values. A comparison of core permeability from UPH 3 and other deep drillholes shows that stress relief damage can often dominate laboratory permeability measurements. We conclude that it may be difficult to make meaningful estimates of in situ permeability based on either borehole samples (possible damage during retrieval) or surface-derived analogs (altered by weathering). Volumetric strain determined from porosity measurements was compared with differential strain analysis (DSA) data reported by other investigators on samples from the same depths in the drillhole. Our strain measurements (0.002 to 0.005 at 100 MPa) were nearly twice as large as the DSA values, probably because of the crack-enhancing effects of fluids present in our samples that are absent in the dry DSA cores, as well as other time-dependent deformation effects. This difference in observed strain magnitudes between the two measurement methods may be an important consideration if strain and/or porosity data from deep core samples are used in models of stress, fluid circulation, and excess fluid pressure generation in the midcrust.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Acceleration spectra for subduction zone earthquakes

We estimate the source spectra of shallow earthquakes from digital recordings of teleseismic P wave groups, that is, P + p P + s P , by making frequency dependent corrections for the attenuation and for the interference of the free surface. The correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths. We apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from M S = 6.2 to 8.1, obtaining corrected P wave acceleration spectra on the frequency band from 0.01 to 2.0 Hz. Seismic moment estimates from surface waves and normal modes are used to extend these P wave spectra to the frequency band from 0.001 to 0.01 Hz. The acceleration spectra of moderate subduction zone earthquakes, that is, earthquakes whose seismic moments are less than 10 27 dyn cm, exhibit ω-square or Brune-type spectra, while the acceleration spectra of large subduction zone earthquakes, that is, earthquakes whose seismic moments are greater than 10 27 dyn cm, exhibit intermediate slopes where ü(ω) ∝ ω 5/4 for frequencies from 0.005 to 0.05 Hz. For this set of earthquakes, spectral shape appears to be a discontinuous function of seismic moment. Using reasonable assumptions for the phase characteristics, we transform the spectral shape observed for large earthquakes into the time domain to fit Ekström's (1987) moment rate functions for the M S = 8.1 Michoacan earthquake of September 19, 1985, and the M S = 7.6 Michoacan aftershock of September 21, 1985.

Journal of Geophysical Research Solid Earth↗

USGS1996

The Border Ranges fault system in Glacier Bay National Park, Alaska: Evidence for major early Cenozoic dextral strike-slip motion

The Border Ranges fault system of southern Alaska, the fundamental break between the arc basement and the forearc accretionary complex, is the boundary between the Peninsular–Alexander–Wrangellia terrane and the Chugach terrane. The fault system separates crystalline rocks of the Alexander terrane from metamorphic rocks of the Chugach terrane in Glacier Bay National Park. Mylonitic rocks in the zone record abundant evidence for dextral strike-slip motion along north-northwest-striking subvertical surfaces. Geochronologic data together with regional correlations of Chugach terrane rocks involved in the deformation constrain this movement between latest Cretaceous and Early Eocene (~50 Ma). These findings are in agreement with studies to the northwest and southeast along the Border Ranges fault system which show dextral strike-slip motion occurring between 58 and 50 Ma. Correlations between Glacier Bay plutons and rocks of similar ages elsewhere along the Border Ranges fault system suggest that as much as 700 km of dextral motion may have been accommodated by this structure. These observations are consistent with oblique convergence of the Kula plate during early Cenozoic and forearc slivering above an ancient subduction zone following late Mesozoic accretion of the Peninsular–Alexander–Wrangellia terrane to North America.

Alaska↗

USGS2019

Hydroacoustic, meteorologic and seismic observations of the 2016 Nansen Ice Shelf calving event and iceberg formation

On 7 April 2016, the Nansen Ice Shelf (NIS) front calved into two icebergs, the first large-scale calving event in >30 years. Three hydrophone moorings were deployed seaward of the NIS in December 2015 and over the following months recorded hundreds of short duration, broadband (10-400 Hz) cryogenic signals, likely caused by fracturing of the ice-shelf. The majority of these icequakes occur between January and early March 2016, several weeks prior to the calving observed by satellite on 7 April. Barometric pressure and wind speed records show the day the icebergs drifted from the NIS coincided with the largest low-pressure storm system recorded in the previous 7 months. A nearby seismic station also shows an increase in low-frequency energy, harmonic tremor, and microseisms on 7 April. Our interpretation is the northern segment of the NIS leading edge broke free during mid-January to February, producing high acoustic energy, but the icebergs remained stationary until the combination of a strong low-pressure system, with high winds freed the icebergs. As the unpinning of Antarctic ice shelves is not a well-documented process, our observations show that storm systems may play an under-appreciated role in Antarctic ice shelf break-up.

Frontiers in Earth Science↗

USGS2007

Statistical analysis of water-quality data containing multiple detection limits II: S-language software for nonparametric distribution modeling and hypothesis testing

Analysis of low concentrations of trace contaminants in environmental media often results in left-censored data that are below some limit of analytical precision. Interpretation of values becomes complicated when there are multiple detection limits in the data-perhaps as a result of changing analytical precision over time. Parametric and semi-parametric methods, such as maximum likelihood estimation and robust regression on order statistics, can be employed to model distributions of multiply censored data and provide estimates of summary statistics. However, these methods are based on assumptions about the underlying distribution of data. Nonparametric methods provide an alternative that does not require such assumptions. A standard nonparametric method for estimating summary statistics of multiply-censored data is the Kaplan-Meier (K-M) method. This method has seen widespread usage in the medical sciences within a general framework termed "survival analysis" where it is employed with right-censored time-to-failure data. However, K-M methods are equally valid for the left-censored data common in the geosciences. Our S-language software provides an analytical framework based on K-M methods that is tailored to the needs of the earth and environmental sciences community. This includes routines for the generation of empirical cumulative distribution functions, prediction or exceedance probabilities, and related confidence limits computation. Additionally, our software contains K-M-based routines for nonparametric hypothesis testing among an unlimited number of grouping variables. A primary characteristic of K-M methods is that they do not perform extrapolation and interpolation. Thus, these routines cannot be used to model statistics beyond the observed data range or when linear interpolation is desired. For such applications, the aforementioned parametric and semi-parametric methods must be used.

Computers & Geosciences↗

USGS2022

The consequences of climate change for dryland biogeochemistry

Drylands, which cover more than 40% of Earth’s terrestrial surface, are dominant drivers of global biogeochemical cycling and home to more than one third of the human population. Climate projections predict warming, drought frequency and severity, and evaporative demand will increase in drylands at faster rates than global means. Due to extreme temperatures and high biological dependency on limited water availability, drylands are predicted to be exceptionally sensitive to climate change and, indeed, significant climate impacts are already being observed. Yet our understanding and ability to forecast climate change effects on dryland biogeochemistry and ecosystem functions lag behind many mesic systems. To improve our capacity to forecast ecosystem change, we propose focusing on the controls and consequences of two key characteristics affecting dryland biogeochemistry: i) high spatial and temporal heterogeneity in environmental conditions and ii) generalized resource scarcity. In addition to climate change, drylands are experiencing accelerating land use change. Building our understanding of dryland biogeochemistry in both intact and disturbed systems will better equip us to address the interacting effects of climate change and landscape degradation. Responding to these challenges will require a diverse, globally distributed, and interdisciplinary community of dryland experts united towards better understanding these vast and important ecosystems.

New Phytologist Foundation↗

USGS2010

Evolution of Ore Deposits and Technology Transfer Project: Isotope and Chemical Methods in Support of the U.S. Geological Survey Science Strategy, 2003-2008

Principal functions of the U.S. Geological Survey (USGS) Mineral Resources Program are providing assessments of the location, quantity, and quality of undiscovered mineral deposits, and predicting the environmental impacts of exploration and mine development. The mineral and environmental assessments of domestic deposits are used by planners and decisionmakers to improve the stewardship of public lands and public resources. Assessments of undiscovered mineral deposits on a global scale reveal the potential availability of minerals to the United States and other countries that manufacture goods imported to the United States. These resources are of fundamental relevance to national and international economic and security policy in our globalized world economy. Performing mineral and environmental assessments requires that predictions be made of the likelihood of undiscovered deposits. The predictions are based on geologic and geoenvironmental models that are constructed for the diverse types of mineral deposits from detailed descriptions of actual deposits and detailed understanding of the processes that formed them. Over the past three decades the understanding of ore-forming processes has benefited greatly from the integration of laboratory-based geochemical tools with field observations and other data sources. Under the aegis of the Evolution of Ore Deposits and Technology Transfer Project (referred to hereinafter as the Project), a 5-year effort that terminated in 2008, the Mineral Resources Program provided state-of-the-art analytical capabilities to support applications of several related geochemical tools to ore-deposit-related studies. The analytical capabilities and scientific approaches developed within the Project have wide applicability within Earth-system science. For this reason the Project Laboratories represent a valuable catalyst for interdisciplinary collaborations of the type that should be formed in the coming years for the United States to meet its natural-resources and natural-science needs. This circular presents an overview of the Project. Descriptions of the Project laboratories are given first including descriptions of the types of chemical or isotopic analyses that are made and the utility of the measurements. This is followed by summaries of select measurements that were carried out by the Project scientists. The studies are grouped by science direction. Virtually all of them were collaborations with USGS colleagues or with scientists from other governmental agencies, academia, or the private sector.

Circular↗

USGS1999

P wave crustal velocity structure in the greater Mount Rainier area from local earthquake tomography

We present results from a local earthquake tomographic imaging experiment in the greater Mount Rainier area. We inverted P wave arrival times from local earthquakes recorded at permanent and temporary Pacific Northwest Seismograph Network seismographs between 1980 and 1996. We used a method similar to that described by Lees and Crosson [1989], modified to incorporate the parameter separation method for decoupling the hypocenter and velocity problems. In the upper 7 km of the resulting model there is good correlation between velocity anomalies and surface geology. Many focal mechanisms within the St. Helens seismic zone have nodal planes parallel to the epicentral trend as well as to a north-south trending low-velocity trough, leading us to speculate that the trough represents a zone of structural weakness in which a moderate ( M 6.5–7.0) earthquake could occur. In contrast, the western Rainier seismic zone does not correlate in any simple way with anomaly patterns or focal mechanism fault planes, leading us to infer that it is less likely to experience a moderate earthquake. A ∼10 km-wide low-velocity anomaly occurs 5 to 18 km beneath the summit of Mount Rainier, which we interpret to be a signal of a region composed of hot, fractured rock with possible small amounts of melt or fluid. No systematic velocity pattern is observed in association with the southern Washington Cascades conductor. A midcrustal anomaly parallels the Olympic-Wallowa lineament as well as several other geophysical trends, indicating that it may play an important role in regional tectonics.

Journal of Geophysical Research B: Solid Earth↗

USGS2021

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗

USGS1998

Episodic fluid flow in the Nankai accretionary complex: Timescale, geochemistry, flow rates, and fluid budget

Down-hole geochemical anomalies encountered in active accretionary systems can be used to constrain the timing, rates, and localization of fluid flow. Here we combine a coupled flow and solute transport model with a kinetic model for smectite dehydration to better understand and quantify fluid flow in the Nankai accretionary complex offshore of Japan. Compaction of sediments and clay dehydration provide fluid sources which drive the model flow system. We explicitly include the consolidation rate of underthrust sediments in our calculations to evaluate the impact that variations in this unknown quantity have on pressure and chloride distribution. Sensitivity analysis of steady state pressure solutions constrains bulk and flow conduit permeabilities. Steady state simulations with 30% smectite in the incoming sedimentary sequence result in minimum chloride concentrations at site 808 of 550 m M , but measured chlorinity is as low as 447 m M . We simulate the transient effects of hydrofracture or a strain event by assuming an instantaneous permeability increase of 3–4 orders of magnitude along a flow conduit (in this case the décollement), using steady state results as initial conditions. Transient results with an increase in décollement permeability from 10 −16 m 2 to 10 −13 m 2 and 20% smectite reproduce the observed chloride profile at site 808 after 80–160 kyr. Modeled chloride concentrations are highly sensitive to the consolidation rate of underthrust sediments, such that rapid compaction of underthrust material leads to increased freshening. Pressures within the décollement during transient simulations rise rapidly to a significant fraction of lithostatic and remain high for at least 160 kyr, providing a mechanism for maintaining high permeability. Flow rates at the deformation front for transient simulations are in good agreement with direct measurements, but steady state flow rates are 2–3 orders of magnitude smaller than observed. Fluid budget calculations indicate that nearly 71% of the incoming water in the sediments leaves the accretionary wedge via diffuse flow out the seafloor, 0–5% escapes by focused flow along the décollement, and roughly 1% is subducted.

Journal of Geophysical Research B: Solid Earth↗

USGS1998

Trace element abundances of high-MgO glasses from Kilauea, Mauna Loa and Haleakala volcanoes, Hawaii

We performed an ion-microprobe study of eleven high-MgO (6.7-14.8 wt%) tholeiite glasses from the Hawaiian volcanoes Kilauea, Mauna Loa and Haleakala. We determined the rare earth (RE), high field strength, and other selected trace element abundances of these glasses, and used the data to establish their relationship to typical Hawaiian shield tholeiite and to infer characteristics of their source. The glasses have trace element abundance characteristics generally similar to those of typical shield tholeiites, e.g. L(light)REE/H(heavy)REE(C1) > 1. The Kilauea and Mauna Loa glasses, however, display trace and major element characteristics that cross geochemical discriminants observed between Kilauea and Mauna Loa shield lavas. The glasses contain a blend of these discriminating chemical characteristics, and are not exactly like the typical shield lavas from either volcano. The production of these hybrid magmas likely requires a complexly zoned source, rather than two unique sources. When corrected for olivine fractionation, the glass data show correlations between CaO concentration and incompatible trace element abundances, indicating that CaO may behave incompatibly during melting of the tholeiite source. Furthermore, the tholeiite source must contain residual garnet and clinopyroxene to account for the variation in trace element abundances of the Kilauea glasses. Inversion modeling indicates that the Kilauea source is flat relative to C1 chondrites, and has a higher bulk distribution coefficient for the HREE than the LREE.

Hawaii↗

USGS1987

Nucleation and triggering of earthquake slip: Effect of periodic stresses

Results of stability analyses for spring and slider systems, with state variable constitutive properties, are applied to slip on embedded fault patches. Unstable slip may nucleate only if the slipping patch exceeds some minimum size. Subsequent to the onset of instability the earthquake slip may propagate well beyond the patch. It is proposed that the seismicity of a volume of the earth's crust is determined by the distribution of initial conditions on the population of fault patches that nucleate earthquake slip, and the loading history acting upon the volume. Patches with constitutive properties inferred from laboratory experiments are characterized by an interval of self-driven accelerating slip prior to instability, if initial stress exceeds a minimum threshold. This delayed instability of the patches provides an explanation for the occurrence of aftershocks and foreshocks including decay of earthquake rates by time −1 . A population of patches subjected to loading with a periodic component results in periodic variation of the rate of occurrence of instabilities. The change of the rate of seismicity for a sinusoidal load is proportional to the amplitude of the periodic stress component and inversely proportional to both the normal stress acting on the fault patches and the constitutive parameter, A 1 , that controls the direct velocity dependence of fault slip. Values of A 1 representative of laboratory experiments indicate that in a homogeneous crust, correlation of earthquake rates with earth tides should not be detectable at normal stresses in excess of about 8 MPa. Correlation of earthquakes with tides at higher normal stresses can be explained if there exist inhomogeneities that locally amplify the magnitude of the tidal stresses. Such amplification might occur near magma chambers or other soft inclusions in the crust and possibly near the ends of creeping fault segments if the creep or afterslip rates vary in response to tides. Observations of seismicity rate variations associated with seasonal fluctuations of reservoir levels appear to be consistent with the model.

Tectonophysics↗

USGS2012

Ash3d: A finite-volume, conservative numerical model for ash transport and tephra deposition

We develop a transient, 3-D Eulerian model (Ash3d) to predict airborne volcanic ash concentration and tephra deposition during volcanic eruptions. This model simulates downwind advection, turbulent diffusion, and settling of ash injected into the atmosphere by a volcanic eruption column. Ash advection is calculated using time-varying pre-existing wind data and a robust, high-order, finite-volume method. Our routine is mass-conservative and uses the coordinate system of the wind data, either a Cartesian system local to the volcano or a global spherical system for the Earth. Volcanic ash is specified with an arbitrary number of grain sizes, which affects the fall velocity, distribution and duration of transport. Above the source volcano, the vertical mass distribution with elevation is calculated using a Suzuki distribution for a given plume height, eruptive volume, and eruption duration. Multiple eruptions separated in time may be included in a single simulation. We test the model using analytical solutions for transport. Comparisons of the predicted and observed ash distributions for the 18 August 1992 eruption of Mt. Spurr in Alaska demonstrate to the efficacy and efficiency of the routine.

Journal of Geophysical Research↗

USGS2021

Three-dimensional electrical resistivity characterization of Mountain Pass, California and surrounding region

The Sulphide Queen carbonatite deposit at Mountain Pass in southeast California is a world class rare earth element (REE) resource. This study images electrical resistivity structure of the REE deposit and surrounding area to characterize resources under cover. An east-west elongated grid (35 × 15 km) of 65 wideband magnetotelluric stations spanning from eastern Shadow Valley to eastern Ivanpah Valley were collected and modeled in three-dimensions (3-D). Gravity, aeromagnetic, and geologic data are used to inform interpretation of structures in the resistivity model, including the following observations. Shadow Valley is filled with conductive sediment that locally dips southward to a depth of 1 km. The Kingston Range-Halloran Hills detachment fault dips westward at ∼15 degrees. The REE deposit is a moderate low resistivity zone dipping southwest to a possible depth of ∼1 km, and is bounded by the North and South faults and bisected by the Middle fault. Ivanpah Dry Lake is underlain by a north striking southward dipping sedimentary basin. Two possible zones of mineralization are observed in Ivanpah Valley, one along the western edge of Ivanpah Dry Lake and one on the western edge of valley along a new inferred fault. The brittle-ductile transition is imaged at ∼10 km below mean sea level. No deep electrically conductive structures are imaged to be related to the REE deposit likely due to the complex geologic history of the Mojave terrane. Future studies should regional target Proterozoic rocks and search within for geophysical signatures similar to Mountain Pass.

California↗

USGS2017

Initial pulse of Siberian Traps sills as the trigger of the end-Permian mass extinction

Mass extinction events are short-lived and characterized by catastrophic biosphere collapse and subsequent reorganization. Their abrupt nature necessitates a similarly short-lived trigger, and large igneous province magmatism is often implicated. However, large igneous provinces are long-lived compared to mass extinctions. Therefore, if large igneous provinces are an effective trigger, a subinterval of magmatism must be responsible for driving deleterious environmental effects. The onset of Earth’s most severe extinction, the end-Permian, coincided with an abrupt change in the emplacement style of the contemporaneous Siberian Traps large igneous province, from dominantly flood lavas to sill intrusions. Here we identify the initial emplacement pulse of laterally extensive sills as the critical deadly interval. Heat from these sills exposed untapped volatile-fertile sediments to contact metamorphism, likely liberating the massive greenhouse gas volumes needed to drive extinction. These observations suggest that large igneous provinces characterized by sill complexes are more likely to trigger catastrophic global environmental change than their flood basalt- and/or dike-dominated counterparts.

Siberia, Siberian Traps↗