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,063 records · Page 59Linked to original sources

Simulated juvenile salmon growth and phenology respond to altered thermal regimes and stream network shape

It is generally accepted that climate change will stress coldwater species such as Pacific salmon. However, it is unclear what aspect of altered thermal regimes (e.g., warmer winters, springs, summers, or increased variability) will have the greatest effect, and what role the spatial properties of river networks play. Thermally diverse habitats may afford protection from climate change by providing opportunities for aquatic organisms to find and use habitats with optimal conditions for growth. We hypothesized that climate‐altered thermal regimes will change growth and timing of life history events such as emergence or migration but that changes will be moderated in topologically complex stream networks where opportunities to thermoregulate are more readily available to mobile animals. Because climate change effects on populations are spatially variable and contingent upon physiological optima, assessments of risk must take a spatially explicit approach. We developed a spatially structured individual‐based model for Chinook Salmon ( Oncorhynchus tshawytscha ) in which movement decisions and growth were governed by water temperature and conspecific density. We evaluated growth and phenology (timing of egg emergence and smolting) under a variety of thermal regimes (each having a different minimum, rate of warming, maximum, and variability) and in three network shapes of increasing spatial complexity. Across networks, fish generally grew faster and were capable of smolting earlier in warmer scenarios where water temperatures experienced by fish were closer to optimal; however, growth decreased for some fish. We found that salmon size and smolt date responded more strongly to warmer springs and summers than to warmer winters or increased variability. Fish in the least complex network grew faster and were ready to smolt earlier than fish in the more spatially complex network shapes in the contemporary thermal regime; patterns were similar but less clear in warmer thermal regimes. Our results demonstrate that network topology may influence how fish respond to thermal landscapes, and this information will be useful for incorporating a spatiotemporal context into conservation decisions that promote long‐term viability of salmon in a changing climate.

Ecosphere↗

Clarifying the trophic state concept to advance macroscale freshwater science and management

For over a century, ecologists have used the concept of trophic state (TS) to characterize an aquatic ecosystem's biological productivity. However, multiple TS classification schemes, each relying on a variety of measurable parameters as proxies for productivity, have emerged to meet use-specific needs. Frequently, chlorophyll a, phosphorus, and Secchi depth are used to classify TS based on autotrophic production, whereas phosphorus, dissolved organic carbon, and true color are used to classify TS based on both autotrophic and heterotrophic production. Both classification approaches aim to characterize an ecosystem's function broadly, but with varying degrees of autotrophic and heterotrophic processes considered in those characterizations. Moreover, differing classification schemes can create inconsistent interpretations of ecosystem integrity. For example, the US Clean Water Act focuses exclusively on algal threats to water quality, framed in terms of eutrophication in response to nutrient loading. This usage lacks information about non-algal threats to water quality, such as dystrophication in response to dissolved organic carbon loading. Consequently, the TS classification schemes used to identify eutrophication and dystrophication may refer to ecosystems similarly (e.g., oligotrophic and eutrophic), yet these categories are derived from different proxies. These inconsistencies in TS classification schemes may be compounded when interdisciplinary projects employ varied TS frameworks. Even with these shortcomings, TS can still be used to distill information on complex aquatic ecosystem function into a set of generalizable expectations. The usefulness of distilling complex information into a TS index is substantial such that usage inconsistencies should be explicitly addressed and resolved. To emphasize the consequences of diverging TS classification schemes, we present three case studies for which an improved understanding of the TS concept advances freshwater research, management efforts, and interdisciplinary collaboration. To increase clarity in TS, the aquatic sciences could benefit from including information about the proxy variables, ecosystem type, as well as the spatiotemporal domains used to classify TS. As the field of aquatic sciences expands and climatic irregularity increases, we highlight the importance of re-evaluating fundamental concepts, such as TS, to ensure their compatibility with evolving science.

Ecosphere↗

Interaction strength and harvest intensity mediate predator–prey dynamics on coral reefs

Understanding predator–prey relationships is fundamental to our knowledge of the stability and resilience of ecological systems. These dynamics are shaped by both ecological factors, like interaction strength, and anthropogenic factors, like harvest intensity, which can have large-scale implications for community structure. However, few studies have focused on the combined impact of these effects and their contribution to phenomena like prey release within two-species frameworks. In this study, we investigate the interactive impact of interaction strength and harvest pressure on two trophic levels in a predator–prey system using a mathematical modeling approach. Our results reveal that interaction strength plays a crucial role in shaping population dynamics, with high interaction strength leading to a predator-dominated system and low interaction strength enabling coexistence between species. The addition of predator harvest into the system reveals complex and counterintuitive behavior not seen in unharvested systems, likely due to the destabilizing impacts of harvest at some interaction strengths. Specifically, the inclusion of harvest on the predator can induce a range of behaviors, such as prey release and predator decline, that alter the equilibrium abundance of both predator and prey populations. Interestingly, predator–prey systems with intermediate to high interaction strengths achieve maximum total abundance with low harvest levels rather than in scenarios with no harvest pressure, as prey populations benefit greatly from reduced predation mortality associated with predator harvest. We gain insights into the complex interplay between predator–prey interactions and human activities in shaping community composition and abundances across trophic levels. This study provides potential mechanisms that may explain the observed variation in numerical prey release in trophically complex systems in which predators and prey are both extracted, like coral reef fisheries. Results highlight the need for resource management to consider the wide range of factors that shape ecosystem dynamics to develop effective strategies that safeguard the long-term health of complex ecosystems and the human communities that they support.

Ecosphere↗

Relations between rainfall–runoff-induced erosion and aeolian deposition at archaeological sites in a semi-arid dam-controlled river corridor

Process dynamics in fluvial-based dryland environments are highly complex with fluvial, aeolian, and alluvial processes all contributing to landscape change. When anthropogenic activities such as dam-building affect fluvial processes, the complexity in local response can be further increased by flood- and sediment-limiting flows. Understanding these complexities is key to predicting landscape behavior in drylands and has important scientific and management implications, including for studies related to paleoclimatology, landscape ecology evolution, and archaeological site context and preservation. Here we use multi-temporal LiDAR surveys, local weather data, and geomorphological observations to identify trends in site change throughout the 446-km-long semi-arid Colorado River corridor in Grand Canyon, Arizona, USA, where archaeological site degradation related to the effects of upstream dam operation is a concern. Using several site case studies, we show the range of landscape responses that might be expected from concomitant occurrence of dam-controlled fluvial sand bar deposition, aeolian sand transport, and rainfall-induced erosion. Empirical rainfall-erosion threshold analyses coupled with a numerical rainfall–runoff–soil erosion model indicate that infiltration-excess overland flow and gullying govern large-scale (centimeter- to decimeter-scale) landscape changes, but that aeolian deposition can in some cases mitigate gully erosion. Whereas threshold analyses identify the normalized rainfall intensity (defined as the ratio of rainfall intensity to hydraulic conductivity) as the primary factor governing hydrologic-driven erosion, assessment of false positives and false negatives in the dataset highlight topographic slope as the next most important parameter governing site response. Analysis of 4+ years of high resolution (four-minute) weather data and 75+ years of low resolution (daily) climate records indicates that dryland erosion is dependent on short-term, storm-driven rainfall intensity rather than cumulative rainfall, and that erosion can occur outside of wet seasons and even wet years. These results can apply to other similar semi-arid landscapes where process complexity may not be fully understood.

Arizona↗

The Ebro Deep-Sea Fan system

The Ebro Fan System consists of en echelon channel-levee complexes, 50??20 km in area and 200-m thick. A few strong reflectors in a generally transparent seismic facies identify the sand-rich channel floors and levee crests. Numerous continuous acoustic reflectors characterize overbank turbidites and hemipelagites that blanket abandoned channel-levee complexes. The interlobe areas between channel complexes fill with homogeneous mud and sand from mass flow and overbank deposition; these exhibit a transparent seismic character. The steep continental rise and sediment 'drainage' of Valencia Trough at the end of the channel-levee complexes prevent the development of distributary channels and midfan lobe deposits. ?? 1984 Springer-Verlag New York Inc.

Geo-Marine Letters↗

Acoustical deterrence of Silver Carp ( Hypophthalmichthys molitrix )

The invasive Silver Carp ( Hypophthalmichthys molitrix ) dominate large regions of the Mississippi River drainage and continue to expand their range northward threatening the Laurentian Great Lakes. This study found that complex broadband sound (0–10 kHz) is effective in altering the behavior of Silver Carp with implications for deterrent barriers or potential control measures (e.g., herding fish into nets). The phonotaxic response of Silver Carp was investigated using controlled experiments in outdoor concrete ponds (10 × 4.9 × 1.2 m). Pure tones (500–2000 Hz) and complex sound (underwater field recordings of outboard motors) were broadcast using underwater speakers. Silver Carp always reacted to the complex sounds by exhibiting negative phonotaxis to the sound source and by alternating speaker location, Silver Carp could be directed consistently, up to 37 consecutive times, to opposite ends of the large outdoor pond. However, fish habituated quickly to pure tones, reacting to only approximately 5 % of these presentations and never showed more than two consecutive responses. Previous studies have demonstrated the success of sound barriers in preventing Silver Carp movement using pure tones and this research suggests that a complex sound stimulus would be an even more effective deterrent.

Biological Invasions↗

Rearing environment influences boldness and prey acquisition behavior, and brain and lens development of bull trout

Animals reared in barren captive environments exhibit different developmental trajectories and behaviors than wild counterparts. Hence, the captive phenotypes may influence the success of reintroduction and recovery programs for threatened and endangered species. We collected wild bull trout embryos from the Metolius River Basin, Oregon and reared them in differing environments to better understand how captivity affects the bull trout Salvelinus confluentus phenotype. We compared the boldness and prey acquisition behaviors and development of the brain and eye lens of bull trout reared in conventional barren and more structurally complex captive environments with that of wild fish. Wild fish and captive reared fish from complex habitats exhibited a greater level of boldness and prey acquisition ability, than fish reared in conventional captive environments. In addition, the eye lens of conventionally reared bull trout was larger than complex reared captive fish or same age wild fish. Interestingly, we detected wild fish had a smaller relative cerebellum than either captive reared treatment. Our results suggest that rearing fish in more complex captive environments can create a more wild-like phenotype than conventional rearing practices. A better understanding of the effects of captivity on the development and behavior of bull trout can inform rearing and reintroduction programs though prediction of the performance of released individuals.

Environmental Biology of Fishes↗

Salt marsh habitats and diamondback terrapins in a rapidly changing climate: A review

The impacts associated with global climate change (e.g., sea-level rise, tropical storms, and warming temperatures) are expected to alter predator–prey interactions, foundation species, and plant community structure in coastal ecosystems. While the complex dynamics of these habitats have been examined under future climate predictions, few ecosystem models incorporate influences from fauna, such as the diamondback terrapin, the only estuarine turtle native to North America. This review examines the impacts of climate change on diamondback terrapins ( Malaclemys terrapin ) and the role that terrapins play as higher trophic level predators and keystone species in driving the dynamics of these ecosystems. We also review the potential implications of changes to terrapin populations on coastal ecosystems as a conservation challenge and suggest solutions to advance our understanding of those complex systems. Because of their role as a keystone and area-sensitive species that helps maintain healthy coastal habitats by foraging on herbivorous periwinkle snails, alterations to terrapin life history from climate change are expected, which could have significant impacts to the conservation of coastal habitats. Life history alterations could occur due to individual stressors, such as warming temperatures altering terrapin sex ratios. However, because of the complexity of these coastal systems, these stressors could also act additively or synergistically. Inclusion of faunal taxa such as the diamondback terrapin in modeling efforts examining climate change impacts to coastal ecosystems would better represent the complexity of these habitats thereby providing a more comprehensive evaluation of the entire ecosystem, resulting in more effective conservation strategies.

Estuaries and Coasts↗

Effect of surface area and chemisorbed oxygen on the SO2 adsorption capacity of activated char

The objective of this study was to determine whether activated char produced from Illinois coal could be used effectively to remove sulfur dioxide from coal combustion flue gas. Chars were prepared from a high-volatile Illinois bituminous coal under a wide range of pyrolysis and activation conditions. A novel char preparation technique was developed to prepare chars with SO2 adsorption capacities significantly greater than that of a commercial activated carbon. In general, there was no correlation between SO2 adsorption capacity and surface area. Temperature-programmed desorption (TPD) was used to determine the nature and extent of carbon-oxygen (C-O) complexes formed on the char surface. TPD data revealed that SO2 adsorption was inversely proportional to the amount of C-O complex. The formation of a stable C-O complex during char preparation may have served only to occupy carbon sites that were otherwise reactive towards SO2 adsorption. A fleeting C(O) complex formed during SO2 adsorption is postulated to be the reaction intermediate necessary for conversion of SO2 to H2SO4. Copyright ?? 1996 Elsevier Science Ltd.

Fuel↗

Interactions between ground water and wetlands, southern shore of Lake Michigan, USA

Wetlands between, and within, dune-beach complexes along the south shore of Lake Michigan are strongly affected by ground water. The hydrogeology of the glacial drift aquifer system in a 26 km 2 area was investigated to determine the effects of ground water on the hydrology and hydrochemistry of Cowles Bog and its adjacent wetlands. The investigation showed that ground water from intermediate- and regional-scale flow systems discharges to Cowles Bog from confined aquifers that underlie the wetland. These flow systems are recharged in moraines south of the dune-beach complexes. Water from the confined aquifers discharges into the surficial aquifer mainly by upward leakage through a buried till sheet that serves as the confining layer. However, the till sheet is breached below a raised peat mound in Cowles Bog, allowing direct upward discharge from the confinef aquifer into the surficial sand, marl, and peat. The shallow ground and wetland water in the area influenced by this leakage is a calcium magnesium bicarbonate type, with low tritium concentrations consistent with mixing of older ground water and more recent precipitation. Ground water and wetland water from surrounding areas are less mineralized and have higher tritium concentrations characteristic of precipitation in the late 1970s. The results of this study suggest that wetlands in complex hydrogeologic settings may be influenced by multiple ground-water flow systems that are affected by geomorphic features, stratigraphic discontinuities, and changes in sediment types. Discharge and recharge zones may both occur in the same wetland. Multidisciplinary studies incorporating hydrological, hydrochemical, geophysical, and sedimentological data are necessary to identify such complexities in wetland hydrology.

Journal of Hydrology↗

Role of heat and detachment in continental extension as viewed from the eastern basin and range province in Arizona

The Bill Williams River area of west-central Arizona includes not only the Rawhide-Buckskin metamorphic core complex, which is part of the lower Colorado River highly extended terrane (HET), but also the boundary between the extended terranes of the Basin and Range Province and the less deformed Arizona Transition Zone/Colorado Plateau. This provides important constraints on models that address the mechanisms for the mid- to late Tertiary deformation. Three phases of extension are present. The oldest is the extension associated with core-complex tectonism, which characteristically shows a lower plate composed of lineated mylonitic gneiss overlain by a detachment fault that is regionally nearly horizontal but undulates at the local scale. The fault in turn is overlain by an upper plate that includes Precambrian basement rocks, recrystallized Paleozoic sedimentary rocks, Mesozoic(?) metasedimentary and metavolcanic rocks of greenschist facies, and unaltered to hydrothermally altered syntectonic sedimentary and volcanic rocks of Miocene age. The upper plate is cut by closely spaced faults of modest structural relief that strike northwest and strongly rotate intervening blocks to face southwest. Most of these faults do not penetrate below the detachment fault. Fault spacing increases, and rotation decreases, to the northeast, away from the trace of the detachment. The second phase consists of “classic” Basin-Range high-angle normal faults that strike about north and have wide spacing, high structural relief, and modest rotation of blocks. These faults have no consistent direction of displacement and so produced horst and graben that form the ranges and basins visible today. This phase is locally superposed on Phase I, and also extends in more subdued form into the Transition Zone/Colorado Plateau. The third phase consists of tectonic quiescence and is present everywhere except parts of the Transition Zone that are still active seismically. The first phase occurred in the early and middle Miocene and was accompanied by deposition of syntectonic fluviolacustrine rocks (Suite I); the second (middle to late Miocene) was marked by interior-basin deposits (Suite II); the third (latest Miocene through Quaternary) is characterized by deposits related to through-flowing drainage. The phases grade into each other and thus are likely to be genetically related. Tectonic models must take into account not only the geographic distribution of deformation at any one time but also the time-dependent succession of deformation at any one place. A model proposed in this paper attempts to do this. The model is thermotectonic. A heating event in the lower crust, (basaltic intrusion, asthenospheric upwelling) combined with stretching, causes a sharp thermal front to rise within the crust. Embedded within the front is an “isotherm” that marks the brittle-ductile transition. As the front rises, it leaves behind a trail of shear zones, each marking a locus of preferred failure defined by mechanical or physical properties, or combinations thereof. The highest shear zone, now preserved in fossil form as the “detachment”, occurs where the front impinges on the meteoric groundwater, a few km below the topographic surface. The water steepens the thermal gradient at the front, which it stabilizes. A convective hydrothermal circulation system is established, causing alteration and mineralization above the ductile-brittle transition, as well as pore overpressure that results in hydrofracturing (producing monolithologic breccias) and the sliding of gravity-glide sheets. During these events, extension is taking place by brittle failure in the upper plate and ductile deformation below the detachment. Simultaneously, the hottest areas (core complexes) are updomed, promoting drainage reversals and the sliding of breccias and glide sheets. All this occurred only in the hottest areas or “blisters”, now marked by the core complexes. Distal areas showed less or no deformation at the surface. With time and the waning of the thermal event, the thermal front, and thus the brittle-ductile transition, smoothed out and sank, again leaving a trail of shear zones. Phase 1 deformation ceased and was replaced by Phase 2 deformation that occurred over a much wider area. Eventually, the front sank so deep that surface deformation ceased. This illustrates how the style of deformation at the surface may be a measure of the depth to the brittle-ductile transition. According to the thermotectonic model, extensional strain does not need to be constant along the detachment, in contrast to models involving simple shear through crustal-scale normal faults. On the contrary, one would expect strain to vary geographically as a function of maximum temperature attained, because of the well known relation between temperature and lithospheric strength. The thermotectonic model is also in good accord with geophysical characteristics of the Basin and Range Province, which suggests that extension was accompanied by intrusion of basalt into the lower crust, with consequent heating and anatexis. Many studies in the U.S. and elsewhere support the model by showing that continental extension commonly is accompanied by near-surface temperatures corresponding to the brittle-ductile transition, by steep thermal gradients, and by hydrothemal convective systems. A possible driving mechanism from the thermotectonic processes described by the model is the rise of asthenospheric domes or welts, which thin the lithosphere by subcrustal transfer while heating and stretching it. An asthenospheric welt that migrates northeastward while dying out might explain the encroachment of relatively subdued extension onto the Colorado Plateau, as well as the juxtaposition of compressive stress on the plateau with extensional stress in the adjacent Transition Zone and Basin and Range Province.

Arizona↗

Precambrian ophiolites of Arabia: Geologic settings, U Pb geochronology, Pb-isotope characteristics, and implications for continental accretion

Disrupted ophiolites occur in linear belts up to 900 km long between microplates that collided during the late Proterozoic to form the Arabian Shield. UPb zircon ages and Pb-isotope data from these ophiolitic rocks help constrain the history of accretion of the Arabian Shield and thereby contribute to the definition of its microplates and terranes. Terranes of the central and western Arabian Shield are generally thought to represent intraoceanic island arcs that range in age from about 900 to 640 Ma; however, a region of the eastern Arabian Shield contains rocks of Early Proterozoic age and may represent an exotic continental fragment entrained between the arc complexes. Ophiolites of the Yanbu suture (northwestern shield), dated by UPb (zircon) and SmNd (mineral isochron) methods, yield model ages of 740–780 Ma. These are among the oldest well-dated rocks in the northwestern Arabian Shield. Ages from the Jabal al Wask complex overlap with ages of adjacent arc rocks. This overlap in age supports geologic and geochemical evidence that the Wask complex represents a fragment of back-arc oceanic lithosphere formed during arc magmatism. Older ages of about 780 Ma for gabbro from the Jabal Ess ophiolite suggest that the ophiolite is either a fragment of fore-arc oceanic crust or oceanic basement on which an arc was built. Gabbro samples from ophiolites of the Bir Umq suture (west-central Arabian Shield) yield zircons with ages of 820–870 Ma and $ ̆ 1250 Ma "> 1250Ma . The 820–870 Ma dates overlap with ages of the oldest nearby arc rocks; this favors an intra-arc or near-arc paleotectonic setting. The older zircons suggest that middle or early Proterozoic crustal material, possibly derived from the Mozambique belt of Africa, was present during back- or intra-arc magmatism. Plagiogranite from the Bir Tuluhah ophiolitic complex at the northern end of the 900 km-long Nabitah mobile belt was dated by the zircon U Pb method at ∼ 830 Ma. This date is in the range of the oldest dated arc rocks along the northern and central parts of the Nabitah suture, but is ∼ 100 Ma older than the oldest arc plutons (tonalites) associated with the southern part of the belt. These age relations suggest that the northern part of the Nabitah belt contains an extension of the Bir Umq suture that was transposed parallel to the Nabitah trend during collision of the arc terranes of the northwest Arabian Shield with the Afif plate to the east. Feldspar lead-isotope data from the ophiolites are of three types: (1) lead from the ophiolitic rocks and arc tonalites of the northwestern Shield and ophiolitic rocks of the Nabitah suture is similar to lead in modern mid-ocean ridge basalt, (2) anomalous radiogenic data from the Thurwah ophiolite are from rocks that contain zircons from pre-late Proterozoic continental crust, and (3) feldspar from the Urd ophiolite shows retarded uranogenic lead growth and is related either to an anomalous and perhaps primitive oceanic mantle source, or in an unknown manner to ancient continental mantle or lower crust of the eastern Arabian Shield.

Precambrian Research↗

Extended history of a 3.5 Ga trondhjemitic gneiss, Wyoming Province, USA: Evidence from U-Pb systematics in zircon

The Beartooth-Bighorn magmatic zone (BBMZ) and the Montana metasedimentary province (MMP) are two major subprovinces of the Archean Wyoming province. In the northwestern Beartooth Mountains, these subprovinces are separated by a structurally, lithologically and metamorphically complex assemblage of lithotectonic units that include: (1) a strongly deformed complex of trondhjemitic gneiss and interlayered amphibolites; and (2) an amphibolite facies mafic unit that occurs in a nappe that structurally overlies the gneiss complex. Zircons from a trondhjemitic blastomylonite in the gneiss complex yield concordant U-Pb ages of 3.5 Ga, establishing it as the oldest rock yet documented in the Wyoming province. Two younger events are also recorded by zircons in this rock: (1) an apparently protracted period of high-grade metamorphism and/or intrusion of additional magmas at ∼ 3.25 Ga; and (2) growth of hydrothermal zircon at ∼ 2.55 Ga, apparently associated with ductile deformation that immediately preceded structural emplacement of the gneiss. Although this latter event appears confined to areas along the BBMZ-MMP boundary, evidence of ∼ 3.25 Ga igneous activity is found in the overlying amphibolite (3.24 Ga) and throughout the MMP. These data suggest that this boundary first developed as a major intracratonic zone of displacement at or before 3.25 Ga. The limited occurrences of 2.8 Ga magmatic activity in the MMP suggest that it had a controlling influence on late Archean magmatism as well.

Precambrian Research↗

Enhancing the Simplified Surface Energy Balance (SSEB) Approach for Estimating Landscape ET: Validation with the METRIC model

Evapotranspiration (ET) can be derived from satellite data using surface energy balance principles. METRIC (Mapping EvapoTranspiration at high Resolution with Internalized Calibration) is one of the most widely used models available in the literature to estimate ET from satellite imagery. The Simplified Surface Energy Balance (SSEB) model is much easier and less expensive to implement. The main purpose of this research was to present an enhanced version of the Simplified Surface Energy Balance (SSEB) model and to evaluate its performance using the established METRIC model. In this study, SSEB and METRIC ET fractions were compared using 7 Landsat images acquired for south central Idaho during the 2003 growing season. The enhanced SSEB model compared well with the METRIC model output exhibiting an r2 improvement from 0.83 to 0.90 in less complex topography (elevation less than 2000 m) and with an improvement of r2 from 0.27 to 0.38 in more complex (mountain) areas with elevation greater than 2000 m. Independent evaluation showed that both models exhibited higher variation in complex topographic regions, although more with SSEB than with METRIC. The higher ET fraction variation in the complex mountainous regions highlighted the difficulty of capturing the radiation and heat transfer physics on steep slopes having variable aspect with the simple index model, and the need to conduct more research. However, the temporal consistency of the results suggests that the SSEB model can be used on a wide range of elevation (more successfully up 2000 m) to detect anomalies in space and time for water resources management and monitoring such as for drought early warning systems in data scarce regions. SSEB has a potential for operational agro-hydrologic applications to estimate ET with inputs of surface temperature, NDVI, DEM and reference ET.

Agricultural Water Management↗

Modules based on the geochemical model PHREEQC for use in scripting and programming languages

The geochemical model PHREEQC is capable of simulating a wide range of equilibrium reactions between water and minerals, ion exchangers, surface complexes, solid solutions, and gases. It also has a general kinetic formulation that allows modeling of nonequilibrium mineral dissolution and precipitation, microbial reactions, decomposition of organic compounds, and other kinetic reactions. To facilitate use of these reaction capabilities in scripting languages and other models, PHREEQC has been implemented in modules that easily interface with other software. A Microsoft COM (component object model) has been implemented, which allows PHREEQC to be used by any software that can interface with a COM server—for example, Excel®, Visual Basic®, Python, or MATLAB". PHREEQC has been converted to a C++ class, which can be included in programs written in C++. The class also has been compiled in libraries for Linux and Windows that allow PHREEQC to be called from C++, C, and Fortran. A limited set of methods implements the full reaction capabilities of PHREEQC for each module. Input methods use strings or files to define reaction calculations in exactly the same formats used by PHREEQC. Output methods provide a table of user-selected model results, such as concentrations, activities, saturation indices, and densities. The PHREEQC module can add geochemical reaction capabilities to surface-water, groundwater, and watershed transport models. It is possible to store and manipulate solution compositions and reaction information for many cells within the module. In addition, the object-oriented nature of the PHREEQC modules simplifies implementation of parallel processing for reactive-transport models. The PHREEQC COM module may be used in scripting languages to fit parameters; to plot PHREEQC results for field, laboratory, or theoretical investigations; or to develop new models that include simple or complex geochemical calculations.

Computers & Geosciences↗

Fluoride geochemistry of thermal waters in Yellowstone National Park: I. Aqueous fluoride speciation

Thermal water samples from Yellowstone National Park (YNP) have a wide range of pH (1–10), temperature, and high concentrations of fluoride (up to 50 mg/l). High fluoride concentrations are found in waters with field pH higher than 6 (except those in Crater Hills) and temperatures higher than 50 °C based on data from more than 750 water samples covering most thermal areas in YNP from 1975 to 2008. In this study, more than 140 water samples from YNP collected in 2006–2009 were analyzed for free-fluoride activity by ion-selective electrode (ISE) method as an independent check on the reliability of fluoride speciation calculations. The free to total fluoride concentration ratio ranged from <1% at low pH values to >99% at high pH. The wide range in fluoride activity can be explained by strong complexing with H + and Al 3+ under acidic conditions and lack of complexing under basic conditions. Differences between the free-fluoride activities calculated with the WATEQ4F code and those measured by ISE were within 0.3–30% for more than 90% of samples at or above 10 −6 molar, providing corroboration for chemical speciation models for a wide range of pH and chemistry of YNP thermal waters. Calculated speciation results show that free fluoride, F − , and major complexes ( HF ( aq ) 0 "> HF(aq)0 , AlF 2+ , AlF 2 + "> AlF2+ and AlF 3 0 "> AlF30 ) account for more than 95% of total fluoride. Occasionally, some complex species like AlF 4 - "> AlF4- , FeF 2+ , FeF 2 + "> FeF2+ , MgF + and BF 2 ( OH ) 2 - "> BF2(OH)2- may comprise 1–10% when the concentrations of the appropriate components are high. According to the simulation results by PHREEQC and calculated results, the ratio of main fluoride species to total fluoride varies as a function of pH and the concentrations and ratios of F and Al.

Yellowstone National Park↗

Geochemistry of glacial sediments in the area of the Bend massive sulfide deposit, north-central Wisconsin

Geochemical exploration in northern Wisconsin has been problematic because of thick glacial overburden and complex stratigraphic record of glacial history. To assess till geochemical exploration in an area of thick glacial cover and complex stratigraphy samples of glacial materials were collected from cores from five rotasonic boreholes near a known massive sulfide deposit, the Bend deposit in north-central Wisconsin. Diamond drilling in the Bend area has defined a long, thin zone of mineralization at least partly intersected at the bedrock surface beneath 30-40 m of unconsolidated glacial sediments. The bedrock surface has remnant regolith and saprolite resulting from pre-Pleistocene weathering. Massive sulfide and mineralized rock collected from diamond drill core from the deposit contain high (10s to 10,000s ppm) concentrations of Ag, As, Au, Bi, Cu, Hg, Se, Te, and Tl. Geochemical properties of the glacial stratigraphic units helped clarify the sequence and source areas of several glacial ice advances preserved in the section. At least two till sheets are recognized. Over the zone of mineralization, saprolite and preglacial alluvial and lacustrine samples are preserved on the bedrock surface in a paleoriver valley. The overlying till sheet is a gray, silty carbonate till with a source hundreds of kilometers to the northwest of the study area. This gray till is overlain by red, sandy till with a source to the north in Proterozoic rocks of the Lake Superior area. The complex glacial stratigraphy confounds down-ice geochemical till exploration. The presence of remnant saprolite, preglacial sediment, and far-traveled carbonate till minimized glacial erosion of mineralized material. As a result, little evidence of down-ice glacial dispersion of lithologic or mineralogic indicators of Bend massive sulfide mineralization was found in the samples from the rotasonic cores. This study points out the importance of determining glacial stratigraphy and history, and identifying favorable lithologies required for geochemical exploration. Drift prospecting in Wisconsin and other areas near the outer limits of the Pleistocene ice sheets may not be unsuccessful, in part, because of complex stratigraphic sequences of multiple glaciations where deposition dominates over erosion. ?? 2004 Elsevier B.V. All rights reserved.

Journal of Geochemical Exploration↗

Forecasting water levels using the ConvLSTM algorithm in the Everglades, USA

Forecasting water levels in complex ecosystems like wetlands can support effective water resource management, ecological conservation, and understanding surface and groundwater hydrology. Predictive models can be used to simulate the complex interactions among natural processes, hydrometeorological factors, and human activities. The Greater Everglades in the USA is a well-known example of an ecosystem where complexity has motivated adoption of machine learning algorithms in water level prediction studies. This paper aims to contribute to extending existing machine learning algorithms by integrating spatiotemporal data with deep-learning algorithms in the forecasting process. In this study, a deep-learning model is developed to predict water levels on a regional scale, covering a large area of approximately 9,138 square kilometers in the Everglades ecosystem. This model has the architecture of Convolutional Long Short-Term Memory which can deal with spatiotemporal data by capturing both spatial and temporal dependencies in the training data. The forecasting capabilities of this model (referred to as the global model) are assessed by comparing the global model to two Artificial Neural Networks developed at two different gaging stations, referred to here as local models. One local model is developed at a gaging station directly influenced by nearby water control structures, whereas the other is developed at a gaging station located farther away from these structures. By leveraging data from the Everglades Depth Estimation Network spanning from January 2002 to May 2023, the global and local models were trained to forecast water levels with a two-day lead time. Our findings suggest that both the global and local models perform with approximately the same level of accuracy, with Mean Absolute Relative Error values ranging from 0.38% to 1.4% at the selected stations. The developed global model has demonstrated strong potential as a standalone forecasting tool for the entire study area in the Everglades and could eliminate the need for developing multiple local models. This finding also highlights how machine learning can capture complex spatial and temporal relationships to generate accurate water level predictions on a regional scale.

Florida↗