Search USGS⌕ Search

SEARCH · Search USGS

Results for “Complexity”

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

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

At least 1,009 records · Page 56Linked to original sources

Genesis of a zoned granite stock, Seward Peninsula, Alaska

A composite epizonal stock of biotite granite has intruded a diverse assemblage of metamorphic rocks in the Serpentine Hot Springs area of north-central Seward Peninsula, Alaska. The metamorphic rocks include amphibolite-facies orthogneiss and paragneiss, greenschist-facies fine-grained siliceous and graphitic metasediments, and a variety of carbonate rocks. Lithologic units within the metamorphic terrane trend generally north-northeast and dip moderately toward the southeast. Thrust faults locally juxtapose lithologic units in the metamorphic assemblage, and normal faults displace both the metamorphic rocks and some parts of the granite stock. The gneisses and graphitic metasediments are believed to be late Precambrian in age, but the carbonate rocks are in part Paleozoic. Dating by the potassium-argon method indicates that the granite stock is Late Cretaceous. The stock has sharp discordant contacts, beyond which is a well-developed thermal aureole with rocks of hornblende hornfels facies. The average mode of the granite is 29 percent plagioclase, 31 percent quartz, 36 percent K-feldspar, and 4 percent biotite. Accessory minerals include apatite, magnetite, sphene, allanite, and zircon. Late-stage or deuteric minerals include muscovite, fluorite, tourmaline, quartz, and albite. The stock is a zoned complex containing rocks with several textural facies that are present in four partly concentric zones. Zone 1 is a discontinuous border unit, containing fine- to coarse-grained biotite granite, that grades inward into zone 2. Zone 2 consists of porphyritic biotite granite with oriented phenocrysts of pinkish-gray microcline in a coarse-grained equigranular groundmass of plagioclase, quartz, and biotite. It is in sharp, concordant to discordant contact with rocks of zone 3. Zone 3 consists of seriate-textured biotite granite that has been intruded by bodies of porphyritic biotite granite containing phenocrysts of plagioclase, K-feldspar, quartz, and biotite in an aplitic groundmass. Flow structures, pegmatite and aplite segregations, and miarolitic cavities are common in the seriate-textured granite. Zone 4, which forms the central part of the complex, consists of fine- to medium-grained biotite granite and locally developed leucogranite. Small miarolitic cavities are common within it. Eight textural facies have been defined within the complex, and mineralogic, petrographic, modal, and chemical variations are broadly systematic within the facies sequence. Study of these variations shows that the gradational facies of zones l and 2 systematically shift toward more mafic compositions inward within the complex. Seriate-textured rocks of zone 3 are similar in composition to those of zone 2, but porphyritic rocks of zone 3 and rocks of zone 4 mark shifts to more felsic compositions. These late-crystallizing felsic rocks are products of an interior residual magma system. This system was enriched in water and certain trace elements including tin, lithium, niobium, lead, and zinc. The complex as a whole has higher concentrations of these elements than many other granites. The nature of this geochemical specialization is particularly well demonstrated by the trace-element composition of biotite. The crystallization history of the pluton was complex. The available data suggest that this history could have included: (1) chilling and metasomatic alteration adjacent to the contact, (2) in-situ crystallization in several marginal facies accompanied by some transfer of residual constituents toward interior parts of the pluton, (3) slight upward displacement of magma that was subjacent to the crystallized walls, accompanied by disequilibrium crystallization and local vapor saturation, (4) upward displacement of part of the residual water-rich interior magma, accompanied by rapid loss of a separated vapor phase, and (5) displacement of the margins of the pluton by normal faults, accompanied by loss of an exsolved vapor phase from th

Open-File Report↗

Alkalic rocks and resources of thorium and associated elements in the Powderhorn District, Gunnison County, Colorado

Alkalic igneous rocks and related concentrations of thorium, niobium, rare-earth elements, titanium, and other elements have long been known in the Powderhorn mining district and have been explored intermittently for several decades. The deposits formed chiefly about 570 m.y. (million years) ago in latest Precambrian or Early Cambrian time. They were emplaced in lower Proterozoic (Proterozoic X) metasedimentary, metavolcanic, and plutonic rocks. The complex of alkalic rocks of Iron Hill occupies 31 km 2 (square kilometers) and is composed of pyroxenite, uncompahgrite, ijolite, nepheline syenite, and carbonatite, in order of generally decreasing age. Fenite occurs in a zone, in places more than 0.6 km (kilometer) wide, around a large part of the margin of the complex and adjacent to alkalic dikes intruding Precambrian host rock. The alkalic rocks have a radioactivity, chiefly due to thorium, greater than that of the surrounding Powderhorn Granite (Proterozoic X) and metamorphic rocks. The pyroxenite, uncompahgrite, ijolite, and nepheline syenite, which form more than 80 percent of the complex, have fairly uniform radioactivity. Radioactivity in the carbonatite stock, carbonatite dikes, and the carbonatite-pyroxenite mixed rock zone, however, generally exceeds that in the other rocks of the complex. The thorium concentrations in the Powderhorn district occur in six types of deposits: thorite veins, a large massive carbonatite body, carbonatite dikes, trachyte dikes, magnetite-ilmeniteperovskite dikes or segregations, and disseminations in small, anomalously radioactive plutons chiefly of granite or quartz syenite that are older than rocks of the alkalic complex. The highest grade thorium concentrations in the district are in veins that commonly occur in steeply dipping, crosscutting shear or breccia zones in the Precambrian rocks. They range in thickness from a centimeter or less to 5 m (meters) and are as much as 1 km long. The thorite veins are composed chiefly of potassic feldspar, white to smoky quartz, calcite, barite, goethite, and hematite, and also contain thorite, jasper, magnetite, pyrite, galena, chalcopyrite, sphalerite, synchysite, apatite, fluorite, biotite, sodic amphibole, rutile, monazite, bastnaesite, and vanadinite. The Th0 2 content of the thorite veins ranges from less than 0.01 percent to as much as 4.9 percent in high-grade samples. The Th0 2 content is generally less than 1 percent, however, and is only 0.05 to 0.1 percent in many of the veins examined in the district. Samples of the dolomitic carbonatite of Iron Hill mostly range from 3 to 145 ppm (parts per million) thorium. Thirty samples of the carbonatite dikes, the most radioactive rocks within the complex of Iron Hill, contain about 30 to 3,200 ppm thorium and a trace to about 1.5 percent rare-earth oxides. The magnetite-ilmenite-perovskite rocks have a radioactivity of 2 to 12 times the background of Precambrian granite that is attributable chiefly to thorium substitution for calcium in the perovskite. In two analyses the perovskite contains 0.12 and 0.15 percent Th0 2 . Trachyte dikes as much as 25 m thick cut the Precambrian rocks; their radioactivity is generally about two to four times the background of typical Precambrian granite, is locally higher, but is low relative to other types of thorium concentrations. A finegrained granite that is anomalously radioactive occurs in thick, dikelike plutons as much as 1.2 km wide, or more. The thorium content varies widely within the granite bodies. Eight samples of the granite contain 32 to 281 ppm thorium (averaging 115 ppm). The economic potential of thorium in the Powderhorn district is related in part to other elements such as niobium, titanium, iron, and rare earths. The proportions of niobium and rare earths to thorium vary in different parts of the district. Within the carbonatite body of Iron Hill, the Nb 2 0 5 content greatly exceeds Th0 2 , but the Th0 2 -Nb 2 0 5

Professional Paper↗

Reconnaissance geologic map of the Kuskokwim Bay region, southwest Alaska

The rocks of the map area range from Proterozoic age metamorphic rocks of the Kanektok metamorphic complex (Kilbuck terrane) to Quaternary age mafic volcanic rocks of Nunivak Island. The map area encompasses much of the type area of the Togiak-Tikchik Complex. The geologic maps used to construct this compilation were, for the most part, reconnaissance studies done in the time period from the 1950s to 1990s. Pioneering work in the map area by J.M. Hoare and W.L. Coonrad forms the basis for much of this map, either directly or as the stepping off point for later studies compiled here. Physiographically, the map area ranges from glaciated mountains, as much as 1,500 m high, in the Ahklun Mountains to the coastal lowlands of northern Bristol Bay and the Kuskokwim River delta. The mountains and the finger lakes (drowned fiords) on the east have been strongly affected by Pleistocene and Holocene glaciation. Within the map area are a number of major faults. The Togiak-Tikchik Fault and its extension to the northeast, the Holitna Fault, are considered extensions of the Denali fault system of central Alaska. Other sub-parallel faults include the Golden Gate, Sawpit, Goodnews, and East Kulukak Faults. Northwest-trending strike-slip faults crosscut and offset northeast-trending fault systems. Rocks of the area are assigned to a number of distinctive lithologic packages. Most distinctive among these packages are the high-grade metamorphic rocks of the Kanektok metamorphic complex or Kilbuck terrane, composed of a high-grade metamorphic orthogneiss core surrounded by greenschist and amphibolite facies schist, gneiss, and rare marble and quartzite. These rocks have yielded radiometric ages strongly suggestive of a 2.05 Ga emplacement age. Poorly known Paleozoic rocks, including Ordovician to Devonian and Permian limestone, are found east of the Kanektok metamorphic complex. A Triassic(?) ophiolite complex is on the southeast side of Kuskokwim Bay; otherwise only minor Triassic rock units are known. The most widespread rocks of the area are Jurassic and Early Cretaceous(?) volcanic and volcaniclastic rocks. The Kuskokwim Group flysch is restricted largely to the northeast part of the map area. It consists primarily of shelf and minor nearshore facies rocks. Primarily exposed in the lowlands west of the Ahklun Mountains, extensive latest Tertiary and Quaternary alkalic basalt flows and lesser pyroclastic rocks form much of the bedrock of the remaining area. On Saint Matthew Island, Cretaceous volcanic and pyroclastic rocks occur that are not found elsewhere within the map area. The Kuskokwim Group and older rocks, including on Saint Matthew Island, but not the Kanektok metamorphic complex, are intruded by widely dispersed Late Cretaceous and (or) Early Tertiary granitic rocks. Much of the lowland area is mantled by unconsolidated deposits that include glacial, alluvial and fluvial, marine, estuarine, and eolian deposits. These formed during several episodes of Quaternary glaciation.

Alaska↗

Geologic map and structure sections along the southern part of the Bartlett Springs Fault Zone and adjacent areas from Cache Creek to Lake Berryessa, northern Coast Ranges, California

Introduction Located in the Coast Ranges of northern California, the Bartlett Springs Fault Zone is the easternmost fault in the San Andreas Fault system in northern California. The fault is a right-lateral, strike-slip structure considered capable of producing an earthquake of moment magnitude 7. The purpose of this mapping is to better characterize the geology and earthquake hazards associated with the southern part of the Bartlett Springs Fault Zone and to help identify any evidence of active uplift on the faults bounding the Coast Ranges. Although the area immediately surrounding the Bartlett Springs Fault Zone is sparsely populated, its southern segment presents a potential seismic hazard to northern California communities as far away as the San Francisco Bay region and Sacramento. There are also nearby water resources, mineral resources, and public lands used for public recreation. The Coast Ranges of northern California are a series of northwest-southeast-oriented mountain ranges and valleys located north of the San Francisco Bay region, between the Pacific Ocean to the west and the Sacramento Valley to the east. The region has rugged terrain, high mountain peaks that reach more than 2,400 meters above sea level, isolated and narrow valley bottoms on which most human settlements are located, and large drainage systems that tend to follow the northwest-southeast-oriented topographic grain. The physiographic character of the region is shaped by its bedrock geology, deformational history, and active faulting. The basement rocks of the northern Coast Ranges consist of the Franciscan Complex and the Great Valley complex, the latter of which consists of two informal units, the Coast Range ophiolite and the Great Valley sequence. The Franciscan Complex and the Great Valley complex are in structural contact along the Coast Range Fault, a regional-scale structure and fundamental crustal boundary. The Franciscan Complex and the Great Valley complex are superposed by active, northwest-southeast-striking strike-slip faults that are associated with seismicity swarms. These active strike-slip faults can produce moderate to large earthquakes that have moment magnitudes of 7–8. In places, these active structures bound large ranges and valleys, suggesting that much of the modern topographic expression is the result of active deformation processes. This report contains new 1:24,000-scale geologic mapping along the southern part of the Bartlett Springs Fault Zone between Clear Lake and Lake Berryessa. The map area spans 738 square kilometers in northern Napa County, southern Lake County, and parts of Yolo and Colusa Counties. The south and east borders of the map are 90 kilometers north of San Francisco and 70 kilometers west of Sacramento, respectively. The map area is within the Knoxville mining district, which has a history of mercury and gold mining dating back to the mid-19th century. The two main towns in the region, Lower Lake and Clearlake, California, are west-northwest of the map area. Approximately 71,000 people live in the cities and rural communities located within a 40-kilometer radius of the center of the map area. The bedrock geology, cross sections, and structural data presented herein are critical for evaluating the long-term evolution of the Bartlett Springs Fault Zone. This work will supplement studies on local seismic hazards, liquefaction potential, landslide hazards, earthquake geology, natural resources, groundwater resources, engineering geology, and tectonic history by providing the background information for site-specific investigations on these subjects.

California↗

Vegetation changes associated with a population irruption by Roosevelt elk

Interactions between large herbivores and their food supply are central to the study of population dynamics. We assessed temporal and spatial patterns in meadow plant biomass over a 23-year period for meadow complexes that were spatially linked to three distinct populations of Roosevelt elk (Cervus elaphus roosevelti) in northwestern California. Our objectives were to determine whether the plant community exhibited a tolerant or resistant response when elk population growth became irruptive. Plant biomass for the three meadow complexes inhabited by the elk populations was measured using Normalized Difference Vegetation Index (NDVI), which was derived from Landsat 5 Thematic Mapper imagery. Elk populations exhibited different patterns of growth through the time series, whereby one population underwent a complete four-stage irruptive growth pattern while the other two did not. Temporal changes in NDVI for the meadow complex used by the irruptive population suggested a decline in forage biomass during the end of the dry season and a temporal decline in spatial variation of NDVI at the peak of plant biomass in May. Conversely, no such patterns were detected in the meadow complexes inhabited by the nonirruptive populations. Our findings suggest that the meadow complex used by the irruptive elk population may have undergone changes in plant community composition favoring plants that were resistant to elk grazing.

Ecology and Evolution↗

Environmentally relevant chemical mixtures of concern in waters of United States tributaries to the Great Lakes

The North American Great Lakes are a vital natural resource that provide fish and wildlife habitat, as well as drinking water and waste assimilation services for millions of people. Tributaries to the Great Lakes receive chemical inputs from various point and nonpoint sources, and thus are expected to have complex mixtures of chemicals. However, our understanding of the co‐occurrence of specific chemicals in complex mixtures is limited. To better understand the occurrence of specific chemical mixtures in the US Great Lakes Basin, surface water from 24 US tributaries to the Laurentian Great Lakes was collected and analyzed for diverse suites of organic chemicals, primarily focused on chemicals of concern (e.g., pharmaceuticals, personal care products, fragrances). A total of 181 samples and 21 chemical classes were assessed for mixture compositions. Basin wide, 1664 mixtures occurred in at least 25% of sites. The most complex mixtures identified comprised 9 chemical classes and occurred in 58% of sampled tributaries. Pharmaceuticals typically occurred in complex mixtures, reflecting pharmaceutical‐use patterns and wastewater facility outfall influences. Fewer mixtures were identified at lake or lake‐influenced sites than at riverine sites. As mixture complexity increased, the probability of a specific mixture occurring more often than by chance greatly increased, highlighting the importance of understanding source contributions to the environment. This empirically based analysis of mixture composition and occurrence may be used to focus future sampling efforts or mixture toxicity assessments.

Integrated Environmental Assessment and Management↗

Bridging the gap between salmon spawner abundance and marine nutrient assimilation by juvenile salmon: Seasonal cycles and landscape effects at the watershed scale

Anadromous Pacific salmon are semelparous, and resource subsidies from spawning adults (marine-derived nutrients, or MDN) benefit juvenile salmonids rearing in freshwater. However, it remains unclear how MDN assimilation relates to spawner abundance within a watershed. To address this, we examined seasonal, watershed-scale patterns of MDN assimilation in rearing coho ( Oncorhynchus kisutch ) and Chinook ( O. tshawytscha ) salmon and compared it with spawner biomass and landscape features in a western Alaska watershed with contrasting structural complexity in two sub-drainages. Adult salmon biomass density was estimated from escapement and spawner distribution data, and MDN assimilation in juvenile salmon was estimated via stable isotopes. In the North River, MDN assimilation was lowest in early summer, prior to annual spawning migrations, increased after spawning, and peaked in late winter. In the more complex mainstem Unalakleet River, MDN assimilation was higher but varied minimally from summer through fall before increasing in late fall and winter. Summer MDN assimilation, prior to salmon spawning, was primarily a function of habitat complexity, where MDN was highest in sloughs and the more complex mainstem river. After salmon spawned, fall MDN assimilation was a function of adult pink and Chinook salmon biomass as well as MDN assimilation that occurred prior to spawning (that is, summer MDN), but unrelated to total summer biomass (all salmon species biomass combined). Thus, MDN assimilation by juvenile salmon in the fall was a function of species-specific adult spawner abundance but seasonal patterns of MDN assimilation were masked in complex habitat where summer MDN assimilation remained high.

Alaska↗

Testing the suitability of geologic frameworks for extrapolating hydraulic properties across regional scales

The suitability of geologic frameworks for extrapolating hydraulic conductivity ( K ) to length scales commensurate with hydraulic data is difficult to assess. A novel method is presented for evaluating assumed relations between K and geologic interpretations for regional-scale groundwater modeling. The approach relies on simultaneous interpretation of multiple aquifer tests using alternative geologic frameworks of variable complexity, where each framework is incorporated as prior information that assumes homogeneous K within each model unit. This approach is tested at Pahute Mesa within the Nevada National Security Site (USA), where observed drawdowns from eight aquifer tests in complex, highly faulted volcanic rocks provide the necessary hydraulic constraints. The investigated volume encompasses 40 mi 3 (167 km 3 ) where drawdowns traversed major fault structures and were detected more than 2 mi (3.2 km) from pumping wells. Complexity of the five frameworks assessed ranges from an undifferentiated mass of rock with a single unit to 14 distinct geologic units. Results show that only four geologic units can be justified as hydraulically unique for this location. The approach qualitatively evaluates the consistency of hydraulic property estimates within extents of investigation and effects of geologic frameworks on extrapolation. Distributions of transmissivity are similar within the investigated extents irrespective of the geologic framework. In contrast, the extrapolation of hydraulic properties beyond the volume investigated with interfering aquifer tests is strongly affected by the complexity of a given framework. Testing at Pahute Mesa illustrates how this method can be employed to determine the appropriate level of geologic complexity for large-scale groundwater modeling.

Nevada↗

Ion association in natural brines

Natural brines, both surface and subsurface, are highly associated aqueous solutions. Ion complexes in brines may be ion pairs in which the cation remains fully hydrated and the bond between the ions is essentially electrostatic, or coordination complexes in which one or more of the hydration water molecules are replaced by covalent bonds to the anion. Except for Cl − , the major simple ions in natural brines form ion pairs; trace and minor metals in brines form mainly coordination complexes. Limitations of the Debye-Hückel relations for activity coefficients and lack of data on definition and stability of all associated species in concentrated solutions tend to produce underestimates of the degree of ion association, except where the brines contain a very high proportion of Cl − . Data and calculations on closed basin brines of highly varied composition have been coupled with electrode measurements of single-ion activities in an attempt to quantify the degree of ion association. Such data emphasize the role of magnesium complexes. Trace metal contents of closed basin brines are related to complexes formed with major anions. Alkaline sulfo- or chlorocarbonate brines (western Great Basin) carry significant trace metal contents apparently as hydroxides or hydroxy polyions. Neutral high chloride brines (Bonneville Basin) are generally deficient in trace metals. With a knowledge of the thermodynamic properties of a natural water, many possible reactions with other phases (solids, gases, other liquids) may be predicted. A knowledge of these reactions is particularly important in the study of natural brines which may be saturated with many solid phases (silicates, carbonates, sulfates, etc.), which may have a high pH and bring about dissolution of other phases (silica, amphoteric hydroxides, CO 2 , etc.), and which because of their high density may form relatively stable interfaces with dilute waters.

Chemical Geology↗

The dissolution of quartz in dilute aqueous solutions of organic acids at 25°C

The dissolution of quartz in dilute aqueous solutions of organic acids at 25° and standard pressure was investigated by the batch dissolution method. The bulk dissolution rate of quartz in 20 mmole/Kg citrate solutions at pH 7 was 8 to 10 times faster than that in pure water. After 1750 hours the concentration of dissolved silica in the citrate solution was 167 μmole/Kg compared to 50 μmole/Kg in water and a 20 mmole/Kg solution of acetate at pH 7. Solutions of salicylic, oxalic, and humic acids also accelerated the dissolution of quartz in aqueous solution at pH 7. The rate of dissolution in organic acids decreased sharply with decreasing pH. The possibility of a silica-organic acid complex was investigated using UV-difference spectroscopy. Results suggest that dissolved silica is complexed by citrate, oxalate and pyruvate at pH 7 by an electron-donor acceptor complex, whereas no complexation occurs between silica and acetate, lactate, malonate, or succinate. Three models are proposed for the solution and surface complexation of silica by organic acid anions which result in the accelerated dissolution and increased solubility of quartz in organic rich water.

Geochimica et Cosmochimica Acta↗

Enrichment of trace elements in garnet amphibolites from a paleo-subduction zone: Catalina Schist, southern California

The abundance, P - T stability, solubility, and element-partitioning behavior of minerals such as rutile, garnet, sphene, apatite, zircon, zoisite, and allanite are critical variables in models for mass transfer from the slab to the mantle wedge in deep regions of subduction zones. The influence of these minerals on the composition of subduction-related magmas has been inferred (and disputed) from inverse modelling of the geochemistry of island-arc basalt, or by experiment. Although direct samples of the dehydration + partial-melting region of a mature subduction zone have not been reported from subduction complexes, garnet amphibolites from melanges of circumpacific and Caribbean blueschist terranes reflect high T (>600°C) conditions in shallower regions. Such rocks record geochemical processes that affected deep-seated, high- T portions of paleo-subduction zones. In the Catalina Schist, a subduction-zone metamorphic terrane of southern California, metasomatized and migmatitic garnet amphibolites occur as blocks in a matrix of meta-ultramafic rocks. This mafic and ultramafic complex may represent either slab-derived material accreted to the mantle wedge of a nascent subduction zone or a portion of a shear zone closely related to the slab-mantle wedge contact, or both. The trace-element geochemistry of the complex and the distribution of trace elements among the minerals of garnet amphibolites were studied by INAA, XRF, electron microprobe, and SEM. In order of increasing alteration from a probable metabasalt protolith, three common types of garnet amphibolite blocks in the Catalina Schist are: (1) non-migmatitic, clinopyroxene-bearing blocks, which are compositionally similar to MORB that has lost an albite component; (2) garnet-amphibolite blocks, which have rinds that reflect local interaction between metabasite, metaperidotite, and fluid; and (3) migmatites that are extremely enriched in Th, HFSE, LREE, and other trace elements. These trace-element enrichments are mineralogically controlled by rutile, garnet, sphene, apatite, zircon, zoisite, and allanite. Alkali and alkaline earth elements are much less enriched in the solid assemblage, and thus appear to be decoupled from the other elements in the inferred metasomatic process(es). The compositions of migmatitic garnet amphibolite blocks seem to complement that of “average” island-arc tholeiite. Trace-element metasomatism reflects fluid-solid, rather than melt-solid, interaction. The metasomatic effects indicate that H 2 O-rich fluid, perhaps with a significant component of Na-Al silicate and alkalis, carried Th, U, Sr, REE, and HFSE. Fractionations of LREE in migmatites resemble those of migmatitic metasedimentary rocks underlying the mafic and ultramafic complex. “Exotic” LREE deposited in allanite in migmatites could have been derived from fluids in equilibrium with subducted sediment. If the paleo-subduction zone represented by the mafic and ultramafic complex of the Catalina Schist had continued its thermal and fluid evolution, a selvage of similarly enriched rocks might have been generated along the slab-mantle wedge contact between ~30 and 85 km depth. Rocks affected by “subduction-zone metasomatism,” although rarely recognized at the surface, could be volumetrically significant products of the initiation of subduction and may prove to be geochemical probes of convergent margins that approach the significance of xenoliths in the study of other magmatic environments.

Geochimica et Cosmochimica Acta↗

Comparison of the partitioning behaviours of yttrium, rare earth elements, and titanium between hydrogenetic marine ferromanganese crusts and seawater

In order to evaluate details of the partitioning behaviours of Y, rare earth elements (REEs), and Ti between inorganic metal oxide surfaces and seawater, we studied the distribution of these elements in hydrogenetic marine ferromanganese (Fe-Mn) crusts from the Central Pacific Ocean. Nonphosphatized Fe-Mn crusts display shale-normalized rare earths and yttrium (REY SN ) patterns (Y inserted between Dy and Ho) that are depleted in light REEs (LREEs) and which show negative anomalies for Y sn , and positive anomalies for La SN , Eu SN , Gd SN , and in most cases, Ce sn . They show considerably smaller Y/ Ho ratios than seawater or common igneous and clastic rocks, indicating that Y and Ho are fractionated in the marine environment. Compared to P-poor crusts, REY SN patterns of phosphatized Fe-Mn crusts are similar, but yield pronounced positive Y sn anomalies, stronger positive La SN anomalies, and enrichment of the HREEs relative to the MREEs. The data suggest modification of REY during phosphatization and indicate that studies requiring primary REY distributions or isotopic ratios should be restricted to non-phosphatized (layers of) Fe-Mn crusts. Apparent bulk coefficients, K d m , describing trace metal partitioning between nonphosphatized hydrogenetic Fe-Mn crusts and seawater, are similar for Pr to Eu and decrease for Eu to Yb. Exceptionally high values of K D Ce , which are similar to those of Ti, result from oxidative scavenging of Ce and support previous suggestions that Ce (IV) is a hydroxide-dominated element in seawater. Yttrium and Gd show lower K D values than their respective neighbours in the REY series. Results of modelling the exchange equilibrium between REY dissolved in seawater and REY sorbed on hydrous Fe-Mn oxides corroborate previous studies that suggested the surface complexation of REY can be approximated by their first hydroxide binding constant. Negative “anomalies” occur for stabilities of bulk surface complexes of Gd, La, and particularly Y. The differences in inorganic surface complex stability between Y and Ho and between Gd and its REE neighbours are similar to those shown by the stabilities of complexes with aminocarboxylic acids and are significantly larger than those shown by stabilities of complexes with carboxylic acids. Hence, sorption of Y and REEs onto hydrous Fe-Mn oxides may contribute significantly to the positive Y SN and Gd SN anomalies in seawater.

Geochimica et Cosmochimica Acta↗

Machine learning for predicting soil classes in three semi-arid landscapes

Mapping the spatial distribution of soil taxonomic classes is important for informing soil use and management decisions. Digital soil mapping (DSM) can quantitatively predict the spatial distribution of soil taxonomic classes. Key components of DSM are the method and the set of environmental covariates used to predict soil classes. Machine learning is a general term for a broad set of statistical modeling techniques. Many different machine learning models have been applied in the literature and there are different approaches for selecting covariates for DSM. However, there is little guidance as to which, if any, machine learning model and covariate set might be optimal for predicting soil classes across different landscapes. Our objective was to compare multiple machine learning models and covariate sets for predicting soil taxonomic classes at three geographically distinct areas in the semi-arid western United States of America (southern New Mexico, southwestern Utah, and northeastern Wyoming). All three areas were the focus of digital soil mapping studies. Sampling sites at each study area were selected using conditioned Latin hypercube sampling (cLHS). We compared models that had been used in other DSM studies, including clustering algorithms, discriminant analysis, multinomial logistic regression, neural networks, tree based methods, and support vector machine classifiers. Tested machine learning models were divided into three groups based on model complexity: simple, moderate, and complex. We also compared environmental covariates derived from digital elevation models and Landsat imagery that were divided into three different sets: 1) covariates selected a priori by soil scientists familiar with each area and used as input into cLHS, 2) the covariates in set 1 plus 113 additional covariates, and 3) covariates selected using recursive feature elimination. Overall, complex models were consistently more accurate than simple or moderately complex models. Random forests (RF) using covariates selected via recursive feature elimination was consistently the most accurate, or was among the most accurate, classifiers between study areas and between covariate sets within each study area. We recommend that for soil taxonomic class prediction, complex models and covariates selected by recursive feature elimination be used. Overall classification accuracy in each study area was largely dependent upon the number of soil taxonomic classes and the frequency distribution of pedon observations between taxonomic classes. Individual subgroup class accuracy was generally dependent upon the number of soil pedon observations in each taxonomic class. The number of soil classes is related to the inherent variability of a given area. The imbalance of soil pedon observations between classes is likely related to cLHS. Imbalanced frequency distributions of soil pedon observations between classes must be addressed to improve model accuracy. Solutions include increasing the number of soil pedon observations in classes with few observations or decreasing the number of classes. Spatial predictions using the most accurate models generally agree with expected soil–landscape relationships. Spatial prediction uncertainty was lowest in areas of relatively low relief for each study area.

New Mexico, Utah, Wyoming↗

Coproduction and modeling spatial contact networks prevent bias about infectious hematopoietic necrosis virus transmission for Snake River Basin salmonids

Much remains unknown about variation in pathogen transmission across the geographic range of a free-ranging fish or animal species and about the influence of movement (associated with husbandry practices or animal behavior) on pathogen transmission. Salmonid hatcheries are an ideal system in which to study these processes. Salmonid hatcheries are managed for endangered species recovery, supplementation of threatened or at-risk fish stocks, support of fisheries, and ecosystem stability. Infectious hematopoietic necrosis virus (IHNV) is a rhabdovirus of significant concern to salmon aquaculture. Landscape IHNV transmission dynamics previously had been estimated only for salmonid hatcheries in the Lower Columbia River Basin (LCRB). The objectives of this study were to estimate IHNV transmission dynamics in a unique geographic region, the Snake River Basin (SRB), and to quantitatively estimate the effect of model coproduction on inference because previous assessments of coproduction have been qualitative. In contrast to the LCRB, the SRB has hatchery complexes consisting of a main hatchery and ≥1 satellite facility. Knowledge about hatchery complexes was held by a subset of project researchers but would not have been available to project modelers without coproduction. Project modelers generated and tested multiple versions of Bayesian susceptible-exposedinfected models to realistically represent the SRB and estimate the effect of coproduction. Models estimated the frequency of transmission routes, route-specific infection probabilities, and infection probabilities for combinations of salmonid hosts and IHNV lineages. Model results indicated that in the SRB, avoiding exposure to IHNV-positive adult salmonids is the most important action to prevent juvenile infections. Migrating adult salmonids exposed juvenile cohort-sites most frequently, and the infection probability was greatest following exposure to migrating adults. Without coproduction, the frequency of exposure by migrating adults would have been overestimated by 70 cohort-sites, and the infection probability following exposure to migrating adults would have been underestimated by∼0.09. The coproduced model had less uncertainty in the infection probability if no transmission route could be identified (Bayesian credible interval (BCI) width = 0.12) compared to the model without coproduction (BCI width = 0.34). Evidence for virus lineage MD specialization on steelhead and rainbow trout (both Oncorhynchus mykiss) was apparent without model coproduction. In the SRB, we found a greater probability of virus lineage UC infection in Chinook salmon (Oncorhynchus tshawytscha) compared to in O. mykiss, whereas in the LCRB, UC more clearly exhibited a generalist approach. Coproduction influenced estimates that depended on transmission routes, which operated differently at main hatcheries and satellite sites within hatchery complexes. Hatchery complexes are found outside of the SRB and are not specific to salmonid hatcheries alone. There is great potential for coproduction and modeling spatial contact networks to advance understanding about infectious disease transmission in complex production systems and surrounding free-ranging animal populations.

Idaho, Oregon, Washington↗

Mapping the 3-D extent of the Northern Lobe of the Bushveld layered mafic intrusion from geophysical data

Geophysical models image the 3D geometry of the mafic portion of the Bushveld Complex north of the Thabazimbi-Murchison Lineament (TML), critical for understanding the origin of the world's largest layered mafic intrusion and platinum group element deposits. The combination of the gravity and magnetic data with recent seismic, MT, borehole and rock property measurements powerfully constrains the models. The intrusion north of the TML is generally shallowly buried (generally <1500 m) with a modeled area of &sim;160 km &times; &sim;125 km. The modeled thicknesses are not well constrained but vary from &sim;<1000 to >12,000 m, averaging &sim;4000 m. A feeder, suggested by a large modeled thickness (>10,000 m) and funnel shape, for Lower Zone magmas could have originated near the intersection of NS and NE trending TML faults under Mokopane. The TML has been thought to be the feeder zone for the entire Bushveld Complex but the identification of local feeders and/or dikes in the TML in the models is complicated by uncertainties on the syn- and post-Bushveld deformation history. However, modeled moderately thick high density material near the intersection of faults within the central and western TML may represent feeders for parts of the Bushveld Complex if deformation was minimal. The correspondence of flat, high resistivity and density regions reflect the sill-like geometry of the Bushveld Complex without evidence for feeders north of Mokopane. Magnetotelluric models indicate that the Transvaal sedimentary basin underlies much of the Bushveld Complex north of the TML, further than previously thought and important because the degree of reaction and assimilation of the Transvaal rocks with the mafic magmas resulted in a variety of mineralization zones.

Precambrian Research↗

Towards globally customizable ecosystem service models

Scientists, stakeholders and decision makers face trade-offs between adopting simple or complex approaches when modeling ecosystem services (ES). Complex approaches may be time- and data-intensive, making them more challenging to implement and difficult to scale, but can produce more accurate and locally specific results. In contrast, simple approaches allow for faster assessments but may sacrifice accuracy and credibility. The ARtificial Intelligence for Ecosystem Services (ARIES) modeling platform has endeavored to provide a spectrum of simple to complex ES models that are readily accessible to a broad range of users. In this paper, we describe a series of five “Tier 1” ES models that users can run anywhere in the world with no user input, while offering the option to easily customize models with context-specific data and parameters. This approach enables rapid ES quantification, as models are automatically adapted to the application context. We provide examples of customized ES assessments at three locations on different continents and demonstrate the use of ARIES' spatial multi-criteria analysis module, which enables spatial prioritization of ES for different beneficiary groups. The models described here use publicly available global- and continental-scale data as defaults. Advanced users can modify data input requirements, model parameters or entire model structures to capitalize on high-resolution data and context-specific model formulations. Data and methods contributed by the research community become part of a growing knowledge base, enabling faster and better ES assessment for users worldwide. By engaging with the ES modeling community to further develop and customize these models based on user needs, spatiotemporal contexts, and scale(s) of analysis, we aim to cover the full arc from simple to complex assessments, minimizing the additional cost to the user when increased complexity and accuracy are needed.

Science of the Total Environment↗

Geologic structures related to New Madrid earthquakes near Memphis, Tennessee, based on gravity and magnetic interpretations

New inversions of gravity and magnetic data in the region north of memphis. Tennessee, and south of latitude 36?? define boundaries of regional structures and igneous complexes in the upper crust. Microseismicity patterns near interpreted boundaries suggest that igneous complexes influence the locations of microseismicity. A weak seismicity cluster occurs near one intrusion (Covington pluton), at the intersection of the southwest margin of the Missouri batholith and the southeast margin of the Reelfoot rift. A narrow seismicity trend along the Reelfoot rift axis becomes diffuse near a second intrusion (Osceola intrusive complex) and changes direction to an area along the northwest flank of the intrusion. The axial seismicity trend also contains a tight cluster of earthquakes located just outside the Osceola intrusive complex. The mechanical explanation of the two seismicity patterns is uncertain, but the first cluster may be caused by stress concentration due to the high elastic stiffness and strength of the Covington intrusion. The spatially changing seismicity pattern near the Osceola complex may be caused by the preceding factors plus interaction with faulting along the rift axis. The axial seismicity strand itself is one of several connected and interacting active strands that may produce stress concentrations at strand ends and junctions. The microseismicity clusters at the peripheries of the two intrusions lead us to conclude that these stress concentrations or stressed volumes may be locations of future moderate to large earthquakes near Memphis. Published by Elsevier Science B.V.

Engineering Geology↗

Constants for mercury binding by organic matter isolates from the Florida Everglades

Dissolved organic matter (DOM) has been implicated as an important complexing agent for Hg that can affect its mobility and bioavailability in aquatic ecosystems. However, binding constants for natural Hg-DOM complexes are not well known. We employed a competitive ligand approach to estimate conditional stability constants for Hg complexes with DOM isolates collected from Florida Everglades surface waters. The isolates examined were the hydrophobic fraction of DOM from a eutrophic, sulfidic site (F1-HPoA) and the hydrophilic fraction from an oligotrophic, low-sulfide site (2BS-HPiA). Our experimental determinations utilized overall octanol-water partitioning coefficients (D ow ) for 203 Hg at 0.01 M chloride and across pH and DOM concentration gradients. Use of this radioisotope allowed rapid determinations of Hg concentrations in both water and octanol phases without problems of matrix interference. Conditional stability constants (I = 0.06, 23°C) were log K′ = 11.8 for F1-HPoA and log K′ = 10.6 for 2BS-HPiA. These are similar to previously published stability constants for Hg binding to low-molecular-weight thiols. Further, F1-HPoA showed a pH-dependent decline in D ow that was consistent with models of Hg complexation with thiol groups as the dominant Hg binding sites in DOM. These experiments demonstrate that the DOM isolates are stronger ligands for Hg than chloride ion or ethylenediamine-tetraacetic acid. Speciation calculations indicate that at the DOM concentrations frequently measured in Everglades, 20 to 40 μM, significant complexation of Hg by DOM would be expected in aerobic (sulfide-free) surface waters.

Florida↗