USGS ScienceSearch

Geology topics

J. R. Budahn

Publications and source records attributed to J. R. Budahn.

17 recordsLinked to original sources

Pillow basalts of the Angayucham terrane: Oceanic plateau and island crust accreted to the Brooks Range

The Angayucham Mountains (north margin of the Yukon-Koyukuk province) are made up of an imbricate stack of four to eight east-west trending, steeply dipping, fault slabs composed of Paleozoic (Devonian to Mississippean), Middle to Late Triassic, and Early Jurassic oceanic upper crustal rocks (pillow basalt, subordinate diabase, basaltic tuff, and radiolarian chert). Field relations and geochemical characteristics of the basaltic rocks suggest that the fault slabs were derived from an oceanic plateau or island setting and were emplaced onto the Brooks Range continental margin. The basalts are variably metamorphosed to prehnite-pumpellyite and low-greenschist facies. Major element analyses suggest that many are hypersthene-normative olivine tholeiites. Classification based on immobile trace elements confirms the tholeiitic character of most of the basalts but suggests that some had primary compositions transitional to alkali basalt. Although field and petrographic features of the basalts are similar, trace element characteristics allow definition of geographically distinct suites. A central outcrop belt along the crest of the mountains is made up of basalt with relatively flat rare earth element (REE) patterns. This belt is flanked to the north and south by LREE (light rare earth element)-enriched basalts. Radiolarian and conodont ages from interpillow and interlayered chert and limestone indicate that the central belt of basalts is Triassic in age, the southern belt is Jurassic in age, and the northern belt contains a mixture of Paleozoic and Mesozoic ages. Data for most of the basalts cluster in the “within-plate basalt” fields of trace element discriminant diagrams; none have trace-element characteristics of island arc basalt. The Triassic and Jurassic basalts are geochemically most akin to modern oceanic plateau and island basalts. Field evidence also favors an oceanic plateau or island setting. The great composite thickness of pillow basalt probably resulted from obduction faulting, but the lack of fault slabs of gabbro or peridotite suggests that obduction faults did not penetrate below oceanic layer 2, a likely occurrence if layer 2 were anomalously thick, as in the vicinity of an oceanic island. The presence of basaltic tuff interbeds indicates proximity to an explosive basaltic eruptive center. The juxtaposition of submarine basalts of differing chemical affinity and age, adjacent to higher-grade Paleozoic metamorphic rocks of the Brooks Range to the north, may be explained by obduction of internally complex (thickened) oceanic crust formed in an ocean plateau setting. Emplacement and rotation of thrust plates to steep attitudes occurred during accretion of the Brooks Range passive margin, probably beginning in the Late to Middle Jurassic.

Journal of Geophysical Research Solid Earth

Radionuclides, trace elements, and radium residence in phosphogypsum of Jordan

Voluminous stockpiles of phosphogypsum (PG) generated during the wet process production of phosphoric acid are stored at many sites around the world and pose problems for their safe storage, disposal, or utilization. A major concern is the elevated concentration of long-lived 226Ra (half-life = 1,600 years) inherited from the processed phosphate rock. Knowledge of the abundance and mode-of-occurrence of radium (Ra) in PG is critical for accurate prediction of Ra leachability and radon (Rn) emanation, and for prediction of radiation-exposure pathways to workers and to the public. The mean (??SD) of 226Ra concentrations in ten samples of Jordan PG is 601 ?? 98 Bq/kg, which falls near the midrange of values reported for PG samples collected worldwide. Jordan PG generally shows no analytically significant enrichment (< 10%) of 226Ra in the finer (< 53 ??m) grain size fraction. Phosphogypsum samples collected from two industrial sites with different sources of phosphate rock feedstock show consistent differences in concentration of 226Ra and rare earth elements, and also consistent trends of enrichment in these elements with increasing age of PG. Water-insoluble residues from Jordan PG constitute <10% of PG mass but contain 30-65% of the 226Ra. 226Ra correlates closely with Ba in the water-insoluble residues. Uniformly tiny (< 10 ??m) grains of barite (barium sulfate) observed with scanning electron microscopy have crystal morphologies that indicate their formation during the wet process. Barite is a well-documented and efficient scavenger of Ra from solution and is also very insoluble in water and mineral acids. Radium-bearing barite in PG influences the environmental mobility of radium and the radiation-exposure pathways near PG stockpiles. ?? 2010 US Government.

Environmental Geochemistry and Health

Compositions of modern dust and surface sediments in the Desert Southwest, United States

Modern dusts across southwestern United States deserts are compositionally similar to dust-rich Av soil horizons (depths of 0-0.5 cm and 1-4 cm at 35 sites) for common crustal elements but distinctly different for some trace elements. Chemical compositions and magnetic properties of the soil samples are similar among sites relative to dust sources, geographic areas, and lithologic substrates. Exceptions are Li, U, and W, enriched in Owens Valley, California, and Mg and Sr, enriched in soils formed on calcareous fan gravel in southeast Nevada. The Av horizons are dominated by dust and reflect limited mixing with substrate sediments. Modern dust samples are also similar across the region, except that Owens Valley dusts are higher in Mg, Ba, and Li and dusts both there and at sites to the north on volcanic substrates are higher in Sb and W. Thus, dust and Av horizons consist of contributions from many different sources that are well mixed before deposition. Modern dusts contain significantly greater amounts of As, Cd, Cr, Cu, Ni, Pb, and Sb than do Av horizons, which record dust additions over hundreds to thousands of years. These results suggest that modern dust compositions are influenced by anthropogenic sources and emissions from Owens (dry) Lake after its artificial desiccation in 1926. Both modern dusts and Av horizons are enriched in As, Ba, Cu, Li, Sb, Th, U, and W relative to average crustal composition, which we interpret to indicate that the geologic sources of dust in the southwestern United States are geochemically distinctive.

Journal of Geophysical Research F: Earth Surface

Geochemical evidence for African dust inputs to soils of western Atlantic islands: Barbados, the Bahamas, and Florida

We studied soils on high-purity limestones of Quaternary age on the western Atlantic Ocean islands of Barbados, the Florida Keys, and the Bahamas. Potential soil parent materials in this region, external to the carbonate substrate, include volcanic ash from the island of St. Vincent (near Barbados), volcanic ash from the islands of Dominica and St. Lucia (somewhat farther from Barbados), the fine-grained component of distal loess from the lower Mississippi River Valley, and wind-transported dust from Africa. These four parent materials can be differentiated using trace elements (Sc, Cr, Th, and Zr) and rare earth elements that have minimal mobility in the soil-forming environment. Barbados soils have compositions that indicate a complex derivation. Volcanic ash from the island of St. Vincent appears to have been the most important influence, but African dust is a significant contributor, and even Mississippi River valley loess may be a very minor contributor to Barbados soils. Soils on the Florida Keys and islands in the Bahamas appear to have developed mostly from African dust, but Mississippi River valley loess may be a significant contributor. Our results indicate that inputs of African dust are more important to the genesis of soils on islands in the western Atlantic Ocean than previously supposed. We hypothesize that African dust may also be a major contributor to soils on other islands of the Caribbean and to soils in northern South America, central America, Mexico, and the southeastern United States. Dust inputs to subtropical and tropical soils in this region increase both nutrient-holding capacity and nutrient status and thus may be critical in sustaining vegetation. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface

Geochemical evidence for the origin of late Quaternary loess in central Alaska

Loess is extensive in central Alaska, but there are uncertainties about its source and the direction of paleo-winds that deposited it. Both northerly and southerly winds have been inferred. The most likely sources of loess are the Tanana River (south), the Nenana River (southeast), and the Yukon River (north). Late Quaternary loess in central Alaska has immobile trace-element compositions (Cr/Sc, Th/Ta, Th/ Sc, Th/U, Eu/Eu*, GdN/YbN) that indicate derivation mostly from the Tanana River. However, other ratios (As/Sb, Zr/Hf, LaN/YbN) and quantitative modeling indicate that the Yukon River was also a source. During the last glacial period, there may have been a longer residence time of the Siberian and Canadian high-pressure cells, along with a strengthened Aleutian low-pressure cell. This would have generated regional-scale northeasterly winds and explains derivation of loess from the Yukon River. However, superim-posed upon this synoptic-scale circulation, there may have been strong, southerly katabatic winds from expanded glaciers on the northern flank of the Alaska Range. These winds could have provided eolian silt from the Tanana River. Yukon River and Tanana River sediments are highly calcareous, whereas Fairbanks-area loess is not. This suggests that carbonate leaching in loess kept ahead of sedimentation and that late Quaternary loess in central Alaska was deposited relatively slowly. ?? 2006 NRC Canada.

Canadian Journal of Earth Sciences

The 2003 phreatomagmatic eruptions of Anatahan volcano - Textural and petrologic features of deposits at an emergent island volcano

Stratigraphic and field data are used in conjunction with textural and chemical evidence (including data from scanning electron microscope, electron microprobe, X-ray fluorescence, X-ray diffraction, and instrumental neutron activation analysis) to establish that the 2003 eruption of Anatahan volcano was mainly phreatomagmatic, dominated by explosive interaction of homogeneous composition low-viscosity crystal-poor andesite magma with water. The hydromagmatic mode of eruption contributed to the significant height of initial eruptive columns and to the excavation and eruption of altered rock debris from the sub-volcanic hydrothermal system. Volatile contents of glass inclusions in equilibrium phenocrysts less abundances of these constituents in matrix glass times the estimated mass of juvenile magma indicate minimum emissions of 19 kt SO2 and 13 kt Cl. This petrologic estimate of SO2 emission is an order-of-magnitude less than an estimate from TOMS. Similarly, inferred magma volumes from the petrologic data are an order of magnitude greater than those modeled from deformation data. Both discrepancies indicate additional sources of volatiles, likely derived from a separate fluid phase in the magma. The paucity of near-source volcanic-tectonic earthquakes preceding the eruption, and the dominance of sustained long-period tremor are attributed to the ease of ascent of the hot low-viscosity andesite, followed by a shallow phreatomagmatic mode of eruption. Phreatomagmatic eruptions are probably more common at emergent tropical island volcanoes, where shallow fresh-water lenses occur at near-sea-level vents. These relations suggest that phreatomagmatic explosions contributed to the formation of many of the near-sea-level craters and possibly even to the small calderas at the other Mariana islands.

Northern Mariana Islands

Rear-arc vs. arc-front volcanoes in the Katmai reach of the Alaska Peninsula: A critical appraisal of across-arc compositional variation

Physical and compositional data and K-Ar ages are reported for 14 rear-arc volcanoes that lic 11-22 km behind the narrowly linear volcanic front defined by the Mount Katmai-to-Devils Desk chain on the Alaska Peninsula. One is a 30-km3 stratocone (Mount Griggs; 51-63% SiO2) active intermittently from 292 ka to Holocene. The others are monogenetic cones, domes, lava flows, plugs, and maars, of which 12 were previously unnamed and unstudied; they include seven basalts (48-52% SiO2), four mafic andesites (53-55% SiO2), and three andesite-dacite units. Six erupted in the interval 500-88 ka, one historically in 1977, and five in the interval 3-2 Ma. No migration of the volcanic front is discernible since the late Miocene, so even the older units erupted well behind the front. Discussion explores the significance of the volcanic front and the processes that influence compositional overlaps and differences among mafic products of the rear-arc volcanoes and of the several arc-front edifices nearby. The latter have together erupted a magma volume of about 200 km3, at least four times that of all rear-arc products combined. Correlation of Sr-isotope ratios with indices of fractionation indicates crustal contributions in volcanic-front magmas (0.7033-0.7038), but lack of such trends among the rear-arc units (0.70298-0.70356) suggests weaker and less systematic crustal influence. Slab contributions and mantle partial-melt fractions both appear to decline behind the front, but neither trend is crisp and unambiguous. No intraplate mantle contribution is recognized nor is any systematic across-arc difference in intrinsic mantle-wedge source fertility discerned. Both rear-arc and arc-front basalts apparently issued from fluxing of typically fertile NMORB-source mantle beneath the Peninsular terrane, which docked here in the Mesozoic. Springer-Verlag 2004.

Alaska

Isotopic, geochemical, and temporal characterization of Proterozoic basement rocks in the Quitovac region, northwestern Sonora, Mexico: Implications for the reconstruction of the southwestern margin of Laurentia

A detailed geochemical characterization of 19 representative Proterozoic basement rocks in the Quitovac region in northwestern Sonora, Mexico, has identified two distinct Paleoproterozoic basement blocks that coincide spatially with the previously proposed Caborca and "North America" blocks. New U-Pb zircon geochronology revises their age ranges, the Caborca (1.78-1.69 Ga) and "North America" (1.71-1.66 Ga) blocks at Quitovac, and precludes a simple age differentiation between them. In addition, Grenvillian-age granitoids (ca. 1.1 Ga), spatially associated with the Caborca block have been identified at Quitovac. Nd isotopes and major- and trace-element geochemistry support the distinction of these Paleoproterozoic blocks. Granitoids of the "North America" block are characterized by depleted εNd values (3.4-3.9) and younger Nd model ages (1800-1740 Ma) and have lower K2O, Y, Rb, Ba, Th, REE, and Fe/Mg values than coeval rocks of the Caborca block. The Caborca block granitoids are likewise characterized by slightly less depleted εNd (0.6-2.6) and older Nd model ages (2070-1880 Ma). Despite the subtle differences, granitoids from both the Caborca and "North America" blocks exhibit island arc-like affinities. We propose that the Proterozoic basement rocks from the Quitovac region are an extension of the Proterozoic crustal provinces in the southwestern United States. Specifically, rocks of the Caborca block exhibit an affinity to rocks of either the Yavapai province or the Mojave-Yavapai transition zone, whereas rocks of the "North America" block have signatures similar to those of the Mazatzal province or possibly the Yavapai province of Arizona. The new isotopic ages and geochemical data do not support the existence of the Late Jurassic Mojave-Sonora megashear at Quitovac, as originally proposed. However, the Quitovac region accounts only for a small fraction of the Proterozoic basement in Sonora, so these findings do not eliminate the possibility of a megashear elsewhere in northern Sonora. Our new data create the possibility of alternative hypotheses for the distribution of Paleoproterozoic crustal provinces in southwestern North America that affect reconstructions of the original southwestern margin of Laurentia, and reduce uncertainties in the configuration, timing, and existence of the Proterozoic supercontinent, Rodinia.

Quitovac

Kilauea east rift zone magmatism: An episode 54 perspective

On January 29 30, 1997, prolonged steady-state effusion of lava from Pu'u'O'o was briefly disrupted by shallow extension beneath Napau Crater, 1 4 km uprift of the active Kilauea vent. A 23-h-long eruption ( episode 54 ) ensued from fissures that were overlapping or en echelon with eruptive fissures formed during episode 1 in 1983 and those of earlier rift zone eruptions in 1963 and 1968. Combined geophysical and petrologic data for the 1994 1999 eruptive interval, including episode 54 , reveal a variety of shallow magmatic conditions that persist in association with prolonged rift zone eruption. Near-vent lava samples document a significant range in composition, temperature and crystallinity of pre-eruptive magma. As supported by phenocryst liquid relations and Kilauea mineral thermometers established herein, the rift zone extension that led to episode 54 resulted in mixture of near-cotectic magma with discrete magma bodies cooled to ≤1100°C. Mixing models indicate that magmas isolated beneath Napau Crater since 1963 and 1968 constituted 32 65% of the hybrid mixtures erupted during episode 54 . Geophysical measurements support passive displacement of open-system magma along the active east rift conduit into closed-system rift -reservoirs along a shallow zone of extension. Geophysical and petrologic data for early episode 55 document the gradual flushing of episode 54 related magma during magmatic recharge of the edifice.

Hawaii

Use of radium isotopes to determine the age and origin of radioactive barite at oil-field production sites

Radium-bearing barite (radiobarite) is a common constituent of scale and sludge deposits that form in oil-field production equipment. The barite forms as a precipitate from radium-bearing, saline formation water that is pumped to the surface along with oil. Radioactivity levels in some oil-field equipment and in soils contaminated by scale and sludge can be sufficiently high to pose a potential health threat. Accurate determinations of radium isotopes (226Ra+228Ra) in soils are required to establish the level of soil contamination and the volume of soil that may exceed regulatory limits for total radium content. In this study the radium isotopic data are used to provide estimates of the age of formation of the radiobarite contaminant. Age estimates require that highly insoluble radiobarite approximates a chemically closed system from the time of its formation. Age estimates are based on the decay of short-lived 228Ra (half-life=5.76 years) compared to 226Ra (half-life=1600 years). Present activity ratios of 228Ra/226Ra in radiobarite-rich scale or highly contaminated soil are compared to initial ratios at the time of radiobarite precipitation. Initial ratios are estimated by measurements of saline water or recent barite precipitates at the site or by considering a range of probable initial ratios based on reported values in modern oil-field brines. At sites that contain two distinct radiobarite sources of different age, the soils containing mixtures of sources can be identified, and mixing proportions quantified using radium concentration and isotopic data. These uses of radium isotope data provide more description of contamination history and can possibly address liability issues. Copyright ?? 2000 .

Environmental Pollution

Late Cenozoic regional collapse due to evaporite flow and Dissolution in the Carbondale Collapse Center, West-Central Colorado

Dissolution and flow of Pennsylvanian evaporitic rocks in west-central Colorado created the Carbondale Collapse Center, a 450 mi2 structural depression with about 4,000 ft of vertical collapse during the late Cenozoic. This paper describes evidence of collapse in the lower Roaring Fork River valley. Both the lateral extent and amount of vertical collapse is constrained by deformed upper Cenozoic volcanic rocks that have been correlated using field mapping, 40Ar/39Ar geochronology, geochemistry, and paleomagnetism. The Carbondale Collapse Center is one of at least two contiguous areas that have experienced major evaporite tectonism during the late Cenozoic. Historic sinkholes, deformed Holocene deposits, and modern high-salinity loads in the rivers and thermal springs indicate the collapse process continues today. Flow of evaporitic rocks is an important element in the collapse process, and during initial stages of collapse it was probably the primary causative mechanism. Dissolution, however, is the ultimate means by which evaporite is removed from the collapse area. As the Roaring Fork River began to rapidly down-cut through a broad volcanic plateau during the late Miocene, the underlying evaporite beds were subjected to differential overburden pressures. The evaporitic rocks flowed from beneath the upland areas where overburden pressures remained high, toward the Roaring Fork River Valley where the pressures were much lower. Along the valley the evaporitic rocks rose upward, sometimes as diapirs, forming or enhancing a valley anticline in bedrock and locally upwarping Pleistocene terraces. Wherever the evaporites encountered relatively fresh ground water, they were dissolved, forming underground voids into which overlying bedrock and surficial deposits subsided. The saline ground water eventually discharged to streams and rivers through thermal springs and by seepage into alluvial aquifers.

Mountain Geologist

Elemental analyses of modern dust in southern Nevada and California

Selected samples of modern dust collected in marble traps at sites in southern Nevada and California (Reheis and Kihl, 1995; Reheis, 1997) have been analyzed for elemental composition using instrumental neutron activation analysis (INAA) and inductively coupled plasma atomic emission spectroscopy (ICP-AES) and inductively coupled plasma mass spectroscopy (ICP-MS). For information on these analytical techniques and their levels of precision and accuracy, refer to Baedecker and McKown (1987) for INAA, to Briggs (1996) for ICP-AES, and to Briggs and Meier (1999) for ICP-MS. This report presents the elemental compositions obtained using these techniques on dust samples collected from 1991 through 1997. The dust-trap sites were established at varying times; some have been maintained since 1984, others since 1991. For details on site location, dust-trap construction, and collection techniques, see Reheis and Kihl (1995) and Reheis (1997). Briefly, the trap consists of a coated angel-food cake pan painted black on the outside and mounted on a post about 2 m above the ground. Glass marbles rest on a circular piece of galvanized hardware cloth (now replaced by stainless-steel mesh), which is fitted into the pan so that it rests 3-4 cm below the rim. The 2-m height eliminates most saltating sand-sized particles. The marbles simulate the effect of a gravelly fan surface and prevent dust that has filtered or washed into the bottom of the pan from being blown back out. The dust traps are fitted with two metal straps looped in an inverted basket shape; the top surfaces of the straps are coated with a sticky material that effectively discourages birds from roosting.

California, Nevada

Actualistic models of mantle metasomatism documented in a composite xenolith from Dish Hill, California

Major and trace-element whole rock and mineral variations in composite hornblendite-peridotite xenolith Ba-2-1, from Dish Hill, CA, are due to a single event of metasomatism in the mantle. The hornblendite is the crystallized selvage of a dike conduit charged with incompatible-element-enriched hydrous mafic magma. The magma infiltrated the refractory peridotite wallrock, reacted with its constituent minerals, and simultaneously deposited amphibole. The systematic data from this study show considerable variation in isotopic values and trace elements. These data provide insight into a mantle process that was defined previously from samples without context, lacking evidence about the number or source of metasomatic events. In the contact zone of Ba-2-1, peridotite is enriched in Fe, Ti, CO 2 ) and H 2 O; clinopyroxene and amphibole also are enriched in Fe and Ti, but clinopyroxene appears slightly depleted in CaO. Compared to chondrites, peridotite, clinopyroxene, and probably amphibole are enriched in light rare earth (LREE cn ) and other incompatible trace elements. Values of 87 Sr 86 Sr "> 87 Sr 86 Sr and 143 Nd 144 Nd "> 143 Nd 144 Nd in the contact zone are close to isotopic equilibrium with the dike. Whole rock and constituent clinopyroxene compositions change to those of refractory peridotite with distance from the contact. These compositional variations were modelled using Gresens' equation for whole-rock major and minor elements, and calculations for isotopic ratios and REEs, which emulate the effects of Chromatographic fractionation. The choice of endmembers was restricted to compositions actually present in mantle samples from Dish Hill. Model results indicate that: 1. the variations can be explained as the result of a single metasomatic event, probably a single pulse of previously fractionated liquid; 2. the ratio of total interacting liquid to peridotite was at least 1:3 by weight in the contact zone; and 3. the composition of the metasomatic liquid changed progressively as it infiltrated beyond that zone. The small distance over which variations occur is due to the small amount of liquid that infiltrated. Only in the contact zone was peridotite wallrock saturated by a liquid composition similar to the dike. Comparison of the Ba-2-1 data with those of another xenolith from Dish Hill suggests that the compositional variations of mantle metasomatism result from both the compositional contrast between the metasomatizing liquid and wallrock and the relative abundances of each. Compositional and volumetric variations of mantle partial melts and their fractionates, and repeated events of melting and reaction in contiguous mantle, can create broad ranges of metasomatic “signatures” from the same process.

California

Back-arc with frontal-arc component origin of Triassic Karmutsen basalt, British Columbia, Canada

The largely basaltic, ???4.5-6.2-km-thick, Middle to Upper Triassic Karmutsen Formation is a prominent part of the Wrangellian sequence. Twelve analyses of major and minor elements of representative samples of pillowed and massive basalt flows and sills from Queen Charlotte and Vancouver Islands are ferrotholeiites that show a range of 10.2-3.8% MgO (as normalized, H2O- and CO2-free) and related increases in TiO2 (1.0-2.5%), Zr (43-147 ppm) and Nb (5-16 ppm). Other elemental abundances are not related simply to MgO: distinct groupings are evident in Al2O3, Na2O and Cr, but considerable scatter is present in FeO* (FeO + 0.9Fe2O3) and CaO. Some of the variation is attributed to alteration during low-rank metamorphism or by seawater - including variation of Ba, Rb, Sr and Cu, but high-field-strength elements (Sc, Ti, Y, Zr and Nb) as well as Cr, Ni, Cu and rare-earth elements (REE's) were relatively immobile. REE's show chondrite-normalized patterns ranging from light-REE depleted to moderately light-REE enriched. On eleven discriminant plots these analyses fall largely into or across fields of within-plate basalt (WIP), normal or enriched mid-ocean-ridge tholeiite (MORB) and island-arc tholeiite (IAT). Karmutsen basalts are chemically identical to the stratigraphically equivalent Nikolai Greenstone of southern Alaska and Yukon Territory. These data and the fact that the Karmutsen rests on Sicker Group island-arc rocks of Paleozoic age suggest to us that: 1. (1) the basal arc, after minor carbonate-shale deposition, underwent near-axial back-arc rifting (as, e.g., the Mariana arc rifted at different times); 2. (2) the Karmutsen basalts were erupted along this rift or basin as "arc-rift" tholeiitite; and 3. (3) after subsequent deposition of carbonates and other rocks, and Jurassic magmatism, a large fragment of this basalt-sediment-covered island arc was accreted to North America as Wrangellia. The major- and minor-elemental abundances of Karmutsen basalt is modeled by first mixing primitive arc magma with enriched basaltic liquid derived either from garnet peridotite or metasomatized mantle, followed by fractionation of olivine, pyroxenes, plagioclase and spinel. ?? 1989.

Chemical Geology