USGS ScienceSearch

SEARCH · USGS Science

Results for “Marine Ornithology: Journal of Seabird Research and Conservation”

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.

1,688 records · Page 45Linked to original sources

Assessment of western Oregon debris-flow hazards in burned and unburned environments

In the steep and mountainous environment of western Oregon, debris flows pose a considerable threat to property, infrastructure and life. Wildfire is commonly known to increase the susceptibility of steep slopes to debris flows, but the extent of this process in the western Cascades is not well understood. The US Geological Survey (USGS) currently estimates postfire debris-flow likelihood and triggering rainfall thresholds using a model calibrated to a southern California inventory of debris flows generated by excess runoff within the first year after fire. Because of a lack of available data, this model has not been tested in western Oregon, or in locations where postfire debris flows initiate via other mechanisms (e.g., shallow landslides or in-channel failures). Using repeat field observations and aerial imagery, we developed two new debris-flow inventories within and adjacent to the perimeters of five 2020 wildfires in western Oregon: Archie Creek, Holiday Farm, Beachie Creek, Lionshead and Riverside. The first inventory focuses on postfire debris flows (2020–2022); the second focuses on debris flows prior to fires (1995–2020). Our inventories of prefire and postfire debris flows were used to document initiation mechanisms in Oregon's western Cascades and to evaluate the effects of wildfire. We found that wildfire changed the distribution of debris-flow initiation mechanisms in the western Cascades. After the wildfires, annual rates of runoff-generated debris flows increased by 22% and the number of shallow landslide-initiated debris flows decreased by 17% relative to before the wildfires. Despite this shift, shallow landsliding was the dominant debris-flow initiation mechanism in both unburned and burned environments. We found the performance of the current USGS debris-flow likelihood model was degraded relative to other previously tested locations across the intermountain western United States. Our results highlight the need for improved postfire hazard assessment in western Oregon based on regional model calibration that is tuned to the dominant debris-flow initiation mechanisms.

Oregon

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California

In-situ geochronology using LA-ICP-MS/MS: Application of the Lu-Hf system in carbonate, apatite and fluorite

Geochronology is a fundamental tool for reconstructing earth history and constraining the timing of ore deposit formation. Traditionally, U-rich minerals like zircon, titanite, monazite and apatite have been analysed to determine the timing of mineralization and alteration events. However, not all deposits or alteration assemblages contain these minerals, and the U-Pb system can be disturbed due to open system behaviour. We present a methodology for obtaining in-situ Lu-Hf dates from various minerals to demonstrate the potential for dating ore deposits that lack U-rich minerals. Carbonate, apatite and fluorite minerals from different tectonic settings and geological environments were analysed. We acquired Lu-Hf and REE (rare earth element) data from carbonate minerals from a carbonatite orebody and dike samples from the Mountain Pass REE carbonatite deposit in California, USA. We also tested a ‘double-dating’ method, where U-Th-Pb dates and Lu-Hf dates were obtained from the same analysis of apatite from rocks in the Michigan’s Great Lakes Tectonic Zone. Preliminary results for Lu-Hf dating of fluorite from carbonatite rocks in the Lemitar Mountains of New Mexico are also included. The analytical methodology shows promising results in each of these cases where there is sufficient parent Lu and U and measurable daughter Hf and Pb respectively.

Conference Paper

Geologic map of the northwest flank of Mauna Loa volcano, Island of Hawai‘i, Hawaii

Mauna Loa, the largest active volcano on Earth, has erupted 34 times since written descriptions became available in A.D. 1832. The most recent eruption of Mauna Loa occurred on November 27, 2022, after a 38 year hiatus; it lasted for 12 days. Some eruptions began with only brief seismic unrest, whereas others followed several months to a year of increased seismicity. Once underway, Mauna Loa’s eruptions can produce lava flows that may reach the sea in less than 24 hours, severing roads and utilities. For example, lava flows that erupted from the Southwest Rift Zone in 1950 advanced at an average rate of 9.3 kilometers per hour (5.8 miles per hour); all three lobes reached the ocean within ~24 hours. Near the eruptive vents, the flows likely traveled even faster. In terms of eruption frequency, pre-eruption warning, and rapid flow emplacement, Mauna Loa has great volcanic-hazard potential for the Island of Hawai‘i. Volcanic hazards on Mauna Loa can be anticipated, and risk substantially mitigated, by documenting its past activity to refine our knowledge of the hazards, and by alerting the public and local government officials of our findings and their implications for hazards assessments and risk. The map of the north and west flanks of Mauna Loa shows the distribution and relation of volcanic and surficial sedimentary deposits. It incorporates previously reported work published as generalized small-scale maps and a more detailed map. Within the mapped area, lava has flowed from three different source regions: the Northeast Rift Zone (22 percent), the summit (64 percent), and radial vents (14 percent). All three have different points of origin which, in turn, affect the flow characteristics and periodicity of activity. The map area includes the uppermost part of the NERZ and extends from the highest elevation––13,040 feet at the south end of the Kokoolau quadrangle, just below the summit caldera––to the sea northwest and west of the summit. Lava that erupts from the north and west flanks typically flows to the west, northwest, or north, depending on the vent location. Both morphologic lava flow types—‘a‘ā and pāhoehoe—are present. Pāhoehoe units tend to spread out or widen in low-slope regions, such as in the saddle regions between Mauna Loa and Mauna Kea or between Mauna Loa and Hualālai. In comparison, ʻaʻā flows generally produce narrower flow lobes that have higher relief. This map is the fifth in a series of five maps that will cover Mauna Loa volcano. NOTE: Map sheet 1 contains lines and type with overprint. This feature may be turned on or off in the Adobe Acrobat page display preferences.

Hawaii

Geophysical characterization of an alkaline‑carbonatite complex using gravity and magnetic methods at Magnet Cove, Arkansas, USA

The Magnet Cove alkaline‑carbonatite complex (MCC), located in the Ouachita Mountains of south-central Arkansas in the United States, hosts an extensive variety of rare rock types and critical mineral resources with physical properties (density and magnetization) that contrast significantly with the sedimentary rocks into which they have intruded. Newly acquired ground-based gravity and magnetic data were used to develop two-dimensional and three-dimensional geophysical models of the Cretaceous-aged Magnet Cove intrusive complex. The models reveal that the MCC: (1) widens out at middle crustal depths to as much 22 km across, and may reach a depth of 20 km; (2) has a total volume (exposed and subsurface) that may be over 800 km 3 ; (3) is likely connected at depth to other intrusions in the Arkansas alkaline province; and (4) has a geometry that is aligned with pre-existing structures such as the Reelfoot rift and the Ouachita orogenic belt, some of which were likely structurally controlled by the Precambrian crystalline basement and the continent-ocean transition zone buried beneath the Ouachita orogen. For the first time, the magnetic models of the MCC account for the presence of strong remanent magnetization. This results in a geophysical workflow necessary to accurately interpret magnetic anomalies over the much larger Arkansas alkaline province, its geologic and structural framework, and critical mineral potential.

Arkansas

Precariously balanced rocks in northern New York and Vermont, U.S.A.: Ground-motion constraints and implications for fault sources

Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake‐source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground‐based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time‐independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site‐specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground‐motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M ∼6.5 earthquake‐source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake‐source faults in regions of cryptic faults.

Bulletin of the Seismological Society of America

Hydrogeologic framework and conceptual model of the Red River alluvial aquifer east of Lake Texoma, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statutes §82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. At present (2025), the Oklahoma Water Resources Board has not established a maximum annual yield for the Red River alluvial aquifer east of Lake Texoma. To support the evaluation and determination of a maximum annual yield, a hydrogeologic framework and conceptual groundwater-flow model were developed to assess groundwater availability in the Red River alluvial aquifer east of Lake Texoma. The scope of this hydrologic investigation is the alluvium and terrace containing the Red River alluvial aquifer in Oklahoma between Lake Texoma, the Texas State line, and the Arkansas State line, an extent referred to in this report as “the eastern part of the Red River alluvial aquifer.” Parts of the alluvium and terrace extent in Arkansas and Texas are included in some analyses to address hydrologic influences from outside the aquifer’s boundaries in Oklahoma. The eastern part of the Red River alluvial aquifer in southeastern Oklahoma consists of approximately 401,280 acres of Quaternary alluvium and terrace deposits associated with the Red River and its major tributaries. Mean annual recharge to the aquifer for the 1980–2022 study period was estimated to be 8.62 inches per year, or 17.98 percent of the mean annual precipitation over the same period (47.94 inches). This mean annual recharge rate is equivalent to an inflow of approximately 288,250 acre-feet per year for the eastern part of the Red River alluvial aquifer. Recharge estimated using the Soil-Water-Balance code accounts for 98.7 percent of the conceptual-model inflows to the eastern part of the Red River alluvial aquifer. Saturated-zone evapotranspiration accounts for 11.9 percent and net streambed seepage accounts for 87.4 percent of the outflows in the conceptual model.

Arkansas, Oklahoma, Texas

Virginia and Landsat

From the shores of Jamestown and spreading north, south, and west, the lands that became the State of Virginia were some of the first in North America top experience rapid landscape change from European settlement. Imagery and data from the USGS Landsat series of satellites offer an unparalleled resource for the study, understanding, and preservation of Virginia’s land and water resources. From monitoring the health of water bodies to managing invasive species to planning for a range of climate change effects, the USGS National Land Imaging Program’s stewardship and public delivery of Landsat data have benefitted Virginians in myriad ways.

Virginia

Development of vegetation based soil quality indices for mineralized terrane in arid and semi-arid regions

Soil quality indices (SQIs) are often management driven and attempt to describe key relationships between above- and below-ground parameters. In terrestrial systems, indices that were initially developed and modified for agroecosystems have been applied to non-agricultural systems in increasing number. We develop an SQI in arid and semi-arid ecosystems of the Western US impacted by different types of geologic mineralization using the relationship between vegetation community parameters and soil abiotic and biotic properties. We analyze these relations in soils associated with three different mineralization types: podiform chromite, Cu/Mo porphyry, and acid-sulfate gold vein systems at four different sites in California and Nevada. Soil samples were collected from undisturbed soils in both mineralized and nearby unmineralized substrates as well as from waste rock and tailings. Aboveground net primary productivity (ANPP), canopy cover and shrub density were measured for the vegetative communities. Minimum data sets were developed based on correlations between the soil and vegetation parameters, refined using principal components analysis, scored using non-linear functions, and combined into an overall SQI. The indices are comprised of one or two microbial parameters and three to six abiotic parameters, the latter consisting of nutrients and metals. Given the preliminary development of this approach, the parameters and combinations to arrive at an SQI for a given site cannot at this time be correlated or compared with that of another site. This SQI approach provides a means of quantifying disturbed ecosystem recovery resulting from mining, and could be applied to other disturbances in a way that readily distills the information for potential use by land managers. However, severely disturbed areas with little to no aboveground biomass, such as unreclaimed tailings, have likely crossed an ecological threshold that precludes the use of this type of monitoring tool.

California

StreamStats—A quarter century of delivering web-based geospatial and hydrologic information to the public, and lessons learned

StreamStats is a U.S. Geological Survey (USGS) web application that provides streamflow statistics, such as the 1-percent annual exceedance probability peak flow, the mean flow, and the 7-day, 10-year low flow, to the public through a map-based user interface. These statistics are used in many ways, such as in the design of roads, bridges, and other structures; in delineation of floodplains for land-use zoning and setting of insurance rates; for regulatory purposes, such as the permitting of wastewater discharges; and for hydrologic and climate change studies. StreamStats was first developed for Massachusetts and released in 2001. The application provided users with the ability to obtain streamflow statistics computed from data collected at USGS streamgages and to obtain estimates of streamflow statistics for user-selected ungaged sites. Massachusetts StreamStats used geographic information system software and digital mapping to compute drainage-basin characteristics, which were then used in statistical models to estimate streamflow statistics for the user-selected sites. The statistical models were in the form of equations that were developed through a process known as regression analysis. StreamStats was the first known web application with the ability to do interactive geoprocessing. The utility of Massachusetts StreamStats was instantly apparent, leading the USGS to develop a version of StreamStats that could be implemented nationally. USGS State offices normally were required to develop custom regression equations and prepare local digital mapping data needed for implementing StreamStats for their States. Funding needed to complete this work usually was provided through cooperative agreements between the USGS and State agencies. In 2004, Idaho became the first to be released in the national version of StreamStats. By 2023, 44 States were fully implemented and six were undergoing implementation. StreamStats has undergone many modifications over the years to keep up with changes to the underlying software and to add functionality. Customized functionality and separate linked StreamStats applications were developed for several States. Meeting the high demand for additions and improvements to StreamStats while also adhering to budgetary constraints has, at times, been challenging. The StreamStats development team has identified numerous additional improvements that could be made to provide better performance and more functionality. The lessons learned from the experience of building and operating StreamStats for nearly 25 years could be relevant to others interested in pursuing efforts of a similar scale.

Circular

Magnitude and frequency of floods in the Suwannee River Water Management District, Florida

Flood-frequency statistics for 2-, 5-, 10-, 25-, 50-. 100-, 200-, and 500-year recurrence intervals, based on three methods of analysis, are presented for 25 continuous-record and seven peak flow partial-record gaging stations in the Suwannee River Water Management District. The first method, for gaged stations, utilizes station records; the second method, for ungaged sites, utilizes regional regression analysis; and the third method uses a weighted combination of the station and regional values. Because the weighted values utilize two more or less independent estimates of the peak flow statistic, they are considered more accurate than the station estimates or the regression estimates alone. Also, the use of another weighting scheme to improve estimates of flood frequency statistics at ungaged sites is demonstrated. The karstic nature of much of the Suwannee River Water Management District significantly attenuates flood peaks in some streams by providing substantial subsurface storage when river stages are high. At such times, springs discharging into rivers may reverse flow temporarily and become sinks.

Florida

Biogeochemical and community ecology responses to the wetting of non-perennial streams

Transitions between dry and wet hydrologic states are the defining characteristic of non-perennial rivers and streams, which constitute the majority of the global river network. Although past work has focused on stream drying characteristics, there has been less focus on how hydrology, ecology and biogeochemistry respond and interact during stream wetting. Wetting mechanisms are highly variable and can range from dramatic floods and debris flows to gradual saturation by upwelling groundwater. This variation in wetting affects ecological and biogeochemical functions, including nutrient processing, sediment transport and the assembly of biotic communities. Here we synthesize evidence describing the hydrological mechanisms underpinning different types of wetting regimes, the associated biogeochemical and organismal responses, and the potential scientific and management implications for downstream ecosystems. This combined multidisciplinary understanding of wetting dynamics in non-perennial streams will be key to predicting and managing for the effects of climate change on non-perennial ecosystems.

Nature Water

Using a time-of-travel sampling approach to quantify per- and polyfluoroalkyl substances (PFAS) stream loading and source inputs in a mixed-source, urban catchment

Understanding per- and polyfluoroalkyl substances (PFAS) mass distribution in surface and groundwater systems can support source prioritization, load reduction, and water management. Thirteen sites within an urban catchment were sampled utilizing a time-of-travel sampling approach to minimize the influence of subdaily fluctuations in mass from PFAS point sources and to quantify PFAS and ancillary chemical loads from various PFAS sources. A larger increase in perfluoroalkyl sulfonate (PFSA) loads (8 to 11 μg/s, up to 618%) than in perfluoroalkyl carboxylate (PFCA) loads (no change to 3.4 μg/s, up to 122%) was observed at sites below tributaries influenced by military bases with known groundwater discharge. Point discharges from two sewage treatment plants (STPs) resulted in increases in PFCA and PFSA loads that were similar (6 and 10 μg/s respectively) below the first STP and greater for PFCA compared to PFSA loads (23 and 13 μg/s respectively) below the second STP. Overall, percent increases in total PFAS load ranged from 20 to 277% for military base inputs and 44 to 77% for STP inputs. A focus catchment that represents only 14% (76.9 km 2 ) of the drainage area at the most downstream site (544 km 2 ) accounted for about 70% of PFSA and 40% of PFCA loads observed at the most downstream site. Results show that by using a time-of-travel sampling approach in mixed, urban settings with several PFAS sources, it is possible to quantify stream loads from individual PFAS sources, thereby improving source attribution and providing actionable data for water-resource managers.

ACS ES&T Water

Hydrogeologic characterization of the Cahuilla Valley and Terwilliger Valley Groundwater Basins, Riverside County, California

The relation between the groundwater and the amount of natural recharge to the Cahuilla Valley and Terwilliger Valley groundwater basins is not well understood. During the 20th century, the reliance on groundwater near Anza, California, used for agricultural, domestic, and municipal reasons has increased, and there is the potential for changes in groundwater availability related to climate change. Several types of existing data were evaluated, and new data were collected for this study, with the goal of characterizing the region’s hydrogeology. The study’s scope included constructing a geologic framework model to show where the groundwater-bearing units are present and their relation to each other, estimating the major components of the groundwater budget, and understanding local short-term and regional long-term groundwater flow and how that has changed since the early 1900s. Two electrical resistivity tomography surveys were done in the Durasno Valley about 2,150 feet apart to identify the thickness of the alluvium, its horizontal extent, and the depth-to-basement along two profiles perpendicular to Cahuilla Creek. The subsurface sediments were mostly horizontally layered and the transitional boundary between the alluvium and basement was thinner and shallower along the upgradient profile where the depth-to-basement was about 70 feet below land surface; the depth-to-basement at the downgradient profile was more than about 140 feet below land surface. The results from the surveys were used to place four monitoring wells at two sites along the survey profiles. Artesian flow from the deepest well at the downgradient site indicated that the decomposed and competent basement likely contributed some groundwater to the overlying alluvium, laterally, from below, or both. A digital three-dimensional geologic framework model was constructed using EarthVision software to represent the subsurface geometry of the alluvium, decomposed basement, and competent basement. Maps and cross sections of the modeled thicknesses of the alluvium and decomposed basement, and the modeled elevation of the top of the competent basement, were made to show the subsurface geometry of vertical faults, selected wells, and the groundwater-bearing units. Because natural recharge is related to the variable cycles of precipitation, estimates are difficult to quantify. Recharge and runoff have extreme interannual variability in the study area; recharge and runoff can be sporadic, and a substantive amount may not occur in some years. Estimates of recharge from a previous study and the regional-scale Basin Characterization Model for California for four different periods ranged from 3,800 acre-feet/year for 1897–1947 to 5,900 acre-feet/year for 1971–2000. Potential recharge from the disposal of domestic septic systems may have been as much as 500 acre-feet in 2020. It was estimated that between about 400 and 2,400 acre-feet/year of groundwater is lost through evapotranspiration by vegetation and evaporation from open water bodies, but the main source of discharge is through pumpage, mainly used for agriculture from the alluvium in the Cahuilla Valley and Terwilliger Valley groundwater basins. The estimated total pumpage for 1991–2021 ranged from about 1,140 acre-feet in 2019 to about 3,450 acre-feet in 1994. When summed, the cumulative amount of estimated pumpage between 1991 and 2021 was about 81,400 acre-feet. The general direction of groundwater flow is from the northeast along the San Jacinto fault zone at the headwaters of Cahuilla and Hamilton Creeks, to the surface-water outlets at the west and southeast parts of the study area. Groundwater-level data from the 1950s and earlier indicate that there was a natural groundwater divide between the Cahuilla Valley and Terwilliger Valley groundwater basins, but the changing magnitude and extent of the groundwater depressions caused by pumping since about 1950 indicate that the location of the natural groundwater boundary between the Cahuilla Valley and Terwilliger Valley groundwater basins has migrated over time. Flow from the upper to the lower parts of the Cahuilla Valley groundwater basin roughly follows the course of Cahuilla Creek through the narrow Durasno Valley where an estimated volume of flow in April 2019 was about 10–150 acre-feet/year. Short-term trends in groundwater levels, particularly in wells where groundwater is shallow and in the basement unit, show how some areas respond quickly to recharge and discharge. Wells located further to the east within the Cahuilla Valley groundwater basin in the alluvium show much less of a response to recharge events; areas of sustained pumpage from the alluvium, primarily for agriculture, show long-term declines in groundwater levels and generally do not show the effects of storm events or recent runoff. Groundwater levels in wells that are farthest from where most of the recharge occurs and where pumping has been the greatest, had some of the largest long-term groundwater-level declines at a rate of about 0.8 foot/year between 1971 and 2021.

California

Seismic velocity changes from repetitive seismicity at Mauna Loa prior to and during its 2022 eruption

Mauna Loa’s short-lived eruption from late November to early December 2022 marked the culmination of nearly a decade of elevated seismic activity and geodetic inflation. The volcano has been monitored by a network of permanent, short period and broadband seismometers. I used the continuous waveform data from that network starting in 2012 to generate a catalog of seismicity that enhances the US Geological Survey Hawaiian Volcano Observatory’s public seismic catalog with four times the number of earthquakes, which were then grouped by waveform similarity. Analysis of subtle delays in the timing of arrivals of scattered waves between pairs of earthquakes in this catalog yields a history of small changes in the shallow seismic velocity structure of the volcano. Seismic velocities have been shown at other volcanoes to change during unrest and eruption. My results show a decrease in seismic velocity centered on the summit beginning in September 2022, corresponding to the onset of a vigorous precursory swarm of seismic activity and shallow inflation. During the eruption itself, I observe large changes due likely to dike opening along the northeast rift zone and deflation of the summit reservoir. However, seismic velocity changes associated with non-volcanic sources such as ground shaking from large earthquakes and meteorological influences at seasonal and diurnal time scales are also observed, and these dominate the velocity changes prior to the eruption. Proper accounting of these effects will be a requirement for use in real-time monitoring, and this work serves as a starting point in that endeavor for Mauna Loa.

Hawaii

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, situated in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the impact of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study we constructed the Strawberry Creek integrated hydrological model (SCIHM) to increase this understanding. The SCIHM is an integrated surface runoff and groundwater model that uses the coupled groundwater and surface-water flow model (GSFLOW), which is based on the integration of the precipitation-runoff modeling system (PRMS) and the modular groundwater flow model commonly called MODFLOW, version MODFLOW-2005 software to simulate surface runoff and infiltration and groundwater flow. The model has three layers, 263 rows, and 176 columns. The model area includes the Strawberry Creek and four adjacent watersheds. The PRMS model was calibrated using two streamflow gaging stations and the GSFLOW model was calibrated to reported spring diversion discharge and a sparse number of groundwater-level measurements. The SCIHM was run with and without diversions active and simulated streamflow was compared, finding that in the headwaters of Strawberry Creek about 35 percent of the diversion flow was captured from base flow.

EarthArXiv

Spatial variation in effectiveness and costs of sagebrush restoration treatments across the western United States

Natural resource managers strive to improve restoration efficiency across the vast network of public lands of the United States, seeking the best return on investment through innovation, geographic prioritization, and adaptive management. This challenge is exemplified by ecosystems dominated by sagebrush ( Artemisia spp.), which occur across the western United States and have been degraded by increasing wildfire frequency driven by invasive plant spread, among other factors. However, both costs and effectiveness of sagebrush restoration can vary spatially due to management practices in addition to biotic and abiotic factors, and characterizing this variation could inform broad-scale planning. We leveraged published models of sagebrush recovery and treatment costs to predict and evaluate cost-effectiveness of aerial and drill seeding Artemisia spp. across 429 718 km 2 that overlapped 12 US states. Compared to natural recovery, effectiveness for both methods was generally low, yet effectiveness was often greater, and costs were always greater, for drill than aerial seeding, resulting in slightly higher cost-effectiveness on average for drill seeding. Cost-effectiveness for both treatments increased near major roads and in areas with more repeated burns. Cost-effectiveness also increased with predisturbance cover of sagebrush but decreased with herbaceous cover. However, we also identified areas where aerial seeding was more cost-effective than drill seeding, particularly with greater slope. Our results provide spatially explicit estimates and potential mechanisms of cost-effectiveness for two common seeding methods for sagebrush, which can help prioritize limited resources, guide land use, and improve restoration efficiency and effectiveness across public rangelands of the western United States.

western United States

Factors affecting the distribution of water-bearing fractures in the bedrock aquifers of West Virginia

Bedrock aquifers cover 23,601 square miles within the State of West Virginia and comprise 97.4 percent of the surficial area within the State; the remaining 2.6 percent (621 square miles) consists of alluvial sand-and-gravel and glacial outwash aquifers bordering the State’s major rivers. While West Virginia’s alluvial aquifers have been studied extensively, bedrock aquifers have only been characterized for studies completed in a few areas in Jefferson, McDowell, and Monroe Counties. Bedrock aquifers are water supplies for public supply, agriculture, industry, and residential homeowner use. In this study, the U.S. Geological Survey, in cooperation with the West Virginia Department of Environmental Protection Division of Water and Waste Management, provides a statewide assessment of the occurrence and distribution of fractures within bedrock aquifers of the State and the various topographic, physiographic, and lithologic influences controlling the occurrence and distribution of bedrock fractures. The results of this study provide an increased understanding of the distribution of fractures in bedrock aquifers in West Virginia and help to verify trends that have been suspected for many years but were never well documented or verified by data. The results confirmed that the density of fractures and those that were determined to be water bearing decrease significantly with depth. A statistically significant difference in the density of fractures was observed at a depth of 215 feet for wells in the Appalachian Plateaus Physiographic Province’s and in the Valley and Ridge Physiographic Province’s aquifers; a higher density of fractures and water-bearing fractures were above a depth of 215 feet than below that depth. This is an important consideration when drilling wells for residential, commercial, industrial, or agricultural water supply. Abandoned underground coal mines are commonly believed to form large pools of water in the interconnected mine entries in abandoned room and pillar coal mines. Such pools of water can and do exist in abandoned underground coal mines, but many mines lack open entries and are held up by overburden strata and pillars that can collapse and form aquifers comprised of vast interconnected rubble zones (gob), especially in older mines. Data assessed for this study showed that shale-corrected values of effective porosity for limestone aquifers in West Virginia had a median value of 2 percent and an average value of 4 percent and generally are mineralized with low effective porosity. Argillaceous or sandy limestone has a median shale-corrected porosity of 4 percent and an average shale-corrected porosity of 5 percent. The median and average shale-corrected porosity of sandstone aquifers was estimated to be 14 percent, but the median shale-corrected porosity for argillaceous or calcareous sandstone was 5 percent and the average shale-corrected porosity for argillaceous or calcareous sandstone was 6 percent. Even though shale has a relatively high total sonic porosity compared to other lithologies, shale and siltstone had relatively low shale-corrected porosity, ranging from 0 to 2 percent. Well yields were previously documented to be highest in valley settings, lowest on hilltops, and intermediate on hillsides. Transmissivity data provided by this study confirm this general pattern within the Appalachian Plateaus Province; however, the Valley and Ridge Province does not follow this pattern. While still lowest on hilltop settings, the highest well yields were in hillside settings. The trend for the Valley and Ridge Province was likely skewed because of 9 high-yield wells specifically targeting deeper thin limestone units, such as the Tonoloway and Helderberg Limestones, at depths with transmissivity in excess of 2,000 feet squared per day in Mineral County, West Virginia, or targeting karst aquifers in Berkeley, Jefferson, or Greenbrier Counties, West Virginia. Finally, water-bearing fractures have been hypothesized to comprise a small number of all fractures within a typical bedrock well in West Virginia. Data collected for this study support this theory. A total of 3,403 fractures were identified during this study; 3,151 (92.6 percent) of those fractures are low-transmissive, and only 252 (7.4 percent) fractures are water-bearing. Even though a well may contain many fractures, less than 8 percent are considered water-bearing fractures.

West Virginia