SEARCH · Geology Reports

Results for “Environmental Modeling and Software”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,663 records · Page 86Linked to original sources

USGS1976

Results of phytoplankton sampling at National Stream Quality Accounting Network stations in Montana, 1975 water year

Twelve National Stream Quality Accounting Network stations were operated in Montana during the 1975 water year. The network was established for the purpose of acquiring a base of hydrologic data for use by agencies engaged in water-resources planning on a national or regional scale. Among the characteristics analyzed were phytoplankton identification and cell counts. Samples consisted of composites of equal aliquots, collected at the center of each quartile of flow, using modified suspended-sediment samplers and sediment collection techniques. Identification and counting were done using the Sedgwick-Rafter cell method. Cell counts ranged from 21 cells per millilitre at Flathead River at Flathead, British Columbia to 27,000 cells per millilitre at Yellowstone River near Sidney. The class Bacillariophyceae was most abundant in both number and variety at all sampling sites. Anabaena and Aphanizomenon of the phyllum Cyanophyta were found at six stations and two stations respectively. These two genera of blue-green algae often become abundant in enriched waters resulting in nuisance conditions.

Montana↗

USGS2024

Correction to A regime shift in sediment export from a coastal watershed during a record wet winter, California: Implications for landscape response to hydroclimatic extremes

In the referenced article, the authors would like to correct text in the first paragraph on page 2571, Figure 9 and its caption. The changes reflect an error made in the processing of the rainfall intensity-duration data used to compare storms to published debris flow triggering thresholds. The correctly processed data does not change the interpretations made in the paper but does correctly indicate that the investigated storms did not exceed the rainfall intensity – duration threshold of Cannon (1988) but did significantly exceed the debris flow triggering threshold of Wieczorek (1987).

Earth Surface Processes and Landforms↗

USGS2024

Unscrambling the Proterozoic supercontinent record of northeastern Washington State, USA

The time interval from Supercontinent Nuna assembly in the late Paleoproterozoic to Supercontinent Rodinia breakup in the Neoproterozoic is considered by some geologists to comprise the “Boring Billion,” an interval possibly marked by a slowdown in plate tectonic processes. In northeastern Washington State, USA, similar to much of western Laurentia, early workers generally thought the tectonostratigraphic framework of this interval of geologic time consisted of two major sequences, the (ca. 1480–1380 Ma) Mesoproterozoic Belt Supergroup and unconformably overlying (<720 Ma) Neoproterozoic Windermere Supergroup. However, recent research indicates that strata considered by early workers as Belt Supergroup equivalents are actually younger, and a post-Belt, pre-Windermere record is present within the <1360 Ma Deer Trail Group and <760 Ma Buffalo Hump Formation. Thus, the northeastern Washington region perhaps comprises the most complete stratigraphic record of the “Boring Billion” time interval in the northwestern United States and holds important insights into global Proterozoic supercontinent tectonic processes. In light of these exciting developments, this field guide will address the early historic economic geology and original mapping of these Proterozoic sequences in the northeastern Washington region, and from that foundation explore more recent isotopic provenance data and their regional to global context. Finally, the guide will end with a discussion of remaining questions with a goal of stimulating interest in these relatively understudied, yet important, rocks.

Washington↗

USGS2024

Aboveground carbon stocks across a hydrological gradient: Ghost forests to non-tidal freshwater forested wetlands

Upper estuarine forested wetlands (UEFWs) play an important role in the sequestration of atmospheric carbon (C), which is facilitated by their position at the boundary of terrestrial and maritime environments but threatened by sea level rise. This study assessed the change in aboveground C stocks along the estuarine–riverine hydrogeomorphic gradient spanning salt-impacted freshwater tidal forested wetlands to freshwater forested wetlands in seasonally tidal and nontidal landscape positions. Standing stocks of C in forested wetlands were measured along two major coastal river systems, the Winyah Bay in South Carolina and the Savannah River in Georgia (USA), replicating and expanding a previous study to allow the assessment of change over time. Aboveground C stocks on these systems averaged 172.9 Mg C ha −1 , comparable to those found in UEFWs across the globe and distinct from the terrestrial forested ecosystems they are often considered to be a part of during large-scale C inventory efforts. Groundwater salinity conditions as low as 1.3 ppt were observed in conjunction with losses of aboveground C. When viewed in context alongside expected sea level rise and corresponding saltwater intrusion estimates, these data suggest a marked decrease in aboveground C stocks in forested wetlands situated in and around tidal estuaries.

Georgia, South Carolina↗

USGS2026

Top-down targeted network analysis of critical mineral commodities applied to international geochemistry database

The global demand for critical mineral commodities is rapidly increasing, making domestic production an important factor in supporting the economy and national security. Large scale, publicly available geochemical databases allow for the application of data informatics methods to interrogate critical mineral commodities data for correlations in deposit formation and distribution, particularly for identifying enrichment of multiple critical mineral commodities at the same deposit. In this study, we applied network analysis to the Critical Minerals Mapping Initiative (CMMI) ore geochemistry (Critical Minerals in Ores, CMiO) database to identify the high concentration (defined as 10× bulk crustal abundance) co-occurrence of different critical mineral commodities across a mineral system hierarchy from deposit environments to individual deposits. Identifying patterns or unique outliers in enrichment in network communities will allow for the location of secondary critical mineral commodity resources from under-utilized deposits. We find trends in the enrichment of critical mineral commodities in network-communities between the elements praseodymium (Pr), neodymium (Nd), terbium (Tb), and dysprosium (Dy) across multiple CMiO database deposit environments and groups down to specific deposit types and sites. A separate trend in network community deposition is observed as well between iridium (Ir) and platinum (Pt) in deposit environments, groups, types, and sites. Network analysis focused on critical minerals in magmatic-hydrothermal deposits identified multiple deposit sites from different deposit types within the CMiO database with concentrations of Dy, Nd, Tb, Pr, Ir, and Pt that are at least ten times greater than the crustal average. This approach can be applied to any target element(s) or deposit(s) of interest, allowing broad investigation of co-enriched critical mineral commodities.

Journal of Geochemical Exploration↗

USGS2023

Thickness of the Saudi Arabian crust

As part of a joint Saudi Geological Survey (SGS) and U.S. Geological Survey (USGS) project, we analyzed P-wave receiver functions from seismic stations covering most of the Kingdom of Saudi Arabia to map the thickness of the crust across the Arabia Plate. We present an update of crustal-thickness estimates and fill in gaps for the western Arabian Shield and the rifted margin at the Red Sea (the coastal plain), as well as the eastern Arabian Platform. We applied a conventional H-k stacking algorithm and included careful attention to stacking weights, two forms of sedimentary corrections for stations located on the Arabian Platform, and additional processing for noisy stations. We obtained useful results at 154 stations from 898 teleseismic events over a 2-year period from 1995–1997 (for non-SGS stations) and a 6-year period from 2008–2014 (for SGS stations). Average crustal thickness (that is, depth to the Mohorovičić discontinuity [Moho] below the surface) beneath the Red Sea coastal plain (the rift margin) is 29 kilometers (km), beneath the volcanic fields (known in Arabic as harra [plural] or harrat [singular]) is 35 km, beneath the Arabian Shield (excluding harrats) is 37 km, and beneath the Arabian Platform is 38 km. Crustal thinning appears not to extend east of the rift escarpment, suggesting uniform extension that is no broader at depth than at the surface. In contrast to some previous claims that the Arabian Platform crust is thicker than that of the Arabian Shield, we find no statistically significant difference between their whole crustal thicknesses. However, the average subsedimentary crustal thickness (that is, the crystalline crust) for stations on the Arabian Platform is 34 km, 3 km thinner than the crust of the Arabian Shield. Individual station P-wave (pressure) velocity and S-wave (shear) velocity ratios ( V P / V S ) are highly variable for the Arabia Plate, ranging from 1.60 to 1.97 and averaging 1.75, with a standard deviation of 0.07. There are no statistically significant differences between V P / V S ratios of the different geologic regions of Saudi Arabia. Similar V P / V S ratios, coupled with similar crustal thicknesses for harrats and the Arabian Shield, indicate that Cenozoic magmatism has contributed negligibly to crustal growth.

Professional Paper↗

USGS2025

The influence of postfire root strength decay on shallow landslide susceptibility in western Oregon

Severe wildfires can increase the potential for postfire landslides, partly due to the loss of vegetation and root reinforcement. This study investigates the multi-year dynamics of root strength following wildfire in coniferous forests of western Oregon and evaluates implications for postfire slope stability. Root tensile strength was measured through laboratory testing of 969 roots collected from 26 test pits over four postfire timeframes: unburned conditions, 4–12 months postfire (Cedar Creek, 2022), 24–35 months postfire (Holiday Farm, 2020), and 59–67 months postfire (Eagle Creek, 2017). Variables analyzed include root diameter, tensile thread strength, stiffness, progressive tensile strength with displacement, time-dependent strength and ductility, root area ratio, and cohesion. Results indicate a 50% reduction in root strength after wildfire, with the lowest values occurring around 4 years postfire. Projections indicate it could take 10–22 years following fire for root strength to recover to 70% and 90% of the prefire original root strength, respectively. This timeline indicates a multi-year window of vulnerability for shallow landslides that reaches a maximum approximately 4 years following fire. Brittle failure behavior was observed in burned root systems compared to unburned root systems, signaling a more abrupt loss of strength at yield. Simple slope stability analyses show that reduced root strength can lead to instability on more gentle slopes and under less saturation in comparison to unburned conditions. These findings highlight the critical role of root reinforcement in postfire slope stability and the long-term implications of wildfire disturbance on landslide susceptibility.

Oregon↗

USGS2025

The U.S. Geological Survey National Water Quality Network—Groundwater, 2024

The U.S. Geological Survey (USGS) operates a National Water Quality Network (NWQN) to monitor trends in groundwater quality and assess emerging contaminants of concern. It is a “network of networks” with 81 subnetworks being sampled on a decadal time scale. Each year, 8 of the subnetworks are sampled. Subnetworks have 20–30 wells each and include studies of domestic supply wells or shallow groundwater (20–50 feet deep) underlying urban land use or agricultural land use. Currently there are 2,089 wells in the network. All wells are sampled for physical properties, nutrients, major ions, trace elements, per- and polyfluoroalkyl substances (PFAS), and a subset of wells are sampled for pesticides, volatile organic compounds, radionuclides, and microbiological contaminants.

Contiguous United States↗

USGS2026

Rapid hydrothermal triggering of induced seismicity at the Coso Geothermal Field

The long-term producing Coso Geothermal Field (CGF) in California operates over 100 wells tapping into a reservoir characterized by an extensive fracture network, complex fluid pathways, and regular seismic activity. Understanding the interaction between seismicity and injection can shed important light on the hydrothermal characteristics of the field. Here, we analyze 15 years of local seismic and daily operational data from the CGF, identifying a strong correlation between short-term increase in seismicity rate and seasonal volumetric and temperature variations in the reinjected fluid. Furthermore, the seismic footprint during peak injection of colder fluids reveals a near-instantaneous response up to 2 km away from the injection well, too rapid for pore pressure diffusion alone. This short-term and distant response is observed to have directional preference, indicating structural or permeability anisotropy within the reservoir. Additionally, the seismic response correlates with the initial volumetric increase of colder fluids, but also with temperature decrease during stable injection periods, suggesting thermal effects alone can play an important role in triggering distant seismicity.

California↗

USGS2026

Towards a refined definition of “tidal influence” for the coastal wetland sciences

The astronomical tide relates to lunar (and to a lesser extent solar) influences on rising and falling water bodies, and is a key control of many coastal processes. However, tides manifest along coastal margins in different ways to muddle application beyond Sir Isaac Newton’s first descriptions. Here, we argue that tide, or specifically “tidal,” has not been adequately defined for the coastal wetland sciences to facilitate contemporary interdisciplinary distinction for rapidly changing coastlines and future climate adaptation planning. We provide a brief history of the tidal concept, and we describe tides in a diversity of uncommonly recognized marginal coastal environments. Common astronomical tide manifestations are described through hydrographic examples, then expanded to include those tidal environments exposed additionally to fluvial, meteorological (wind), or groundwater influences. Our expanded definition of “tidal” is intended to be a forward-looking construct as environments will inevitably change with sea-level rise, water extraction, and anthropogenic flow alteration. Our definition of wetland tidal influence, which we broaden, includes those “wetlands affected by astronomic water level fluctuations within coastally restricted environments causing frequent or infrequent surface inundation or groundwater tidal variation that affects the biogeochemistry of the soil (e.g., ion exchange capacity, oxygen state, mineralization processes)”. Astronomic water level fluctuations are often amplified or dampened by non-astronomical events, such as wind. Expanding the definition of what is considered tidal by the scientific community is critical for improved understanding of coastal wetland function and the suite of ecosystem services they provide people and society, including climate mitigation opportunity, management, and adaptation.

Wetlands↗

USGS2024

Projections of multiple climate-related coastal hazards for the US Southeast Atlantic

Faced with accelerating sea level rise and changing ocean storm conditions, coastal communities require comprehensive assessments of climate-driven hazard impacts to inform adaptation measures. Previous studies have focused on flooding but rarely on other climate-related coastal hazards, such as subsidence, beach erosion and groundwater. Here, we project societal exposure to multiple hazards along the Southeast Atlantic coast of the United States. Assuming 1 m of sea level rise, more than 70% of the coastal residents and US \ $ 1 trillion in property are in areas projected to experience shallow and emerging groundwater, 15 times higher than daily flooding. Storms increase flooding exposure by an order of magnitude over daily flooding, which could impact up to ~50% of all coastal residents and US \ $ 770 billion in property value. The loss of up to ~80% of present-day beaches and high subsidence rates that currently affect over 1 million residents will exacerbate flooding and groundwater hazard risks.

Florida, Georgia, North Carolina, South Carolina, ↗

USGS2025

3D Converted wave reverse time migration imaging

We describe a newly developed method for recovering high-resolution images of seismic discontinuities, such as subducting slabs, in 3D. Our method makes use of converted S → P or P → S waves observed by dense arrays of seismometers to infer the locations and relative strengths of seismic discontinuities at depth in a target region. Observed direct and converted waves are backpropagated to their times of origin. The time-reversed wavefield is then separated into its constituent P and S components via the Helmholtz decomposition, and those separated wavefields are used to compute imaging functions that characterize the locations and relative strengths of seismic discontinuities. Imaging functions may be designed to use either S → P or P → S waves, so that users can target those arrivals expected to be most dominant in a given dataset. We have previously demonstrated the efficacy of our method in two dimensions, and we now present a 3D implementation of our technique which addresses the significant computational challenges posed by the size of volumetric wavefield data in three dimensions. Through a series of synthetic examples, we demonstrate that our method is capable of recovering the fine scale structure of a subducting slab given realistic station coverage and earthquake sources. We investigate optimal seismic station geometries for our technique and explore image interpretability in regions with poor data coverage. We find that linear station geometries yield more optimal, interpretable imaging functions than collections of small arrays can. We also show that our method can successfully recover bothS → P or P → S images when realistic shear earthquake sources are used, and we explore the additional computational challenges presented by the high frequency content of S waves. Our results demonstrate the potential for our technique to recover high-resolution information about subducting slabs in real-world regions, given that relatively sparse seismic arrays with only approximately 100 stations are capable of recovering interpretable imaging functions from just a few realistic earthquake sources for multiple discontinuities at significant depth in an area of approximately 400~sq~km.

Seismica↗

USGS2024

On the provenance of field reports of the 1886 Charleston, South Carolina, earthquake: A seismo-historical whodunnit

Much of what is known about the effects of the 1886 Charleston, South Carolina, earthquake throughout the epicentral region can be attributed to meticulous field investigations by an individual with training in geology and engineering, Earle Sloan ( Clendenin, 1926 ). In a recent study, Bilham and Hough (2024) undertook a detailed analysis of the effects of the earthquake on railroads in the Charleston region, drawing heavily from Sloan’s reports. This exercise identified several inconsistencies in Sloan’s field reports, including understandable measurement imprecision, inferred data entry mistakes, and transcription errors. The study also begged the question, where was Sloan at the time of the mainshock and over the following week? And to what extent did he draw from secondhand information in compiling his reports? On this question Sloan’s reports were sometimes enigmatic, lending themselves to misinterpretation in contemporaneous as well as modern interpretations. Beyond the details that were germane for, and briefly summarized by, the studies of Bilham and Hough (2023 , 2024) , in this report we don our historical seismologist caps to chronicle Sloan’s activities following the earthquake. We summarize our inferences here for the benefit of future scholars who might attempt to retrace either Sloan’s footsteps or our own. This study also serves to highlight Sloan’s singular contributions to earthquake science, which were never published separately.

South Carolina↗

USGS1999

Effects of snowmobile use on snowpack chemistry in Yellowstone National Park, 1998

Snowmobile use in Yellowstone National Park has increased substantially in the past three decades. In areas of greatest snowmobile use, elevated levels of by-products of gasoline combustion such as ammonium and benzene have been detected in snowpack samples. Annual snowpacks and snow-covered roadways trap deposition from local and regional atmospheric emissions. Snowpack samples representing most of the winter precipitation were collected at about the time of maximum annual snow accumulation at a variety of locations in the park to observe effects of a range of snowmobile traffic levels. Concentrations of organic and inorganic compounds in snow samples from pairs of sites located directly in and off snow-packed roadways used by snowmobiles were compared to concentrations in samples collected at nearby off-road sites. Concentrations of ammonium were 2 to 5 times higher for the in-road snow compared to off-road snow for each pair of sites. Thus, concentrations decreased rapidly with distance from roadways. In addition, concentrations of ammonium, nitrate, sulfate, benzene, and toluene in snow were positively correlated with snowmobile use.

Wyoming↗

USGS2026

Conventional hydrothermal power-producing systems of the Great Basin, USA

As part of the update to the electric-grade conventional hydrothermal assessment of the Great Basin, USA, Monte Carlo analyses of identified resources within explored regions will be performed to make estimates of discovered resources and associated uncertainty. Analyses use conditional statistics where estimates are conditioned upon a hydrothermal favorability map, allowing for the likelihood that more resources exist in regions of higher hydrothermal favorability. For these analyses, a dataset of identified hydrothermal systems is compiled, and the new compilation is described herein. Recognizing that a single hydrothermal system may be developed with multiple power plants, and that the hydrothermal upflow zone may be several kilometers across with many measurements characterizing a single hydrothermal system, a procedure was developed and employed to create clusters of points (power plants, measurements, etc.) that are associated with a single system, and a new central point was defined as the best estimator of the center of the hydrothermal system. Hydrothermal systems were uniquely identified by grouping electric-grade hydrothermal measurements and operating power plants within a distance of 10 km. Groups that are >10 km apart are assumed to be different electric-grade hydrothermal systems. While 10 km was used as the threshold, most systems were significantly further apart, and most points within groups were typically within 5 km of each other. A well measurement was considered an electric-grade measurement of a hydrothermal system if it had two properties: a measured temperature of >85 °C and evidence of hydrothermal convection. Other points that were added to the dataset are locations of operating powerplants or locations that have been classified as an electric-grade hydrothermal resource by either the U.S Geological Survey (USGS) or the Great Basin Center for Geothermal Energy. After all points are assigned to systems, new points were computed with the goal of identifying the center of the throat of the hydrothermal upflow zone. If operating powerplants exist for a system, then the arithmetic average of all power plant locations is used. Otherwise, if USGS made an estimate, that location is used. In the absence of both powerplants or USGS estimates, the arithmetic average of all electric-grade measurement locations is used. An example is shown of how these newly compiled locations might be ranked for uncertainty analyses, where higher confidence is assumed if measured temperature is higher and there are many supporting measurements indicating an electric-grade resource. In summary, 28 systems have operating power plants, an additional 78 systems are known identified electric-grade hydrothermal resources, and 99 new systems were identified as probable electric-grade systems with varying levels of confidence. These 205 locations are shown as a function of a recent hydrothermal favorability map, conceptually illustrating the conditional statistics that can be used to make estimates of the undiscovered resources of the Great Basin. An accompanying data release provides summaries of developed capacity by system and USGS estimates of likely total capacity and associated uncertainty.

Arizona, California, Idaho, Nevada, Oregon, Utah↗

USGS1976

Effect of deicing chemicals on the hydrologic environment in Massachusetts; evaluation of surface-water data collection

The objective of this study is to develop predictive relationships that can be used to describe the impact of highway deicing salts on the hydrologic environment. A method, presented in an earlier report, for estimating yearly mean chloride concentrations from estimated or actual runoff and salt-application data was applied to the period 1971 through1974 and the results compared with yearly mean chloride concentrations computed from records of specific conductance. Estimates were found to be from minus 63 percent to plus 56 percent in error when the error was expressed as the difference between estimated and computed values as a percentage of the computed values. Chloride concentrations used in this study are computed from specific conductance/chloride concentration relationships and records of specific conductance. Preliminary results from graphic analyses and least-squares regression analysis show that relationships between measured values of specific conductance and chloride concentration have correlation coefficients ranging from 0.30 to 0.97. Analysis of streamflow for all major dissolved constituents is recommended with the purpose of attempting to describe the variations in the specific conductance/chloride concentration relationships.

Massachusetts↗

USGS2026

Repeated coral bleaching events have eroded reef accretion potential in the Florida Keys

Coral bleaching is now the primary driver of coral-reef degradation globally, threatening the persistence of coral-reef functions and the invaluable ecosystem services they support. The Florida Keys have been impacted by four severe, regional-scale coral-bleaching events, which occurred in 1997–1998, 2005, 2014–2015, and 2023. Our aim in this study was to quantify the effects of those thermal-stress events on the persistence of coral-reef structures. We leveraged long-term coral-reef monitoring data to quantify the impacts of coral bleaching on coral cover, carbonate production, and reef-accretion potential throughout the Florida Keys from 1996–2024. Our study indicates that coral bleaching drove significant declines of reef-building corals, which successively diminished regional-scale reef-accretion potential (from 1.65 ± 0.33 mm y -1 in 1996) and culminated with a regional-scale transition to net reef erosion (-0.12 ± 0.08 mm y -1 ) following the 2023 event. Coral bleaching was the primary driver of declining reef-building capacity, with the cumulative impacts of the four coral-bleaching events accounting for >75% of total declines since 1996. The 1997–1998 event alone drove ~50% of overall declines, largely due to losses of the shallow-water ecosystem engineer, Acropora palmata . More moderate impacts of thermal stress on regional reef-accretion potential in 2005 (6%) and 2014–2015 (9%) could lend some support to the hypothesis that the remaining corals were more thermally tolerant; however, the unprecedented thermal stress event in 2023 caused widespread coral-bleaching and bleaching-related mortality, and caused an additional 12% of the total decline in reef-accretion potential. Stony-coral tissue loss disease accounted for an additional 6% of the decline, with the remaining 19% unexplained. Coral bleaching has already dramatically reshaped the structure and function of coral reefs in the Florida Keys over the last three decades. With the high likelihood that thermal stress will accelerate moving forward, the persistence of Florida’s essential coral-reef functions and services is uncertain.

Florida↗

USGS2025

Igneous differentiation in arcs

This chapter assesses crystallization-differentiation of convergent-margin magmas both from the standpoints of natural rock suites and from laboratory phase-equilibrium experiments. Crystallization-differentiation of dry subalkaline and alkaline basalts are summarized at the beginning for reference and context. Experiments crystallizing arc basalts through andesites at middle to upper crustal pressures (100–400 MPa), saturated either with H 2 O or mixed H 2 O-CO 2 , produce liquid lines of descent similar to arc igneous suites, with higher f O 2 s giving results closer to calc-alkaline magma series. Early enrichment of melts in SiO 2 with little or no enrichment in FeO* (total Fe as FeO) comes about by the combined effects of H 2 O suppressing crystallization of plagioclase relative to olivine, shifting plagioclase to more calcic – therefore lower SiO 2 – compositions, and suppressing silicate minerals to greater extents than Fe-rich spinel. Also contributing to SiO 2 enrichment is crystallization near the liquidus of higher proportions of olivine relative to calcic clinopyroxene, compared with MORB, but it is unclear if this is brought about by H 2 O or elevated pressure. Exhumed arc terranes, however, generally fail to expose evidence for primitive liquids having reached and differentiated at middle to upper crustal depths, at least in any abundance. Instead, most differentiation in arcs appears to take place in the lower crust or near the crust-mantle transition, and although hydrous experiments at such pressures produce liquids broadly similar to arc suites in their SiO 2 enrichment, the experimental liquids differ in becoming peraluminous at early stages of differentiation while melts have basaltic, basaltic andesitic, or andesitic SiO 2 concentrations. Crystallization of abundant calcic clinopyroxene close to the liquidus causes the high-pressure experimental liquids to become peraluminous with only small increases in SiO 2 concentration. In contrast, natural arc volcanic and plutonic suites become gradually less metaluminous with increasing SiO 2 , crossing to peraluminous at dacitic to rhyolitic SiO 2 concentrations. Hybridism, encompassing magma-mixing, assimilation of old crust, and incorporation of near-solidus melts from only shortly older deep, antecedent arc intrusions may accompany natural crystallization-differentiation, preventing melts from attaining peraluminous compositions at low-SiO 2 concentrations. If so, then nearly all arc suites are to substantial extents composites, rather than preserving simple, progressive liquid lines of descent. Some comments and assessments are also made on the participation of garnet in arc differentiation and on the production of arc batholiths.

Book chapter↗