SEARCH · Search USGS

Results for “Journal of Geophysical Research--Solid Earth”

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 289 records · Page 16Linked to original sources

USGS2012

Reflected and mode-converted seismic waves within the shallow Aleutian subduction zone, southern Kenai Peninsula, Alaska

Pronounced secondary phases observed in local recordings of quarry shots and earthquakes on the southern Kenai Peninsula, Alaska, are identified as reflected P and S and converted S-to-P phases originating within four depth ranges: in the upper few kilometers of the Cook Inlet Tertiary basin, at midcrustal depths within the overthrust North American plate, at about 35 km depth near the top of the Wadati-Benioff seismic zone in proximity to the inferred interplate megathrust, and at about 5–10 km below the megathrust in the subducted Pacific plate. Beneath the landward margin of the accreted Chugach terrane, the mid-upper plate (MUP) discontinuity dips as steeply as 20°–30° to the west-northwest over the depth range 12–18 km. At shallower depths it flattens and possibly arches with the crest at about 10 km depth. Similar midcrustal reflectors have been imaged about 125 km to the southwest and about 350 km to the northeast along the structural trend of the Chugach terrane. The extensive reflectors may have a common origin, possible caused by the presence of underplated rocks. Relative amplitudes and polarities of the secondary phases originating at MUP depths provide few constraints on the nature of the discontinuity. It appears that the MUP discontinuity is seismically inactive and does not represent a brittle-ductile transition zone within the upper plate. The two converted S -to- P phases generated near the top of the subducted plate could indicate a low-velocity zone associated with subducted oceanic crust.

Journal of Geophysical Research Solid Earth↗

USGS2012

Slip partitioning of the Calaveras Fault, California, and prospects for future earthquakes

Examination of main shock and microearthquake data from the Calaveras Fault during the last 20 years reveals that main shock hypocenters occur at depths of 8-9 km near the base of the zone of microearthquakes. Microseismicity extends between depths of 4 and 10 km and defines zones of concentrated microseismicity and aseismic zones. Estimates of the fault regions which slipped during the Coyote Lake and Morgan Hill earthquakes as derived from seismic radiation coincide with zones which are otherwise aseismic. We propose that these persistent aseismic zones represent stuck patches which slip only during moderate earthquakes. From the pattern of microearthquake locations we recognize six aseismic zones where we expect future main shocks will rupture the Calaveras Fault. -from Authors

Journal of Geophysical Research Solid Earth↗

USGS2012

A tomographic glimpse of the upper mantle source of magmas of the Jemez lineament, New Mexico

The 800-km-long Jemez lineament is the most active volcanic feature in the southwestern United States. It is the southeastern tectonic boundary of the Colorado Plateau and crosses the Rio Grande rift at the Jemez Mountains. The primary volcanism of the lineament is basaltic and has occurred in the last 4.5 m.y. To infer spatial distributions of partial melt in the upper mantle source zones for the Rio Grande rift and the Jemez lineament, we investigated the lateral variations of P wave velocity in the upper mantle beneath these features. We used teleseismic P wave delays recorded at a 22-station network to perform a damped least squares, three-dimensional inversion for these lateral variations. Our technique employed velocity interpolation within a three-dimensional grid of points, rather than using blocks of constant P wave velocity. This method allows highly realistic computation of seismic ray paths as well as accurate computation of the matrix elements in our system of equations. Determinations of resolution of results were done in two independent ways, both of which gave consistent estimates of resolution. In our best resolved volume the inversion showed no significant concentration of relative low velocity for P waves beneath the Rio Grande rift. However, directly beneath the Jemez lineament there is a ∼100-km-wide, 1–2% low-velocity feature in the depth range of 50–160 km. Because of the association of the low P wave velocity with the Jemez volcanic lineament but not with the Rio Grande rift, because lowered P wave velocity can be associated with increased partial melt, and because the volume of recent volcanism at the lineament greatly exceeds that at the rift, we infer that a large magmatic source zone exists beneath the Jemez lineament but not beneath the Rio Grande rift. This implies that the volcanic potential of the Jemez lineaments continues to greatly exceed that of the Rio Grande rift. The mantle source zones for volcanics of the Jemez lineament are not overridden by, but rather track, the motion of the North American plate; this implies that these sources are within the lithospheric plate, as is clarified in the discussion. The magmatic source zones of the Jemez lineament are modeled as due to clockwise rotation of the Colorado Plateau about a pole in northeastern Colorado. This rotation caused extension of the lithosphere beneath the Jemez lineament, permitting concentration there of partially melted rock in the upper mantle.

Journal of Geophysical Research Solid Earth↗

USGS2012

Midcontinent rift volcanism in the Lake Superior region: Sr, Nd, and Pb isotopic evidence for a mantle plume origin

Between 1091 and 1098 Ma, most of a 15- to 20-km thickness of dominantly tholeiitic basalt erupted in the Midcontinent Rift System of the Lake Superior region, North America. The Portage Lake Volcanics in Michigan, which are the youngest MRS flood basalts, fall into distinctly high- and low-TiO 2 types having different liquid lines of descent. Incompatible trace elements in both types of tholeiites are enriched compared to depleted or primitive mantle (La/Yb = 4.3–5.3; Th/Ta = 2.12–2.16; Zr/Y = 4.3–4.4), and both basalt types are isotopically indistinguishable. Sr, Nd, and Pb isotopic compositions of the Portage Lake tholeiites have 87 Sr/ 86 Sr i ≈ 0.7038, ε Nd (1095 Ma) ≈ 0±2, and μ 1 ≈ 8.2. Model ages with respect to a depleted mantle source ( T DM ) average about 1950–2100 Ma. Portage Lake rhyolites fall into two groups. Type I rhyolites have Nd and Pb isotopic characteristics (ε Nd (1095 Ma) ≈ 0 to −4.7; μ 1 ≈ 8.2–7.8) consistent with contamination of tholeiitic rocks by 5–10% Archean crust. The one type II rhyolite analyzed has Nd and Pb isotopic compositions (ε Nd (1095 Ma) ≈ −13 to −16; μ 1 ≈ 7.6–7.7) which are consistent with partial melting of Archean crust. Early Proterozoic crust was not a major contaminant of MRS rocks in the Lake Superior region. Most reported Nd and Pb isotopic compositions of MRS tholeiites from the main stage of volcanism in the Lake Superior region and of the Duluth Complex are comparable to the Nd and Pb isotopic data for Portage Lake tholeiites. The isotopic enrichment of the MRS source compared to depleted mantle is striking and must have occurred at least 700 m.y. before 1100 Ma. There are two likely sources for such enriched MRS tholeiitic magmatism: subcontinental lithospheric mantle enriched during the early Proterozoic or enriched mantle derived from an upwelling plume. Subcontinental lithospheric mantle alone as a source cannot be ruled out strictly on isotopic grounds, but melting of enriched lithosphere potentially ranging in age from 3.6 to 1.8 Ga would have been less likely to generate either the homogeneity or the volume of MRS tholeiites. Decompression melting of an upwelling enriched mantle plume in a region of lithosphere thinned by extension could have successfully generated the enormous volume (850×10 3 km 3 ) of relatively homogeneous magma in a restricted time interval.

Journal of Geophysical Research Solid Earth↗

USGS2012

Gabbroic xenoliths from the northern Gorda Ridge: Implications for magma chamber processes under slow spreading centers

Abundant gabbroic xenoliths in porphyritic pillow basalt were dredged from the northern Gorda Ridge. The host lava is a moderately fractionated, normal mid-ocean ridge basalt with a heterogeneous glass rind (Mg numbers 56-60). Other lavas in the vicinity range from near primary (Mg number 69) to fractionated (Mg number 56). On the basis of textures and mineral compositions, the xenoliths are divided into five types. The xenoliths are not cognate to the host lava, but they are genetically related. Chemistry of mineral phases in conjunction with textural features suggests that the xenoliths formed in different parts of a convecting magma chamber that underwent a period of closed system fractionation. The chamber was filled with a large proportion of crystalline mush when new, more primitive, and less dense magma was injected and mixed incompletely with the contents in the chamber, forming the hybrid host lava. -from Authors

Journal of Geophysical Research Solid Earth↗

USGS2012

Teleseismic tomography of the compressional wave velocity structure beneath the Long Valley region, California

In 1982 and 1984 the U.S. Geological Survey used several seismic networks, totaling over 90 stations, to record teleseismic P waves and measure travel time residuals in an area centered on the Long Valley caldera. We inverted the travel time residuals to obtain a three-dimensional image of the velocity structure with resolution of 5–6 km to depths of 70 km beneath the array. Direct inversion of these data indicates that the 2- to 4-km-thick low-velocity caldera fill contaminates the signal from any midcrustal velocity anomalies beneath the caldera. Thus two methods were used to strip the effects of the upper crust from the travel time residuals: (1) ray tracing through upper crustal velocity models provided by seismic refraction experiments and gravity surveys, and (2) an iterative stripping scheme using the inversion itself. The methods produce essentially identical results and adequately remove the effects of the shallowest crustal structures, including the caldera fill and hydrothermal alteration effects. The resulting “stripped” models show two well-resolved midcrustal low-velocity bodies in the Long Valley region. The first body is centered between 7 and 20 km depth beneath the resurgent dome of the Long Valley caldera and has a volume of 150–600 km 3 . The second, with a similar volume, is centered between 10 and 20 km depth beneath the Mono Craters, about 10 km north of Long Valley. Velocity contrasts in both of these bodies are about 6–10%, and the features are interpreted as silicic magma chambers. This experiment does not preclude the presence of additional pockets of magma smaller than 5 km across in the upper crust, particularly beneath the resurgent dome of the caldera (which would be removed with the stripping methods). It is likely that the midcrust anomaly beneath the resurgent dome is a remnant of the caldera forming magma chamber. Upper mantle velocities are lower than average beneath Mono Craters and higher than average beneath Long Valley. We tentatively interpret these anomalies as mantle partial melt beneath the Mono Craters and refractory upper mantle residuum beneath Long Valley. The high eruptive rate of the Mono Craters and these upper mantle structures suggest that the focus of volcanism is shifting north from Long Valley to the Mono Craters.

Journal of Geophysical Research Solid Earth↗

USGS2012

Source complexity of the 1987 Whittier Narrows, California, earthquake from the inversion of strong motion records

Strong motion records for the Whittier Narrows earthquake are inverted to obtain the history of slip. Both constant rupture velocity models and variable rupture velocity models are considered. The results show a complex rupture process within a relatively small source volume, with at least four separate concentrations of slip. Two sources are associated with the hypocenter, the larger having a slip of 55–90 cm, depending on the rupture model. These sources have a radius of approximately 2–3 km and are ringed by a region of reduced slip. The aftershocks fall within this low slip annulus. Other sources with slips from 40 to 70 cm each ring the central source region and the aftershock pattern. All the sources are predominantly thrust, although some minor right-lateral strike-slip motion is seen. The overall dimensions of the Whittier earthquake from the strong motion inversions is 10 km long (along the strike) and 6 km wide (down the dip). The preferred dip is 30° and the preferred average rupture velocity is 2.5 km/s. Moment estimates range from 7.4 to 10.0×10 24 dyn cm, depending on the rupture model.

Journal of Geophysical Research Solid Earth↗

USGS2012

Discontinuities in the shallow Martian crust at Lunae, Syria, and Sinai Plana

Detailed photoclinometric profiles across 125 erosional features and 141 grabens in the western equatorial region of Mars indicate the presence of three discontinuities within the shallow crust. Pits, troughs, and wall valleys (tributary canyons) within Noctis Labyrinthus and Valles Marineris and escarpments within the fretted terrain of Sacra Fossae and Kasei Valles show distinct erosional base levels at depths of 0.3–0.6 km, 1 km, and 2–3 km. The shallowest discontinuity corresponds to thickness estimates for the ridged plains unit in this region, and thus the discontinuity probably is the contact between a sequence of layered rock making up this unit and the underlying megaregolith. The 1-km discontinuity is reflected in the base levels of erosion of all the features studied, and it may correspond to the base of the proposed layer of ground ice. Model calculations using photoclinometric profiles of simple grabens (corrected for the effects of mass wasting) and the 60° dip of bounding faults (measured from the exposed traces of faults on trough walls) show that graben-bounding faults consistently intersect at the mechanical discontinuity at about 1 km depth. This discontinuity may represent an interface between ice-laden and dry regolith, ice-laden and water-laden regolith, or pristine and cemented regolith. A correlation between wall valley head depth and local thickness of the faulted layer suggests that the 1-km discontinuity also controlled the depth of the heads of sapping canyons. There is no apparent relation between the pit and trough depths and the local thickness of the faulted layer, which can be explained if pits and troughs developed by subsidence into underlying tension cracks that disrupted the 1-km discontinuity. The data do suggest, however, that erosion of shallow pits and troughs was influenced by the 1-km discontinuity. The third discontinuity, at a depth of 2–3 km, corresponds to the proposed base of the Martian megaregolith and is probably the interface between overlying, ejected breccia and in situ, fractured basement rocks.

Journal of Geophysical Research Solid Earth↗

USGS2012

The relation of catastrophic flooding of Mangala Valles, Mars, to faulting of Memnonia Fossae and Tharsis volcanism

Detailed stratigraphic relations indicate two coeval periods of catastrophic flooding and Tharsis centered faulting (producing Memnonia Fossae) in the Mangala Valles region of Mars. Major sequences of lava flows of the Tharsis Montes Formation and local, lobate plains flows were erupted during and between these channeling and faulting episodes. First, Late Hesperian channel development overlapped in time the Tharsis-centered faulting that trends north 75° to 90°E. Next, Late Hesperian/Early Amazonian flooding was coeval with faulting that trends north 55° to 70°E. In some reaches, resistant lava flows filled the early channels, resulting in inverted channel topography after the later flooding swept through. Both floods likely originated from the same graben, which probably was activated during each episode of faulting. Faulting broke through groundwater barriers and tapped confined aquifers in higher regions west and east of the point of discharge. The minimum volume of water required to erode Mangala Valles (about 5×10 12 m 3 ) may have been released through two floods that drained a few percent pore volume from a relatively permeable aquifer. The peak discharges of the floods may have lasted from days to weeks. The perched water discharged from the aquifer may have been produced by hydrothermal groundwater circulation induced by Tharsis magmatism, tectonic uplift centered at Tharsis Montes, and compaction of saturated crater ejecta due to loading by lava flows.

Journal of Geophysical Research Solid Earth↗

USGS2012

Three-dimensional P and S velocity structure in the Coalinga Region, California

The Coalinga earthquake sequence of 1983 provided a unique opportunity to perform a three-dimensional velocity and hypocenter inversion in an area of complex three-dimensional structure dominated by folding and blind thrusts. Additionally, since other varied geological and geophysical studies have been completed in this area, the three-dimensional inversion solution could be compared to other interpretations. Inversion of 7696 P and 1511 S first arrivals from earthquakes and 696 P first arrivals from refraction shots produced a three-dimensional velocity model with grid spacing of 1–2 km in the hypocentral area. The overall shape and location of velocity features correspond well to the mapped surface geology. The three-dimensional inversion yields details of folds where the resolution is good and the general shape where resolution is lower. The amounts of structural relief inferred for the local folds are similar to values inferred from geologic data. The three-dimensional velocity solution has several distinctive features. There is a linear high-velocity body (6.1–6.5 km/s), about 25 km long, from 6 to 8 km depth, that may represent a fragment of Coast Range ophiolite. A shallow low-velocity zone (LVZ), which extends for 20 km along the fold axis at about 6 km depth and correlates with LVZs observed in both the refraction and the reflection data, may indicate high pore pressure caused by lateral compressive strain. Deeper LVZs occur within inferred Franciscan material and are characterized by horizontal or southwest dipping zones of varied thickness, 4–8 km wide and 5–10 km long. These LVZs may represent multiple thrust faults. Their locations and geometry are consistent with thrust faults inferred with seismic reflection data. The three-dimensional velocity solution compares well to prior two-dimensional seismic reflection and refraction models and observed gravity. The shape of the inferred sedimentary section agrees well with the reflections from the Cenozoic strata. Compared to the refraction model, the three-dimensional solution has similar velocities and similar locations of velocity features but is more detailed in the hypocentral zone where it uses more data. The gravity computed from the three-dimensional velocities is similar to the observed gravity in both shape and amplitude. Both a simple one-dimensional initial model and a complex initial model derived from the refraction interpretation were tried. A simple starting model gave the best results. The S velocity solution has different resolution than the P velocity solution because it uses a different set of stations, and it has lower resolution because it uses fewer arrival times. While the general patterns of velocity variation are similar for both V p and V s , the V s solution tends to have more smearing of velocity features and can have somewhat different locations of velocity features.

Journal of Geophysical Research Solid Earth↗

USGS2012

Some additional observations on inclusions in the granitic rocks of the Sierra Nevada

Microgranular quartz diorite and diorite inclusions are widespread in central Sierra Nevada granitoid rocks and are almost exclusively restricted to hornblende-bearing rocks, most commonly felsic tonalites and mafic granodiorites. Strontium 87/strontium 86 values of the inclusions and host granitoids from individual plutons generally plot on single isochrons that agree closely with previously determined ages. Neodymium 143/neodymium 144 values of five inclusions ranged from 0.51229 to 0.51248 and yield ε Nd values that correlate closely with values of their respective hosts. The Nd-Sm and Rb-Sr systematics indicate that most inclusions were in isotopic equilibrium with enclosing materials at the time of formation. Silica contents of inclusions and granitoids are contiguous, but inclusions generally contain less than, and granitoids more than, 60% SiO 2 . Linear continuity of the femic oxide and H 2 O + trends relative to SiO 2 suggests many inclusions formed as concentrations of hydrous mafic minerals. Variation of other major element oxides and trace elements support this inference. The aforementioned data are generally inconsistent with a direct genetic link between the mafic inclusions and the mantle-derived basalt. However, many initial Sierra Nevada magmas were produced by mixing of mafic, mantle-derived melts and felsic crustal materials, resulting in isotopically varied tonalitic and granodioritic magmas. Most inclusions thus represent fragmented crystal accumulations of early-formed, near-liquidus minerals generated from these previously mixed magmas.

Journal of Geophysical Research Solid Earth↗

USGS2012

Paleomagnetic rotations and the Cenozoic tectonics of the Cascade Arc, Washington, Oregon, and California

Paleomagnetic results from Cenozoic (62–12 Ma) volcanic rocks of the Cascade arc and adjacent areas indicate that moderate to large clockwise rotations are an important component of the tectonic history of the arc. Two mechanisms of rotation are suggested by the regional pattern of paleomagnetic rotations. The progressive increase in rotation toward the coast in arc and forearc rocks results from distributed dextral shear, which is likely driven by oblique subduction of oceanic plates to the west. Simple shear rotation is accommodated in the upper crust by strike-slip faulting. The right-lateral Mount St. Helens seismic zone may be an active manifestation of this process. Dextral shear probably obscures a subequal contribution to arc and forearc rotation that is driven by intraarc or backarc extension. This rotation is suggested by the average southward increase in continental margin rotations into the region outboard of the Basin and Range. The southward increase in rotation parallels a change in the arc tectonic regime from largely compressional in northern Washington to extensional in Oregon. Concomitant with this change is a southward increase in the volume of eruptive rocks and the number of basaltic vents in the arc. A progressive eastward shift of the arc volcanic front with time in the rotated arc terrane is the result of the westward pivoting of the arc block in front of a zone of extension since Eocene time. Westward migration of bimodal Basin and Range volcanism since at least 16 Ma is tracking westward rotation of the frontal arc block and growth of the Basin and Range in its wake.

Journal of Geophysical Research Solid Earth↗

USGS2012

Late Cenozoic volcanism, subduction, and extension in the Lassen region of California, southern Cascade Range

Some 537 volcanic vents younger than 7 Ma are identified and these are classified into five age intervals and five compositional categories based on SiO2 content. Maps of vents by age and composition illustrate regionally representative volcanic trends. By 2 Ma, the eastern limit of volcanism had contracted westward toward the late Quaternary arc. Late Quaternary volcanism is concentrated around and north of the silicic Lassen volcanic center. The belt of most recent volcanism (25-0 ka) has been active since at least 2 Ma. Most mafic volcanism is calcalkaline basalt and basaltic andesite. Normal faults and linear groups of vents are evidence of widespread crustal extension throughout most of the Lassen region. NNW orientation of maximum horizontal stress is indicated. The Lassen volcanic region is thought to occur above the subducting Gorda North plate but also within a broad zone of distributed extension that occurs in the North American lithosphere east and southeast of the present Cascadia subduction zone. The scarcity of volcanic rocks older than 7 Ma suggests that a more compressive lithospheric stress regime prior to the late Miocene extensional episode may have suppressed volcanism. -from Authors

Journal of Geophysical Research Solid Earth↗

USGS2012

Quaternary extrusion rates of the Cascade Range, northwestern United States and southern British Columbia

Quaternary (2–0 Ma) extrusion rates change significantly along the Cascade Range volcanic arc. The extrusion rate north of Mount Rainier is about 0.21 km 3 km −1 m.y. −1 ; the rate in southern Washington and northern Oregon south to Mount Hood is about 1.6 km 3 km −1 m.y. −1 ; in central Oregon the rate is 3–6 km 3 km −1 m.y. −1 ; and in northern California, the rate is 3.2 km 3 km −1 m.y. −1 . Eruption style also changes along the arc but at latitudes different from rate changes. At the ends of the arc, volcanism is focused at isolated intermediate to silicic composite volcanoes. The composite volcanoes represent ∼30% of the total volume of the arc. Mafic volcanic fields partly ring some composite volcanoes, especially in the south. In contrast, volcanism is diffused in the middle of the arc, where numerous overlapping mafic shields and a few composite volcanoes have built a broad ridge. Contrasting eruption style may signify diffuse versus focused heat sources or may reflect changes in permeability to ascending magma along the arc.

Journal of Geophysical Research Solid Earth↗

USGS2012

Geophysical constraints on Washington convergent margin structure

Gravity and magnetic maps of western Washington reveal the lateral structure and fabric of the Washington Coast Range, Puget Basin, and southern Washington Cascade Range. The magnetic and gravity maps show large amplitude positive anomalies associated with the shallow but largely buried section of Washington Coast Range mafic rocks which are separated by negative anomalies over deep sedimentary basins. The positive anomalies indicate that the Coast Range mafic basement extends farther east than previously thought, at least as far east as the longitude of Seattle. Linear and steep gravity and magnetic gradients indicate many unmapped, often buried faults in the Washington Coast Range Province. Magnetic highs are also associated with mapped batholiths in the Cascade arc. Several magnetic highs observed east of the Coast Range rocks and west of these batholiths may be associated with buried Tertiary plutons or ophiolites. Two-dimensional gravity and magnetic modeling constrained with geological and other geophysical data indicate that the Coast Range Province rocks are about 1 km thick at the coast, thickening to as much as 30 km near their postulated eastern edge. A maximum boundary on the average density of the upper 15–20 km of the rocks that compose the Coast Range Province of 2920 kg/m 3 was established by the modeling, suggesting a composition largely of basalt and gabbro with little interbedded sediments. Under these rocks may be mantle or a subduction complex composed of dense mafic, ultramafic, and sedimentary rocks like that proposed to underlie Vancouver Island. Previous gravity models of the Washington margin include lower densities for the proposed subduction complex than for Vancouver Island, suggesting a lower component of mafic and ultramafic rocks than the rocks underlying Vancouver Island. However, my Washington model requires that the proposed subduction complex be more dense than the trench sediments and, therefore, that material denser than sediments be incorporated within it. The absence of continental mantle and the modeled wedge shape of the Coast Range Province upper crust suggest that erosion of the bottom of the overriding plate by subduction processes may have occurred.

Washington↗

USGS2012

Stratigraphic, lithologic, and major element geochemical constraints on magmatic evolution at Lassen volcanic center, California

The Lassen volcanic center is the most recent of several long-lived volcanic centers in the southernmost Cascade Range. These centers have erupted products ranging from basaltic andesite to rhyolite and are superimposed on a background of regional basaltic to andesitic volcanism. The evolution of the Lassen volcanic center is described in three stages. Stages I and II comprise the Brokeoff volcano, and 80 km 3 andesitic stratocone, active from 600 to 400 ka. Brokeoff volcano is compositionally equivalent to the regional basaltic andesite to andesite volcanism in the Lassen region and is the result of structurally controlled focusing of the diffuse regional mafic magmatism. Stage III comprises a silicic dome field and adjacent area of hybrid andesites and has a total volume of about 100 km 3 . Volcanism during stage III was episodic and is subdivided into four sequences of lithologically and temporally distinct lavas. Stage III began at 400 ka with a rhyolitic, caldera-forming pyroclastic eruption and chemically related lavas. Additional sequences of dacite erupted between 250–200 ka and 100–0 ka. Hybrid andesites erupted adjacent to the silicic dome field between 300 and 0 ka. Porphyritic andesite and dacite with high Al 2 O 3 , low TiO 2 , medium K 2 O and FeO/MgO ratios of 1.5–2.0 are the most abundant rock types in the Lassen volcanic center. However, the single most voluminous unit is sparsely phyric rhyolite pumice. In general, the lavas of Lassen volcanic center form a single coherent trend on major element variation diagrams and in pseudo-quaternary phase space, consistant with an origin either by fractional crystallization or magma mixing. In detail, however, the lack of systematic temporal change in silica and subtly crossing trends indicate a complex origin. A variety of statistically successful fractional crystallization models can be constructed that derive Brokeoff andesites from regional magmas. An important conclusion of the modeling is that if fractional crystallization is the process responsible for generation of Brokeoff andesite, then the parent magma must be low to medium K in geochemical affinity in order to explain the variation in K 2 O. However, although major element variation can be modeled by fractional crystallization, petrographic and stratigraphic evidence indicates that magma mixing is an important but subtle process in Brokeoff lavas and suggests that lavas evolved in small independent batches. Lavas erupted during stage III, while predominantly silicic, range from 53 to 75% SiO 2 . Disequilibrium mineral assemblages in the stage III lavas indicate that they are not directly derived from Brokeoff andesite by fractional crystallization. Mixing of silicic magma with regional mafic magma and disaggregation of andesitic quenched magmatic inclusions play dominant roles in the compositional diversity of stage III lavas.

Journal of Geophysical Research Solid Earth↗

USGS2012

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

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.

Journal of Geophysical Research Solid Earth↗

USGS2012

Post-11,000-year volcanism at Medicine Lake Volcano, Cascade Range, northern California

Eruptive activity during the past 11,000 years at Medicine Lake volcano has been episodic. Eight eruptions produced about 5.3 km 3 of basaltic lava during an interval of a few hundred years about 10,500 years B.P. After a hiatus of about 6000 years, eruptive activity resumed with a small andesite eruption at about 4300 years B.P. Approximately 2.5 km 3 of lava with compositions ranging from basalt to rhyolite vented in nine eruptions during an interval of about 3400 years in late Holocene time. The most recent eruption occurred about 900 years B.P. A compositional gap in SiO 2 values of erupted lavas occurs between 58 and 63%. The gap is spanned by chilled magmatic inclusions in late Holocene silicic lavas. Late Holocene andesitic to rhyolitic lavas were probably derived by fractionation, assimilation, and mixing from high-alumina basalt parental magma, possibly from basalt intruded into the volcano during the early mafic episode. Many basaltic to andesitic lavas contain iron-rich crystals and have high FeO*/MgO (characteristics caused by mixing of high-alumina basalt with ferrobasalt liquid produced by fractionation of parental high-alumina basalt). When ferrobasalt and high-alumina basalt are contaminated with a granitic crustal component, a calc-alkaline trend is produced. Some eruptions have produced both tholeiitic and calc-alkaline compositions. The eruptive activity is probably driven by intrusions of basalt that occur during east-west stretching of the crust in an extensional tectonic environment. Vents are typically aligned parallel or subparallel to major structural features, most commonly within 30° of north. Intruded magma should provide adequate heat for commercial geothermal development if sufficient fluids can be found. The nature and timing of future volcanic activity cannot be predicted from the observed pattern, but eruptions high on the edifice could produce high-silica products that might be accompanied by explosive activity, whereas eruptions lower on the flanks are likely to vent more fluid mafic lavas.

Journal of Geophysical Research Solid Earth↗