SEARCH · Geology Reports

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 37 records · Page 2Linked to original sources

USGS1983

Late Proterozoic island-arc complexes and tectonic belts in the southern part of the Arabian Shield, Kingdom of Saudi Arabia

Two main subdivisions of layered rocks are recognized in the southern Arabian Shield south of lat 22? N. These are an older ensimatic-arc complex, which formed 1100-800 m.y. ago, and a younger marginal-arc complex, which formed 800-690 m.y. ago. The older ensimatic-arc complex, located in the southwestern part of the Shield, includes graywacke and mafic to intermediate volcanic rocks of the essentially contemporaneous Baish, Bahah, and Jiddah groups. Although the younger arc complex is also dominantly ensimatic in character, it is also partly superimposed over the older ensimaticarc complex. The superimposed portions of the younger arc complex are represented by the Ablah, Samran, and possibly the Ararat groups. The ensimatic portion of the younger arc group is represented by the Halaban group, which was deposited to the east and northeast of the older ensimatic-arc complex. The Halaban group includes andesitic and dacitic volcanic rocks and associated clastic sedimentary rocks. The layered rocks of both arc complexes are intruded by dioritic (quartz diorite, tonalite, trondhjemite) plutonic rocks. The southern Shield is also subdivided into a number of structurally bounded, north-trending tectonic belts. Within the older ensimatic complex, three belts are recognized. From west to east, these are the Lith, Bidah, and Tayyah belts. Within these three belts, progressive facies changes indicate a gradation from deep-water facies in the south to shallow-water or-terrestrial facies in the north. The distribution of dioritic batholiths, as well as the distribution of layered-rock facies, suggests a northwest-trending axis for the older ensimatic-arc complex. The younger arc complex is present within six belts, the Makkah source papers. In Fleck and others (1980), the term 'quartz diorite' includes both tonalite and quartz diorite as defined in the International Union of Geological Sciences (IUGS) system of plutonic rock classification (Streckeisen, 1973). Initial 87Sr/86Sr ratios are not included in the appendix, but all rocks more than 660 m.y. old have initial ratios in the range 0.7021-0.7035, with only two greater than 0.7030. Thus, nothing in the Rb-Sr data suggests involvement of an older continental crust during the evolution of the southern Shield. A lead isotope study of ore minerals and potassium feldspars of the Arabian Shield by Stacey and others (1980) also suggests that no older (Archean to early Proterozoic) evolved continental-type crust underlies the southern Shield. An early summary of mapping (Schmidt and others, 1973) suggests that older sialic basement underlies the late Proterozoic layered rocks in the southern Shield. However, subsequent-mapping and the isotopic studies cited above have established that all of these rocks are of late Proterozoic age and that all rocks of the southern Shield that are more than 660 m.y. old have ensimatic or mantle isotopic characteristics. Figure 2 shows, with only two exceptions, that rocks more than 800 m.y. old are present west of the boundary separating the Tayyah and Khadra belts. The exceptions are two poorly controlled Rb-Sr ages obtained by Fleck (1980) on two quartz diorite plutons in the Malahah region (appendix 1, localities 26 and 27). Preliminary uranium-thorium zircon data of Stacey now suggest that one of these quartz diorite plutons (locality 26) has an age of approximately 640 m.y. Therefore, we prefer to discount the two dates of Fleck until further information is available. As noted earlier and as described below, most of the rocks of the southern Arabian Shield have characteristics typical of those formed in the island-arc environment by subduction-related processes. We shall refer to the group of rocks in the western part of the southern Shield, which formed from 1100 to 800 m.y. ago, as the 'older ensimatic-arc complex' and those in the eastern and northwestern parts, which formed from 800 to 690 m.y. ago, as the 'younger marginal-arc compl

Open-File Report↗

USGS1991

A lead isotopic study of the Stillwater Complex, Montana: constraints on crustal contamination and source regions

Analyses of the Pb isotopic compositions of plagioclase from 23 samples covering the stratigraphic thickness of the Stillwater Complex indicate a narrow range of apparent initial isotopic compositions (206Pb/ 204Pb=13.95; 207Pb/204Pb=14.95-15.01; 208Pb/204Pb=33.6). The uniformity of our data is in contrast to, but not necessarily contradictory to, other recent investigations which give indications that the complex formed by repeated injection of magmas with at least two distinct compositions that were presumably derived from different source regions. Samples from the Basal series of the complex have consistently higher 207Pb/204Pb ratios, suggesting either minor contamination from adjacent country rocks or a slight distinction between parental magmas. Apparent initial Pb isotopic compositions of the complex are very radiogenic compared to Late Archean model-mantle values, but are nearly identical to initial Pb isotopic compositions found for the the adjacent, slightly older (2.73-2.79 Ga), Late Archean crustal suite in the Beartooth Mountains. Contamination of magmas parental to the Stillwater Complex by the Late Archean crustal suite is rejected for two reasons: (1) Th and U concentrations in Stillwater rocks and plagioclase are very low (about 0.08 and 0.02 ppm respectively), yet Th/U ratios are uniform at about 4, in contrast to the highly variable (2-26) but often high Th/U ratios found for the Late Archean crustal complex; (2) it seems improbable that any contamination process would have adjusted the isotopic compositions of the diverse magmas entering the Stillwater chamber to near-identical values. The preferred hypothesis to explain the Pb isotopic data for the Stillwater Complex and the associated Late Archean crustal suite involves a major Late Archean crust-forming event that resulted in a compositionally complex crust/mantle system with relatively homogeneous and unusual Pb isotopic compositions. The parental magmas of the Stillwater Complex were generated at different levels within this crust/mantle system, before isotopic contrasts could develop by radioactive decay within compositionally discrete reservoirs. This situation limits the utility of all isotopic tracer systems in discriminating among the various mantle and crustal reservoirs that may have affected the final isotopic character of the Stillwater magmas. The late Archean crustal complex and the Stillwater Complex melts were ultimately derived from the same distinct mantle without obvious direct interaction with the Middle to Early Archean crust present in the region. ?? 1991 Springer-Verlag.

Contributions to Mineralogy and Petrology↗

USGS2012

Geophysical study of the San Juan Mountains batholith complex, southwestern Colorado

One of the largest and most pronounced gravity lows over North America is over the rugged San Juan Mountains of southwestern Colorado (USA). The mountain range is coincident with the San Juan volcanic field (SJVF), the largest erosional remnant of a widespread mid-Cenozoic volcanic field that spanned much of the southern Rocky Mountains. A buried, low-density silicic batholith complex related to the volcanic field has been the accepted interpretation of the source of the gravity low since the 1970s. However, this interpretation was based on gravity data processed with standard techniques that are problematic in the SJVF region. The combination of high-relief topography, topography with low densities, and the use of a common reduction density of 2670 kg/m 3 produces spurious large-amplitude gravity lows that may distort the geophysical signature of deeper features such as a batholith complex. We applied an unconventional processing procedure that uses geologically appropriate densities for the uppermost crust and digital topography to mostly remove the effect of the low-density units that underlie the topography associated with the SJVF. This approach resulted in a gravity map that provides an improved representation of deeper sources, including reducing the amplitude of the anomaly attributed to a batholith complex. We also reinterpreted vintage seismic refraction data that indicate the presence of low-velocity zones under the SJVF. Assuming that the source of the gravity low on the improved gravity anomaly map is the same as the source of the low seismic velocities, integrated modeling corroborates the interpretation of a batholith complex and then defines the dimensions and overall density contrast of the complex. Models show that the thickness of the batholith complex varies laterally to a significant degree, with the greatest thickness (∼20 km) under the western SJVF, and lesser thicknesses (<10 km) under the eastern SJVF. The largest group of nested calderas on the surface of the SJVF, the central caldera cluster, is not correlated with the thickest part of the batholith complex. This result is consistent with petrologic interpretations from recent studies that the batholith complex continued to be modified after cessation of volcanism and therefore is not necessarily representative of synvolcanic magma chambers. The total volume of the batholith complex is estimated to be 82,000–130,000 km 3 . The formation of such a large felsic batholith complex would inevitably involve production of a considerably greater volume of residuum, which could be present in the lower crust or uppermost mantle. The interpreted vertically averaged density contrast (–60 to –110 kg/m 3 ), density (2590–2640 kg/m 3 ), and seismic expression of the batholith complex are consistent with results of geophysical studies of other large batholiths in the western United States.

Colorado↗

USGS2005

Contrasting Proterozoic basement complexes near the truncated margin of Laurentia, northwestern Sonora–Arizona international border region

We utilize new geological mapping, conventional isotope dilution–thermal ionization mass spectrometry (ID-TIMS) and sensitive high-resolution ion microprobe (SHRIMP) U-Pb zircon analyses, and whole-rock radiogenic isotope characteristics to distinguish two contrasting Proterozoic basement complexes in the international border region southeast of Yuma, Arizona. Strategically located near the truncated southwest margin of Laurentia, these Proterozoic exposures are separated by a northwest-striking Late Cretaceous batholith. Although both complexes contain strongly deformed Paleoproterozoic granitoids (augen gneisses) intruded into fine-grained host rocks, our work demonstrates marked differences in age, host rock composition, and structure between the two areas. The Western Complex reveals a >5-km-thick tilted section of finely banded felsic, intermediate, and mafic orthogneiss interspersed with tabular intrusive bodies of medium-grained leucocratic biotite granite (1696 ± 11 Ma; deepest level), medium-grained hornblende-biotite granodiorite (1722 ± 12 Ma), and coarse-grained porphyritic biotite granite (1725 ± 19 Ma; shallowest level). Penetrative ductile deformation has converted the granites to augen gneisses and caused isoclinal folding and transposition of primary contacts. Exposed in a belt of northwest-trending folds, these rocks preserve southwest-vergent shear fabric annealed during amphibolite facies metamorphism, when crystalloblastic textures developed. Deformation and regional metamorphism occurred before emplacement of 1.1 Ga(?) mafic dikes. Throughout the Eastern Complex, meta-arkose, quartzite, biotite schist, and possible felsic metavolcanic rocks comprise the country rocks of strongly foliated medium- and coarse-grained biotite granite augen gneisses that yield mean 207 Pb/ 206 Pb ages of 1646 ± 10 Ma, 1642 ± 19 Ma, and 1639 ± 15 Ma. Detrital zircons from four samples of host sandstone are isotopically disturbed; nevertheless, the data indicate a restricted provenance (ca. 1665 Ma to 1650 Ma), with two older grains (1697 and 1681 Ma). The pervasively recrystallized Paleoproterozoic map units strike parallel to foliation and are repeated in south-trending folds that are locally refolded about easterly hinges. Southeasterly lineation developed in augen gneiss and host strata becomes penetrative in local domains of L-tectonite. Regional metamorphism associated with this tectonism persisted until ca. 1590 Ma, as recorded by metamorphic growths within some zircon grains. Mesoproterozoic intrusions that crosscut the Paleoproterozoic metasediments and augen gneisses include coarsely porphyritic biotite granite (1432 ± 6 Ma) and diabase dikes (1.1 Ga?). Emplacement of the granite was accompanied by secondary high-U overgrowths, dated at 1433 ± 8 Ma, on some of the Paleoproterozoic detrital zircons, and apparently was also responsible for resetting the whole-rock Pb isotopic systematics (1441 ± 39 Ma) within these Eastern Complex augen gneisses. Younger plutons emplaced into both Proterozoic basement complexes include medium-grained quartz diorite (73.4 ± 3.3 Ma and 72.8 ± 1.7 Ma), Late Cretaceous hornblende-biotite granodiorite, and Paleogene leucocratic biotite granite. Neogene sedimentary and volcanic strata overlie basement along unconformities that are tilted to the northeast, southeast, or southwest. A brittle normal fault, dipping gently northeast, juxtaposes Tertiary andesite with Paleoproterozoic metasandstone. These relationships suggest that the area shares a common history of mid-Tertiary extension with southwestern Arizona. Later influence of the southern San Andreas fault system is implied by multiple dextral offsets of pre-Tertiary units across northwest-trending valleys. Our structural, geochronologic, and isotopic data provide new information to constrain pre–750 Ma Rodinia reconstructions involving southwestern Laurentia. Whole-rock U-Th-Pb and Rb-Sr isotopic systematics in both Paleoproterozoic gneiss complexes are disturbed, however, well-behaved Sm-Nd analyses preserve depleted initial ε Nd values (+2 to +4) that are distinct from the Mojave crustal province, but overlapping with the Yavapai and Mazatzal Provinces of Arizona. The Eastern Complex has the appropriate age and Nd isotopic signature to be part of the Mazatzal Province, but records major tectonism and metamorphism at ca. 1.6 Ga that postdates the Mazatzal orogeny. Deformed granitoids of the Western Complex have “Yavapai-type” ages and ε Nd but display structures discordant to the southwesterly Yavapai trend in central Arizona. The Western Complex lies along-strike with similar-age rocks (1.77 Ga to 1.69 Ga) of the “Caborca block” that have only been studied in detail near Quitovac and south of Caborca. Collectively, these rocks form a northwest-trending strip of basement situated at the truncated edge of Laurentia. The present-day basement geography may reflect an original oroclinal bend in the Yavapai orogenic belt. Alternatively, the western Proterozoic belt of Sonora may represent displaced fragments of basement juxtaposed against the Yavapai-Mazatzal Provinces along a younger sinistral transform fault (e.g., the Late Jurassic Mojave-Sonora megashear or the Permian Coahuila transform). Crustal blocks with these specific petrologic, geochronologic, and isotopic characteristics can be found in south-central and northeastern portions of the Australian Proterozoic basement, further supporting a connection between the two continents prior to breakup of the Rodinian supercontinent.

Arizona, California, Nevada, Sonora↗

USGS2015

What lies beneath: geophysical mapping of a concealed Precambrian intrusive complex along the Iowa–Minnesota border

Large-amplitude gravity and magnetic highs over northeast Iowa are interpreted to reflect a buried intrusive complex composed of mafic–ultramafic rocks, the northeast Iowa intrusive complex (NEIIC), intruding Yavapai province (1.8–1.72 Ga) rocks. The age of the complex is unproven, although it has been considered to be Keweenawan (∼1.1 Ga). Because only four boreholes reach the complex, which is covered by 200–700 m of Paleozoic sedimentary rocks, geophysical methods are critical to developing a better understanding of the nature and mineral resource potential of the NEIIC. Lithologic and cross-cutting relations interpreted from high-resolution aeromagnetic and airborne gravity gradient data are presented in the form of a preliminary geologic map of the basement Precambrian rocks. Numerous magnetic anomalies are coincident with airborne gravity gradient (AGG) highs, indicating widespread strongly magnetized and dense rocks of likely mafic–ultramafic composition. A Yavapai-age metagabbro unit is interpreted to be part of a layered intrusion with subvertical dip. Another presumed Yavapai unit has low density and weak magnetization, observations consistent with felsic plutons. Northeast-trending, linear magnetic lows are interpreted to reflect reversely magnetized diabase dikes and have properties consistent with Keweenawan rocks. The interpreted dikes are cut in places by normally magnetized mafic–ultramafic rocks, suggesting that the latter represent younger Keweenawan rocks. Distinctive horseshoe-shaped magnetic and AGG highs correspond with a known gabbro, and surround rocks with weaker magnetization and lower density. Here, informally called the Decorah complex, the source body has notable geophysical similarities to Keweenawan alkaline ring complexes, such as the Coldwell and Killala Lake complexes, and Mesoproterozoic anorogenic complexes, such as the Kiglapait, Hettasch, and Voisey’s Bay intrusions in Labrador. Results presented here suggest that much of the NEIIC is composed of such complexes, and broadly speaking, may be a discontinuous group of several intrusive bodies. Most units are cut by suspected northwest-trending faults imaged as magnetic lineaments, and one produces apparent sinistral fault separation of a dike in the eastern part of the survey area. The location, trend, and apparent sinistral sense of motion are consistent with the suspected faults being part of the Belle Plaine fault zone, a complex transform fault zone within the Midcontinent rift system that is here proposed to correspond with a major structural discontinuity.

Iowa, Minnesota↗

USGS2011

Relation of hydrologic processes to groundwater and surface-water levels and flow directions in a dune-beach complex at Indiana Dunes National Lakeshore and Beverly Shores, Indiana

The potential for high groundwater levels to cause wet basements (groundwater flooding) is of concern to residents of communities in northwestern Indiana. Changes in recharge from precipitation increases during 2006-9, water-level changes from restoration of nearby wetlands in the Great Marsh in 1998-2002, and changes in recharge due to the end of groundwater withdrawals for water supply since 2005 in a community at Beverly Shores, Ind., were suspected as factors in increased groundwater levels in an unconfined surficial aquifer beneath nearby parts of a dune-beach complex. Results of this study indicate that increased recharge from precipitation and snowmelt was the principal cause of raised water levels in the dune-beach complex from 2006 to 2009. Annual precipitation totals in 2006-9 ranged from 43.88 to 55.75 inches per year (in/yr) and were substantially greater than the median 1952-2009 precipitation of 36.35 in/yr. Recharge to groundwater from precipitation in 2006-9 ranged from 13.5 to 22 in/yr; it was higher than the typical 11 in/yr because of large precipitation events and precipitation amounts received during non-growing-season months. An estimated increase in net recharge from reduced groundwater use in Beverly Shores since 2005 ranged from 1.6 in/yr in 2006 to 1.9 in/yr in 2009. Surface-water levels in the wetland were as much as about 1.1 feet higher in 2007-9 (after the 1998-2002 wetland restoration) than during seasonally wet periods in 1979-89. Similar surface-water levels and ponded water were likely during winter and spring wet periods before and after wetland restoration. High water levels similar to those in 2009 were measured elsewhere in the dune-beach complex near a natural wetland during the spring months in 1991 and 1993 after receipt of near record precipitation. Recharge from similarly high precipitation amounts in 2008-9 was also a likely cause of high groundwater levels in other parts of the dune-beach complex, such as at Beverly Shores. Perennial mounding of the water table in the surficial aquifer indicates that the recharge that created the water-table mound originates within the dune-beach complex and not through flow from the adjacent hydrologic boundaries: the restored wetland, Lake Michigan, and Derby Ditch. Infiltrating precipitation causes most seasonal and episodic rises in groundwater levels beneath the dune-beach complex. Groundwater-level fluctuations lasting days to weeks in the dune-beach complex in 2008-9 were superimposed on a seasonal high water-table altitude that began with the recharge from snowmelt and rain in February 2009 and maintained through July 2009. Increases in water-table-mound altitude under the dune-beach complex recurred in 2008-9 in response to the largest rain events of 1 inch or more and to snowmelt. Smaller, shorter-term rises in water level after individual rain events persisted over hours to less than 1 week. Groundwater-level fluctuations varied over a relatively narrow range of about 2 to 3 feet, with no net fluctuations greater than 4 feet. Groundwater levels in or near low parts of the dune-beach complex were frequently within 0 to 6 feet of the land surface and indicate the potential for groundwater flooding. Groundwater-level gradients from the water-table mound to wells next to surface-water discharges increase after rainfall and snowmelt events and recede slowly as groundwater discharges from the aquifer. Evapotranspiration is responsible for part of the general pattern of decreasing water-table altitudes observed from May to August 2009. Rapid water-level rises in the restored wetland after precipitation do not likely have an effect on groundwater flooding elsewhere in the dune-beach complex. Surface-water-level fluctuations during this study generally varied over a narrower range, approximately from 1 to 1.5 feet, as compared with groundwater fluctuations, except after a very large, 10.77-inch rainfall. Time-delayed and smaller groundwater-level

Indiana↗

USGS2024

Investigation of land cover within wetland complexes at Dixie Meadows, Churchill County, Nevada, from October 2015 to January 2022

The U.S. Geological Survey investigated land cover at subannual time steps within six wetland areas in Dixie Valley, Churchill County, Nevada, from October 2015 to January 2022. As requested by the U.S. Fish and Wildlife Service, we used aerial photography and satellite remote sensing data to map surface water and other land cover types within the wetland complexes. We identified five land cover classes using the green normalized difference vegetation index (gNDVI) and its inverse relationship to the normalized difference water index (NDWI) within three U.S. Department of Agriculture National Agriculture Imagery Program aerial images (acquired in 2015, 2017, and 2019) and 110 European Space Agency Sentinel-2 satellite images (acquired 2015–2022). The relative wetness of soil conditions within each land cover class is estimated by comparison to previously published observations of relative conductivity measured by 79 field-based sensors within the wetlands from 2019 to 2021. We mapped the areal coverage of the five land cover classes for approximately 385 acres (1,559,000 square meters [m²]) comprising six individual wetland complexes as well as a larger 1,298- acre (5,254,000-m 2 ) area of interest inclusive of the wetland complexes and adjacent landscape. Land cover of open water (Class 5) primarily within ponds at one of the wetland complexes comprised 8,333 m 2 , on average, of the wetland complexes. Land cover of mixed shallow surface water, saturated soil, and vegetation (Class 4) comprised 111,723 m 2 on average of the wetland complexes. Land cover of dense green vegetation canopy cover (Class 3) that often (46 percent of observations) had underlying surface water or saturated soil conditions comprised 592,522 m 2 on average of the wetland complexes. The remaining areas of the wetland complexes not mapped as these three land cover types (Classes 2 and 1) had sparse vegetation or bare soil cover and commonly (greater than or equal to 67 percent of observations) had dry soil conditions. The investigation of land cover detailed in this report could inform future efforts to map land cover more precisely via higher resolution remote sensing or ground-based surveying or could be incorporated with other environmental monitoring data to characterize habitat and hydrology of the wetland complexes at Dixie Meadows.

Nevada↗

USGS2016

An index of floodplain surface complexity

Floodplain surface topography is an important component of floodplain ecosystems. It is the primary physical template upon which ecosystem processes are acted out, and complexity in this template can contribute to the high biodiversity and productivity of floodplain ecosystems. There has been a limited appreciation of floodplain surface complexity because of the traditional focus on temporal variability in floodplains as well as limitations to quantifying spatial complexity. An index of floodplain surface complexity (FSC) is developed in this paper and applied to eight floodplains from different geographic settings. The index is based on two key indicators of complexity, variability in surface geometry (VSG) and the spatial organisation of surface conditions (SPO), and was determined at three sampling scales. FSC, VSG, and SPO varied between the eight floodplains and these differences depended upon sampling scale. Relationships between these measures of spatial complexity and seven geomorphological and hydrological drivers were investigated. There was a significant decline in all complexity measures with increasing floodplain width, which was explained by either a power, logarithmic, or exponential function. There was an initial rapid decline in surface complexity as floodplain width increased from 1.5 to 5 km, followed by little change in floodplains wider than 10 km. VSG also increased significantly with increasing sediment yield. No significant relationships were determined between any of the four hydrological variables and floodplain surface complexity.

Hydrology and Earth System Sciences↗

USGS2019

Fish assemblages in a Mississippi reservoir mudflat with low structural complexity

In shallow reservoirs, seasonal water drawdowns expose littoral areas and over time produce barren mudflats. When flooded, mudflats provide homogeneous substrates, turbid water, and eroding shorelines of limited ecological value. We hypothesized that in mudflats structurally complex habitats are occupied by more fish, smaller fish of a larger range in sizes, more species, and fish assemblages that are different from those in simpler habitats. We tested these hypotheses over two consecutive years with fish collections made in sites with varying structural complexity. Results indicated that structural complexity harbors more fish in transects and enclosures. Structural complexity did not influence median length, but length range increased with structural complexity. Average species richness increased with structural complexity. Fish assemblage composition changed as structural complexity increased. The ability of cover to provide survival, growth, and carrying capacity benefits is fundamental to programs aimed at increasing structural complexity. Results suggest observed effects on fish assemblages can lead to such benefits. Considering mudflats are a major component of reservoirs, expand as reservoirs age, and there is a potential to exert meaningful change on fish assemblages of impounded rivers by managing mudflats, we suggest additional attention is needed to develop practical habitat restoration options.

Mississippi↗

USGS2020

Mapping the 3-D extent of the Stillwater Complex, Montana—Implications for potential platinum group element exploration and development

Geophysical models characterize the exposed and interpreted buried extent of the Stillwater Complex, critical for understanding the origin of the layered mafic intrusion and its associated high-grade platinum group element resources. The 3D models, constrained by gravity, magnetic, xenolith, seismic, borehole, and rock property data indicate that the likely maximum extent of the Stillwater Complex beneath Phanerozoic cover is ~10 times greater than its outcrop, ~2240 km 2 . The thickness values are poorly constrained but vary from ~7000 to 12,000 m, depending on crustal and mantle density variations and depths to the top of the lower crust and mantle. This thickness may include dense metasedimentary units of the basin into which the Stillwater Complex intruded. Using the modeled thickness results in a volume estimate of ~24,700 km 3 , albeit poorly constrained. New analyses of xenoliths from the Cretaceous Sliderock and Suzie Peak intrusions produce ages of 2706–2716 Ma, corresponding to the age of the Stillwater Complex, and 2813 Ma, corresponding to the age of Archean gneissic basement. Seismic reflectors in inferred Archean crystalline basement, possibly including the Stillwater Complex, dip ~25–30° north, with segments dipping as much as 70° north. Layered reflectors beneath the Phanerozoic sedimentary section and above the inferred Archean crystalline basement may represent metasedimentary units, perhaps a southern extension of the Mesoproterozoic Belt Basin. The potential field models and the seismic reflection data suggest that the Stillwater Complex was dipping northward prior to deposition of Cambrian strata, perhaps uplifted in the late Archean or Proterozoic as previously proposed, and that during Laramide times, faulting and intrusions highly disrupted the complex. Temperature measurements from boreholes help constrain the depths of feasible mining of the complex.

Montana↗

USGS2024

Deep magmatic staging chambers for crustal layered mafic intrusions: An example from the Bushveld Complex of southern Africa

The deep mafic magmatic staging chambers of layered mafic intrusions have been conjectured but not imaged. Their existence has long been postulated from geochemical models which require multiple magma injections from staging chambers to account for their multi-scale igneous layering and variations in sources and degrees of crustal contamination. For the Bushveld Complex of southern Africa, the world’s largest layered mafic intrusion, seismic receiver functions identify a diffuse crust-mantle transition beneath the Complex that suggests high velocity lower crust and/or uppermost lithospheric mantle. Here we present 3D gravity modelling of the Bushveld Complex that includes dense material at the crust-mantle boundary, imaging for the first time, remnants of magma staging chambers. They underlie the whole Bushveld Complex and extend westwards to the Molopo Farms Complex in Botswana. Feeders to the Bushveld Complex coincide with intersections of major faults like the Thabazimbi-Murchison Lineament and Sugarbush Fault with the staging chamber. This identification of magmatic staging chambers beneath the Bushveld relates to similar geophysically imaged lower crustal features beneath the Duluth and Stillwater Complexes. Comparison of seismic and gravity data with geochemical models from other complexes aid in development of models for the magmatic architecture of layered mafic intrusions in general.

Precambrian Research↗

USGS2020

Balancing transferability and complexity of species distribution models for rare species conservation

Aim Species distribution models (SDMs) are valuable for rare species conservation and are commonly used to extrapolate predictions of habitat suitability geographically to regions where species occurrence is unknown (i.e., transferability). Spatially structured cross-validation can be used to infer transferability, yet, few studies have evaluated how delineation of cross-validation folds affects model complexity and predictions. We developed SDMs using multiple cross-validation approaches to understand the implications for predicting habitat suitability for northern Idaho ground squirrels, a rare, federally threatened species that has been extensively surveyed in regions where known populations occur, resulting in >8000 presence locations. Location Idaho, USA. Methods We delineated cross-validation folds by mimicking the manner in which predictions would be geographically extrapolated or by using existing dispersal barriers. We varied the distance between, number, and directionality of folds. We conducted a grid search on statistical regularization parameters to optimize model complexity, covering a range of values exceeding that typically implemented. For each cross-validation approach, we selected optimal regularization and model complexity based on out-of-sample predictive ability. Results Delineation of cross-validation folds substantially affected resulting model complexity and extrapolated predictions. All cross-validation approaches resulted in models with apparently high out-of-sample predictive ability, yet optimal model complexity varied substantially among the approaches. Regularization demonstrated a noisy relationship between model complexity and prediction, where local optima in predictive performance were common at small values. Main conclusion Subtle modelling decisions can have large consequences for predictions of habitat suitability and transferability of SDMs. When transferability is the goal, cross-validation approaches should be considered carefully and mimic the manner in which spatial extrapolation will occur, else overly complex models with inflated assessments of predictive accuracy may result. Further, spatially structured cross-validation may not guard against over-parameterization, and assessing a broader range of regularization parameters may be necessary to optimize model complexity for transferability.

Idaho↗

USGS2014

What do we gain from simplicity versus complexity in species distribution models?

Species distribution models (SDMs) are widely used to explain and predict species ranges and environmental niches. They are most commonly constructed by inferring species' occurrence–environment relationships using statistical and machine-learning methods. The variety of methods that can be used to construct SDMs (e.g. generalized linear/additive models, tree-based models, maximum entropy, etc.), and the variety of ways that such models can be implemented, permits substantial flexibility in SDM complexity. Building models with an appropriate amount of complexity for the study objectives is critical for robust inference. We characterize complexity as the shape of the inferred occurrence–environment relationships and the number of parameters used to describe them, and search for insights into whether additional complexity is informative or superfluous. By building ‘under fit’ models, having insufficient flexibility to describe observed occurrence–environment relationships, we risk misunderstanding the factors shaping species distributions. By building ‘over fit’ models, with excessive flexibility, we risk inadvertently ascribing pattern to noise or building opaque models. However, model selection can be challenging, especially when comparing models constructed under different modeling approaches. Here we argue for a more pragmatic approach: researchers should constrain the complexity of their models based on study objective, attributes of the data, and an understanding of how these interact with the underlying biological processes. We discuss guidelines for balancing under fitting with over fitting and consequently how complexity affects decisions made during model building. Although some generalities are possible, our discussion reflects differences in opinions that favor simpler versus more complex models. We conclude that combining insights from both simple and complex SDM building approaches best advances our knowledge of current and future species ranges.

Ecography↗

USGS2006

Deciphering igneous and metamorphic events in high-grade rocks of the Wilmington complex, Delaware: Morphology, cathodoluminescence and backscattered electron zoning, and SHRIMP U-Pb geochronology of zircon and monazite

High-grade rocks of the Wilmington Complex, northern Delaware and adjacent Maryland and Pennsylvania, contain morphologically complex zircons that formed through both igneous and metamorphic processes during the development of an island-arc complex and suturing of the arc to Laurentia. The arc complex has been divided into several members, the protoliths of which include both intrusive and extrusive rocks. Metasedimentary rocks are interlayered with the complex and are believed to be the infrastructure upon which the arc was built. In the Wilmingto n Complex rocks, both igneous and metamorphic zircons occur as elongate and equant forms. Chemical zoning, shown by cathodoluminescence (CL), includes both concentric, oscillatory patterns, indicative of igneous origin, and patchwork and sector patterns, suggestive of metamorphic growth. Metamorphic monazites are chemically homogeneous, or show oscillatory or spotted chemical zoning in backscattered electron images. U-Pb geochronology by sensitive high resolution ion microprobe (SHRIMP) was used to date complexly zoned zircon and monazite. All but one member of the Wilmington Complex crystallized in the Ordovician between ca. 475 and 485 Ma; these rocks were intruded by a suite of gabbro-to-granite plutonic rocks at 434 ?? Ma. Detrital zircons in metavolcanic and metasedimentary units were derived predominantly from 0.9 to 1.4 Ga (Grenvillian) basement, presumably of Laurentian origin. Amphibolite to granulite facies metamorphism of the Wilmington Complex, recorded by ages of metamorphic zircon (428 ?? 4 and 432 ?? 6 Ma) and monazite (429 ?? 2 and 426 ?? 3 Ma), occurred contemporaneously with emplacement of the younger plutonic rocks. On the basis of varying CL zoning patterns and external morphologies, metamorphic zircons formed by different processes (presumably controlled by rock chemistry) at slightly different times and temperatures during prograde metamorphism. In addition, at least three other thermal episodes are recorded by monazite growth at 447 ?? 4, 411 ?? 3, and 398 ?? 3 Ma. ?? 2006 Geological Society of America.

Geological Society of America Bulletin↗

USGS2010

Detrital zircon provenance from three turbidite depocenters of the Middle-Upper Triassic Songpan-Ganzi complex, central China: Record of collisional tectonics, erosional exhumation, and sediment production

To test the idea that the voluminous upper Middle to Upper Triassic turbidite strata in the Songpan-Ganzi complex of central China archive a detrital record of Dabie ultrahigh-pressure (UHP) terrane unroofing, we report 2080 single detrital U-Pb zircon ages by sensitive high-resolution ion microprobe-reverse geometry (SHRIMP-RG) and laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) analysis from 29 eastern Songpan-Ganzi complex sandstone samples. Low (<0.07) Th/U zircons, consistent with crystallization under UHP conditions, are rare in eastern Songpan-Ganzi complex zircon, and U-Pb ages of low Th/U zircons are incompatible with a Dabie terrane source. An unweighted pair group method with arithmetic mean nearest-neighbor analysis of Kolmogorov-Smirnov two-sample test results reveals that the eastern Songpan-Ganzi complex is not a single contiguous turbidite system but is instead composed of three subsidiary depocenters, each associated with distinct sediment sources. The northeastern depocenter contains zircon ages characterized by Paleozoic and bimodally distributed Precambrian zircon populations, which, together with south-to southeast-directed paleocurrent data, indicate derivation from the retro-side of the Qinling-Dabie (Q-D) collisional orogen wedge. In the central depocenter, the dominantly Paleozoic detrital zircon signature and south-to southwest-oriented paleocurrent indicators reflect a profusion of Paleozoic zircon grains. These data are interpreted to reflect an influx of material derived from erosion of Paleozoic supra-UHP rocks of the Dabie terrane in the eastern Qinling-Dabie orogen, which we speculate may have been enhanced by development of a monsoonal climate. This suggests that erosional unroofing played a significant role in the initial phase of UHP exhumation and likely influenced the petrotectonic and structural evolution of the Qinling-Dabie orogen, as evidenced by compressed Triassic isotherms/grads reported in the Huwan shear zone that bounds the Dabie terrane to the north. The central depocenter deposits reflect a later influx of bimodally distributed Precambrian zircon, signifying either a decrease in the influx of Paleozoic zircon grains due to stalled UHP exhumation and/or dilution of the same influx of Paleozoic zircons by spilling of Precambrian zircon from the northeastern depocenter into the central depocenter basin, perhaps due to infilling and bypass of sediment from the northern depocenter or due to initial collapse and constriction of the eastern Songpan-Ganzi complex basin. The southeastern depocenter of the eastern Songpan-Ganzi complex bears significant Paleozoic, Neoproterozoic, and Paleoproterozoic zircon populations derived from the South China block and Yidun arc complex, likely recording nascent uplift of the Longmenshan deformation belt due to impingement of the Yidun arc complex upon the western margin of the South China block. ?? 2010 Geological Society of America.

Geological Society of America Bulletin↗

USGS2011

The effect of structural complexity, prey density, and "predator-free space" on prey survivorship at created oyster reef mesocosms

Interactions between predators and their prey are influenced by the habitat they occupy. Using created oyster ( Crassostrea virginica ) reef mesocosms, we conducted a series of laboratory experiments that created structure and manipulated complexity as well as prey density and “predator-free space” to examine the relationship between structural complexity and prey survivorship. Specifically, volume and spatial arrangement of oysters as well as prey density were manipulated, and the survivorship of prey (grass shrimp, Palaemonetes pugio ) in the presence of a predator (wild red drum, Sciaenops ocellatus ) was quantified. We found that the presence of structure increased prey survivorship, and that increasing complexity of this structure further increased survivorship, but only to a point. This agrees with the theory that structural complexity may influence predator-prey dynamics, but that a threshold exists with diminishing returns. These results held true even when prey density was scaled to structural complexity, or the amount of “predator-free space” was manipulated within our created reef mesocosms. The presence of structure and its complexity (oyster shell volume) were more important in facilitating prey survivorship than perceived refugia or density-dependent prey effects. A more accurate indicator of refugia might require “predator-free space” measures that also account for the available area within the structure itself (i.e., volume) and not just on the surface of a structure. Creating experiments that better mimic natural conditions and test a wider range of “predator-free space” are suggested to better understand the role of structural complexity in oyster reefs and other complex habitats.

PLoS ONE↗

USGS2002

Surface complexation and precipitate geometry for aqueous Zn(II) sorption on ferrihydrite I: X-ray absorption extended fine structure spectroscopy analysis

“Two-line” ferrihydrite samples precipitated and then exposed to a range of aqueous Zn solutions (10 −5 to 10 −3 M), and also coprecipitated in similar Zn solutions (pH 6.5), have been examined by Zn and Fe K-edge X-ray absorption spectroscopy. Typical Zn complexes on the surface have Zn-O distances of 1.97(.02) Å and coordination numbers of about 4.0(0.5), consistent with tetrahedral oxygen coordination. This contrasts with Zn-O distances of 2.11(.02) Å and coordination numbers of 6 to 7 in the aqueous Zn solutions used in sample preparation. X-ray absorption extended fine structure spectroscopy (EXAFS) fits to the second shell of cation neighbors indicate as many as 4 Zn-Fe neighbors at 3.44(.04) Å in coprecipitated samples, and about two Zn-Fe neighbors at the same distance in adsorptionsamples. In both sets of samples, the fitted coordination number of second shell cations decreases as sorption density increases, indicating changes in the number and type of available complexing sites or the onset of competitive precipitation processes. Comparison of our results with the possible geometries for surface complexes and precipitates suggests that the Zn sorption complexes are inner sphere and at lowest adsorption densities are bidentate, sharing apical oxygens with adjacent edge-sharing Fe(O,OH) 6 octahedra. Coprecipitation samples have complexes with similar geometry, but these are polydentate, sharing apices with more than two adjacent edge-sharing Fe(O,OH) 6 polyhedra. The results are inconsistent with Zn entering the ferrihydrite structure (i.e., solid solution formation) or formation of other Zn-Fe precipitates. The fitted Zn-Fe coordination numbers drop with increasing Zn density with a minimum of about 0.8(.2) at Zn/(Zn + Fe) of 0.08 or more. This change appears to be attributable to the onset of precipitation of zinc hydroxide polymers with mainly tetrahedral Zn coordination. At the highest loadings studied, the nature of the complexes changes further, and a second type of precipitate forms. This has a structure based on a brucite layer topology, with mainly octahedral Zn coordination. Amorphous zinc hydroxide samples prepared for comparison had a closely similar local structure. Analysis of the Fe K-edge EXAFS is consistent with surface complexation reactions and surface precipitation at high Zn loadings with little or no Fe-Zn solid solution formation. The formation of Zn-containing precipitates at solution conditions two or more orders of magnitude below their solubility limit is compared with other sorption and spectroscopic studies that describe similar behavior.

Geochimica et Cosmochimica Acta↗

USGS2002

Effect of adsorbed metals ions on the transport of Zn- and Ni-EDTA complexes in a sand and gravel aquifer

Adsorption, complexation, and dissolution reactions strongly influenced the transport of metal ions complexed with ethylenediaminetetraacetic acid (EDTA) in a predominantly quartz-sand aquifer during two tracer tests conducted under mildly reducing conditions at pH 5.8 to 6.1. In tracer test M89, EDTA complexes of zinc (Zn) and nickel (Ni), along with excess free EDTA, were injected such that the lower portion of the tracer cloud traveled through a region with adsorbed manganese (Mn) and the upper portion of the tracer cloud traveled through a region with adsorbed Zn. In tracer test S89, Ni- and Zn-EDTA complexes, along with excess EDTA complexed with calcium (Ca), were injected into a region with adsorbed Mn. The only discernable chemical reaction between Ni-EDTA and the sediments was a small degree of reversible adsorption leading to minor retardation. In the absence of adsorbed Zn, the injected Zn was displaced from EDTA complexes by iron(III) [Fe(III)] dissolved from the sediments. Displacement of Zn by Fe(III) on EDTA became increasingly thermodynamically favorable with decreasing total EDTA concentration. The reaction was slow compared to the time-scale of transport. Free EDTA rapidly dissolved aluminum (Al) from the sediments, which was subsequently displaced slowly by Fe. In the portion of tracer cloud M89 that traveled through the region contaminated with adsorbed Zn, little displacement of Zn complexed with EDTA was observed, and Al was rapidly displaced from EDTA by Zn desorbed from the sediments, in agreement with equilibrium calculations. In tracer test S89, desorption of Mn dominated over the more thermodynamically favorable dissolution of Al oxyhydroxides. Comparison with results from M89 suggests that dissolution of Al oxyhydroxides in coatings on these sediment grains by Ca-EDTA was rate-limited whereas that by free EDTA reached equilibrium on the time-scale of transport. Rates of desorption are much faster than rates of dissolution of Fe oxyhydroxides from sediment-grain surfaces and, therefore, adsorbed metal ions can strongly influence the speciation of ligands like EDTA in soils and sediments, especially over small temporal and spatial scales.

Geochimica et Cosmochimica Acta↗