Search USGS⌕ Search

SEARCH · Search USGS

Results for “One Earth”

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 847 records · Page 47Linked to original sources

Characterizing riverbed sediment using high-frequency acoustics 1: spectral properties of scattering

Bed-sediment classification using high-frequency hydro-acoustic instruments is challenging when sediments are spatially heterogeneous, which is often the case in rivers. The use of acoustic backscatter to classify sediments is an attractive alternative to analysis of topography because it is potentially sensitive to grain-scale roughness. Here, a new method is presented which uses high-frequency acoustic backscatter from multibeam sonar to classify heterogeneous riverbed sediments by type (sand, gravel,rock) continuously in space and at small spatial resolution. In this, the first of a pair of papers that examine the scattering signatures from a heterogeneous riverbed, methods are presented to construct spatially explicit maps of spectral properties from geo-referenced point clouds of geometrically and radiometrically corrected echoes. Backscatter power spectra are computed to produce scale and amplitude metrics that collectively characterize the length scales of stochastic measures of riverbed scattering, termed ‘stochastic geometries’. Backscatter aggregated over small spatial scales have spectra that obey a power-law. This apparently self-affine behavior could instead arise from morphological- and grain-scale roughnesses over multiple overlapping scales, or riverbed scattering being transitional between Rayleigh and geometric regimes. Relationships exist between stochastic geometries of backscatter and areas of rough and smooth sediments. However, no one parameter can uniquely characterize a particular substrate, nor definitively separate the relative contributions of roughness and acoustic impedance (hardness). Combinations of spectral quantities do, however, have the potential to delineate riverbed sediment patchiness, in a data-driven approach comparing backscatter with bed-sediment observations (which is the subject of part two of this manuscript).

Journal of Geophysical Research F: Earth Surface↗

The diffusion of ions in unconsolidated sediments

Diffusion in unconsolidated sediments generally proceeds at rates ranging from half to one twentieth of those applying to diffusion of ions and molecules in free solution. Diffusion rates are predictable with respect to porosity and path tortuosity in host sediments, and can be conveniently measured by determinations of electrical resistivity on bulk sediment samples. Net ion flux is further influenced by reactions of diffusing species with enclosing sediments, but such influences should not be confused with or lumped with diffusion processes.

Earth and Planetary Science Letters↗

Influence of stretching and density contrasts on the chemical evolution of continental magmas: An example from the Ivrea-Verbano Zone

The southern Ivrea-Verbano Zone of the Italian Western Alps contains a huge mafic complex that intruded high-grade metamorphic rocks while they were resident in the lower crust. Geologic mapping and chemical variations of the igneous body were used to study the evolution of underplated crust. Slivers of crustal rocks (septa) interlayered with igneous mafic rocks are concentrated in a narrow zone deep in the complex (Paragneiss-bearing Belt) and show evidence of advanced degrees of partial melting. Variations of rare-earth-element patterns and Sr isotope composition of the igneous rocks across the sequence are consistent with increasing crustal contamination approaching the septa. Therefore, the Paragneiss-bearing Belt is considered representative of an “assimilation region” where in-situ interaction between mantle- and crust-derived magmas resulted in production of hybrid melts. Buoyancy caused upwards migration of the hybrid melts that incorporated the last septa and were stored at higher levels, feeding the Upper Mafic Complex. Synmagmatic stretching of the assimilation region facilitated mixing and homogenization of melts. Chemical variations of granitoids extracted from the septa show that deep septa are more depleted than shallow ones. This suggests that the first incorporated septa were denser than the later ones, as required by the high density of the first-injected mafic magmas. It is inferred that density contrasts between mafic melts and crustal rocks play a crucial role for the processes of contamination of continental magmas. In thick under plated crust, the extraction of early felsic/hybrid melts from the lower crust may be required to increase the density of the lower crust and to allow the later mafic magmas to penetrate higher crustal levels.

Contributions to Mineralogy and Petrology↗

Non-double-couple microearthquakes at Long Valley caldera, California, provide evidence for hydraulic fracturing

Most of 26 small (0.4≲ M ≲3.1) microearthquakes at Long Valley caldera in mid-1997, analyzed using data from a dense temporary network of 69 digital three-component seismometers, have significantly non-double-couple focal mechanisms, inconsistent with simple shear faulting. We determined their mechanisms by inverting P - and S -wave polarities and amplitude ratios using linear-programming methods, and tracing rays through a three-dimensional Earth model derived using tomography. More than 80% of the mechanisms have positive (volume increase) isotropic components and most have compensated linear-vector dipole components with outward-directed major dipoles. The simplest interpretation of these mechanisms is combined shear and extensional faulting with a volume-compensating process, such as rapid flow of water, steam, or CO 2 into opening tensile cracks. Source orientations of earthquakes in the south moat suggest extensional faulting on ESE-striking subvertical planes, an orientation consistent with planes defined by earthquake hypocenters. The focal mechanisms show that clearly defined hypocentral planes in different locations result from different source processes. One such plane in the eastern south moat is consistent with extensional faulting, while one near Casa Diablo Hot Springs reflects en echelon right-lateral shear faulting. Source orientations at Mammoth Mountain vary systematically with location, indicating that the volcano influences the local stress field. Events in a ‘spasmodic burst’ at Mammoth Mountain have practically identical mechanisms that indicate nearly pure compensated tensile failure and high fluid mobility. Five earthquakes had mechanisms involving small volume decreases, but these may not be significant. No mechanisms have volumetric moment fractions larger than that of a force dipole, but the reason for this fact is unknown.

California↗

Functional metagenomic selection of RubisCOs from uncultivated bacteria

Ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) is a critical yet severely inefficient enzyme that catalyses the fixation of virtually all of the carbon found on Earth. Here, we report a functional metagenomic selection that recovers physiologically active RubisCO molecules directly from uncultivated and largely unknown members of natural microbial communities. Selection is based on CO 2 -dependent growth in a host strain capable of expressing environmental deoxyribonucleic acid (DNA), precluding the need for pure cultures or screening of recombinant clones for enzymatic activity. Seventeen functional RubisCO-encoded sequences were selected using DNA extracted from soil and river autotrophic enrichments, a photosynthetic biofilm and a subsurface groundwater aquifer. Notably, three related form II RubisCOs were recovered which share high sequence similarity with metagenomic scaffolds from uncultivated members of the G allionellaceae family. One of the G allionellaceae RubisCOs was purified and shown to possess CO 2 /O 2 specificity typical of form II enzymes. X-ray crystallography determined that this enzyme is a hexamer, only the second form II multimer ever solved and the first RubisCO structure obtained from an uncultivated bacterium. Functional metagenomic selection leverages natural biological diversity and billions of years of evolution inherent in environmental communities, providing a new window into the discovery of CO 2 -fixing enzymes not previously characterized.

Environmental Microbiology↗

Geologic field-trip guide of volcaniclastic sediments from snow- and ice-capped volcanoes—Mount St. Helens, Washington, and Mount Hood, Oregon

This field guide for the International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Scientific Assembly 2017 focuses on volcaniclastic sediments from Mount St. Helens in Washington and Mount Hood in Oregon. The trip spends four days in the field and includes nine stops at each volcano. For completeness, this guidebook also includes sixteen optional stops in the Mount St. Helens area and three in the Mount Hood area. These two volcanoes provide excellent depositional records of the broad spectrum of volcanic hazards that involve the flow or fall of volcaniclastic particles. At the field-trip stops we will contrast and compare the different types of deposits—sediments that can be frustratingly difficult to distinguish from one another. Correct identification of deposit origin leads to a better understanding of hazards that can impact vulnerable communities at particular volcanoes. During this four-day field trip, we will spend two days in the Mount St. Helens area and two days near Mount Hood. On the first day, we will visit debris-avalanche and lahar deposits in the Toutle River valley in Washington. The second day will focus on pyroclastic density current, pyroclastic-fall, lahar, and volcano-fluvial deposits in the Lewis River drainage system near Mount St. Helens. On the third day, we will visit hyaloclastites, glacial tills, and other diamicts in the Hood River and Sandy River basins near Mount Hood. The fourth and final day will focus on distal lahar and volcano-fluvial deposits downstream of Mount Hood. What follows is an in-depth introduction to the deposits studied on this field trip. We provide criteria that are observable in the field to aid in differentiating between pyroclastic density current, pyroclastic-fall, debris-avalanche, lahar, water-flood, and glacial deposits. We also introduce the Holocene eruptive histories of Mount St. Helens and Mount Hood and discuss the processes responsible for deposit emplacement. Field-trip stops and features of interest along the route are described in detail in a road log that provides daily cumulative mileage.

Oregon, Washington↗

An estimate of hydrothermal fluid residence times and vent chimney growth rates based on 210Pb Pb ratios and mineralogic studies of sulfides dredged from the Juan de Fuca Ridge

The 210 PbPb ratios across two sulfide samples dredged from the Juan de Fuca Ridge are used to estimate the growth rate of the sulfide material and the residence time of the hydrothermal fluid within the oceanic crust from the onset of basalt alteration. 210 Pb is added to the hydrothermal fluid by two processes: (1) high-temperature alteration of basalt and (2) if the residence time of the fluid is on the order of the 22.3-year half-life of 210 Pb, by in-situ growth from 222 Rn ( Krishnaswami and Turekian, 1982 ). Stable lead is derived only from the alteration of basalt. The 210 Pb/Pb ratio across one sample was ∼ 0.5 dpm/10 −6 g Pb, and across the other it was ∼ 0.4 dpm/10 −6 g Pb. These values are quite close to the 238 UP/b ratios of basalts from the area, suggesting that the residence time of the hydrothermal fluid from the onset of basalt alteration is appreciably less than the mean life of 210 Pb, i.e., the time required for ingrowth from the radon. An apparent growth rate of 1.2 cm/yr is derived from the slope of the 210 Pb/Pb curve for one of the samples. This is consistent with its mineralogy and texture which suggest an accretionary pattern of development. There is no obvious sequential growth pattern, and virtually no gradient in 210 Pb/Pb across the second sample. This is consistent with alteration of the original 210 Pb/Pb distribution by extensive remobilization reactions which are inferred from the mineralogic and textural relationships of the sample.

Earth and Planetary Science Letters↗

Age and Sr isotopic composition of volcanic rocks in the Maricunga Belt, Chile: Implications for magma sources

K-Ar and 40 Ar/ 39 Ar dating of volcanic rocks from the Maricunga belt of north-central Chile indicate that igneous activity took place throughout most of Miocene time at various places in the 150 by 30 km belt. No migration patterns of volcanism appear in the Miocene rocks of the belt. Volcanic activity ceased by the end of the Miocene. All the Miocene volcanic rocks studied are calcic andesites to dacites with about 62% SiO 2 , 18% A1 2 O 3 , 4% Fe 2 O 3 (total), and 2% K 2 O. Initial 87 Sr/ 86 Sr (Sr i ) values fall into two groups, one of lower values around 0.7050 to the west and the other of higher values around 0.7060 to the east. It is postulated that the two Sr i groups reflect two adjoining coherent lower-crustal magma sources of possibly different age and subtly different composition that form part of the western edge of the South American craton.

Journal of South American Earth Sciences↗

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

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

Alaska↗

Basin and range crustal and upper mantle structure, northwest to central Nevada

We present an interpretation of the crustal and uppermost mantle structure of the Basin and Range of northwestern Nevada based on seismic refraction/wide‐angle reflection, near‐vertical reflection, and gravity data. In comparison to most previous estimates, we find that the crust is somewhat thicker (32–36 km versus 22–30 km), and the uppermost mantle velocity is somewhat higher (8.0 km/s versus 7.3–7.9 km/s). Along our transects, the crust is thinnest (32 km) in the Carson Sink‐Buena Vista Valley region and increases by 2–4 km to the west and east, respectively. There is considerable complexity throughout the crust where velocities range from of 2.5 km/s at the surface to 7.4 km/s in the lowermost crust. Variations in velocity and structure of the upper crustal layers reveal apparent basement velocity depressions (areas of lower velocities extending up to 10 km in depth) that underlie some surface ranges as well as the basins. The middle crust rises from about 20 km beneath central Nevada to within 12 km of the surface beneath the area of thinnest crust and is characterized by a modest (∼0.1 km/s) change in velocity and low‐velocity gradients. These midcrustal layers mark the onset of high crustal reflectivity and the apparent limiting depth to which Basin and Range faults can be traced in near‐vertical reflection profiles, suggesting that these midcrustal layers represent the transition between the brittle and ductile zones of the crust. The lower crust is more structurally complex, with layers thickening and thinning in a systematic manner with the upper crustal layers; generally, where there are velocity depressions in the upper crust, the lower crust is thickest and shallowest. The geometry of these lower crustal layers (derived from refraction modeling) coincides with changes in the crustal reflectivity, determined from the Consortium of Continental Reflection Profiling reflection data. The lower crustal layer is unusually high in velocity (7.4 km/s) and is likely the layer identified as mantle in some previous studies. We do not identify the 7.4 km/s layer as mantle because (1) there is an underlying layer with a velocity (8.0 km/s) that is more consistent with the worldwide average velocity for the upper mantle, and (2) the 7.4 km/s layer does not correspond to the “reflection” Moho. Gravity modeling and comparison to existing seismic models show a general consensus in many aspects with respect to crustal structure. This new model forms the basis for speculation on some of the processes associated with rifting of the Basin and Range Province. One such process, lithospheric magmatism, is inferred from the strong attenuation of transmitted seismic waves, which occurs at the same interface at which high‐amplitude, bright spot reflections originate. Unlike previous models, the overall structure and velocity of the crust and uppermost mantle of our new model are similar to other regions worldwide which have undergone high degrees of extension.

Nevada↗

Estimation of rod scale errors in geodetic leveling

Comparisons among repeated geodetic levelings have often been used for detecting and estimating residual rod scale errors in leveled heights. Individual rod-pair scale errors are estimated by a two-step procedure using a model based on either differences in heights, differences in section height differences, or differences in section tilts. It is shown that the estimated rod-pair scale errors derived from each model are identical only when the data are correctly weighted, and the mathematical correlations are accounted for in the model based on heights. Analyses based on simple regressions of changes in height versus height can easily lead to incorrect conclusions. We also show that the statistically estimated scale errors are not a simple function of height, height difference, or tilt. The models are valid only when terrain slope is constant over adjacent pairs of setups (i.e., smoothly varying terrain). In order to discriminate between rod scale errors and vertical displacements due to crustal motion, the individual rod-pairs should be used in more than one leveling, preferably in areas of contrasting tectonic activity. From an analysis of 37 separately calibrated rod-pairs used in 55 levelings in southern California, we found eight statistically significant coefficients that could be reasonably attributed to rod scale errors, only one of which was larger than the expected random error in the applied calibration-based scale correction. However, significant differences with other independent checks indicate that caution should be exercised before accepting these results as evidence of scale error. Further refinements of the technique are clearly needed if the results are to be routinely applied in practice.

Journal of Geophysical Research B: Solid Earth↗

Depth dependence of earthquake frequency-magnitude distributions in California: Implications for rupture initiation

Statistics of earthquakes in California show linear frequency-magnitude relationships in the range of M 2.0 to M 5.5 for various data sets. Assuming Gutenberg-Richter distributions, there is a systematic decrease in b value with increasing depth of earthquakes. We find consistent results for various data sets from northern and southern California that both include and exclude the larger aftershock sequences. We suggest that at shallow depth (∼0 to 6 km) conditions with more heterogeneous material properties and lower lithospheric stress prevail. Rupture initiations are more likely to stop before growing into large earthquakes, producing relatively more smaller earthquakes and consequently higher b values. These ideas help to explain the depth-dependent observations of foreshocks in the western United States. The higher occurrence rate of foreshocks preceding shallow earthquakes can be interpreted in terms of rupture initiations that are stopped before growing into the mainshock. At greater depth (9–15 km), any rupture initiation is more likely to continue growing into a larger event, so there are fewer foreshocks. If one assumes that frequency-magnitude statistics can be used to estimate probabilities of a small rupture initiation growing into a larger earthquake, then a small ( M 2) rupture initiation at 9 to 12 km depth is 18 times more likely to grow into a M 5.5 or larger event, compared to the same small rupture initiation at 0 to 3 km.

Journal of Geophysical Research B: Solid Earth↗

Populations and home range relationships of the box turtle, Terrapene c. carolina (Linnaeus)

SUMMARY: A population study of the box turtle (Terrapene c. carolina Linnaeus) was made during the years 1944 to 1947 at the Patuxent Research Refuge, Maryland. A thirty acre area in well drained bottomland forest on the flood plain of the Patuxent River was selected for intensive study. Similarly forested land extended in all directions from the study plot. Markers were established at eighty-three foot intervals over the study plot for reference in recording locality data. Individuals were marked by filing notches in the marginal scutes according to a code system. There were 2109 collections of study area turtles. Records of collecting sites and turtle behavior showed that in the bottomlands habitat cover is utilized extensively during the day as well as at night. Turtles not actively moving about are almost always found in or around brush piles, heaps of debris, and tangles of vines and briars. Gully banks and woods openings are used for sunning. Turtles are occasionally found in the mud or water of the gullies. The commonest type of night retreat is a cavity constructed by the turtle in leaves, debris, or earth. These cavities, termed 'forms,' may be used only once, but are sometimes used repeatedly, often at intervals of several days or more. Different turtles sometimes use the same form on successive nights. Weather conditions most favorable to turtle activity are high humidity, warm sunny days, and frequent rains. The most unfavorable influences are low temperatures and drought. On most summer days there are some active turtles but individual turtles are not active every day. Periods of activity are alternated with periods of quiet even in favorable weather. This behavior is most pronounced in early spring and late fall when inactive days are often more numerous than active ones. Adult turtles occupy specific home ranges which they maintain from year to year. The turtles living in the study plot retained their ranges even through a flood that completely covered the area. Maximum home range diameters were determined by measurements of the mapped ranges of individual turtles. The average range of adult males was 330 feet, adult females 370 feet. The difference between male and female ranges was not statistically significant. There was no evidence of defense of territory. Ranges of turtles of all ages and both sexes overlapped grossly. Turtles were frequently found near each other and no antagonistic behavior was observed. A trail-laying device was developed in order to follow individual travel routes. The trailer consists of a light weight housing fastened to the turtle's back. It contains a spool of white thread that unwinds as the turtle moves, thus marking its exact route. Turtles selected for this more detailed study were followed with trailers for a total of 456 turtle days. Maps illustrating their travels are shown. Normal movements within the home range are characterized by, (1) turns, doublings, detours, and criss-crossing paths completely covering the area, (2) interspersion of fairly direct traverses of the home range, (3) frequently repeated travels over certain paths or routes. Trailer records and mapped collection records both show that the maximum limits of the home range are ordinarily reached within a few days or weeks. This general procedure is varied by some turtles to include intensive coverage of only one portion of the range at a time. Some turtles have two home ranges and travel between them at infrequent intervals. One turtle showing this behavior was followed with a trailer for 161 days during 1946 and 1947. Trips outside the home range are made by some turtles. These include egg laying trips by females as well as trips of unexplained nature made by both males and females. Turtles from other areas occasionally occur as transients in the study plot. The size of the population was estimated on the basis of collections during on

Ecological Monographs↗

Site distribution of iron in staurolite

The Mössbauer spectrum of staurolite does not necessarily indicate that ferrous iron is distributed over both the tetrahedral Fe site and one or more of the octahedral sites, as previously thought. Two doublets are present in the spectrum, but their quadrupole spliting and isomer shift are more consistent with both being due to tetrahedrally coordinated ferrous iron. Smith's X-ray work has shown that cations in the Fe site actually assume several slightly different positions, and iron in each of these positions may give rise to a different Mössbauer doublet.

Earth and Planetary Science Letters↗

Simulations of historical and future trends in snowfall and groundwater recharge for basins draining to Long Island Sound

A regional watershed model was developed for watersheds contributing to Long Island Sound, including the Connecticut River basin. The study region covers approximately 40 900 km 2 , extending from a moderate coastal climate zone in the south to a mountainous northern New England climate zone dominated by snowmelt in the north. The input data indicate that precipitation and temperature have been increasing for the last 46 years (1961– 2006) across the region. Minimum temperature has increased more than maximum temperature over the same period (1961–2006). The model simulation indicates that there was an upward trend in groundwater recharge across most of the modeled region. However, trends in increasing precipitation and groundwater recharge are not significant at the 0.05 level if the drought of 1961–67 is removed from the time series. The trend in simulated snowfall is not significant across much of the region, although there is a significant downward trend in southeast Connecticut and in central Massachusetts. To simulate future trends, two input datasets, one assuming high carbon emissions and one assuming low carbon emissions, were developed from GCM forecasts. Under both of the carbon emission scenarios, simulations indicate that historical trends will continue, with increases in groundwater recharge over much of the region and substantial snowfall decreases across Massachusetts, Connecticut, southern Vermont, and southern New Hampshire. The increases in groundwater recharge and decreases in snowfall are most pronounced for the high emission scenario.

Connecticut, Massachusetts, New Hampshire, Rhode I↗

The USGS National Geologic Map Database – Resource for geoscience professionals and the public

Since 1996, the US Geological Survey (USGS) National Geologic Map Database (NGMDB) has served the public and private sectors and is widely recognized as one of the few authoritative resources for the information published by the nation’s geological surveys. On a monthly basis, the NGMDB website receives more than 600,000 visits by about 200,000 users (~20,000 visits per day). In a typical month, between 2 and 3 TB of content are viewed or downloaded. This high level of web traffic results in many user requests for information and assistance. Users vary widely in interest and background and include homeowners, students, local government planners, and professional geologists and engineers. Many professionals use NGMDB information for understanding natural conditions that affect infrastructure, such as drainage, ground stability, properties of near-surface earth materials, and depth to rock.

Conference Paper↗

Lead concentration and isotopic composition in five peridotite inclusions of probable mantle origin

The lead content of five whole-rock peridotite inclusions (four lherzolites and one harzburgite) in alkali basalt ranges from 82 to 570 ppb (parts per billion). Approximately 30–60 ppb of this amount can be accounted for by analyzed major silicate minerals (olivine ≤ 10 ppb; enstatite 5–28 ppb; chrome diopside ∼400 ppb). Through a series of acid leaching experiments, the remainder of the lead is shown to be quite labile and to reside in either glassy or microcrystalline veinlets or accessory mineral phases, such as apatite and mica. The lead isotopic composition of the peridotites ( 206 Pb/ 204 Pb= 18.01–18.90; 207 Pb/ 204 Pb= 15.52–15.61; 208 Pb/ 204 Pb= 37.80–38.86) lies within the range of values defined by many modern volcanic rocks and, in particular, is essentially coextensive with the abyssal tholeiite field. In all but one instance, isotopic differences were found between the peridotite and its host alkali basalt. Two of the peridotites clearly demonstrated internal isotopic heterogeneity between leachable and residual fractions that could not simply be due to contamination by the host basalt. However, there is no evidence that these ultramafic rocks form some layer in the mantle with isotopic characteristics fundamentally different from those of the magma sources of volcanic rocks.

Earth and Planetary Science Letters↗

Partially melted granodiorite and related rocks ejected from Crater Lake caldera, Oregon

Blocks of medium-grained granodiorite to 4 m, and minor diabase, quartz diorite, granite, aplite and granophyre, are common in ejecta of the ∼6,900 yrBP calderaforming eruption of Mount Mazama. The blocks show degrees of melting from 0–50 vol%. Because very few have adhering juvenile magma, it is thought that the blocks are fragments of the Holocene magma chamber's walls. Primary crystallisation of granodiorite produced phenocrystic pl + hyp + aug + mt + il + ap + zc, followed by qz + hb + bt + alkali feldspar (af). Presence of fluid inclusions in all samples implies complete crystallisation before melting. Subsolidus exchange with meteoric hydrothermal fluids before melting is evident in δ 18 O values of −3·4+4·9‰ for quartz and plagioclase in partially melted granodiorites (fresh lavas from the region have δ 18 O values of +5·8−+7·0‰); δ 18 O values of unmelted granodiorites from preclimatic eruptive units suggest hydrothermal exchange began between ∼70 and 24 ka. Before eruption, the granitic rocks equilibrated at temperatures, estimated from Fe-Ti oxide compositions, of up to ∼1000°C for c. 10 2 –10 4 years at a minimum pressure of 100-180 MPa. Heating caused progressive breakdown or dissolution of hb, af, bt, and qz, so that samples with the highest melt fractions have residual pl + qz and new or re-equilibrated af + hyp + aug + mt + il in high-silica rhyolitic glass (75-77% SiO 2 ). Mineral compositions vary systematically with increasing temperature. Hornblende is absent in rocks with Fe-Ti oxide temperatures >870°C, and bt above 970°C. Oxygen isotope fractionation between qz, pl, and glass in partially fused granodiorite also is consistent with equilibration at T≥900°C (Δ 18 O qz.pl = +0·7 ± 0·5‰). Element partitioning between glass and crystals reflects the large fraction of refractory pl, re-equilibration of af and isolation or incomplete dissolution of accessory phases. Ba and REE contents of analysed glass separates can be successfully modelled by observed degrees of partial melting of granodiorite, but Rb, Sr and Sc concentrations cannot. Several samples have veins of microlite-free glass 1–5 mm thick that are compositionally and physically continuous with intergranular melt and which apparently formed after the climactic eruption began. Whole-rock H 2 O content, microprobe glass analysis sums near 100% and evidence for high temperature suggest liquids in the hotter samples were nearly anhydrous. The occurrence of similar granodiorite blocks at all azimuths around the 8 × 10 km caldera implies derivation from one pluton. Compositional similarity between granodiorite and pre-Mazama rhyodacites suggests that the pluton may have crystallised as recently as 0·4 Ma; compositional data preclude crystallisation from the Holocene chamber. The history of crystallisation, hydrothermal alteration, and remelting of the granitic rocks may be characteristic of shallow igneous systems in which the balance between hydrothermal cooling and magmatic input changes repeatedly over intervals of 10 4 -10 6 years.

Oregon↗