Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal Interpretation”

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.

1,223 records · Page 68Linked to original sources
USGS2008

Methanogenic pathways of coal-bed gas in the Powder River Basin, United States: The geologic factor

Coal-bed gas of the Tertiary Fort Union and Wasatch Formations in the Powder River Basin in Wyoming and Montana, U.S. was interpreted as microbial in origin by previous studies based on limited data on the gas and water composition and isotopes associated with the coal beds. To fully evaluate the microbial origin of the gas and mechanisms of methane generation, additional data for 165 gas and water samples from 7 different coal-bed methane-bearing coal-bed reservoirs were collected basinwide and correlated to the coal geology and stratigraphy. The C1/(C2 + C3) ratio and vitrinite reflectance of coal and organic shale permitted differentiation between microbial gas and transitional thermogenic gas in the central part of the basin. Analyses of methane ??13C and ??D, carbon dioxide ??13C, and water ??D values indicate gas was generated primarily from microbial CO2 reduction, but with significant gas generated by microbial methyl-type fermentation (aceticlastic) in some areas of the basin. Microbial CO2 reduction occurs basinwide, but is generally dominant in Paleocene Fort Union Formation coals in the central part of the basin, whereas microbial methyl-type fermentation is common along the northwest and east margins. Isotopically light methane ??13C is distributed along the basin margins where ??D is also depleted, indicating that both CO2-reduction and methyl-type fermentation pathways played major roles in gas generation, but gas from the latter pathway overprinted gas from the former pathway. More specifically, along the northwest basin margin gas generation by methyl-type fermentation may have been stimulated by late-stage infiltration of groundwater recharge from clinker areas, which flowed through highly fractured and faulted coal aquifers. Also, groundwater recharge controlled a change in gas composition in the shallow Eocene Wasatch Formation with the increase of nitrogen and decrease of methane composition of the coal-bed gas. Other geologic factors, such as burial, thermal and maturation history, lateral and vertical continuity, and coalification of the coal beds, also played a significant role in controlling methanogenic pathways and provided new perspectives on gas evolution and emplacement. The early-stage gas produced by CO2 reduction has mixed with transitional thermogenic gas in the deeper, central parts of the Powder River Basin to form 'old' gas, whereas along the basin margins the overprint of gas from methyl-type fermentation represents 'new' gas. Thus, a clear understanding of these geologic factors is necessary to relate the microbiological, biogeochemical, and hydrological processes involved in the generation of coal-bed gas.

International Journal of Coal Geology↗
USGS2017

Magmatic degassing, lava dome extrusion, and explosions from Mount Cleveland volcano, Alaska, 2011–2015: Insight into the continuous nature of volcanic activity over multi-year timescales

Mount Cleveland volcano (1730 m) is one of the most active volcanoes in the Aleutian arc, Alaska, but heightened activity is rarely accompanied by geophysical signals, which makes interpretation of the activity difficult. In this study, we combine volcanic gas emissions measured for the first time in August 2015 with longer-term measurements of thermal output and lava extrusion rates between 2011 and 2015 calculated from MODIS satellite data with the aim to develop a better understanding of the nature of volcanic activity at Mount Cleveland. Degassing measurements were made in the month following two explosive events (21 July and 7 August 2015) and during a period of new dome growth in the summit crater. SO 2 emission rates ranged from 400 to 860 t d − 1 and CO 2 /SO 2 ratios were < 3, consistent with the presence of shallow magma in the conduit and the observed growth of a new lava dome. Thermal anomalies derived from MODIS data from 2011 to 2015 had an average repose time of only 4 days, pointing to the continuous nature of volcanic activity at this volcano. Rapid increases in the cumulative thermal output were often coincident with visual confirmation of dome growth or accumulations of tephra in the crater. The average rate of lava extrusion calculated for 9 periods of rapid increase in thermal output was 0.28 m 3 s − 1 , and the total volume extruded from 2011 to 2015 was 1.9–5.8 Mm 3 . The thermal output from the lava extrusion events only accounts for roughly half of the thermal budget, suggesting a continued presence of shallow magma in the upper conduit, likely driven by convection. Axisymmetric dome morphology and occasional drain back of lava into the conduit suggests low-viscosity magmas drive volcanism at Mount Cleveland. It follows also that only small overpressures can be maintained given the small domes and fluid magmas, which is consistent with the low explosivity of most of Mount Cleveland's eruptions. Changes between phases of dome growth and explosive activity are somewhat unpredictable and likely result from plugs that are related to the dome obtaining a critical dimension, or from small variations in the magma ascent rate that lead to crystallization-induced blockages in the upper conduit, thereby reducing the ability of magma to degas. We suggest the small magma volumes, slow ascent rates, and low magma viscosity lead to the overall lack of anomalous geophysical signals prior to eruptions, and that more continuous volcanic degassing measurements might lead to more successful eruption forecasting at this continuously-active open-vent volcano.

Alaska↗
USGS2010

Depositional setting, petrology and chemistry of Permian coals from the Paraná Basin: 2. South Santa Catarina Coalfield, Brazil

In Brazil economically important coal deposits occur in the southern part of the Paraná Basin, where coal seams occur in the Permian Rio Bonito Formation, with major coal development in the states of Rio Grande de Sul and Santa Catarina. The current paper presents results on sequence stratigraphic interpretation of the coal-bearing strata , and petrological and geochemical coal seam characterization from the South Santa Catarina Coalfield, Paraná Basin. In terms of sequence stratigraphic interpretation the precursor mires of the Santa Catarina coal seams formed in an estuarine-barrier shoreface depositional environment , with major peat accumulation in a high stand systems tract (Pre-Bonito and Bonito seams), a lowstand systems tract (Ponta Alta seam, seam A, seam B) and a transgressive systems tract (Irapuá, Barro Branco and Treviso seams). Seam thicknesses range from 1.70 to 2.39 m, but high proportions of impure coal (coaly shale and shaley coal), carbonaceous shale and partings reduce the net coal thickness significantly. Coal lithoypes are variable, with banded coal predominant in the Barro Branco seam, and banded dull and dull coal predominantly in Bonito and Irapuá seams, respectively. Results from petrographic analyses indicate a vitrinite reflectance range from 0.76 to 1.63 %Rrandom (HVB A to LVB coal). Maceral group distribution varies significantly, with the Barro Branco seam having the highest vitrinite content (mean 67.5 vol%), whereas the Irapuá seam has the highest inertinite content (33.8 vol%). Liptinite mean values range from 7.8 vol% (Barro Branco seam) to 22.5 vol% (Irapuá seam). Results from proximate analyses indicate for the three seams high ash yields (50.2 – 64.2 wt.%). Considering the International Classification of in-Seam Coals, all samples are in fact classified as carbonaceous rocks (> 50 wt.% ash). Sulfur contents range from 3.4 to 7.7 wt.%, of which the major part occurs as pyritic sulfur. Results of X-ray diffraction indicate the predominance of quartz and kaolinite (also pyrite). Gypsum , gibbsite , jarosite and calcite were also identified in some samples. Feldspar was noted but is rare. The major element distribution in the three seams (coal basis) is dominated by SiO 2 (31.3 wt.%, mean value), Al 2 O 3 (14.5 wt.%, mean value) and Fe 2 O 3 (6.9 wt.%, mean value). Considering the concentrations of trace elements that are of potential environmental hazards the Barro Branco, Bonito and Irapuá seams (coal base) are significantly enriched in Co (15.7 ppm), Cr (54.5 ppm), Li (59.3 ppm), Mn (150.4 ppm), Pb (58.0 ppm) and V (99.6 ppm), when compared to average trace elements contents reported for U. S. coals. Hierarchical cluster analysis identified, based on similarity levels, three groups of major elements and seven groups of trace elements. Applying discriminant analyses using trace and major element distribution, it could be demonstrated that the three seams from Santa Catarina show distinct populations in the discriminant analyses plots, and also differ from the coals of Rio Grande do Sul analyzed in a previous study.

Santa Catarina↗
USGS2016

RMT focal plane sensitivity to seismic network geometry and faulting style

Modern tectonic studies often use regional moment tensors (RMTs) to interpret the seismotectonic framework of an earthquake or earthquake sequence; however, despite extensive use, little existing work addresses RMT parameter uncertainty. Here, we quantify how network geometry and faulting style affect RMT sensitivity. We examine how data-model fits change with fault plane geometry (strike and dip) for varying station configurations. We calculate the relative data fit for incrementally varying geometries about a best-fitting solution, applying our workflow to real and synthetic seismograms for both real and hypothetical station distributions and earthquakes. Initially, we conduct purely observational tests, computing RMTs from synthetic seismograms for hypothetical earthquakes and a series of well-behaved network geometries. We then incorporate real data and station distributions from the International Maule Aftershock Deployment (IMAD), which recorded aftershocks of the 2010 M W 8.8 Maule earthquake, and a set of regional stations capturing the ongoing earthquake sequence in Oklahoma and southern Kansas. We consider RMTs computed under three scenarios: (1) real seismic records selected for high data quality; (2) synthetic seismic records with noise computed for the observed source-station pairings and (3) synthetic seismic records with noise computed for all possible station-source pairings. To assess RMT sensitivity for each test, we observe the ‘fit falloff’, which portrays how relative fit changes when strike or dip varies incrementally; we then derive the ranges of acceptable strikes and dips by identifying the span of solutions with relative fits larger than 90 per cent of the best fit. For the azimuthally incomplete IMAD network, Scenario 3 best constrains fault geometry, with average ranges of 45° and 31° for strike and dip, respectively. In Oklahoma, Scenario 3 best constrains fault dip with an average range of 46°; however, strike is best constrained by Scenario 1, with a range of 26°. We draw two main conclusions from this study. (1) Station distribution impacts our ability to constrain RMTs using waveform time-series; however, in some tectonic settings, faulting style also plays a significant role and (2) increasing station density and data quantity (both the number of stations and the number of individual channels) does not necessarily improve RMT constraint. These results may be useful when organizing future seismic deployments (e.g. by concentrating stations in alignment with anticipated nodal planes), and in computing RMTs, either by guiding a more rigorous data selection process for input data or informing variable weighting among the selected data (e.g. by eliminating the transverse component when strike-slip mechanisms are expected).

Geophysical Journal International↗
USGS1989

Multiple hydrothermal and metamorphic events in the Kidd Creek volcanogenic massive sulphide deposit, Timmins, Ontario: evidence from tourmalines and chlorites

Tourmaline and chlorite are the principal ferromagnesian silicate minerals in the Kidd Creek massive sulphide deposit. Tourmaline is most common in sphalerite-rich peripheral margins of the chalcopyrite stringer zone. Within the north orebody, samples typically contain <1% tourmaline, but small areas (hand-specimen scale) may have 10–20%. Chlorite is more widely distributed and in places constitutes 30–50% of rock volumes. Associated assemblages may include quartz, sulphides (principally chalcopyrite, sphalerite, and (or) pyrite), carbonate, albite, sericite, and rare fluorite, allanite, or zoisite(?).The tourmalines and chlorites record a series of multiple hydrothermal and metamorphic events. Paragenetic studies suggest that tourmaline was deposited during several discrete stages of mineralization, as evidenced by brecciation and cross-cutting relationships. Most of the tourmalines have two concentric growth zones defined by different colours (green, brown, blue, yellow). Some tourmalines also display pale discordant rims that cross-cut and embay the inner growth zones and polycrystalline, multiple-extinction domains. Late sulphide veinlets (chalcopyrite, pyrrhotite) transect the inner growth zones and pale discordant rims of many crystals. The concentric growth zones are interpreted as primary features developed by the main ore-forming hydrothermal system, whereas the discordant rims, polycrystalline domains, and cross-cutting sulphide veinlets reflect post-ore metamorphic processes.Detailed electron microprobe analyses of tourmalines show a wide compositional range, from Fe-rich dravite nearly to end-member schorl, with Fe/(Fe + Mg) ratios varying from 0.33 to 0.92; only minor amounts of Ca are present, yielding uniformly high Na/(Na + Ca) ratios of 0.84–0.99. Two sets of chemical zoning trends are identified in the tourmalines, involving systematic changes in Fe/(Fe + Mg), Na/(Na + Ca), Al, and Ti that are believed to reflect internal coupled substitutions (e.g.,  + Ti = Na + Al) and local mineral equilibria (e.g., tourmaline–chlorite). Analyses of the pale discordant reaction rims show consistent depletion of Fe, Ca, and Ti, presumably by fluid–solid reactions during post-ore metamorphism.Chlorites also show an extensive range in composition, from ripidolite nearly to end-member daphnite, with Fe/(Fe + Mg) ratios of 0.43–0.98 and Si cation values of 5.00–5.39. Chlorites from the fringes of the footwall stringer zone have narrow compositional ranges, whereas chlorites near footwall rhyolite sills in the core of the stringer zone display major variations in Fe/(Fe + Mg) ratios, including one sample with a range of 0.68–0.95. The former group of chlorites has Fe/(Fe + Mg) ratios that correlate well with those of coexisting tourmalines (exclusive of late reaction rims). Data for the latter group, in contrast, fall off equilibrium K D curves, indicating that the tourmalines and chlorites within these samples are not in chemical equilibrium. The chlorites are believed to have been altered (overprinted) by Fe-rich hydrothermal fluids apparently generated during intrusion of the rhyolite sills. The tourmalines, however, are unaffected and retain primary chemical signatures.Variations in mineral proportions and mineral chemistry within the deposit mainly depend on fluctuations in temperature, pH, water/rock ratios, and amounts of entrained seawater. The major proposed control is mixing between high-temperature, Fe-rich end-member hydrothermal fluids and cold, Mg-rich entrained seawater. Fe/(Fe + Mg) variations in footwall tourmalines (and equilibrium chlorites) are believed to largely reflect the progressive infiltration of Mg-rich seawater into the margins and top of the hydrothermal system. The more Fe-rich compositions of Kidd Creek tourmalines relative to those from sediment-hosted massive sulphide deposits (e.g., Sullivan, British Columbia) may be related to the preferential generation of end-member hydrothermal fluids in proximal volcanic environments like that at Kidd Creek.

Ontario↗
USGS2013

Merapi 2010 eruption—Chronology and extrusion rates monitored with satellite radar and used in eruption forecasting

Despite dense cloud cover, satellite-borne commercial Synthetic Aperture Radar (SAR) enabled frequent monitoring of Merapi volcano's 2010 eruption. Near-real-time interpretation of images derived from the amplitude of the SAR signals and timely delivery of these interpretations to those responsible for warnings, allowed satellite remote sensing for the first time to play an equal role with in situ seismic, geodetic and gas monitoring in guiding life-saving decisions during a major volcanic crisis. Our remotely sensed data provide an observational chronology for the main phase of the 2010 eruption, which lasted 12 days (26 October–7 November, 2010). Unlike the prolonged low-rate and relatively low explosivity dome-forming and collapse eruptions of recent decades at Merapi, the eruption began with an explosive eruption that produced a new summit crater on 26 October and was accompanied by an ash column and pyroclastic flows that extended 8 km down the flanks. This initial explosive event was followed by smaller explosive eruptions on 29 October–1 November, then by a period of rapid dome growth on 1–4 November, which produced a summit lava dome with a volume of ~ 5 × 10 6 m 3 . A paroxysmal VEI 4 magmatic eruption (with ash column to 17 km altitude) destroyed this dome, greatly enlarged the new summit crater and produced extensive pyroclastic flows (to ~ 16 km radial distance in the Gendol drainage) and surges during the night of 4–5 November. The paroxysmal eruption was followed by a period of jetting of gas and tephra and by a second short period (12 h) of rapid dome growth on 6 November. The eruption ended with low-level ash and steam emissions that buried the 6 November dome with tephra and continued at low levels until seismicity decreased to background levels by about 23 November. Our near-real-time commercial SAR documented the explosive events on 26 October and 4–5 November and high rates of dome growth (> 25 m 3 s − 1 ). An event tree analysis for the previous 2006 Merapi eruption indicated that for lava dome extrusion rates > 1.2 m 3 s − 1 , the probability of a large (1872-scale) eruption was ~ 10%. Consequently, the order-of-magnitude greater rates in 2010, along with the explosive start of the eruption on 26 October, the large volume of lava accumulating at the summit by 4 November, and the rapid and large increases in seismic energy release, deformation and gas emissions were the basis for warnings of an unusually large eruption by the Indonesian Geological Agency's Center for Volcanology and Geologic Hazard Mitigation (CVGHM) and their Volcano Research and Technology Development Center (BPPTK) in Yogyakarta — warnings that saved thousands of lives.

Merapi Volcano↗
USGS2007

Fault locking, block rotation and crustal deformation in the Pacific Northwest

We interpret Global Positioning System (GPS) measurements in the northwestern United States and adjacent parts of western Canada to describe relative motions of crustal blocks, locking on faults and permanent deformation associated with convergence between the Juan de Fuca and North American plates. To estimate angular velocities of the oceanic Juan de Fuca and Explorer plates and several continental crustal blocks, we invert the GPS velocities together with seafloor spreading rates, earthquake slip vector azimuths and fault slip azimuths and rates. We also determine the degree to which faults are either creeping aseismically or, alternatively, locked on the block-bounding faults. The Cascadia subduction thrust is locked mainly offshore, except in central Oregon, where locking extends inland. Most of Oregon and southwest Washington rotate clockwise relative to North America at rates of 0.4–1.0 ° Myr −1 . No shear or extension along the Cascades volcanic arc has occurred at the mm/yr level during the past decade, suggesting that the shear deformation extending northward from the Walker Lane and eastern California shear zone south of Oregon is largely accommodated by block rotation in Oregon. The general agreement of vertical axis rotation rates derived from GPS velocities with those estimated from palaeomagnetic declination anomalies suggests that the rotations have been relatively steady for 10–15 Ma. Additional permanent dextral shear is indicated within the Oregon Coast Range near the coast. Block rotations in the Pacific Northwest do not result in net westward flux of crustal material—the crust is simply spinning and not escaping. On Vancouver Island, where the convergence obliquity is less than in Oregon and Washington, the contractional strain at the coast is more aligned with Juan de Fuca—North America motion. GPS velocities are fit significantly better when Vancouver Island and the southern Coast Mountains move relative to North America in a block-like fashion. The relative motions of the Oregon, western Washington and Vancouver Island crustal blocks indicate that the rate of permanent shortening, the type that causes upper plate earthquakes, across the Puget Sound region is 4.4 ± 0.3 mm yr −1 . This shortening is likely distributed over several faults but GPS data alone cannot determine the partitioning of slip on them. The transition from predominantly shear deformation within the continent south of the Mendocino Triple Junction to predominantly block rotations north of it is similar to changes in tectonic style at other transitions from shear to subduction. This similarity suggests that crustal block rotations are enhanced in the vicinity of subduction zones possibly due to lower resisting stress.

Geophysical Journal International↗
USGS2001

Anomalously high b-values in the South Flank of Kilauea volcano, Hawaii: Evidence for the distribution of magma below Kilauea's East rift zone

The pattern of b-value of the frequency-magnitude relation, or mean magnitude, varies little in the Kaoiki-Hilea area of Hawaii, and the b-values are normal, with b = 0.8 in the top 10 km and somewhat lower values below that depth. We interpret the Kaoiki-Hilea area as relatively stable, normal Hawaiian crust. In contrast, the b-values beneath Kilauea's South Flank are anomalously high (b = 1.3-1.7) at depths between 4 and 8 km, with the highest values near the East Rift zone, but extending 5-8 km away from the rift. Also, the anomalously high b-values vary along strike, parallel to the rift zone. The highest b-values are observed near Hiiaka and Pauahi craters at the bend in the rift, the next highest are near Makaopuhi and also near Puu Kaliu. The mildest anomalies occur adjacent to the central section of the rift. The locations of the three major and two minor b-value anomalies correspond to places where shallow magma reservoirs have been proposed based on analyses of seismicity, geodetic data and differentiated lava chemistry. The existence of the magma reservoirs is also supported by magnetic anomalies, which may be areas of dike concentration, and self-potential anomalies, which are areas of thermal upwelling above a hot source. The simplest explanation of these anomalously high b-values is that they are due to the presence of active magma bodies beneath the East Rift zone at depths down to 8 km. In other volcanoes, anomalously high b-values correlate with volumes adjacent to active magma chambers. This supports a model of a magma body beneath the East Rift zone, which may widen and thin along strike, and which may reach 8 km depth and extend from Kilauea's summit to a distance of at least 40 km down rift. The anomalously high b-values at the center of the South Flank, several kilometers away from the rift, may be explained by unusually high pore pressure throughout the South Flank, or by anomalously strong heterogeneity due to extensive cracking, or by both phenomena. The major b-value anomalies are located SSE of their parent reservoirs, in the direction of motion of the flank, suggesting that magma reservoirs leave an imprint in the mobile flank. We hypothesize that the extensive cracking may have been acquired when the anomalous parts of the South Flank, now several kilometers distant from the rift zone, were generated at the rift zone near persistent reservoirs. Since their generation, these volumes may have moved seaward, away from the rift, but earthquakes occurring in them still use the preexisting complex crack distribution. Along the decollement plane at 10 km depth, the b-values are exceptionally low (b = 0.5), suggesting faulting in a more homogeneous medium. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS1999

Hydrated salt minerals on Europa's surface from the Galileo near-infrared mapping spectrometer (NIMS) investigation

We reported evidence of heavily hydrated salt minerals present over large areas of Europa's surface from analysis of reflectance spectra returned by the Galileo mission near infrared mapping spectrometer (NIMS) [ McCord et al. , 1997a, b, 1998a, b]. Here we elaborate on this earlier evidence, present spatial distributions of these minerals, examine alternate water-ice interpretations, expand on our hydrated-salts interpretation, consider salt mineral stability on Europa, and discuss the implications. Extensive well-defined areas on Europa show distinct, asymmetric water-related absorption bands in the 1 to 2.5-μm region. Radiative transfer modeling of water ice involving different particle sizes and layers at Europa temperatures does not reproduce the distinctive Europa water bands. However, ice near its melting temperature, such as in terrestrial environments, does have some characteristics of the Europa spectrum. Alternatively, some classes of heavily hydrated minerals do exhibit such water bands. Among plausible materials, heavily hydrated salt minerals, such as magnesium and sodium sulfates, sodium carbonate and their mixtures, are preferred. All Europa spectral features are present in some salt minerals and a very good match to the Europa spectrum can be achieved by mixing several salt spectra. However, no single or mix of salt mineral spectra from the limited library available has so far been found to perfectly match the Europa spectrum in every detail. The material is concentrated at the lineaments and in chaotic terrain, which are technically disrupted areas on the trailing side. Since the spectrum of the material on Europa is nearly the same everywhere so-far studied, the salt or salt-mixture composition may be nearly uniform. This suggests similar sources and processes over at least a near-hemispheric scale. This would suggest that an extensive subsurface ocean containing dissolved salts is the source, and several possible mechanisms for deposit emplacement are considered. The hydrogen bonds associated with hydration of these salts are similar or greater in strength and energy to those in pure water ice. Thus, once on the surface, the salt minerals should be as stable to disruption as water ice at the Europa temperatures, and mechanisms are suggested to enhance the stability of both materials. Spectra obtained of MgSO 4 ·6H 2 O at 77 K show only small differences from room temperature spectra. The main difference is the appearance of the individual absorptions composing the broad, composite water features and associated with the several different H 2 O sites in the salt hydrate molecule. This suggests that the Europa absorption bands are also composites. Thus higher spectral resolution may reveal these diagnostic features in Europa's spectrum. The specific salts present and their relative abundances would be indicators of the chemistry and conditions of an ocean environment, and areas of fresh, heavy concentration of these minerals should make ideal lander mission sampling sites.

Journal of Geophysical Research E: Planets↗
USGS2008

Trace fossils and sedimentary facies from a Late Cambrian‐Early Ordovician tide‐dominated shelf (Santa Rosita Formation, northwest Argentina): Implications for ichnofacies models of shallow marine successions

The Santa Rosita Formation is one the most widely distributed lower Paleozoic units of northwest Argentina. At the Quebrada del Salto Alto section, east of Purmamarca, Jujuy Province, it is represented by four sedimentary facies: thick‐bedded planar cross‐stratified quartzose sandstones (A), thin‐bedded planar cross‐stratified quartzose sandstones and mudstones (B), wave‐rippled sandstones and bioturbated mudstones (C), and black and greenish gray shales (D). Paleocurrent data, sandstone architecture, and sedimentary structures from faciès A and B indicate bipolar/bimodal paleoflows, suggesting the action of tidal currents. The succession is interpreted as that of a tide‐dominated shelf, with only secondary influence of wave processes. Trace fossils are restricted to facies B and C. The Cruziana ichnocoenosis is preserved on the soles of thin‐bedded planar cross‐stratified quartzose sandstones (faciès B). This ichnocoenosis consists of Conostichus isp., Cruziana omanica, C. semiplicata, C. cf. tortworthi, Cruziana isp. Helminthopsis abeli, Monomorphichnus bilinearis, M. multilineatus, Palaeophycus tubularis, Rusophycus carbonarias, R. latus , and R. isp. The occurrence of Cruziana semiplicata , C. omanica, C. cf. tortworthi , and Rusophycus latus supports a Late Cambrian‐Tremadoc age. Slabbing of Cruziana shows complex interactions between biologic and sedimentologic processes, and suggests a predominance of exhumed traces, washed out and recast by tractive sand deposition. Sandstone soles are densely packed with biogenic structures and exhibit distinctive clusters of Rusophycus isp. that most likely represent trilobite nesting burrows. The Cruziana ichnocoenosis records the resident fauna of a protected, lower intertidal to subtidal interbar setting. The Skolithos ichnocoenosis is represented by high to low density vertical burrows of Skolithos linearis , which extend downwards to the quartzose sandstone soles of faciès B and cross the Cruziana ichnocoenosis. The Skolithos ichnocoenosis represents colonization by suspension‐feeding organisms following a major change in environmental conditions, related to the migration of lower intertidal to subtidal sandwaves. The Planolites ichnocoenosis consists exclusively of Planolites montanus within mudstones overlying wave‐rippled sandstones (facies C). The Planolites ichnocoenosis records opportunistic colonization by inf aunal deposit feeders that mined the organic‐rich fine‐grained sediment during the waning phase of storms that scoured organic detritus from the sea bottom. The section records, from base to top, a Cruziana‐Skolithos ichnofacies zone, a Skolithos ichnofacies zone and an unbioturbated zone typified by the thick‐bedded cross‐stratified quartzose sandstone (fades A). This trend reflects progressively higher energy conditions linked to the establishment of a large sand wave complex. The presence of a mixed Cruziana‐Skolithos ichnofacies in the lower interval reflects changes in substrate and energy levels, rather than water depth. Accordingly, contrasting ichnocoenoses from interbars (Cruziana) and sandwaves (Skolithos) must be considered an example of ichnofacies controlled by local parameters instead of general bathymétrie trends. Conversely, the vertical replacement of the Cruziana ichnofacies by the Skolithos ichnofacies towards the middle interval of the section reflects the environmental changes associated with the transition between the intertidal and subtidal zones. As overall tidal energy increases from supratidal to subtidal settings, the Skolithos ichnofacies tends to occur seaward of the Cruziana ichnofacies in tide‐dominated shallow marine environments. Therefore, onshore‐offshore ichnofacies replacement in tide‐dominated shallow seas is opposite to that in wave‐dominated marine settings.

Ichnos: An International Journal for Plant and Ani↗
USGS2012

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↗
USGS2016

Erosional and depositional history of the Atlantic passive margin as recorded in detrital zircon fission-track ages and lithic detritus in Atlantic Coastal plain sediments

Comparison of fission-track (FT) ages of detrital zircons recovered from Atlantic Coastal Plain sediments to FT ages of zircons from bedrock in source terranes in the Appalachians provides a key to understanding the provenance of the sediments and, in turn, the erosional and depositional history of the Atlantic passive margin. In Appalachian source terranes, the oldest zircon fission-track (ZFT) ages from bedrock in the western Appalachians (defined for this paper as the Appalachian Plateau, Valley and Ridge, and far western Blue Ridge) are notably older than the oldest ages from bedrock in the eastern Appalachians (Piedmont and main part of the Blue Ridge). The age difference is seen both in ZFT sample ages and in individual zircon grain ages and reflects differences in the thermotectonic history of the rocks. In the east, ZFT data indicate that the rocks cooled from temperatures high enough to partially or totally reset ZFT ages during the Paleozoic and (or) Mesozoic. The majority of the rocks are interpreted to have cooled through the ZFT closure temperature (∼235 °C) at various times during the late Paleozoic Alleghanian orogeny. In contrast, most of the rocks sampled in the western Appalachians have never been heated to temperatures high enough to totally reset their ZFT ages. Reflecting their contrasting thermotectonic histories, nearly 80 percent of the sampled western rocks yield one or more zircon grains with very old FT ages, in excess of 800 Ma; zircon grains yielding FT ages this old have not been found in rocks in the Piedmont and main part of the Blue Ridge. The ZFT data suggest that the asymmetry of zircon ages of exposed bedrock in the eastern and western Appalachians was in evidence by no later than the Early Cretaceous and probably by the Late Triassic. Detrital zircon suites from sands collected in the Atlantic Coastal Plain provide a record of detritus eroded from source terranes in the Appalachians during the Mesozoic and Cenozoic. In Virginia and Maryland, sands of Early Cretaceous through late early Oligocene age do not yield any old zircons comparable in age to the old zircons found in bedrock in the western Appalachians. Very old zircons yielding FT ages >800 Ma are only encountered in Coastal Plain sands of middle early Miocene and younger age. Miocene and younger fluvial-deltaic deposits associated with the major mid-Atlantic Coastal Plain rivers that now head in the western Appalachians (the Hudson, Delaware, Susquehanna, Potomac, James, and Roanoke) contain abundant clasts of fossiliferous chert and quartzite and other distinctive rock types derived from Paleozoic rocks of the western Appalachians. These distinctive clasts have not been reported in older Coastal Plain sediments. The ZFT and lithic detritus data indicate that the drainage divide for one or more east-flowing mid-Atlantic rivers migrated west into the western Appalachians, and the river(s) began transporting western Appalachian detritus to the Atlantic Coastal Plain, sometime between the late early Oligocene and middle early Miocene. By no later than late middle Miocene most if not all of the major rivers that now head west of the Blue Ridge were transporting western Appalachian detritus to the Coastal Plain. Prior to the drainage divide migrating into the western Appalachians, the ZFT data are consistent with the dominant source of Atlantic Coastal Plain sediments being detritus from the Piedmont and main part of the Blue Ridge, with possible input from distant volcanic sources. The ZFT data suggest that the rapid increase in the rate of siliciclastic sediment accumulation in middle Atlantic margin offshore basins that peaked in the middle Miocene and produced almost 30 percent of the total volume of post-rift siliciclastic sediments in the offshore basins began in the early Miocene when Atlantic river(s) gained access to the relatively easily eroded Paleozoic sedimentary rocks of the western Appalachians.

American Journal of Science↗
USGS1982

Application of automated image analysis to coal petrography

The coal petrologist seeks to determine the petrographic characteristics of organic and inorganic coal constituents and their lateral and vertical variations within a single coal bed or different coal beds of a particular coal field. Definitive descriptions of coal characteristics and coal facies provide the basis for interpretation of depositional environments, diagenetic changes, and burial history and determination of the degree of coalification or metamorphism. Numerous coal core or columnar samples must be studied in detail in order to adequately describe and define coal microlithotypes, lithotypes, and lithologic facies and their variations. The large amount of petrographic information required can be obtained rapidly and quantitatively by use of an automated image-analysis system (AIAS). An AIAS can be used to generate quantitative megascopic and microscopic modal analyses for the lithologic units of an entire columnar section of a coal bed. In our scheme for megascopic analysis, distinctive bands 2 mm or more thick are first demarcated by visual inspection. These bands consist of either nearly pure microlithotypes or lithotypes such as vitrite/vitrain or fusite/fusain, or assemblages of microlithotypes. Megascopic analysis with the aid of the AIAS is next performed to determine volume percentages of vitrite, inertite, minerals, and microlithotype mixtures in bands 0.5 to 2 mm thick. The microlithotype mixtures are analyzed microscopically by use of the AIAS to determine their modal composition in terms of maceral and optically observable mineral components. Megascopic and microscopic data are combined to describe the coal unit quantitatively in terms of (V) for vitrite, (E) for liptite, (I) for inertite or fusite, (M) for mineral components other than iron sulfide, (S) for iron sulfide, and (VEIM) for the composition of the mixed phases (X i ) i = 1,2, etc. in terms of the maceral groups vitrinite V, exinite E, inertinite I, and optically observable mineral content M. The volume percentage of each component present is indicated by a subscript. For example, a lithologic unit was determined megascopically to have the composition (V) 13 (I) 1 (S) 1 (X 1 ) 83 (X 2 ) 2 . After microscopic analysis of the mixed phases, this composition was expressed as (V) 13 (I) 1 (S) 1 (V 63 E 19 I 14 M 4 ) 83 (V 67 E 11 I 13 M 9 ) 2 . Finally, these data were combined in a description of the bulk composition as V 67 E 16 I 13 M 3 S 1 . An AIAS can also analyze textural characteristics and can be used for quick and reliable determination of rank (reflectance). Our AIAS is completely software based and incorporates a television (TV) camera that has optimum response characteristics in the range of reflectance less than 5%, making it particularly suitable for coal studies. Analysis of the digitized signal from the TV camera is controlled by a microprocessor having a resolution of 64 gray levels between full illumination and dark current. The processed image is reconverted for display on a TV monitor screen, on which selection of phases or features to be analyzed is readily controlled and edited by the operator through use of a lightpen. We expect that automated image analysis, because it can rapidly provide a large amount of pertinent information, will play a major role in the advancement of coal petrography.

International Journal of Coal Geology↗
USGS2008

The evolution of volcano-hosted geothermal systems based on deep wells from Karaha-Telaga Bodas, Indonesia

Temperature and pressure surveys, fluid samples, and petrologic analyses of rock samples from deep drill holes at the Karaha - Telaga Bodas geothermal field on the volcanic ridge extending northward from Galunggung Volcano, West Java, have provided a unique opportunity to characterize the evolution of an active volcano-hosted geothermal system. Wells up to 3 km in depth have encountered temperatures as high as 353??C and a weakly altered granodiorite that intruded to within 2 to 3 km of the surface. The intrusion is shallowest beneath the southern end of the field where an acid lake overlies a nearly vertical low resistivity structure (<10 ohm-m) defined by magnetotelluric measurements. This structure is interpreted to represent a vapor-dominated chimney that provides a pathway to the surface for magmatic gases. Four distinct hydrothermal mineral assemblages document the evolution of the geothermal system and the transition from liquid- to vapor-dominated conditions. The earliest assemblage represents the initial liquid-dominated system generated during emplacement of the granodiorite between 5910 ?? 76 and 4200 ?? 150 y BP. Tourmaline, biotite, actinolite, epidote and clay minerals were deposited contemporaneously at progressively greater distances from the intrusive contact (assemblage 1). At 4200 ?? 150 y BP, flank collapse and the formation of the volcano's crater, Kawah Galunggung, resulted in catastrophic decompression and boiling of the hydrothermal fluids. This event initiated development of the modern vapor-dominated regime. Chalcedony and then quartz were deposited as the early low salinity liquids boiled (assemblage 2). Both vapor- and liquid-rich fluid inclusions were trapped in the quartz crystals. Liquid-rich fluid inclusions from the southern part of the field record salinities ranging from 0 to 26 weight percent NaCl- CaCl2 equivalent and locally contain fluorite daughter crystals. We suggest, based on temperature-salinity relationships and evidence of boiling, that these fluids were progressively concentrated as steam was lost from the system. However, mixing with fluids derived from the underlying intrusion or generated during the formation of acid SO4 water on the vapor-dominated chimney margins could have contributed to the observed salinities. As pressures declined, CO2- and SO4-rich steam-heated water drained downward, depositing anhydrite and calcite (assemblage 3) in the fractures, limiting further recharge. Fluid inclusions with salinities up to 31 weight percent NaCl equivalent were trapped in these minerals as the descending water vaporized. The final assemblage is represented by precipitates of NaCl, KCl and FeClx deposited on rock surfaces in portions of the vapor-dominated zone that boiled dry. Vapor-dominated conditions extend over a distance of at least 10 km and to depths below sea level. Deep wells drilled into the underlying liquid-dominated reservoir in the northern and central part of the volcanic ridge produce low salinity fluids representing recent recharge of meteoric and steam-heated water. The evolution of volcanic-hosted vapor-dominated geothermal systems can be described by a five stage model. Stage 1 involves the formation of an over-pressured liquid-dominated geothermal system soon after magmatic intrusion. In Stages 2 and 3, pressures progressively decrease, and a curtain of steam-heated water surrounding a magmatic vapor-dominated chimney at 350??C and 14 ?? 2 MPa develops. The relatively low pressure near the base of the chimney causes liquid inflow adjacent to the intrusion and the development of a secondary marginal vapor-dominated zone. In Stage 4, the magmatic vapor discharge from the intrusion becomes small, vapor pressure declines, and the secondary vapor-dominated zone expands above the intrusion. In Stage 5, the vapor-dominated zone floods because heat from the intrusion is insufficient to boil all liquid inflow. A more common, liquid-dominated volcanic-hosted system the

American Journal of Science↗
USGS2012

The mechanics and three-dimensional internal structure of active magmatic systems: Kilauea volcano, Hawaii

Interpretation of abundant seismic data suggests that Kilauea's primary conduit within the upper mantle is concentrically zoned to about 34-km depth. This zoned structure is inferred to contain a central core region of relatively higher permeability, surrounded by numerous dikes that are in intermittent hydraulic communication with each other and with the central core. During periods of relatively high magma transport, the entire cross section of the conduit is utilized. During periods of relatively low to moderate transport, however, only the central core is active. As the conduit penetrates the oceanic crust and enters the volcanic shield, it simultaneously supplies the deeper sections of the rift zones (6-to 10-km depth) and the roots of the summit reservoir with picritic magma. The rift zones at depth are inferred to be almost wholly molten and to possess a high degree of fluid continuity from Heiheiahulu in the East Rift Zone, 45 km westward through the roots of the summit magma reservoir, and well into the Southwest Rift Zone. Higher in the shield, the subcaldera magma reservoir and the shallow rift zones occupy the 2-to 4-km depth interval. Summit-differentiated olivine tholeiite (ρ ≈ 2.62 g cm −3 ) is periodically injected laterally along a horizon of neutral buoyancy within the rift zones, where the density of the magma is just balanced by the in situ density of the shield (Ryan, 1987 a , b ). Deep rift zone intrusions push seaward the deep tectonic blocks of the volcano's south flank. Shallow rift intrusions build a sheeted dike complex, inferred to be in isostatic equilibrium with the higher-density deep rift cores below. General finite element analyses are presented for the deformation and stress fields surrounding such dikes in the horizontal and vertical planes. The dike tip in two and three dimensions is surrounded by a tubular core of tensile (σ 1 , σ 2 ) and shear stress (τ max ). The displacement field is characterized by counterrotating cells on either side of the dike tip which, in vertical orientation, produce the characteristic subsidence above the dike complex, with uplift on either side, forming a ridge-trough-ridge structure. A finite element model of Kilauea's shield computes the displacement fields and principal stress (σ 1 ) distributions resulting from intrusive activity on each or both of the rift zones. Within the summit region, tensile stress lobes produced by the three-dimensional upward extension of the intrusions superpose constructively to produce calderawide regimes of tensile stress, conducive to caldera development. Parametric studies of (1) intrusion in the East Rift Zone only, (2) intrusion in the Southwest Rift Zone only, and (3) intrusion in both rift zones demonstrate their unique kinematic contributions. For case 1, the caldera undergoes a counterclockwise rotation (torque up state) conducive to the development of rightstepping en echelon eruptive fissures, as exemplified by the August 14, 1971, eruption. For case 2, the caldera undergoes a clockwise rotation (torque down state) conducive to the development of left-stepping eruptive fissures, as occurred during the December 31, 1974, eruption. For case 3, the caldera substructure is driven due southward, producing the southward migration of the upper portions of the summit magma reservoir.

Journal of Geophysical Research Solid Earth↗
USGS2004

Shear wave velocity, seismic attenuation, and thermal structure of the continental upper mantle

Seismic velocity and attenuation anomalies in the mantle are commonly interpreted in terms of temperature variations on the basis of laboratory studies of elastic and anelastic properties of rocks. In order to evaluate the relative contributions of thermal and non-thermal effects on anomalies of attenuation of seismic shear waves, Q −1 s , and seismic velocity, V s , we compare global maps of the thermal structure of the continental upper mantle with global Q −1 s and V s maps as determined from Rayleigh waves at periods between 40 and 150 s. We limit the comparison to three continental mantle depths (50, 100 and 150 km), where model resolution is relatively high. The available data set does not indicate that, at a global scale, seismic anomalies in the upper mantle are controlled solely by temperature variations. Continental maps have correlation coefficients of <0.56 between V s and T and of <0.47 between Q s and T at any depth. Such low correlation coefficients can partially be attributed to modelling artefacts; however, they also suggest that not all of the V s and Q s anomalies in the continental upper mantle can be explained by T variations. Global maps show that, by the sign of the anomaly, V s and Q s usually inversely correlate with lithospheric temperatures: most cratonic regions show high V s and Q s and low T , while most active regions have seismic and thermal anomalies of the opposite sign. The strongest inverse correlation is found at a depth of 100 km, where the attenuation model is best resolved. Significantly, at this depth, the contours of near-zero Q s anomalies approximately correspond to the 1000 °C isotherm, in agreement with laboratory measurements that show a pronounced increase in seismic attenuation in upper mantle rocks at 1000–1100 °C. East–west profiles of V s , Q s and T where continental data coverage is best (50°N latitude for North America and 60°N latitude for Eurasia) further demonstrate that temperature plays a dominant, but non-unique, role in determining the value of lithospheric V s and Q s . At 100 km depth, where the resolution of seismic models is the highest, we compare observed seismic V s and Q s with theoretical V T s and Q T s values, respectively, that are calculated solely from temperature anomalies and constrained by experimental data on temperature dependencies of velocity and attenuation. This comparison shows that temperature variations alone are sufficient to explain seismic V s and Q s in ca 50 per cent of continental regions. We hypothesize that compositional anomalies resulting from Fe depletion can explain the misfit between seismic and theoretical V s in cratonic lithosphere. In regions of active tectonics, temperature effects alone cannot explain seismic V s and Q s in the lithosphere. It is likely that partial melts and/or fluids may affect seismic parameters in these regions. This study demonstrates that lithospheric temperature plays the dominant role in controlling V s and Q s anomalies, but other physical parameters, such as compositional variations, fluids, partial melting and scattering, may also play a significant role in determining V s and Q s variations in the continental mantle.

Geophysical Journal International↗
USGS2012

Lead and strontium isotopes and related trace elements as genetic tracers in the Upper Cenozoic rhyolite-basalt association of the Yellowstone Plateau volcanic field

Supported by various field geologic and petrologic data, the contents of Pb, U, Th, Rb, and Sr and the isotopic compositions of Pb and Sr for upper Cenozoic volcanic rocks of the Yellowstone Plateau volcanic field are consistent with the hypothesis of derivation of the basaltic and rhyolitic magmas by partial melting of distinct source regions in the upper mantle and lower crust, respectively. All the basalt samples analyzed but one have systematically lower values of 207 Pb/ 204 Pb and 87 Sr/ 86 Sr than the rhyolites. The values of 206 Pb/ 204 Pb are smaller, and 87 Sr/ 86 Sr are mostly larger than known values in oceanic basalts. In all but one case, the values of 207 Pb/ 204 Pb are higher than expected from an extrapolation of known values in oceanic basalts to less radiogenic values of 206 Pb/ 204 Pb. Because there are no xenoliths, phenocrysts are only moderate to sparse in abundance, REE patterns are low and flat at the radiogenic end of lead isotopic compositions, several values of Rb/Sr are low, and 80% of the basalt samples form a well-developed secondary isochron separate from the rhyolites, we favor an interpretation for basalt genesis wherein isotopic signatures of most mafic magmas were attained in a continental ‘keel’ of mantlelike character about 2.6 b.y. old or somewhat older attached to the crust, and these signatures were unaltered by magma passage through the crust. At the very least, the current data continue to cast serious doubt as to the inevitability of crustal contamination for basaltic magma intruding the continental environment and postulate that much can be learned about the mantle under continents through the study of continental basalts. One basalt unit with an unusually low value of 207 Pb/ 204 Pb and an 87 Ar/ 86 Ar less than 0.704 may represent subcontinental ‘keel’-derived magma that rose unaltered to the surface. Our data also are not consistent with formation of this rhyolite-basalt association primarily by such processes as crystal fractionation, separation of immiscible silicate liquids from a common parental magma, or fractional melting of a homogeneous source. Rather as a conceptual model, we envision large mafic intrusions to have been injected into the lower crust resulting in rhyolite generation through partial anatexis of the adjacent wall rocks which probably had a 206 Pb/ 204 Pb < 17 and 87 Sr/ 86 Sr > 0.709; a model that has much in common with that proposed by Holmes (1931). All the other hypotheses listed have the necessary added complication that either the basalt or the rhyolite or both become contaminated after the two magma types separated, have problems accounting for the lack of igneous rocks of intermediate compositions or production of such large volumes of rhyolitic material (∼5000 km 3 ), and fail to explain why rhyolitic magma is not a more common occurrence in the ocean basin. We appeal to bouyancy of rhyolites to generate a barrier for basalt magma migration and account for the great preponderance of rhyolite relative to basalt at the surface. Furthermore, the complex isotopic picture in the rhyolites indicates that many of these magmas interacted with the upper crustal geologic units that they traversed. The interactions involved diverse processes, probably including reacton with hydrothermal fluids or hydrothermally altered rocks at high levels as well as by contamination with Phanerozoic sedimentary and Precambrian crystalline rocks at deeper levels. At the very least, we feel our study adds a cautionary note to the currently increasingly popular hypothesis that differentiation of basalt or gabbro magmas to rhyolite or granite (as distinct from tonalite or dacite) is a common occurrence and is therefore an important continential building process. Models for formation of rhyolite and granite predominantly by reworking of crust (anatexis) must still be considered. The primitive Archean mantle of the region was characterized by higher Rb/Sr, U/Pb, and Th/U values than are typical of modern suboceanic mantle. The mantle residuum within the continental subcrustal lithosperic ‘keel’ that resulted from the Archean crustal differentiation event probably was depleted in Rb/Sr and U/Pb, and the crust was correspondingly enriched in these ratios. The crust probably was further differentiated by an Archean high-grade metamorphism, during or after the primary event, into a granulitic lower crust depleted in U/Pb and Rb/Sr and a lower-grade upper crust enriched in these ratios.

Journal of Geophysical Research Solid Earth↗