Search USGS⌕ Search

USGS · 70015903

Trace element and isotopic constraints on magmatic evolution at Lassen volcanic center

Abstract

Magmatic evolution at the Lassen volcanic center (LVC) is characterized by a transition from predominantly andesitic to predominantly silicic volcanism with time. Magmas of the andesitic, or “Brokeoff phase” of volcanism range in composition from basaltic andesite io dacite, whereas those of silicic, or “Lassen phase” range in composition from basaltic andesite to rhyolite. The compositions of magmas from each phase define well organized but distinct variation trends. Compared with Brokeoff-phase magmas of similar SiO 2 content, most Lassen-phase magmas contain lower concentrations of most incompatible minor and trace elements. Based on the behavior of both incompatible and compatible trace elements, the geochemical trends defined by the Brokeoff-phase magmas cannot be ascribed merely to fractional crystallization from a single or multiple mafic parental magmas, Moreover, the Lassen-phase magmas cannot be derived from the Brokeoff-phase magmas by fractional crystallization. Rather, the geochemical trends that characterize each volcanic phase define arrays that primarily indicate mixing between well-homogenized silicic and heterogeneous mafic magmas. The distinctive mixing-dominated arrays for each volcanic phase manifest the generation and evolution of two physically distinct, but genetically related magma systems. The LVC magmas have Sr, Nd, and Pb isotopic characteristics that approximate two-component mixing arrays. One isotopic component is similar in composition to that of NE Pacific Ocean ridge and seamount basalts (“MORB” component), the other to mafic Mesozoic granitoids sampled from the neighboring KSamath and Sierra Nevada provinces (“KSN” component). The isotopic compositions of the most silicic LVC magmas lie within the ranges defined by the more mafic LVC magmas, which in turn lie within broad ranges defined by primitive mafic lavas sampled from the Lassen region. The lack of a correlation between the major element and isotopic compositions of LVC magmas seriously limits any model for magmatic evolution that relies on assimilation of old middle to upper crust by isotopically homogeneous mafic magmas during their ascent through the crust. Alternatively, the isotopic and geochemical uniformity of the most silicic magmas of the Brokeoff and Lassen phases suggests that they are well-homogenized partial melts. The likely source region for these silicic melts is the lower crust, which we envision to consist primarily of mafic igneous rocks that are similar in geochemical and isotopic diversity to the regional mafic lavas. Magmatic evolution at LVC can be viewed in terms of a series of mantle melting events that subsequently stimulated meiting in a progressively increasing volume of the lower crust. In general, the LVC magmas represent slightly fractionated mixtures of the mantle-derived mafic magmas and silicic partial melts of the lower crust, the latter melts increasing in relative proportion over the history of the volcanic center. The voluminous rhyolitic lavas and pyroclastic materials erupted during the early Lassen phase represent lower crustal melts thai pooled into sufficient volumes to avoid significant blending and dilution with intruding mafic magmas. The geochemical and isotopic heterogeneity of both Brokeoff-phase andésites and mafic magmatic inclusions in the silicic Lassen-phase magmas must in part record the corresponding diversity of the mantle-derived contributions to arc magmagenesis in this region.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

T.D. Bullen, M.A. Clynne. 2012-09-20. Trace element and isotopic constraints on magmatic evolution at Lassen volcanic center. https://doi.org/10.1029/jb095ib12p19671

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗

Exploring the dynamic interactions between the Southern San Andreas Fault and a normal fault under the Salton Sea

We investigate the dynamic interactions between the Southern San Andreas Fault (SSAF) and a proximal normal fault (NF) beneath the Salton Sea in southern California. The NF, positioned near the SSAF terminus at Bombay Beach, exhibits 11–15 displacement events across 14 stratigraphic sequences, with a range of 0.2–1.4 m of vertical offset since ∼2–3 ka. Notably, four of these events may align temporally with SSAF earthquakes, raising questions about the possible interplay between the two faults. Utilizing dynamic rupture models, we analyze the coseismic interactions between the SSAF and NF, addressing under what conditions the SSAF induces slip on the NF. Our findings reveal that a suite of SSAF ruptures, particularly those propagating from north to south, can trigger slip on the normal fault and replicate observed vertical offsets. If the SSAF extends beneath the Salton Sea, earthquakes originating south of the NF intersection are less likely to trigger normal fault slip, although we cannot exclude this possibility. Some SSAF ruptures do not trigger discernible slip on the NF, rendering such events undetectable in the stratigraphic record. Our research contributes toward discussions regarding the seismic hazard in southern California, shedding light on the interplay between the SSAF and NF.

California↗

Quantifying site effects and their influence on earthquake source parameter estimations using a dense array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

Oklahoma↗