SEARCH · USGS Science

Results for “Complexity”

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

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18Linked to original sources

USGS2011

Evidence from lava flows for complex polarity transitions: The new composite Steens Mountain reversal record

Geomagnetic polarity transitions may be significantly more complex than are currently depicted in many sedimentary and lava-flow records. By splicing together paleomagnetic results from earlier studies at Steens Mountain with those from three newly studied sections of Oregon Plateau flood basalts at Catlow Peak and Poker Jim Ridge 70–90 km to the southeast and west, respectively, we provide support for this interpretation with the most detailed account of a magnetic field reversal yet observed in volcanic rocks. Forty-five new distinguishable transitional (T) directions together with 30 earlier ones reveal a much more complex and detailed record of the 16.7 Ma reversed (R)-to-normal (N) polarity transition that marks the end of Chron C5Cr. Compared to the earlier R-T-N-T-N reversal record, the new record can be described as R-T-N-T-N-T-R-T-N. The composite record confirms earlier features, adds new west and up directions and an entire large N-T-R-T segment to the path, and fills in directions on the path between earlier directional jumps. Persistent virtual geomagnetic pole (VGP) clusters and separate VGPs have a preference for previously described longitudinal bands from transition study compilations, which suggests the presence of features at the core–mantle boundary that influence the flow of core fluid and distribution of magnetic flux. Overall the record is consistent with the generalization that VGP paths vary greatly from reversal to reversal and depend on the location of the observer. Rates of secular variation confirm that the flows comprising these sections were erupted rapidly, with maximum rates estimated to be 85–120 m ka −1 at Catlow and 130–195 m ka −1 at Poker Jim South. Paleomagnetic poles from other studies are combined with 32 non-transitional poles found here to give a clockwise rotation of the Oregon Plateau of 11.4°± 5.6° with respect to the younger Columbia River Basalt Group flows to the north and 14.5°± 4.6° with respect to cratonic North America (95 per cent confidence interval).

Geophysical Journal International↗

USGS2011

Refugial isolation and divergence in the Narrowheaded Gartersnake species complex (Thamnophis rufipunctatus) as revealed by multilocus DNA sequence data

Glacial–interglacial cycles of the Pleistocene are hypothesized as one of the foremost contributors to biological diversification. This is especially true for cold‐adapted montane species, where range shifts have had a pronounced effect on population‐level divergence. Gartersnakes of the Thamnophis rufipunctatus species complex are restricted to cold headwater streams in the highlands of the Sierra Madre Occidental and southwestern USA. We used coalescent and multilocus phylogenetic approaches to test whether genetic diversification of this montane‐restricted species complex is consistent with two prevailing models of range fluctuation for species affected by Pleistocene climate changes. Our concatenated nuDNA and multilocus species analyses recovered evidence for the persistence of multiple lineages that are restricted geographically, despite a mtDNA signature consistent with either more recent connectivity (and introgression) or recent expansion (and incomplete lineage sorting). Divergence times estimated using a relaxed molecular clock and fossil calibrations fall within the Late Pleistocene, and zero gene flow scenarios among current geographically isolated lineages could not be rejected. These results suggest that increased climate shifts in the Late Pleistocene have driven diversification and current range retraction patterns and that the differences between markers reflect the stochasticity of gene lineages (i.e. ancestral polymorphism) rather than gene flow and introgression. These results have important implications for the conservation of T. rufipunctatus (sensu novo), which is restricted to two drainage systems in the southwestern US and has undergone a recent and dramatic decline.

Molecular Ecology↗

USGS2009

Constraints on the depth of generation and emplacement of a magmatic epidote-bearing quartz diorite pluton in the Coast Plutonic Complex, British Columbia

Petrology and P–T estimates indicate that a magmatic epidote-bearing quartz diorite pluton from Mt. Gamsby, Coast Plutonic Complex, British Columbia, was sourced at pressures below ∼1.4 GPa and cooled nearly isobarically at ∼0.9 GPa. The P–T path indicates that the magma was within the stability field of magmatic epidote early and remained there upon final crystallization. The pluton formed and crystallized at depths greater than ∼30 km. REE data indicate that garnet was absent in the melting region and did not fractionate during crystallization. This suggests that the crust was less than or equal to ∼55 km thick at 188 Ma during the early phases of magmatism in the Coast Plutonic Complex. Late Cretaceous contractional deformation and early Tertiary extension exhumed the rocks to upper crustal levels. Textures of magmatic epidote and other magmatic phases, combined with REE data, can be important for constraining the P–T path followed by magmas.

British Columbia↗

USGS2010

Stage structure alters how complexity affects stability of ecological networks

Resolving how complexity affects stability of natural communities is of key importance for predicting the consequences of biodiversity loss. Central to previous stability analysis has been the assumption that the resources of a consumer are substitutable. However, during their development, most species change diets; for instance, adults often use different resources than larvae or juveniles. Here, we show that such ontogenetic niche shifts are common in real ecological networks and that consideration of these shifts can alter which species are predicted to be at risk of extinction. Furthermore, niche shifts reduce and can even reverse the otherwise stabilizing effect of complexity. This pattern arises because species with several specialized life stages appear to be generalists at the species level but act as sequential specialists that are hypersensitive to resource loss. These results suggest that natural communities are more vulnerable to biodiversity loss than indicated by previous analyses.

Ecology Letters↗

USGS2005

Hydrostratigraphic modeling of a complex, glacial-drift aquifer system for importation into MODFLOW

Deposition from at least three episodes of glaciation left a complex glacial-drift aquifer system in central Illinois. The deepest and largest of these aquifers, the Sankoty-Mahomet Aquifer, occupies the lower part of a buried bedrock valley and supplies water to communities throughout central Illinois. Thin, discontinuous aquifers are present within glacial drift overlying the Sankoty-Mahomet Aquifer. This study was commissioned by local governments to identify possible areas where a regional water supply could be obtained from the aquifer with minimal adverse impacts on existing users. Geologic information from more than 2200 existing water well logs was supplemented with new data from 28 test borings, water level measurements in 430 wells, and 35 km of surface geophysical profiles. A three-dimensional (3-D) hydrostratigraphic model was developed using a contouring software package, a geographic information system (GIS), and the 3-D geologic modeling package, EarthVision??. The hydrostratigraphy of the glacial-drift sequence was depicted as seven uneven and discontinuous layers, which could be viewed from an infinite number of horizontal and vertical slices and as solid models of any layer. Several iterations were required before the 3-D model presented a reasonable depiction of the aquifer system. Layers from the resultant hydrostratigraphic model were imported into MODFLOW, where they were modified into continuous layers. This approach of developing a 3-D hydrostratigraphic model can be applied to other areas where complex aquifer systems are to be modeled and is also useful in helping lay audiences visualize aquifer systems.

Ground Water↗

USGS2021

Gulf Coast vicariance shapes phylogeographic history of a North American freshwater mussel species complex

Aim Freshwater mussels share habitat and are parasites of freshwater fishes during the larval life stage. Therefore, models of fish biogeography may also explain the historical biogeography of freshwater mussels. We tested this assumption using predictions of three biogeographic models constructed for northern Gulf of Mexico drainages on a freshwater mussel species complex. Specifically, we tested (1) if speciation was due to vicariant events of fluctuating sea levels that separated lineages east‐west of the Mobile Basin (Central Gulf Coast speciation hypothesis), (2) if the timing of divergences occurred 8.5–3.5 MYA (Gulf Coast allopatric speciation model) and (3) if diversification in Mississippi River populations was recent and for evidence of population increase consistent with range expansion into northern deglaciated regions (Pleistocene glaciation model). Location Eastern North America. Taxon Freshwater mussels (Bivalvia: Unionidae), Lampsilis teres and L. floridensis . Methods We collected 249 specimens from 73 localities across the group's distribution. We used three molecular markers (COI, NDI & ITSI) to conduct time calibrated Bayesian phylogenetic analyses, phylogeographic analyses (AMOVA & SAMOVA) and demographic analyses including Bayesian skyline plots. Results Lampsilis teres and L. floridensis are allopatric species whose distributions meet at the eastern edge of the Mobile Basin. Speciation was estimated to occur in the late Miocene. Populations from isolated river systems surrounding the Gulf of Mexico are almost all monophyletic. Mississippi drainage samples formed a shallow clade with recent diversification and showed evidence of recent population expansion. Main conclusions The historical biogeography of the L. teres species complex is broadly consistent with tested ichthyofaunal models. The timing of speciation and intraspecific divergences correspond to low sea‐level events suggesting that Gulf Coast sea‐level fluctuations are responsible for dispersal (sea‐level recession) and subsequent cladogenesis (sea‐level inundation). Consistent with numerous other freshwater studies, we found the Mobile Basin to be a suture zone, which may be due to the narrow, offshore continental shelf.

Journal of Biogeography↗

USGS1985

Origin of the Sudbury Complex by meteoritic impact: Neodymium isotopic evidence

Samarium-neodymium isotopic data on whole rocks and minerals of the Sudbury Complex in Canada gave an igneous crystallization age of 1840 ± 21 × 10 6 years. The initial epsilon neodymium values for 15 whole rocks are similar to those for average upper continental crust, falling on the crustal trend of neodymium isotopic evolution as defined by shales. The rare earth element concentration patterns of Sudbury rocks are also similar to upper crustal averages. These data suggest that the Sudbury Complex formed from melts generated in the upper crust and are consistent with a meteoritic impact.

Ontario↗

USGS2000

Detrital zircon geochronology of the Shoo Fly Complex, northern Sierra terrane, northeastern California

U-Pb analyses have been conducted on 92 individual detrital zircon grains from 4 of the main thrust sheets of the Shoo Fly Complex . Samples from the Culbertson Lake allochthon, Duncan Peak allochthon, and Lang sequence yield mainly 1.80-2.10, 2.20-2.45, and 2.55-2.70 Ga ages, which suggests that sediments in these units originated in a cratonal region containing Paleoproterozoic and Archean igneous rocks. These ages match those of basement provinces from the northwestern Canadian shield, suggesting a provenance link with northwestern North America during early Paleozoic time. The Sierra City melange, however, has significantly different zircon ages of 551-635 and 1170-1319 Ma, with only a subordinate population of >1.8 Ga grains. These grains apparently were derived originally from an outboard Neoproterozoic-Cambrian(?) volcanic arc and from 1.0-1.7 Ga basement rocks of south- western North America . The occurrence of all three sets of ages in a sandstone that accumulated outboard of the Lang, Culbertson Lake, and Duncan Peak thrust sheets indicates that most rocks of the Shoo Fly Complex formed inboard of a volcanic arc located in proximity to the southern portion of the Cordilleran margin.

California↗

USGS1997

Mid-crustal flow during Tertiary extension in the Ruby Mountains core complex, Nevada

Structural analysis and geochronologic data indicate a nearly orthogonal, late Eocene–Oligocene flow pattern in migmatitic infrastructure immediately beneath the kilometer-thick, extensional, mylonitic shear zone of the Ruby Mountains metamorphic core complex, Nevada. New U-Pb radiometric dating indicates that the development of a northward-trending lineation in the infrastructure is partly coeval with the development of a pervasive, west-northwest–trending lineation in the mylonitic shear zone. U-Pb monazite data from the leucogranite orthogneiss of Thorpe Creek indicate a crystallization age of ca. 36–39 Ma. Zircon fractions from a biotite monzogranite dike yield an age of ca. 29 Ma. The three dated samples from these units exhibit a penetrative, approximately north-south–trending elongation lineation. This lineation is commonly defined by oriented bundles of sillimanite and/or elongated aggregates of quartz and feldspar, indicating a synmetamorphic and syndeformational origin. The elongation lineation can be interpreted as a slip line in the flow plane of the migmatitic, nonmylonitic infrastructural core of the northern Ruby Mountains. A portion of this midcrustal flow is coeval with the well-documented, west-northwest sense of slip in the structurally overlying kilometer-thick, mid-Tertiary mylonitic shear zone. Lineations in the mylonitic zone are orthogonal to those in the deeper infrastructure, suggesting fundamental plastic decoupling between structural levels in this core complex. Furthermore, the infrastructure is characterized by overlapping, oppositely verging fold nappes, which are rooted to the east and west. One of the nappes may be synkinematic with the intrusion of the late Eocene orthogneiss of Thorpe Creek. In addition, the penetrative, elongation lineation in the infrastructure is subparallel to hinge lines of parasitic folds developed synchronous with the fold nappes, suggesting a kinematically related evolution.

Nevada, Utah↗

USGS1998

Nature and regional significance of unconformities associated with the Middle Ordovician Hagan K-bentonite complex in the North American midcontinent

Stratal patterns of the Middle Ordovician Hagan K-bentonite complex and associated rocks show that the Black River–Trenton unconformity in the North American midcontinent formed through the complex interplay of eustasy, sediment accumulation rates, siliciclastic influx, bathymetry, seawater chemistry, and perhaps local tectonic uplift. The unconformity is diachronous and is an amalgamated surface that resulted from local late Turinian lowstand exposure followed by regional early Chatfieldian transgressive drowning and sediment starvation. The duration of the unconformity is greatest in southern Wisconsin, northern Illinois, and northern Indiana, where the Deicke and Millbrig K-bentonite Beds converge at the unconformity. On the basis of published isotopic ages for the Deicke and Millbrig beds, it is possible that in these regions erosion and nondeposition spanned a period of as much as 3.2 m.y.

Geological Society of America Bulletin↗

USGS1992

Synmagmatic deformation in the underplated igneous complex of the Ivrea-Verbano zone

The Ivrea-Verbano zone, northern Italy, contains an igneous complex up to 10 km thick that is thought to have been intruded near the interface between the continental crust and mantle during the late Paleozoic. New data indicate that this complex is pervasively deformed and concentrically foliated. Widespread deformation under hypersolidus conditions is indicated by growth of undeformed poikilitic phases across the foliation and segregation of late-stage melts into high-temperature faults and pressure shadows of boudins. The presence of analogous features in ophiolitic gabbros suggests that emplacement of the Ivrea- Verbano zone plutonic rocks involved large-scale flow of crystal mush in a dynamic, and possibly extensional, tectonic environment.

Geology↗

USGS1993

Massive sulfide metallogenesis at a late Mesozoic sediment-covered spreading axis: Evidence from the Franciscan complex and contemporary analogues

The Island Mountain deposit, an anomalous massive sulfide in the Central belt of the Franciscan subduction complex, northern California Coast Ranges, formed during hydrothermal activity in a sediment-dominated paleo-sea-floor environment. Although the base of the massive sulfide is juxtaposed against a 500-m-wide melange band, its gradational upper contact within a coherent sequence of sandstone, siltstone, and mudstone indicates that hydrothermal activity was concurrent with turbidite deposition. Accumulations of sulfide breccia and clastic sulfide were produced by mass wasting of the sulfide mound prior to burial by turbidites. The bulk composition of sulfide samples (pyrrhotite rich; high Cu, As, and Au contents; radiogenic Pb isotope ratios) is consistent with a hydrothermal system dominated by fluid-sediment interaction. On the basis of a comparison with possible contemporary tectonic analogues at the southern Gorda Ridge and the Chile margin triple junction, we propose that massive sulfide mineralization in the Central belt of the Franciscan complex resulted from hydrothermal activity at a late Mesozoic sediment-covered ridge axis prior to collision with the North American plate.

California↗

USGS2019

Eruption age and duration of the ~9 km3 Burney Mountain dacite dome complex, northern California

At ~9 km3, the six dacite domes of Burney Mountain (db1–db6) constitute the most voluminous Quaternary dome complex in the Cascades volcanic arc. Whole-rock geochemistry, electron microprobe, and petrographic data indicate that the domes are magmatically related, which, when integrated with geomorphology and stratigraphy, indicate early (db1, db2, db3) and late (db4, db5, db6) erupted groups. We present 40Ar/39Ar ages of 271.9±4.6 ka (db1), 280.8±8.2 and 281.7±6.8 ka (db2), and 290.2±6.0 ka (db3) along with a previous age of 280±12 ka (db1). These ages scatter over 20 kyr, whereas remanent magnetic directions are similar between 53.3–59.0° inclination and 352.7–355.9° declination. The latter dataset indicates that the dacite domes were emplaced over a geologically brief time interval; not thousands of years. Crystal-size distribution patterns of plagioclase were used to calculate residence times, which we use to infer the duration over which the eruptions likely occurred. Three slopes represent three populations of plagioclase crystals (fine-grained groundmass, coarse-grained groundmass, phenocrysts). A commonly used growth rate for plagioclase in dacitic magmas (10-10 mm/s) yields 9–10 years of growth for the coarse-grained groundmass (early erupted domes of db1, db2, db3), whereas plagioclase in the fine-grained groundmass (late erupted domes of db4, db5, db6) grew over 4–5 years. All plagioclase phenocrysts have apparent residence times of 26–36 years; however, they contain high An>70 resorbed cores with sieve textures, which have euhedral, lower An<65 overgrowth rims. Similarities in chemistry between groundmass plagioclase and phenocryst overgrowth rims indicate that they grew concurrently, and we therefore propose that both have similar residence times. Thus, the Burney Mountain dacite dome complex was emplaced during a single eruptive episode over the course of years to decades at 281.1±4.8 ka (weighted mean age).

California↗

USGS2006

Detrital zircon provenance of the Late Triassic Songpan-Ganzi complex: Sedimentary record of collision of the North and South China blocks

Using detrital zircon geochronology, turbidite deposystems fed from distinct sediment sources can be distinguished within the Songpan-Ganzi complex, a collapsed Middle to Late Triassic turbidite basin of central China. A southern Songpan-Ganzi deposystem initially was sourced solely by erosion of the Qinling-Dabie orogen during early Late Triassic time, then by Qinling-Dabie orogen, North China block, and South China block sources during middle to late Late Triassic time. A northern Songpan-Ganzi system was sourced by erosion of the Qinling-Dabie orogen and the North China block throughout its deposition. These separate deposystems were later tectonically amalgamated to form one complex and then uplifted as the eastern Tibet Plateau. ?? 2006 Geological Society of America.

Geology↗

USGS2013

Previously unrecognized regional structure of the Coastal Belt of the Franciscan Complex, northern California, revealed by magnetic data

Magnetic anomalies provide surprising structural detail within the previously undivided Coastal Belt, the westernmost, youngest, and least-metamorphosed part of the Franciscan Complex of northern California. Although the Coastal Belt consists almost entirely of arkosic graywacke and shale of mainly Eocene age, new detailed aeromagnetic data show that it is pervasively marked by long, narrow, and regularly spaced anomalies. These anomalies arise from relatively simple tabular bodies composed principally of magnetic basalt or graywacke confi ned mainly to the top couple of kilometers, even though metamorphic grade indicates that these rocks have been more deeply buried, at depths of 5–8 km. If true, this implies surprisingly uniform uplift of these rocks. The basalt (and associated Cretaceous limestone) occurs largely in the northern part of the Coastal Belt; the graywacke is recognized only in the southern Coastal Belt and is magnetic because it contains andesitic grains. The magnetic grains were not derived from the basalt, and thus require a separate source. The anomalies defi ne simple patterns that can be related to folding and faulting within the Coastal Belt. This apparent simplicity belies complex structure mapped at outcrop scale, which can be explained if the relatively simple tabular bodies are internally deformed, fault-bounded slabs. One mechanism that can explain the widespread lateral extent of the thin layers of basalt is peeling up of the uppermost part of the oceanic crust into the accretionary prism, controlled by porosity and permeability contrasts caused by alteration in the upper part of the subducting slab. It is not clear, however, how this mechanism might generate fault-bounded layers containing magnetic graywacke. We propose that structural domains defined by anomaly trend, wavelength, and source reflect imbrication and folding during the accretion process and local plate interactions as the Mendocino triple junction migrated north, a hypothesis that should be tested by more detailed structural studies.

California↗

USGS1993

Extension and contraction within an evolving divergent strike-slip fault complex: The San Andreas and San Jacinto fault zones at their convergence in southern California

A variety of extensional and contractional structures is produced by strike slip faulting. The variety and extent of the structures are directly related to the kind and extent of geometric complexities of the fault zone or system. The area of convergence of the San Andreas fault zone and the much younger San Jacinto fault zone in the eastern Transverse Ranges is exquisitely complex. We propose that the San Jacinto fault zone formed in response to a structural knot in San Gorgonio Pass probably within the past 1.5 Ma. In the area of their convergence we propose that slip is transferred both east and west from the San Jacinto fault zone northward to the San Andreas fault zone over a 60-to 70-km band that extends northwestward from the south end of the San Bernardino basin to the east end of the San Gabriel Mountains. We further propose several structural adjustments as a consequence of onset or acceleration of lateral movement on the San Jacinto fault zone: accelerated uplift of the eastern San Gabriel Mountains, development or accentuation of an arcuate schuppen-like structure in the eastern San Gabriel Mountains, inception of the San Bernardino basin, cessation of deposition in the present-day San Timoteo badlands area, inception of the San Jacinto basin, and an increase in compression and uplift in the San Gorgonio Pass area. We interpret the uplift and compression in San Gorgonio Pass to result from two formerly disparate structural blocks—the eastern San Bernardino and San Jacinto blocks—becoming a relatively coherent block, and the San Gorgonio Pass area constituting a left step between the San Andreas fault zone in the Coachella Valley area and the San Jacinto fault zone in the San Jacinto Valley area. The compression and uplift led to the formation of the San Gorgonio Pass thrust faults and disruption of any through-going San Andreas strands, at least at the surface. In partitioning slip between the San Andreas and San Jacinto fault zones, consideration should be given to the bandwidth over which horizontal strain has accumulated. The average slip rate of the northern part of the San Jacinto fault zone during the past 1.5 m.y. may have been about 20 mm/yr and about 15 mm/yr on the San Andreas. South of the San Bernardino basin, current strain accumulation based on repeated geodetic surveys is nearly equally divided between the San Jacinto and San Andreas fault zones.

California↗

USGS2023

Incorporating wave climate complexity into modeling lower shoreface morphology and transport

The lower shoreface, a transitional subaqueous region extending from the seaward limit of the surf zone to beyond the closure depth, serves as a sediment reservoir and pathway in sandy beach environments over annual to millennial time scales. Despite the important role this region plays in shoreline dynamics, the morphodynamics of the lower shoreface remain poorly quantified and understood. To better understand controls on shoreface morphology, here we combine energetics-based suspended sediment transport formulae (Ortiz & Aston 2016) with empirical wave climate data to incorporate temporal complexity in modeled equilibrium profiles and sediment flux rates. The equilibrium shoreface shape computed using a full wave climate is steeper in shallower water and less steep in the deeper reaches compared to profiles computed using single wave characteristics. Using a full wave climate to simulate steady-state morphology will yield steeper profiles in shallow water. Suspended sediment transport rates also vary in direction and magnitude at different equilibrium profile depths and can potentially inform the location of morphodynamic boundaries in the shoreface. Our results reveal how infrequent storm waves affect shoreface slopes, with large events tending to drive sediment onshore in the deeper portions of the profile. This work explores a few ways to add complexity to simple energetics-based frameworks to reproduce empirical bathymetric data more accurately and provides insight toward refining coastal source-to-sink models.

Conference Paper↗

USGS2009

Arc-continent collision and the formation of continental crust: A new geochemical and isotopic record from the Ordovician Tyrone Igneous Complex, Ireland

Collisions between oceanic island-arc terranes and passive continental margins are thought to have been important in the formation of continental crust throughout much of Earth's history. Magmatic evolution during this stage of the plate-tectonic cycle is evident in several areas of the Ordovician Grampian-Taconic orogen, as we demonstrate in the first detailed geochemical study of the Tyrone Igneous Complex, Ireland. New U-Pb zircon dating yields ages of 493 2 Ma from a primitive mafic intrusion, indicating intra-oceanic subduction in Tremadoc time, and 475 10 Ma from a light rare earth element (LREE)-enriched tonalite intrusion that incorporated Laurentian continental material by early Arenig time (Early Ordovician, Stage 2) during arc-continent collision. Notably, LREE enrichment in volcanism and silicic intrusions of the Tyrone Igneous Complex exceeds that of average Dalradian (Laurentian) continental material that would have been thrust under the colliding forearc and potentially recycled into arc magmatism. This implies that crystal fractionation, in addition to magmatic mixing and assimilation, was important to the formation of new crust in the Grampian-Taconic orogeny. Because similar super-enrichment of orogenic melts occurred elsewhere in the Caledonides in the British Isles and Newfoundland, the addition of new, highly enriched melt to this accreted arc terrane was apparently widespread spatially and temporally. Such super-enrichment of magmatism, especially if accompanied by loss of corresponding lower crustal residues, supports the theory that arc-continent collision plays an important role in altering bulk crustal composition toward typical values for ancient continental crust. ?? 2009 Geological Society of London.

Journal of the Geological Society↗