Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Solutions Journal”

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 541 records · Page 30Linked to original sources
USGS2015

Peclet number as affected by molecular diffusion controls transient anomalous transport in alluvial aquifer-aquitard complexes

This study evaluates the role of the Peclet number as affected by molecular diffusion in transient anomalous transport, which is one of the major knowledge gaps in anomalous transport, by combining Monte Carlo simulations and stochastic model analysis. Two alluvial settings containing either short- or long-connected hydrofacies are generated and used as media for flow and transport modeling. Numerical experiments show that 1) the Peclet number affects both the duration of the power-law segment of tracer breakthrough curves (BTCs) and the transition rate from anomalous to Fickian transport by determining the solute residence time for a given low-permeability layer, 2) mechanical dispersion has a limited contribution to the anomalous characteristics of late-time transport as compared to molecular diffusion due to an almost negligible velocity in floodplain deposits, and 3) the initial source dimensions only enhance the power-law tail of the BTCs at short travel distances. A tempered stable stochastic (TSS) model is then applied to analyze the modeled transport. Applications show that the time-nonlocal parameters in the TSS model relate to the Peclet number, P e . In particular, the truncation parameter in the TSS model increases nonlinearly with a decrease in P e due to the decrease of the mean residence time, and the capacity coefficient increases with an increase in molecular diffusion which is probably due to the increase in the number of immobile particles. The above numerical experiments and stochastic analysis therefore reveal that the Peclet number as affected by molecular diffusion controls transient anomalous transport in alluvial aquifer–aquitard complexes.

Journal of Contaminant Hydrology↗
USGS1996

Chemical and isotopic characteristics of fluids along the cameroon volcanic line, cameroon

Results of the chemical and isotopic analysis of the water and gases discharged from volcanic crater lakes and soda springs located along the Cameroon Volcanic Line were used to characterize and infer their genetic relationships. Variations in the solute compositions of the waters indicate the dominant influence of silicate hydrolysis. Na+ (40-95%) constitutes the major cation in the springs while Fe2+ + Mg2+ (70%) dominate in the CO2-rich lakes. The principal anion is HCO3 (>90%), except in the coastal springs where Cl predominates. Lakes Nyos and Monoun have Fe-Mg-Ca-HCO3 type signatures; the soda springs are essentially Na-HCO3 type, while all other lakes show similar ionic compositions to dilute surface waters. Dissolved gases show essentially CO2 (>90%), with small amounts of Ar and N2, while CH4 constitutes the principal component in the non-gassy lakes. Active volcanic gases are generally absent, except in the Lobe spring with detectable H2S. Stable isotope ratio evidence indicates that the bicarbonate waters are essentially of meteoric origin. CO2 (??13C = -2 to -8???) and He (3He/4He = 1 to 5.6Ra) infer a mantle contribution to the total CO2. CH4 has a biogenic source, while Ar and N2 are essentially atmospheric in origin, but mixing is quite common.

Journal of African Earth Sciences↗
USGS2013

Source mechanisms of persistent shallow earthquakes during eruptive and non-eruptive periods between 1981 and 2011 at Mount St. Helens, Washington

Shallow seismicity between 0 and 3-km depth has persisted at Mount St. Helens, Washington (MSH) during both eruptive and non-eruptive periods for at least the past thirty years. In this study we investigate the source mechanisms of shallow volcano-tectonic (VT) earthquakes at MSH by calculating high-quality hypocenter locations and fault plane solutions (FPS) for all VT events recorded during two eruptive periods (1981–1986 and 2004–2008) and two non-eruptive periods (1987–2004 and 2008–2011). FPS show a mixture of normal, reverse, and strike-slip faulting during all periods, with a sharp increase in strike-slip faulting observed in 1987–1997 and an increase in normal faulting in 1998–2004. FPS P-axis orientations show a ~ 90° rotation with respect to regional σ1 (N23°E) during 1981–1986 and 2004–2008, bimodal orientations (~ N-S and ~ E-W) during 1987–2004, and bimodal orientations at ~ N-E and ~ S-W from 2008–2011. We interpret these orientations to likely be due to pressurization accompanying the shallow intrusion and subsequent eruption of magma as domes during 1981–1986 and 2004–2008 and the buildup of pore pressure beneath a seismogenic volume (located at 0–1 km) with a smaller component due to the buildup of tectonic forces during 1987–2004 and 2008–2011.

Washington↗
USGS2000

Estimating formation properties from early-time oscillatory water levels in a pumped well

Hydrologists often attempt to estimate formation properties from aquifer tests for which only the hydraulic responses in a pumped well are available. Borehole storage, turbulent head losses, and borehole skin, however, can mask the hydraulic behavior of the formation inferred from the water level in the pumped well. Also, in highly permeable formations or in formations at significant depth below land surface, where there is a long column of water in the well casing, oscillatory water levels may arise during the onset of pumping to further mask formation responses in the pumped well. Usually borehole phenomena are confined to the early stages of pumping or recovery, and late-time hydraulic data can be used to estimate formation properties. In many instances, however, early-time hydraulic data provide valuable information about the formation, especially if there are interferences in the late-time data. A mathematical model and its Laplace transform solution that account for inertial influences and turbulent head losses during pumping is developed for the coupled response between the pumped borehole and the formation. The formation is assumed to be homogeneous, isotropic, of infinite areal extent, and uniform thickness, with leakage from an overlying aquifer, and the screened or open interval of the pumped well is assumed to fully penetrate the pumped aquifer. Other mathematical models of aquifer flow can also be coupled with the equations describing turbulent head losses and the inertial effects on the water column in the pumped well. The mathematical model developed in this paper is sufficiently general to consider both underdamped conditions for which oscillations arise, and overdamped conditions for which there are no oscillations. Through numerical inversion of the Laplace transform solution, type curves from the mathematical model are developed to estimate formation properties through comparison with the measured hydraulic response in the pumped well. The mathematical model is applied to estimate formation properties from a singlewell test conducted near Waialua, Oahu, Hawaii. At this site, both the drawdown and recovery showed oscillatory water levels in the pumped well, and a step-drawdown test showed that approximately 86% of the drawdown is attributed to turbulent head losses. Analyses at this site using late-time drawdown data were confounded by the noise present in the measured water levels due primarily to nearby irrigation wells and ocean tides. By analyzing the early-time oscillatory recovery data at the Waialua site, upper and lower bounds were placed on the transmissivity, T, storage coefficient, S, and the leakance of the confining unit, K′/B′. The upper and lower bounds on T differ by a factor of 2. Upper and lower bounds on S and K′/B′ are much larger, because drawdown stabilized relatively quickly after the onset of pumping.

Journal of Hydrology↗
USGS2025

Mitigating flood risks in urban estuaries: Tidal dynamics, shoreline hardening, nature-based solutions, and floodgates in San Francisco Bay

Hydrodynamic models are valuable tools for understanding the primary factors influencing daily and peak water levels and for guiding discussions on potential adaptation strategies for managing flood risk in coastal areas. This analysis uses the Delft3D San Francisco Bay-Delta Community Model to simulate water levels and incorporates the effects of a number of adaptation measures in the urban San Francisco Bay estuary, California. In particular, we examine the influence of shoreline hardening, nature-based solutions, and subregional floodgates on regional water levels. The result shows that under present conditions, tidal amplification is responsible for generating a wide distribution of extreme water levels across San Francisco Bay. Tidal amplification is found to decrease under sea level rise, thereby producing a relative damping effect on extremes. A comparison of different shoreline scenarios demonstrates that hard frontal shorelines result in higher tidal amplification, whereas restored (soft) shorelines lower amplification. The current shoreline configuration has both hard and soft characteristics and results in an intermediate tidal response. In some areas, wetland restoration reduces extreme water levels by as much as 20 cm, whereas hard-shoreline addition elevates them by as much as 10 cm for 1.5 m of sea level rise. Furthermore, local floodgates can significantly reduce high water levels without major adverse effects elsewhere in San Francisco Bay. These findings point toward the justification for a range of adaptive measures across political boundaries, weighing hard and soft options in addressing the mounting danger of sea level rise.

California↗
USGS2017

Optimization of human, animal, and environmental health by using the One Health approach

Emerging diseases are increasing burdens on public health, negatively affecting the world economy, causing extinction of species, and disrupting ecological integrity. One Health recognizes that human, domestic animal, and wildlife health are interconnected within ecosystem health and provides a framework for the development of multidisciplinary solutions to global health challenges. To date, most health-promoting interventions have focused largely on single-sector outcomes. For example, risk for transmission of zoonotic pathogens from bush-meat hunting is primarily focused on human hygiene and personal protection. However, bush-meat hunting is a complex issue promoting the need for holistic strategies to reduce transmission of zoonotic disease while addressing food security and wildlife conservation issues. Temporal and spatial separation of humans and wildlife, risk communication, and other preventative strategies should allow wildlife and humans to co-exist. Upstream surveillance, vaccination, and other tools to prevent pathogen spillover are also needed. Clear multi-sector outcomes should be defined, and a systems-based approach is needed to develop interventions that reduce risks and balance the needs of humans, wildlife, and the environment. The ultimate goal is long-term action to reduce forces driving emerging diseases and provide interdisciplinary scientific approaches to management of risks, thereby achieving optimal outcomes for human, animal, and environmental health.

Journal of Veterinary Science↗
USGS2001

Travel times of P and S from the global digital seismic networks: Implications for the relative variation of P and S velocity in the mantle

We present new data sets of P and S arrival times which have been handpicked from long-period vertical and transverse component recordings of the various global seismic networks. Using events which occurred from 1976 to 1994 results in ∼38,000 globally well-distributed measurements of teleseismic P and ∼41,000 measurements of S . These data are particularly useful for looking at the relative variation of S and P velocities in the lower mantle. We describe both the measurement techniques and the gross characteristics of the data sets. The size of our data sets allows us to exploit the internal consistency of the data to identify outliers using a summary ray analysis. Since the polarity of each arrival is also known, we can construct fault plane solutions and/or compare with polarities predicted by the Harvard centroid moment tensor solutions to further diagnose phase misidentification. This analysis results in ∼5% of the data being identified as outliers. An analysis of variance indicates that the S residual travel times are dominated by the effects of three-dimensional structure but the P data have comparable contributions from noise and source mislocation effects. The summary ray analysis reveals the basic character of lower mantle structure, and there are large-scale patterns in both the S and P data sets that correlate quite well with each other. This analysis suggests that on average, d ln v s / d ln v p is an increasing function of depth in the mantle going from a value of ∼1.7 at the top of the lower mantle to an apparent value of 4 near the base of the mantle. This latter extreme value of R seems to result mainly from data which sample one region in the lowermost mantle under the central Pacific, where large positive S residuals are associated with very small P residuals. Such an anomaly cannot be thermal in origin.

Journal of Geophysical Research B: Solid Earth↗
USGS1973

Karst hydrology: A review

Karst regions of the world are characterized by limestones and other soluble rocks at or near land surface that have been modified by solutional erosion. Such surface features as sinks, long dry valleys, sparse streams, and bare rock and such subsurface features as caverns, arterial solution openings leading to large springs, and a deep water table are typical of karst terranes. These features result in an uneven distribution of permeability in karst systems and surface and subsurface hydrologic conditions that require special hydrogeologic studies. Local high permeability at shallow depth in mature karst regions leads to an ecology associated with a soilless and water-scarce surface environment. Many practical problems result from this high permeability, including: (1) scarcity and poor predictability of groundwater supplies; (2) scarcity of surface-water supplies; (3) instability of the ground; (4) leakage of surface reservoirs; and (5) an unreliable waste-disposal environment. Interest in karst hydrology has increased greatly in the past decade; this interest has resulted in the international exchange of numerous published reports on local areas and on special topical karst problems. Many of these reports have been used by the authors in preparing this paper, which synthesizes results of many workers and focuses attention on: (1) the development of karst features through hydrologic processes; and (2) hydrologic systems of karst terranes.

Journal of Hydrology↗
USGS2012

Uniformity and diversity in the composition of mineralizing fluids from hydrothermal vents on the southern Juan de Fuca Ridge

Abundances of Li, Na, K, Rb, Ca, Sr, Ba, Mn, Fe, Zn, and Si have been determined in fluid samples from seven vents located in three areas on the southern Juan de Fuca Ridge. The hydrothermal component estimated from the Mg contents of the samples ranges from 7% to 76%. Concentrations of Fe and Si, among other elements, in acid-stabilized solutions appear to be generally representative of the parental hydrothermal fluids, but some Zn determinations and most Ba values appear to be too low. Thermodynamic calculations indicate that the acidified samples remain supersaturated with respect to silica, barite, and pyrite; unacidified samples are supersaturated, in addition, with respect to ZnS, FeS, and many silicate phases. Within the constraints of limited sampling there appear to be differences in fluid compositions both within and between the three vent areas. Some uniform differences in the elemental abundances predicted for hydrothermal end-member fluids might be due to inmixing of fresh seawater at depth in the hydrothermal system. The Juan de Fuca hydrothermal fluids contain more Fe but otherwise have relative elemental abundances fairly similar to those in 13°N (East Pacific Rise) fluids, albeit at higher levels. In contrast, fluids from 21°N (East Pacific Rise) and Galapagos have lower K/Rb and much lower Sr and Na abundances; these compositional features probably result from interaction of these fluids with a different mineral assemblage, possibly more mature greenstone.

Journal of Geophysical Research Solid Earth↗
USGS2010

Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration

Shallow rainfall-induced landslides commonly occur under conditions of transient infiltration into initially unsaturated soils. In an effort to predict the timing and location of such landslides, we developed a model of the infiltration process using a two-layer system that consists of an unsaturated zone above a saturated zone and implemented this model in a geographic information system (GIS) framework. The model links analytical solutions for transient, unsaturated, vertical infiltration above the water table to pressure-diffusion solutions for pressure changes below the water table. The solutions are coupled through a transient water table that rises as water accumulates at the base of the unsaturated zone. This scheme, though limited to simplified soil-water characteristics and moist initial conditions, greatly improves computational efficiency over numerical models in spatially distributed modeling applications. Pore pressures computed by these coupled models are subsequently used in one-dimensional slope-stability computations to estimate the timing and locations of slope failures. Applied over a digital landscape near Seattle, Washington, for an hourly rainfall history known to trigger shallow landslides, the model computes a factor of safety for each grid cell at any time during a rainstorm. The unsaturated layer attenuates and delays the rainfall-induced pore-pressure response of the model at depth, consistent with observations at an instrumented hillside near Edmonds, Washington. This attenuation results in realistic estimates of timing for the onset of slope instability (7 h earlier than observed landslides, on average). By considering the spatial distribution of physical properties, the model predicts the primary source areas of landslides.

Journal of Geophysical Research F: Earth Surface↗
USGS1996

Crustal structure of the Colorado Plateau, Arizona: Application of new long-offset seismic data analysis techniques

The Colorado Plateau is a large crustal block in the southwestern United States that has been raised intact nearly 2 km above sea level since Cretaceous marine sediments were deposited on its surface. Controversy exists concerning the thickness of the plateau crust and the source of its buoyancy. Interpretations of seismic data collected on the plateau vary as to whether the crust is closer to 40 or 50 km thick. A thick crust could support the observed topography of the Colorado Plateau isostatically, while a thinner crust would indicate the presence of an underlying low-density mantle. This paper reports results on long-offset seismic data collected during the 1989 segment of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that extended from the Transition Zone into the Colorado Plateau in northwest Arizona. We apply two new methods to analyze long-offset data that employ finite difference travel time calculations: (1) a first-arrival time inverter to find upper crustal velocity structure and (2) a forward-modeling technique that allows the direct use of the inverted upper crustal solution in modeling secondary reflected arrivals. We find that the crustal thickness increases from 30 km beneath the metamorphic core complexes in the southern Basin and Range province to about 42 km beneath the northern Transition Zone and southern Colorado Plateau margin. We observe some crustal thinning (to ∼37 km thick) and slightly higher lower crustal velocities farther inboard; beneath the Kaibab uplift on the north rim of the Grand Canyon the crust thickens to a maximum of 48 km. We observe a nonuniform crustal thickness beneath the Colorado Plateau that varies by ∼15% and corresponds approximately to variations in topography with the thickest crust underlying the highest elevations. Crustal compositions (as inferred from seismic velocities) appear to be the same beneath the Colorado Plateau as those in the Basin and Range province to the southwest, implying that the plateau crust represents an unextended version of the Basin and Range. Some of the variability in crustal structure appears to correspond to preserved lithospheric discontinuities that date back to the Proterozoic Era.

Journal of Geophysical Research B: Solid Earth↗
USGS2019

The expectations and challenges of wildlife disease research in the era of genomics: Forecasting with a horizon scan-like exercise

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented ‘Big Data’ tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural ‘Genomics of Disease in Wildlife’ workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed five themes considered to be the most pressing to the application of genomics in wildlife disease research: i) “Improving Communication”, ii) “Methodological and Analytical Advancements”, iii) “Translation into Practice”, iv) “Integrating Landscape Ecology and Genomics”, and v) “Emerging New Questions”. Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

Journal of Heredity↗
USGS2009

Monitoring suspended sediments and associated chemical constituents in urban environments: Lessons from the city of Atlanta, Georgia, USA Water Quality Monitoring Program

Background, aim, and scope The City of Atlanta, Georgia (COA) is part of the ninth largest metropolitan area in the USA and one of the fastest growing (e.g., >24% between 2000 and 2007). Since 2003, the US Geological Survey has been operating an extensive long-term water-quantity and water-quality monitoring network for the COA. The experience gained in operating this network has provided insights into the challenges as well as some solutions associated with determining urban effects on water quality, especially in terms of estimating the annual fluxes of suspended sediment, trace/major elements, and nutrients. Discussion and findings The majority (>90%) of the annual fluxes of suspended sediment and discharge (>60%) from the COA occur in conjunction with stormflow. Typically, stormflow averages ≤20% of the year. Normally, annual flux calculations employ a daily time-step; however, due to the “flashy” nature of the COA’s streams, this approach can produce substantial underestimates (from 25% to 64%). Greater accuracy requires time-steps as short as every 2 to 3 h. The annual fluxes of ≥75% of trace elements (e.g., Cu, Pb, Zn), major elements (e.g., Fe, Al), and total P occur in association with suspended sediment; in turn, ≥90% of the transport of these constituents occurs in conjunction with stormflow. With the possible exception of nitrogen, baseflow sediment-associated and both baseflow and stormflow dissolved contributions represent relatively insignificant portions of the total annual load; hence, nonpoint (diffuse) sources are the dominant contributors to the fluxes of almost all of these constituents.

Georgia↗
USGS2020

The influence of snow cover, air temperature, and groundwater flow on the active-layer thermal regime of Arctic hillslopes drained by water tracks

Permafrost in Arctic watersheds limits soil biological activity to a thin, seasonally thawed active layer that contributes water to streams. In many hillslopes, relatively wet drainage features called water tracks have distinct freeze-thaw patterns that affect groundwater flow and storage, and thus the export of heat and solutes to Arctic streams. This study uses groundwater flow and energy transport models to examine potential controls on the timing and duration of freeze–thaw conditions and the magnitude of temperature fluctuations within water tracks and their adjacent hillslopes. The simulated length of the active-layer thaw season varies by 1 month over the range of snow-cover and mean annual air-temperature scenarios simulated. The timing and duration of freezing is particularly sensitive to depth and duration of snow cover. Thus, the deeper snowpack covers that can accumulate in water tracks contribute to their more persistent thaw conditions and their ability to conduct groundwater downslope. A three-dimensional simulation shows that during the summer thaw season, the water track captures groundwater laterally from half way across the hillslope. The models presented here elucidate key mechanisms driving small-scale variation in the active-layer thermal regime of tundra hillslopes, which may be responsible for changes in drainage-network geometry and Arctic biogeochemical fluxes under a warming climate.

Alaska↗
USGS2023

Detection and quantification of preferential flow using artificial rainfall with multiple experimental approaches

Preferential flow in the unsaturated zone strongly influences important hydrologic processes, such as infiltration, contaminant transport, and aquifer recharge. Because it entails various combinations of physical processes arising from the interactions of water, air, and solid particles in a porous medium, preferential flow is highly complex. Major research is needed to improve the ability to understand, quantify, model, and predict preferential flow. Toward a solution, a combination of diverse experimental measurements at multiple scales, from laboratory scale to mesoscale, has been implemented to detect and quantify preferential paths in carbonate and karstic unsaturated zones. This involves integration of information from (1) core samples, by means of mercury intrusion porosimeter, evaporation, quasi-steady centrifuge and dewpoint potentiometer laboratory methods, to investigate the effect of pore-size distribution on hydraulic characteristics and the potential activation of preferential flow, (2) field plot experiments with artificial sprinkling, to visualize preferential pathways related to secondary porosity, through use of geophysical measurements, and (3) mesoscale evaluation of field data through episodic master recession modeling of episodic recharge. This study demonstrates that preferential flow processes operate from core scale to two different field scales and impact on the qualitative and quantitative groundwater status, by entailing fast flow with subsequent effects on recharge rate and contaminant mobilizing. The presented results represent a rare example of preferential flow detection and numerical modeling by reducing underestimation of the recharge and contamination risks.

Bari↗
USGS2016

A microbial arsenic cycle in sediments of an acidic mine impoundment: Herman Pit, Clear Lake, California

The involvement of prokaryotes in the redox reactions of arsenic occurring between its +5 [arsenate; As(V)] and +3 [arsenite; As(III)] oxidation states has been well established. Most research to date has focused upon circum-neutral pH environments (e.g., freshwater or estuarine sediments) or arsenic-rich “extreme” environments like hot springs and soda lakes. In contrast, relatively little work has been conducted in acidic environments. With this in mind we conducted experiments with sediments taken from the Herman Pit, an acid mine drainage impoundment of a former mercury (cinnabar) mine. Due to the large adsorptive capacity of the abundant Fe(III)-rich minerals, we were unable to initially detect in solution either As(V) or As(III) added to the aqueous phase of live sediment slurries or autoclaved controls, although the former consumed added electron donors (i.e., lactate, acetate, hydrogen), while the latter did not. This prompted us to conduct further experiments with diluted slurries using the live materials from the first incubation as inoculum. In these experiments we observed reduction of As(V) to As(III) under anoxic conditions and reduction rates were enhanced by addition of electron donors. We also observed oxidation of As(III) to As(V) in oxic slurries as well as in anoxic slurries amended with nitrate. We noted an acid-tolerant trend for sediment slurries in the cases of As(III) oxidation (aerobic and anaerobic) as well as for anaerobic As(V) reduction. These observations indicate the presence of a viable microbial arsenic redox cycle in the sediments of this extreme environment, a result reinforced by the successful amplification of arsenic functional genes ( aioA , and arrA ) from these materials.

California↗
USGS2006

Multiphase, multicomponent parameter estimation for liquid and vapor fluxes in deep arid systems using hydrologic data and natural environmental tracers

Multiphase, multicomponent numerical models of long-term unsaturated-zone liquid and vapor movement were created for a thick alluvial basin at the Nevada Test Site to predict present-day liquid and vapor fluxes. The numerical models are based on recently developed conceptual models of unsaturated-zone moisture movement in thick alluvium that explain present-day water potential and tracer profiles in terms of major climate and vegetation transitions that have occurred during the past 10 000 yr or more. The numerical models were calibrated using borehole hydrologic and environmental tracer data available from a low-level radioactive waste management site located in a former nuclear weapons testing area. The environmental tracer data used in the model calibration includes tracers that migrate in both the liquid and vapor phases (δD, δ 18 O) and tracers that migrate solely as dissolved solutes (Cl), thus enabling the estimation of some gas-phase as well as liquid-phase transport parameters. Parameter uncertainties and correlations identified during model calibration were used to generate parameter combinations for a set of Monte Carlo simulations to more fully characterize the uncertainty in liquid and vapor fluxes. The calculated background liquid and vapor fluxes decrease as the estimated time since the transition to the present-day arid climate increases. However, on the whole, the estimated fluxes display relatively little variability because correlations among parameters tend to create parameter sets for which changes in some parameters offset the effects of others in the set. Independent estimates on the timing since the climate transition established from packrat midden data were essential for constraining the model calibration results. The study demonstrates the utility of environmental tracer data in developing numerical models of liquid- and gas-phase moisture movement and the importance of considering parameter correlations when using Monte Carlo analysis to characterize the uncertainty in moisture fluxes.

Vadose Zone Journal↗
USGS2021

Clothianidin decomposition in Missouri wetland soils

Neonicotinoid pesticides can persist in soils for extended time periods; however, they also have a high potential to contaminate ground and surface waters. Studies have reported negative effects associated with neonicotinoids and nontarget taxa, including aquatic invertebrates, pollinating insect species, and insectivorous birds. This study evaluated factors associated with clothianidin (CTN) degradation and sorption in Missouri wetland soils to assess the potential for wetland soils to mitigate potential environmental risks associated with neonicotinoids. Solid-to-solution partition coefficients ( K d ) for CTN sorption to eight wetland soils were determined via single-point sorption experiments, and sorption isotherm experiments were conducted using the two most contrasting soils. Clothianidin degradation was determined under oxic and anoxic conditions over 60 d. Degradation data were fit to zero- and first-order kinetic decay models to determine CTN half-life ( t 0.5 ). Sorption results indicated CTN sorption to wetland soil was relatively weak (average K d , 3.58 L kg –1 ); thus, CTN has the potential to be mobile and bioavailable within wetland soils. However, incubation results showed anoxic conditions significantly increased CTN degradation rates in wetland soils (anoxic average t 0.5 , 27.2 d; oxic average t 0.5 , 149.1 d). A significant negative correlation was observed between anoxic half-life values and soil organic C content ( r 2 = .782; p = .046). Greater CTN degradation rates in wetland soils under anoxic conditions suggest that managing wetlands to facilitate anoxic conditions could mitigate CTN presence in the environment and reduce exposure to nontarget organisms.

Missouri↗