SEARCH · Geology Reports

Results for “Applied Sciences”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,135 records · Page 63Linked to original sources

USGS2022

Predicting regional fluoride concentrations at public and domestic supply depths in basin-fill aquifers of the western United States using a random forest model

A random forest regression (RFR) model was applied to over 12,000 wells with measured fluoride (F) concentrations in untreated groundwater to predict F concentrations at depths used for domestic and public supply in basin-fill aquifers of the western United States. The model relied on twenty-two regional-scale environmental and surficial predictor variables selected to represent factors known to control F concentrations in groundwater. The testing model fit R 2 and RMSE were 0.52 and 0.78 mg/L. Comparisons of measured to predicted proportions of four F-concentrations categories (<0.7 mg/L, 0.7–2 mg/L, >2 mg/L – 4 mg/L, and > 4 mg/L) indicate that the model performed well at making regional-scale predictions. Differences between measured and predicted proportions indicate underprediction of measured F at values by between 4 and 20 mg/L, representing less than 1% of the regional scale predicted values. These residuals most often map to geographic regions where local-scale processes including evaporative discharge in closed basins or intermittent streams concentrate fluoride in shallow groundwater. Despite this, the RFR model provides spatially continuous F predictions across the basin-fill aquifers where discrete samples are missing. Further, the predictions capture documented areas that exceed the F maximum contaminant level for drinking water of 4 mg/L and areas that are below the oral-health benchmark of 0.7 mg/L. These predictions can be used to estimate fluoride concentrations in unmonitored areas and to aid in identifying geographic areas that may require further investigation at localized scales.

Arizona, California, Colorado, New Mexico, Nevada,↗

USGS2019

Patterns of resource allocation in a coastal marsh plant (Schoenoplectus americanus) along a sediment-addition gradient

Reductions in sediment delivery to coastal marshes increase their vulnerability to relative sea-level rise. Sediment pulses from storm events and commercial dredge-spray operations (e.g., beneficial use) represent increasingly important sediment sources to otherwise sediment-deprived marshes. These sediments can stimulate plant growth by providing nutrient and elevation subsidies, with plant growth predicted to peak at intermediate sediment depths. Most sediment subsidy studies have compared plant responses among discrete sediment depths applied across a uniform marsh platform, and often assessed impacts in the years following sediment deposition. Less is known about the immediate effects of sediment additions on plant growth or the mechanisms driving differences in these responses along a sediment-addition gradient. To investigate immediate plant responses to sediment additions (i.e., within 6 months), we added sediment to Schoenoplectus americanus -dominated brackish marsh sods in a greenhouse setting. Sediment was added incrementally along a sediment addition gradient (0–20 cm) to 18 sods, which were placed in clear planting enclosures to permit monitoring of above- and below-ground responses and hydro-edaphic properties. Shoot production, which increased linearly along the sediment addition gradient, was the first observable response to sediment addition and was likely a result of apical dominance interruption. This trend continued throughout the study even as survival of those shoots declined with increasing sediment depth. At the end of the experiment, however, shoot biomass production was biphasic, responding favorably to sediment addition up to intermediate depths before collapsing at higher depths. While fine root production in sods was maximal at intermediate sediment addition levels, root production in new sediment layers was limited. There were limited differences observed for hydro-edaphic properties along the sediment addition gradient. Thus, plant resources are allocated to shoot production immediately following sediment addition at the expense of root colonization of new sediment layers, suggesting that, in certain conditions, sediment deposition can represent a stress rather than a subsidy.

Louisiana↗

USGS2025

Assessing nonpoint-source uranium pollution in an irrigated stream-aquifer system

Uranium (U) in rocks and soils of arid and semi-arid environments can be mobilized by irrigation and fertilization, posing environmental and health risks. Elevated U, along with selenium (Se) and nitrate (NO 3 ) co-constituents, necessitates careful monitoring and management. We developed a distributed-parameter numerical model to assess U pollution in an irrigated stream-aquifer system, applying it to a 552 km 2 region in Colorado's Lower Arkansas River Valley (LARV) over 14 years. A MODFLOW model, describing groundwater and stream flow, was coupled with an RT3D-OTIS model to portray reactive U transport. Calibration using the PESTPP-iES iterative ensemble smoother (iES) software indicated good agreement with observed U concentrations. The model revealed substantial and variable U levels across the LARV, highlighting potential hotspots and possible contributing factors, such as geological composition of the bedrock and near-surface shale and aquifer sediments derived from them, irrigation practices, and riparian landscape. U levels exceed the chronic standard (85th percentile = 30 μg/L, set by the US Environmental Protection Agency), which is the permissible regulatory threshold, in groundwater across 44 % of the region and along the river by an average factor of 2.9. Simulated average U concentrations in the non-riparian aquifer and river are 124 μg/L and 60 μg/L, respectively, compared with 112 μg/L and 62 μg/L for measured values. The average 85th percentile U concentration is 222 μg/L in the aquifer and 82 μg/L in the river. Average simulated U mass loading to the river is 0.17 kg/day per km, compared to an estimated 0.23 kg/day per km. Findings provide a baseline for comparing future simulated outcomes of alternative best management practices (BMPs) for U pollution mitigation and offer a methodology applicable to other irrigated regions.

Colorado↗

USGS2013

Economic impacts of the SAFRR tsunami scenario in California: Chapter H in The SAFRR (Science Application for Risk Reduction) Tsunami Scenario

This study evaluates the hypothetical economic impacts of the SAFRR (Science Application for Risk Reduction) tsunami scenario to the California economy. The SAFRR scenario simulates a tsunami generated by a hypothetical magnitude 9.1 earthquake that occurs offshore of the Alaska Peninsula (Kirby and others, 2013). Economic impacts are measured by the estimated reduction in California’s gross domestic product (GDP), the standard economic measure of the total value of goods and services produced. Economic impacts are derived from the physical damages from the tsunami as described by Porter and others (2013). The principal physical damages that result in disruption of the California economy are (1) about $100 million in damages to the twin Ports of Los Angeles (POLA) and Long Beach (POLB), (2) about $700 million in damages to marinas, and (3) about $2.5 billion in damages to buildings and contents (properties) in the tsunami inundation zone on the California coast. The study of economic impacts does not include the impacts from damages to roads, bridges, railroads, and agricultural production or fires in fuel storage facilities because these damages will be minimal with respect to the California economy. The economic impacts of damage to other California ports are not included in this study because detailed evaluation of the physical damage to these ports was not available in time for this report. The analysis of economic impacts is accomplished in several steps. First, estimates are made for the direct economic impacts that result in immediate business interruption losses in individual sectors of the economy due to physical damage to facilities or to disruption of the flow of production units (commodities necessary for production). Second, the total economic impacts (consisting of both direct and indirect effects) are measured by including the general equilibrium (essentially quantity and price multiplier effects) of lost production in other sectors by ripple effects upstream and downstream along the supply chain. An appropriate measure of the economic impacts on the California economy for the SAFRR tsunami scenario is the reduction in GDP. The economic impacts are first calculated without resilience, the ability of the economy to adjust to disruptions in ways that mute potential negative impacts. There are many types of resilience, including using existing inventories of materials, using unused capacity, conserving inputs, substituting for disrupted supplies, recapturing production after the disruption is restored, and many others. A method for estimating resilience, identified in the port system and sectors affected by property damages, is applied to indicate potential reductions of direct and total economic impacts. In this SAFRR tsunami scenario analysis of economic impacts to California, we implement established techniques used to model the economic impacts for two previous U.S. Geological Survey (USGS) scenarios: the southern California Shakeout earthquake (Rose and others, 2011) and the California ARkStorm severe winter storm (Sue Wing and others, written commun., 2013). For the SAFRR tsunami scenario, we reviewed the relevant studies that assess economic impacts from previous tsunami events affecting California and elsewhere and estimate the economic impacts of potential tsunami and other threats to POLA and POLB. To our knowledge, assessment of impacts to the California economy from distant source tsunamis does not exist. Previous tsunamis, including those from the 1960 Chile earthquake, the 1964 Alaska earthquake, the 2008 Chile earthquake and the 2011 Japan earthquake, had only relatively minor or very localized severe damage (such as that in Crescent City in 1964), and no studies of the economic impacts were completed. A rare study of the economic impacts of a tsunami event has recently been produced for the Tohoku earthquake and tsunami (Kajitani and others, 2013). Quarterly declines in Japan’s GDP are observed to peak at ‒1.63 percent in the second quarter after the event and stagnate for the rest of the year. The majority of the economic impacts are attributed to the tsunami rather than the earthquake. The hardest hit sectors are identified as agriculture, fisheries, manufacturing, retail, and tourism. Other relevant studies have focused on the economic impacts of threats that close POLA and POLB. We find one analysis of a potential tsunami scenario affecting the California economy through disruption of port operations. Borrero and others (2005) estimated economic impacts to the southern California economy of $7 to $40 billion from a locally generated tsunami that closes POLA and POLB for as much as 1 year. There have also been several studies of the economic impacts of non-tsunami events affecting POLA and POLB. Analyses of an 11-day labor lockout produced a range of estimated national impacts of as much as $1.94 billion/day (Park and others 2008, Martin Associates 2001). Examination of a potential terrorist attack that closes the San Pedro port for 1 month yielded a $29 billion impact to the California economy (Park, 2008). These studies have reinforced the importance of recognizing economic resilience in economic impact analyses. Hall (2004) criticized the upper-end estimate of national economic impacts from the labor lockout based on model shortcomings that neglected short-run substitution behavior and fixed the long-run economic behaviors. Following the 2011 Japanese tsunami, resilience was observed in the forms of rapid recovery of manufacturing sectors, energy conservation, and insurance (Kajitani and others, 2013).

California↗

USGS2022

Seed treatments containing neonicotinoids and fungicides reduce aquatic insect richness and abundance in midwestern USA–managed floodplain wetlands

Agrochemicals including neonicotinoid insecticides and fungicides are frequently applied as seed treatments on corn, soybeans, and other common row crops. Crops grown from pesticide-treated seed are often directly planted in managed floodplain wetlands and used as a soil disturbance or food resource for wildlife. We quantified invertebrate communities within mid-latitude floodplain wetlands and assessed their response to use of pesticide-treated seeds within the floodplain. We collected and tested aqueous and sediment samples for pesticides in addition to sampling aquatic invertebrates from 22 paired wetlands. Samples were collected twice in 2016 (spring [pre-water level drawdown] and autumn [post-water level flood-up]) followed by a third sampling period (spring 2017). Meanwhile, during the summer of 2016, a portion of study wetlands were planted with either pesticide-treated or untreated corn seed. Neonicotinoid toxic equivalencies (NI-EQs) for sediment ( X̅ = 0.58 μg/kg), water ( X̅ = 0.02 μg/L), and sediment fungicide concentrations ( X̅ = 0.10 μg/kg) were used to assess potential effects on wetland invertebrates. An overall decrease in aquatic insect richness and abundance was associated with greater NI-EQs in wetland water and sediments, as well as with sediment fungicide concentration. Post-treatment, treated wetlands displayed a decrease in insect taxa-richness and abundance before recovering by the spring of 2017. Information on timing and magnitude of aquatic insect declines will be useful when considering the use of seed treatments for wildlife management. More broadly, this study brings attention to how agriculture is used in wetland management and conservation planning.

Missouri↗

USGS2013

A critique of the use of indicator-species scores for identifying thresholds in species responses

Identification of ecological thresholds is important both for theoretical and applied ecology. Recently, Baker and King (2010, King and Baker 2010) proposed a method, threshold indicator analysis (TITAN), to calculate species and community thresholds based on indicator species scores adapted from Dufrêne and Legendre (1997). We tested the ability of TITAN to detect thresholds using models with (broken-stick, disjointed broken-stick, dose-response, step-function, Gaussian) and without (linear) definitive thresholds. TITAN accurately and consistently detected thresholds in step-function models, but not in models characterized by abrupt changes in response slopes or response direction. Threshold detection in TITAN was very sensitive to the distribution of 0 values, which caused TITAN to identify thresholds associated with relatively small differences in the distribution of 0 values while ignoring thresholds associated with large changes in abundance. Threshold identification and tests of statistical significance were based on the same data permutations resulting in inflated estimates of statistical significance. Application of bootstrapping to the split-point problem that underlies TITAN led to underestimates of the confidence intervals of thresholds. Bias in the derivation of the z-scores used to identify TITAN thresholds and skewedness in the distribution of data along the gradient produced TITAN thresholds that were much more similar than the actual thresholds. This tendency may account for the synchronicity of thresholds reported in TITAN analyses. The thresholds identified by TITAN represented disparate characteristics of species responses that, when coupled with the inability of TITAN to identify thresholds accurately and consistently, does not support the aggregation of individual species thresholds into a community threshold.

Freshwater Science↗

USGS2016

Western Lake Erie Basin: Soft-data-constrained, NHDPlus resolution watershed modeling and exploration of applicable conservation scenarios

Complex watershed simulation models are powerful tools that can help scientists and policy-makers address challenging topics, such as land use management and water security. In the Western Lake Erie Basin (WLEB), complex hydrological models have been applied at various scales to help describe relationships between land use and water, nutrient, and sediment dynamics. This manuscript evaluated the capacity of the current Soil and Water Assessment Tool (SWAT2012) to predict hydrological and water quality processes within WLEB at the finest resolution watershed boundary unit (NHDPlus) along with the current conditions and conservation scenarios. The process based SWAT model was capable of the fine-scale computation and complex routing used in this project, as indicated by measured data at five gaging stations. The level of detail required for fine-scale spatial simulation made the use of both hard and soft data necessary in model calibration, alongside other model adaptations. Limitations to the model's predictive capacity were due to a paucity of data in the region at the NHDPlus scale rather than due to SWAT functionality. Results of treatment scenarios demonstrate variable effects of structural practices and nutrient management on sediment and nutrient loss dynamics. Targeting treatment to acres with critical outstanding conservation needs provides the largest return on investment in terms of nutrient loss reduction per dollar spent, relative to treating acres with lower inherent nutrient loss vulnerabilities. Importantly, this research raises considerations about use of models to guide land management decisions at very fine spatial scales. Decision makers using these results should be aware of data limitations that hinder fine-scale model interpretation.

Lake Erie Basin↗

USGS2015

A field comparison of multiple techniques to quantify groundwater - surface-water interactions

Groundwater–surface-water (GW-SW) interactions in streams are difficult to quantify because of heterogeneity in hydraulic and reactive processes across a range of spatial and temporal scales. The challenge of quantifying these interactions has led to the development of several techniques, from centimeter-scale probes to whole-system tracers, including chemical, thermal, and electrical methods. We co-applied conservative and smart reactive solute-tracer tests, measurement of hydraulic heads, distributed temperature sensing, vertical profiles of solute tracer and temperature in the stream bed, and electrical resistivity imaging in a 450-m reach of a 3 rd -order stream. GW-SW interactions were not spatially expansive, but were high in flux through a shallow hyporheic zone surrounding the reach. NaCl and resazurin tracers suggested different surface–subsurface exchange patterns in the upper ⅔ and lower ⅓ of the reach. Subsurface sampling of tracers and vertical thermal profiles quantified relatively high fluxes through a 10- to 20-cm deep hyporheic zone with chemical reactivity of the resazurin tracer indicated at 3-, 6-, and 9-cm sampling depths. Monitoring of hydraulic gradients along transects with MINIPOINT streambed samplers starting ∼40 m from the stream indicated that groundwater discharge prevented development of a larger hyporheic zone, which progressively decreased from the stream thalweg toward the banks. Distributed temperature sensing did not detect extensive inflow of ground water to the stream, and electrical resistivity imaging showed limited large-scale hyporheic exchange. We recommend choosing technique(s) based on: 1) clear definition of the questions to be addressed (physical, biological, or chemical processes), 2) explicit identification of the spatial and temporal scales to be covered and those required to provide an appropriate context for interpretation, and 3) maximizing generation of mechanistic understanding and reducing costs of implementing multiple techniques through collaborative research.

Pennsylvania↗

USGS2003

Effects of a fire-retardant chemical to fathead minnows in experimental streams

Background Each year millions of liters of fire-retardant chemicals are applied to wildfires across the nation. Recent laboratory studies with long-term fire-retardant chemicals indicate a significant photoenhanced toxicity of products containing sodium ferrocyanide corrosion inhibitors. Our objective of this study was to determine the toxicity of fire-retardant chemicals to fathead minnows during exposure in experimental outdoor streams. Methods Stream tests were conducted to determine the potential toxicity of a pulse of exposure as might occur when fire retardant chemical is rinsed from the watershed by rainfall. Two artificial 55-meter experimental streams were dosed with different concentrations of Fire-Trol® GTS-R, or uncontaminated for a control. Replicate groups of fathead minnows were added to screened containers (10 fish per container) and exposed to retardant chemicals in the recirculating flow of the stream for up to 6 hours. Results and Discussion Under field conditions toxicity of GTS-R only occurred in the presence of sunlight. When GTS-R was tested on sunny days, 100% mortality occurred. However, when tested during heavily overcast conditions, no mortality occurred. Conclusions Lethal concentrations of cyanide were measured when GTS-R with YPS exposures were conducted under sunny conditions, but not under cloudy conditions, indicating that a minimum UV level is necessary to induce toxicity as well as the release of cyanide from YPS. The toxicity observed with GTS-R was likely associated with lethal concentrations of cyanide. Rainwater runoff following applications of this fire-retardant at the recommended rate could result in lethal concentrations in small ponds and streams receiving limited water flow under sunny conditions.

Environmental Science and Pollution Research↗

USGS2012

A brief history and summary of the effects of river engineering and dams on the Mississippi River system and delta

The U.S. Geological Survey Forecast Mekong project is providing technical assistance and information to aid management decisions and build science capacity of institutions in the Mekong River Basin. A component of this effort is to produce a synthesis of the effects of dams and other engineering structures on large-river hydrology, sediment transport, geomorphology, ecology, water quality, and deltaic systems. The Mississippi River Basin (MRB) of the United States was used as the backdrop and context for this synthesis because it is a continental scale river system with a total annual water discharge proportional to the Mekong River, has been highly engineered over the past two centuries, and the effects of engineering have been widely studied and documented by scientists and engineers. The MRB is controlled and regulated by dams and river-engineering structures. These modifications have resulted in multiple benefits including navigation, flood control, hydropower, bank stabilization, and recreation. Dams and other river-engineering structures in the MRB have afforded the United States substantial socioeconomic benefits; however, these benefits also have transformed the hydrologic, sediment transport, geomorphic, water-quality, and ecologic characteristics of the river and its delta. Large dams on the middle Missouri River have substantially reduced the magnitude of peak floods, increased base discharges, and reduced the overall variability of intraannual discharges. The extensive system of levees and wing dikes throughout the MRB, although providing protection from intermediate magnitude floods, have reduced overall channel capacity and increased flood stage by up to 4 meters for higher magnitude floods. Prior to major river engineering, the estimated average annual sediment yield of the Mississippi River Basin was approximately 400 million metric tons. The construction of large main-channel reservoirs on the Missouri and Arkansas Rivers, sedimentation in dike fields, and protection of channel banks by revetments throughout the basin, have reduced the overall sediment yield of the MRB by more than 60 percent. The primary alterations to channel morphology by dams and other engineering projects have been (1) channel simplification and reduced dynamism; (2) lowering of channel-bed elevation; and (3) disconnection of the river channel from the flood plain, except during extreme flood events. Freshwater discharge from the Mississippi River and its associated sediment and nutrient loads strongly influence the physical and biological components in the northern Gulf of Mexico. Ninety percent of the nitrogen load reaching the Gulf of Mexico is from nonpoint sources with about 60 percent coming from fertilizer and mineralized soil nitrogen. Much of the phosphorus is from animal manure from pasture and rangelands followed by fertilizer applied to corn and soybeans. Increased nutrient enrichment in the northern Gulf of Mexico has resulted in the degradation of water quality as more phytoplankton grow, which increases turbidity and depletes oxygen in the lower depths creating what is known as the "dead zone." In 2002, the dead zone was 22,000 square kilometers (km2), an area similar to the size of the State of Massachusetts. Changes in the flow regime from engineered structures have had direct and indirect effects on the fish communities. The navigation pools in the upper Mississippi River have aged, and these overwintering habitats, which were created when the pools filled, have declined as sedimentation reduces water depth. Reproduction of paddlefish may have been adversely affected by dams, which impede access to suitable spawning habitats. Fishes that inhabit swift-current habitats in the unimpounded lower Mississippi River have not declined as much as in the upper Mississippi River. The decline of the pallid sturgeon may be attributable to channelization of the Missouri River above St. Louis, Missouri. The Missouri River supports a rich fish community and remains relatively intact. Nevertheless, the widespread and long history of human intervention in river discharge has contributed to the declines of about 25 percent of the species. The Mississippi River Delta Plain is built from six delta complexes composed of a massive area of coastal wetlands that support the largest commercial fishery in the conterminous United States. Since the early 20th century, approximately 4,900 km2 of coastal lands have been lost in Louisiana. One of the primary mechanisms of wetland loss on the Plaquemines-Balize complex is believed to be the disconnection of the river distributary network from the delta plain by the massive system of levees on the delta top, which prevent overbank flooding and replenishment of the delta top by sediment and nutrient deliveries. Efforts by Federal and State agencies to conserve and restore the Mississippi River Delta Plain began over three decades ago and have accelerated over the past decade. Regardless of these efforts, however, land losses are expected to continue because the reduced upstream sediment supplies are not sufficient to keep up with the projected depositional space being created by the combined forces of delta plain subsidence and global sea-level rise.

Alabama;Alberta;Arkansas;Colorado;Georgia;Illinois↗

USGS2015

Exposure to runoff from coal-tar-sealed pavement induces genotoxicity and impairment of DNA repair capacity in the RTL-W1 fish liver cell line

Coal-tar-based (CTB) sealcoat, frequently applied to parking lots and driveways in North America, contains elevated concentrations of polycyclic aromatic hydrocarbons (PAHs) and related compounds. The RTL-W1 fish liver cell line was used to investigate two endpoints (genotoxicity and DNA-repair-capacity impairment) associated with exposure to runoff from asphalt pavement with CTB sealcoat or with an asphalt-based sealcoat hypothesized to contain about 7% CTB sealcoat (AS-blend). Genotoxic potential was assessed by the Formamido pyrimidine glycosylase (Fpg)-modified comet assay for 1:10 and 1:100 dilutions of runoff samples collected from 5 h to 36 d following sealcoat application. DNA-repair capacity was assessed by the base excision repair comet assay for 1:10 dilution of samples collected 26 h and 36 d following application. Both assays were run with and without co-exposure to ultraviolet-A radiation (UVA). With co-exposure to UVA, genotoxic effects were significant for both dilutions of CTB runoff for three of four sample times, and for some samples of AS-blend runoff. Base excision repair was significantly impaired for CTB runoff both with and without UVA exposure, and for AS-blend runoff only in the absence of UVA. This study is the first to investigate the effects of exposure to the complex mixture of chemicals in coal tar on DNA repair capacity. The results indicate that co-exposure to runoff from CT-sealcoated pavement and UVA as much as a month after sealcoat application has the potential to cause genotoxicity and impair DNA repair capacity.

Science of the Total Environment↗

USGS2023

Dissolved organic carbon dynamics and fluxes in Mississippi-Atchafalaya deltaic system impacted by an extreme flood event and hurricanes: A multi-satellite approach using Sentinel-2/3 and Landsat-8/9 data

Transport of riverine and wetland-derived dissolved organic carbon (DOC) spanning tidal wetlands, estuaries, and continental shelf waters functionally connects terrestrial and aquatic carbon reservoirs, yet the magnitude and ecological significance of this variable and its spatiotemporal linkage remains uncertain for coastal deltaic regions, such as Mississippi River Delta Plain, which includes Mississippi (MR) and Atchafalaya (AR) rivers and estuaries with vast expanses of wetlands and coastal forests. We examined DOC dynamics and fluxes in this large river-dominated wetland-estuarine system for the period between 2019 and 2021 that included an extreme river flood event in 2019, two major hurricanes (Barry in 2019 and Ida in 2021), and cold front passage using an improved adaptive quasi-analytical algorithm (QAA-AD) applied to multi-satellite sensors (Sentinel 3A/B OLCI, Landsat-8/OLI and Sentinel-2A/B MSI) with varying spectral and spatial (10/30/300 m) resolutions. The DOC estimates from multi-satellite sensors in combination with water fluxes were used to assess DOC fluxes from two large rivers (MR and AR) and small channels across the delta plain. Overall, this system delivered a total of 6.7 Tg C yr -1 (1 Tg = 10 12 g) into the estuarine zone and the northern Gulf of Mexico (nGoM) during 2019. High DOC fluxes from the AR (1.3 Tg C yr -1 ) and MR (4.5 Tg C yr -1 ) were associated with the extreme flood event in 2019. Hurricanes that occurred in the study period also contributed to the wetland and estuarine DOC fluxes into continental shelf waters; for example, the passage of Hurricane Barry in July 2019, delivered over a 3-day period ~1.33 ×10 9 g DOC from Barataria Basin into the nGoM. Sentinel 2-MSI land and water classification revealed that Hurricane Ida eroded a total of 1.34×10 8 m 2 of marshes in middle Barataria Basin, converting those habitats into open water with 3.0 m inundation depth and high DOC concentrations (16.4 mg L -1 ), a potentially large DOC source to the coastal waters. Overall, storms and flood events are major sources of DOC flux that facilitate transport of upstream carbon as well as transformation of carbon in the wetlands, through the conversion of vegetated wetland to open water.

Louisiana↗

USGS1996

Chemical evaluation of soil-solution in acid forest soils

Soil-solution chemistry is commonly studied in forests through the use of soil lysimeters.This approach is impractical for regional survey studies, however, because lysimeter installation and operation is expensive and time consuming. To address these problems, a new technique was developed to compare soil-solution chemistry among red spruce stands in New York, Vermont, New Hampshire, Maine. Soil solutions were expelled by positive air pressure from soil that had been placed in a sealed cylinder. Before the air pressure was applied, a solution chemically similar to throughfall was added to the soil to bring it to approximate field capacity. After the solution sample was expelled, the soil was removed from the cylinder and chemically analyzed. The method was tested with homogenized Oa and Bs horizon soils collected from a red spruce stand in the Adirondack Mountains of New York, a red spruce stand in east-central Vermont, and a mixed hardwood stand in the Catskill Mountains of New York. Reproducibility, effects of varying the reaction time between adding throughfall and expelling soil solution (5-65 minutes) and effects of varying the chemical composition of added throughfall, were evaluated. In general, results showed that (i) the method was reproducible (coefficients of variation were generally < 15%), (ii) variations in the length of reaction-time did not affect expelled solution concentrations, and (iii) adding and expelling solution did not cause detectable changes in soil exchange chemistry. Concentrations of expelled solutions varied with the concentrations of added throughfall; the lower the CEC, the more sensitive expelled solution concentrations were to the chemical concentrations of added throughfall. Addition of a tracer (NaBr) showed that the expelled solution was a mixture of added solution and solution that preexisted in the soil. Comparisons of expelled solution concentrations with concentrations of soil solutions collected by zero-tension and tension lysimetry indicated that expelled solution concentrations were higher than those obtained with either type of lysimeter, although there was less difference with tension lysimeters than zero-tension lysimeters. The method used for collection of soil solution should be taken into consideration whenever soil solution data are being interpreted.

New York, Vermont↗

USGS2016

Bayesian estimation of magma supply, storage, and eruption rates using a multiphysical volcano model: Kīlauea Volcano, 2000–2012

Estimating rates of magma supply to the world's volcanoes remains one of the most fundamental aims of volcanology. Yet, supply rates can be difficult to estimate even at well-monitored volcanoes, in part because observations are noisy and are usually considered independently rather than as part of a holistic system. In this work we demonstrate a technique for probabilistically estimating time-variable rates of magma supply to a volcano through probabilistic constraint on storage and eruption rates. This approach utilizes Bayesian joint inversion of diverse datasets using predictions from a multiphysical volcano model, and independent prior information derived from previous geophysical, geochemical, and geological studies. The solution to the inverse problem takes the form of a probability density function which takes into account uncertainties in observations and prior information, and which we sample using a Markov chain Monte Carlo algorithm. Applying the technique to Kīlauea Volcano, we develop a model which relates magma flow rates with deformation of the volcano's surface, sulfur dioxide emission rates, lava flow field volumes, and composition of the volcano's basaltic magma. This model accounts for effects and processes mostly neglected in previous supply rate estimates at Kīlauea, including magma compressibility, loss of sulfur to the hydrothermal system, and potential magma storage in the volcano's deep rift zones. We jointly invert data and prior information to estimate rates of supply, storage, and eruption during three recent quasi-steady-state periods at the volcano. Results shed new light on the time-variability of magma supply to Kīlauea, which we find to have increased by 35–100% between 2001 and 2006 (from 0.11–0.17 to 0.18–0.28 km 3 /yr), before subsequently decreasing to 0.08–0.12 km 3 /yr by 2012. Changes in supply rate directly impact hazard at the volcano, and were largely responsible for an increase in eruption rate of 60–150% between 2001 and 2006, and subsequent decline by as much as 60% by 2012. We also demonstrate the occurrence of temporal changes in the proportion of Kīlauea's magma supply that is stored versus erupted, with the supply “surge” in 2006 associated with increased accumulation of magma at the summit. Finally, we are able to place some constraints on sulfur concentrations in Kīlauea magma and the scrubbing of sulfur by the volcano's hydrothermal system. Multiphysical, Bayesian constraint on magma flow rates may be used to monitor evolving volcanic hazard not just at Kīlauea but at other volcanoes around the world.

Hawai'i↗

USGS2009

Geometry of the Nojima fault at Nojima-Hirabayashi, Japan - I. A simple damage structure inferred from borehole core permeability

The 1995 Kobe (Hyogo-ken Nanbu) earthquake, M = 7.2, ruptured the Nojima fault in southwest Japan. We have studied core samples taken from two scientific drillholes that crossed the fault zone SW of the epicentral region on Awaji Island. The shallower hole, drilled by the Geological Survey of Japan (GSJ), was started 75 m to the SE of the surface trace of the Nojima fault and crossed the fault at a depth of 624 m. A deeper hole, drilled by the National Research Institute for Earth Science and Disaster Prevention (NIED) was started 302 m to the SE of the fault and crossed fault strands below a depth of 1140 m. We have measured strength and matrix permeability of core samples taken from these two drillholes. We find a strong correlation between permeability and proximity to the fault zone shear axes. The half-width of the high permeability zone (approximately 15 to 25 m) is in good agreement with the fault zone width inferred from trapped seismic wave analysis and other evidence. The fault zone core or shear axis contains clays with permeabilities of approximately 0.1 to 1 microdarcy at 50 MPa effective confining pressure (10 to 30 microdarcy at in situ pressures). Within a few meters of the fault zone core, the rock is highly fractured but has sustained little net shear. Matrix permeability of this zone is approximately 30 to 60 microdarcy at 50 MPa effective confining pressure (300 to 1000 microdarcy at in situ pressures). Outside this damage zone, matrix permeability drops below 0.01 microdarcy. The clay-rich core material has the lowest strength with a coefficient of friction of approximately 0.55. Shear strength increases with distance from the shear axis. These permeability and strength observations reveal a simple fault zone structure with a relatively weak fine-grained core surrounded by a damage zone of fractured rock. In this case, the damage zone will act as a high-permeability conduit for vertical and horizontal flow in the plane of the fault. The fine-grained core region, however, will impede fluid flow across the fault.

Pure and Applied Geophysics↗

USGS2022

Urban proximity while breeding is not a predictor of perfluoroalkyl substance contamination in the eggs of brown pelicans

Identifying sources of exposure to chemical stressors is difficult when both target organisms and stressors are highly mobile. While previous studies have demonstrated that populations of some organisms proximal to urban centers may display increased burdens of human-created chemicals compared to more distal populations, this relationship may not be universal when applied to organisms and stressors capable of transboundary movements. We examined eggs of brown pelicans ( Pelecanus occidentalis ), a nearshore seabird with daily movements ranging from local to 50 km and annual migrations ranging from year-round residency to 1500 km. Thirty-six eggs from three breeding colonies located at increasing distances to a major urban center (Charleston, South Carolina, USA) were analyzed for concentrations of per - and polyfluoroalkyl substances (PFAS). Areas of high use for each colony during the breeding season were also assessed via the tracking of adult pelicans from each colony using GPS-PTT satellite transmitters and overlapped with measures of relative urbanization via land cover data. We report potentially significant ∑PFAS concentrations in the eggs of pelicans (175.4 ± 120.1 ng/g w wt. SD), driven largely by linear perfluorooctane sulfonate (n-PFOS) (48–546 ng/g w wt.). Residues of the precursor compound perfluorooctane sulfonamide (FOSA) were also present in pelican eggs, suggesting continued exposure of local wildlife beyond implemented phaseouts of some PFAS. For most analytes, egg concentrations did not exhibit a significant spatial structure despite some differentiation in high-use areas unlike similar data for another regional apex predator, the bottlenose dolphin ( Tursiops truncatus ). We suggest that the partially migratory nature of brown pelicans during the non-breeding season, combined with daily ranges that may extend to 50 km from local point sources, may have homogenized exposure across individuals. Charleston likely remains a major source for PFAS in the overall region, however, given the high concentrations observed as well as known releases of PFAS in the nearshore environment.

South Carolina↗

USGS2020

Using advanced population genomics to better understand the relationship between offshore and spawning habitat use for Atlantic Sturgeon

Atlantic Sturgeon ( Acipenser oxyrinchus oxyrinchus ) are a large-bodied anadromous fish that historically supported important fisheries along the east coast of the United States. Following years of overharvest and habitat degradation, populations experienced severe declines. In 2012, the National Marine Fisheries Service listed Atlantic Sturgeon under the Endangered Species Act (ESA; 61 FR 4722). Their listing named five Distinct Population Segments (DPSs), predicated on genetic groups composed of geographically proximate populations. Federal management of Atlantic Sturgeon presents challenges, as sturgeon from each of the five DPSs mix extensively in coastal and marine habitats yet take and recovery progress must be evaluated separately for each unit. Genetic assignment testing based on mitochondrial and microsatellite markers allows individuals to be assigned back to their natal river and DPS. However, this approach is not perfect and some individuals may be incorrectly assigned. Recent advances in genomics offer the potential of a higher resolution approach to genetic assignment testing, and thus may reduce uncertainty associated with assignment testing. In addition, genomics allows a greater number of markers to be examined from across a broader portion of the sturgeon genome, thus may provide an enhanced perspective of population structure for the species, and potentially allow other previously intractable questions to be addressed (Bernatchez et al. 2017, Supple and Shapiro 2018). We used next-generation sequencing to develop a draft genome for Atlantic Sturgeon and identify single nucleotide polymorphisms (SNPs) that could be used to resolve the natal river and DPS of individual Atlantic Sturgeon. We identified 1,210 candidate SNPs within the nuclear genome as well as 49 SNPs within the mitochondrial genome. After filtering and review, we selected 161 nuclear SNPs and 39 mitochondrial SNPs for further testing and evaluation. We used genotyping-in-thousands by sequencing (GT-seq) to simultaneously sequence nuclear SNP loci, mitochondrial SNP loci, and the existing panel of twelve microsatellite loci. This effort required a pilot sequencing run on a single sturgeon sample to test marker amplification and refine primer strengths, followed by a series of sequencing runs to generate baseline data for 288 individuals representing nine populations of Atlantic Sturgeon in four DPSs. Using baseline data from the nine populations, we ran a series of genomic analyses to characterize diversity within and among populations, providing a benchmark for this species using the new SNP markers. Allelic richness was similar for all populations, although there was a general trend of more northern population containing greater levels of allelic richness. Interestingly, we observed linkage disequilibrium among many pairs of loci within many populations. This might be the result of physical linkage but could also suggest these populations are recovering from genetic bottlenecks and/or are effectively small, leading to specific haplotypes to be favored by chance. Pairwise differentiation among populations varied among the populations ( F ST range: 0.010-0.098) and was significantly correlated ( r = 0.771; P < 0.001) to pairwise F ST observed using microsatellite markers). Population clustering and ordination techniques using the new genomic data both support an overall population structure that is similar to the current DPS management units (which were developed primarily based on microsatellite genetic data). Overall, this suggests that existing microsatellite markers and the panel of SNP markers developed in this study provide similar information about the populations structure and ecology of Atlantic Sturgeon. Given the observed differences in allele frequencies among populations, our genomic baseline supports previous assertations that Atlantic Sturgeon show natal homing, despite mixing extensively in marine waters during non-breeding periods. Lower levels of differentiation between populations in the South Atlantic DPS suggest that populations in this region may have greater levels of gene flow relative to their more northerly conspecifics, which has also previously been suggested based on microsatellite data. The observed differentiation among populations provides the necessary foundation for determining the natal river and DPS of Atlantic Sturgeon using assignment testing. We tested the utility of our new genomic baseline for resolving the population and DPS of Atlantic Sturgeon. Our nuclear SNP markers showed utility for identifying the origin of unknown Atlantic Sturgeon samples, as 86.5% were assigned to the correct DPS and 66.3% were assigned to the correct natal river. However, since this study was funded the Conservation Genetics and Genomics Laboratory at Leetown Science Center has made significant improvements to their microsatellite genetic baseline, which now performs more effectively than our new genomic approach (the genetic baseline includes 12 populations and 5 DPSs, and correctly assigns 95.8% of individuals to DPS and 84.9% of individuals to their natal population using 12 microsatellite loci). We conducted an ad hoc exploration of how additional microsatellite or nuclear SNP loci may further improve the accuracy of assignment testing. We found that additional microsatellite markers are likely to result in greater improvements in assignment efficiency than additional nuclear SNPs. However, a much larger number of SNP loci (which if identified could be sequenced using other methods that are now available; e.g., the RAD-capture approach published by Ali et al. 2016) could produce assignment efficiencies that are greater than what is currently feasible using microsatellites. In the absence of further research and development of additional SNP markers for Atlantic Sturgeon (possibly using an approach other than GT-seq), the existing microsatellite loci are the most effective means available to determine the natal river and DPS of Atlantic Sturgeon encountered in offshore waters. Because our new genomic markers were less effective than the existing panel of 12 microsatellite markers, we chose to use the existing microsatellite markers to assign Atlantic Sturgeon captured in another BOEM-funded study (cooperative agreement M16AC00003; Monitoring endangered Atlantic Sturgeon and commercial finfish habitat use offshore New York) following consultation with our project officer. Using this approach, we genotyped and assigned 186 Atlantic Sturgeon captured in coastal waters off the Rockaway Peninsula, New York. The vast majority of these sturgeon were assigned to the New York Bight DPS (94.62%), and most appear to belong to the Hudson River population (87.10%) with smaller contributions from the Delaware River population (7.53%). Smaller contributions (2.15%) were observed from six other populations, including those from the James, York, Kennebec, Ogeechee, and Edisto rivers. Although most of the fish we assigned were assigned to the nearest spawning rivers (Hudson and Delaware), the contributions from distant rivers is consistent with the propensity of this species to move long distances and form mixed stock aggregations along the continental shelf. This finding indicates that spawning populations (and their corresponding DPS) from distant locations may potentially be impacted by offshore activities. In fact, activities in this region of the New York Bight could negatively impact Atlantic Sturgeon population from at least four different DPSs. Genetic or genomic assignment testing remains an essential tool to characterize potential impacts to Atlantic Sturgeon populations and should be applied more broadly to better characterize potential impacts of activities in other locations.

Atlantic Coast↗

USGS2021

Improved wetland soil organic carbon stocks of the conterminous U.S. through data harmonization

Wetland soil stocks are important global repositories of carbon (C) but are difficult to quantify and model due to varying sampling protocols, and geomorphic/spatio-temporal discontinuity. Merging scales of soil-survey spatial extents with wetland-specific point-based data offers an explicit, empirical and updatable improvement for regional and continental scale soil C stock assessments. Agency-collected (U.S. Department of Agriculture, U.S. Environmental Protection Agency) and community-contributed soil datasets were compared for representativeness and bias, with the goal of producing a harmonized national map of wetland soil C stocks with error quantification for wetland areas of the conterminous United States (CONUS) identified by the USGS National Landcover Change Dataset (NLCD). This allowed application of an empirical predictive model of SOC density to be applied across the entire CONUS using relational %OC distribution alone. A broken-stick quantile-regression model identified %OC with its relatively high analytical confidence as a key predictor of SOC density in soil segments; soils less than 6%OC (hereafter, mineral wetland soils, 85% of the dataset) had a strong linear relationship of %OC to SOC density (RMSE = 0.0059, ~4% mean RMSE) and soils greater than 6%OC (organic wetland soils, 15% of the dataset) had virtually no predictive relationship of %OC to SOC density (RMSE = 0.0348 g C cm-3, ~56% mean RMSE). Disaggregation by vegetation type (woody v. emergent herbaceous), or region did not alter the breakpoint significantly (6% OC) nor improve model accuracies for inland and tidal wetlands. Similarly, SOC stocks in tidal wetlands were related to %OC, but without a mappable product for disaggregation to improve accuracy by soil class, region or depth. Our layered, harmonized CONUS wetland soil maps have now revised wetland SOC stock estimates downward by 24% (9.5 vs. 12.5Pg C) with the overestimation being entirely an issue of inland, organic wetland soils, (35% lower than SSURGO-derived SOC stocks). Further, SSURGO underestimated soil carbon stocks at depth, as modeled wetland SOC stocks for organic-rich soils showed significant preservation downcore in the NWCA dataset (<3% loss between 0-30 cm and 30-100 cm depths) in contrast to mineral-rich soils (37% downcore stock loss). Future CONUS wetland soil C assessments will benefit from focused attention on improved organic wetland soil measurements, land history, and spatial representativeness.

Frontiers in Soil Science↗