USGS Science⌕ Search

SEARCH · USGS Science

Results for “Soil Systems”

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

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

At least 1,693 records · Page 94Linked to original sources

Simulation of rainfall-runoff response in mined and unmined watersheds in coal areas of West Virginia

Meteorologic and hydrologic data from five small watersheds in the coal areas of West Virginia were used to calibrate and test the U.S. Geological Survey Precipitation-Runoff Modeling System for simulating streamflow under various climatic and land-use conditions. Three of the basins--Horsecamp Run, Gilmer Run, and Collison Creek--are primarily forested and relatively undisturbed. The remaining basins--Drawdy Creek and Brier Creek-are extensively mined, both surface and underground above stream drainage level. Low-flow measurements at numerous synoptic sites in the mined basins indicate that coal mining has substantially altered the hydrologic system of each basin. The effects of mining on streamflow that were identified are (1) reduced base flow in stream segments underlain by underground mines, (2) increased base flow in streams that are downdip and stratigraphically below the elevation of the mined coal beds, and (3) interbasin transfer of ground water through underground mines. These changes probably reflect increased permeability of surface rocks caused by subsidence fractures associated with collapsed underground mines in the basin. Such fractures would increase downward percolation of precipitation, surface and subsurface flow, and ground-water flow to deeper rocks or to underground mine workings. Model simulations of the water budgets for the unmined basins during the 1972-73 water years indicate that total annual runoff averaged 60 percent of average annual precipitation; annual evapotranspiration losses averaged 40 percent of average annual precipitation. Of the total annual runoff, approximately 91 percent was surface and subsurface runoff and 9 percent was groundwater discharge. Changes in storage in the soil zone and in the subsurface and ground-water reservoirs in the basins were negligible. In contrast, water-budget simulations for the mined basins indicate significant differences in annual recharge and in total annual runoff. Model simulations of the water budget for Drawdy Creek basin indicate that total annual runoff during 1972-73 averaged only 43 percent of average annual precipitation--the lowest of all study basins; annual evapotranspiration losses averaged 49 percent, and interbasin transfer of ground-water losses averaged about 8 percent. Of the total annual runoff, approximately 74 percent was surface and subsurface flow and 26 percent was ground-water discharge. The low total annual runoff at Drawdy Creek probably reflects increased recharge of precipitation and surface and subsurface flow losses to ground water. Most of the increase in ground-water storage is, in turn, lost to a ground-water sink--namely, interbasin transfer of ground water by gravity drainage and (or) mine pumpage from underground mines that extend to adjacent basins. Hypothetical mining situations were posed for model analysis to determine the effects of increased mining on streamflow in the mined basins. Results of model simulations indicate that streamflow characteristics, the water budget, and the seasonal distribution of streamflow would be significantly modified in response to an increase in mining in the basins. Simulations indicate that (1) total annual runoff in the basins would decrease because of increased surface- and subsurface-flow losses and increased recharge of precipitation to ground water (these losses would tend to reduce medium to high flows mainly during winter and spring when losses would be greatest), (2) extreme high flows in response to intense rainstorms would be negligibly affected, regardless of the magnitude of mining in the basins, (3) ground-water discharge also would decrease during winter and spring, but the amount and duration of low flows during summer and fall would substantially increase in response to increased ground-water storage in rocks and in underground mines, and (4) the increase in ground-water storage in the basins would be depleted, mostly by increased losses to a grou

Water Supply Paper↗

Conceptual model of sediment processes in the upper Yuba River watershed, Sierra Nevada, CA

This study examines the development of a conceptual model of sediment processes in the upper Yuba River watershed; and we hypothesize how components of the conceptual model may be spatially distributed using a geographical information system (GIS). The conceptual model illustrates key processes controlling sediment dynamics in the upper Yuba River watershed and was tested and revised using field measurements, aerial photography, and low elevation videography. Field reconnaissance included mass wasting and channel storage inventories, assessment of annual channel change in upland tributaries, and evaluation of the relative importance of sediment sources and transport processes. Hillslope erosion rates throughout the study area are relatively low when compared to more rapidly eroding landscapes such as the Pacific Northwest and notable hillslope sediment sources include highly erodible andesitic mudflows, serpentinized ultramafics, and unvegetated hydraulic mine pits. Mass wasting dominates surface erosion on the hillslopes; however, erosion of stored channel sediment is the primary contributor to annual sediment yield. We used GIS to spatially distribute the components of the conceptual model and created hillslope erosion potential and channel storage models. The GIS models exemplify the conceptual model in that landscapes with low potential evapotranspiration, sparse vegetation, steep slopes, erodible geology and soils, and high road densities display the greatest hillslope erosion potential and channel storage increases with increasing stream order. In-channel storage in upland tributaries impacted by hydraulic mining is an exception. Reworking of stored hydraulic mining sediment in low-order tributaries continues to elevate upper Yuba River sediment yields. Finally, we propose that spatially distributing the components of a conceptual model in a GIS framework provides a guide for developing more detailed sediment budgets or numerical models making it an inexpensive way to develop a roadmap for understanding sediment dynamics at a watershed scale.

Geomorphology↗

Sensitivity of high-elevation streams in the Southern Blue Ridge Province to acidic deposition

The Southern Blue Ridge Province, which encompasses parts of northern Georgia, eastern Tennessee, and western North Carolina, has been predicted to be sensitive to impacts from acidic deposition, owing to the chemical composition of the bedrock geology and soils. This study confirms the predicted potential sensitivity, quantifies the level of total alkalinity and describes the chemical characteristics of 30 headwater streams of this area. Water chemistry was measured five times between April 1983 and June 1984 at first and third order reaches of each stream during baseflow conditions. Sensitivity based on total alkalinity and the Calcite Saturation Index indicates that the headwater streams of the Province are vulnerable to acidification. Total alkalinity and p11 were generally higher in third order reaches (mean, 72 ?eq/? and 6.7) than in first order reaches (64 ?eq/? and 6.4). Ionic concentrations were low, averaging 310 and 340 ?eq/? in first and third order reaches, respectively. A single sampling appears adequate for evaluating sensitivity based on total alkalinity, but large temporal variability requires multiple sampling for the detection of changes in pH and alkalinity over time. Monitoring of stream water should continue in order to detect any subtle effects of acidic deposition on these unique resource systems.

Journal of the American Water Resources Associatio↗

Porewater chemistry of Louisiana marshes with contrasting salinities and its implications for coastal acidification

Dissolved inorganic carbon (DIC) and total alkalinity (TA) are fundamental components of carbonate systems that control pH and buffering capacity of the receiving water body. Three coastal marshes with contrasting salinities in Barataria Basin, Louisiana, USA, were sampled to understand seasonal changes in porewater carbonate chemistry and its impact on surrounding water bodies. Each marsh was sampled five times between December 2018 and October 2019. Porewater DIC and TA increased with depth irrespective of marsh type and ranged from 4.47 to 31.61 mmol/kg and from 1.78 to 28.56 mmol/kg, respectively. The salt marsh had higher porewater DIC and TA compared to the lower salinity intermediate and brackish marshes, probably due to sulfate reduction in the salt marsh. However, it is likely that denitrification is the dominant anaerobic process in these marshes because of low porewater TA/DIC ratios in all three marshes. Porewater TA and DIC concentrations were generally higher during warmer months than colder months. However, the marsh flooding regime had a profound influence on TA and DIC concentrations by changing the redox potential of the marsh soil. Porewater TA/DIC ratios in all three marshes were generally less than 1, while surface water TA/DIC ratios were around 1, suggesting that export of DIC and TA from coastal marshes have the potential to contribute to coastal acidification.

Louisiana↗

Digital method for regional mapping of surficial basin deposits in arid regions, example from central Death Valley, Inyo County, California

Derivative maps generated from DEM's and panchromatic remote sensing data (TM Landsat 7 or SPOT) can be used to characterize surficial basin deposits in arid regions dominated by basin and range topography. Results indicate the technique is useful for rapid digital mapping of surficial deposits where a first order, systematic subdivision of bedrock, alluvial fan units, and playas is unavailable at regional scales. Digital mapping can provide information about relative age and material properties of units that in part can be derived from the position of units within the basin. This automated mapping, implemented in a GIS system, involves an iterative process applied to a combination of digital elevation models (DEM) and satellite image data, such as SPOT or the high-resolution panchromatic Band 8 of Landsat 7 scenes. The method first discriminates the region into first-order terrains consisting of bedrock mountain highlands, basin piedmonts, and playa-basin interiors based on user-defined slope cutoffs applied to DEM data. The basin areas are subsequently classified into surficial map units such as active channels, ground-water discharge zones, and multiple age alluvial-fan piedmont units based on reflective properties of the associated surfaces in the satellite imagery. The surficial units are differentiated through systematic classification based on specific user-defined ranges of spectral values for each unit. The spectral ranges used in the classification are largely dependent on the composite effects of surface characteristics and material properties, including depositional morphology and texture, pavement development, degree of surface clast varnishing, and (or) properties of exposed soils of the alluvial fan units. We have used the slope-curvature properties derived from the DEM data to discriminate the bajada areas that exhibit non-unique spectral characteristics. Slope curvature is particularly effective at differentiating young undissected surfaces from older dissected piedmont units. Available geologic maps and field observations may be used both to iteratively calibrate the spectral classification scheme and to provide additional verification of the digital map output. Digital mapping combined with detailed field studies in selected areas provides useful regional maps of surficial units until time and funding is available for more field intensive studies. In addition, anomalous areas on the thematic maps indicate where more detailed field or air photo work is warranted. The technique successfully distinguishes between bedrock, alluvial fans (generally multiple fan units), active washes, playas, playa rimming marshes and seeps and other active and inactive discharge zones in arid basin and mountain regions. Limitations occur in the subdivision of some fan units where the dominant detrital clast lithologies are not susceptible to varnish development.

California↗

Adaptive management framework and decision support tool for invasive annual bromes in seven Northern Great Plains National Park Service units

National Park Service (NPS) units in the northern Great Plains (NGP) were established to preserve and interpret the history of the United States, protect and showcase unusual geology and paleontology, and provide a home for vanishing large wildlife. A unifying feature among these national parks, monuments, and historic sites is northern mixed-grass prairie, which not only provides background scenery and habitat but is the foundation of many park missions. As recognition of the prairie’s importance to park fundamental resources and values has grown, so too has the realization that invasive plants threaten these values by reducing native species diversity, altering food webs, and marring the visitor experience. Cheatgrass ( Bromus tectorum ) and Japanese brome ( Bromus japonicus )—collectively referred to as “annual bromes”—are of particular concern because of their documented increase through time, and their association with lower native plant diversity, in NGP parks. A variety of grazing, herbicide-application, and prescribed-fire experiments have shown promising short-term results for controlling annual bromes in research-scale plots in the NGP, but it is unclear whether these management actions will be as effective at the larger spatial and longer temporal scales relevant to park management. When uncertainties about the effectiveness of different management actions cannot be answered with traditional research approaches in time to prevent resource degradation, yet recurrent management decisions must be made, an adaptive management approach may be appropriate. Thus, in 2017, we began to develop the ABAM—Annual Brome Adaptive Management—framework. The aim of this framework is to reduce uncertainties about methods for controlling annual bromes in seven NGP parks through a formal process of learning from the application of on-going management. A uniform framework across seven parks provides greater opportunities for reducing these uncertainties compared to a single park acting alone or to multiple parks using different adaptive management frameworks. This technical report details the development and expected implementation of the ABAM framework. After briefly introducing the issue (Section 1) and describing the context in which the framework was developed (Section 2), the report describes how a structured decision-making process was used to frame the problem, determine concrete objectives, and decide the alternative actions for achieving those objectives that the framework would be designed around (Section 3). Then the report describes the process used to develop the ABAM decision support system (Section 4). At the core of this system is the ABAM decision support tool, a Bayesian decision network built on nearly two decades of vegetation monitoring data from NGP parks, as well as current literature and ABAMspecific experiments. This tool, referred to as the ABAM model by its intended users, is built to work with the existing vegetation monitoring, prescribed fire, and invasive plant management programs that support the seven ABAM parks. In Section 5, the report describes how output from the ABAM model is produced and used in annual vegetation management decision making. It describes the ABAM R package (Baldwin et al. 2021) and an example R script that leads a user through an annual workflow using the model and the package. This workflow updates the model with information from new monitoring events following management actions of prescribed fire, herbicide application, or a combination thereof. With the updated model, data describing the current condition of vegetation in park management units, and current data for environmental factors included in the model (soil texture, slope, weather, and grazing), the user then runs the model to predict future vegetation conditions—and managers’ happiness with the outcome—in response to each of 10 management actions for each management unit in each park. These predictions inform managers’ decisions regarding locations and types of management actions to apply in the upcoming year. The ABAM framework is in its infancy, and the report concludes (Section 6) with a discussion of its longer-term viability. Successful adaptive management requires commitment for the long term, likely decades. Currently, the predictions of the decision support tool are not expected to be highly accurate, but they ideally will improve over time as more management actions are applied and their outcomes are captured by monitoring. We designed the ABAM decision support tool to work with the existing management and monitoring resources in ABAM parks to maximize the sustainability of the model’s use, but the ABAM framework requires more than the model. Because this application of an adaptive management framework supported by a quantitative decision support tool to guide vegetation management is unique within the NPS (to our knowledge), institutional knowledge and mechanisms for long-term implementation of the ABAM framework do not exist within the agency. Additionally, the ABAM model and the data that inform it could be improved in a variety of ways. Thus, this report concludes with a discussion of ways to both sustain and improve upon the work completed so far.

Montana, Nebraska, South Dakota, Wyoming↗

Vegetation of wetlands of the prairie pothole region

Five themes dominate the literature dealing with the vegetation of palustrine and lacustrine wetlands of the prairie pothole region: environmental conditions (water or moisture regime, salinity), agricultural disturbances (draining, grazing, burning, sedimentation, etc.), vegetation dynamics, zonation patterns, and classification of the wetlands.The flora of a prairie wetland is a function of its water regime, salinity, and disturbance by man. Within a pothole, water depth and duration determines distribution of species. In potholes deep enough to have standing water even during droughts, the central zone will be dominated by submersed species (open water). In wetlands that go dry during periods of drought or annually, the central zone will be dominated by either tall emergent species (deep marsh) or midheight emergents (shallow marsh), respectively. Potholes that are only flooded briefly in the spring are dominated by grasses, sedges, and forbs (wet meadow). Within a pothole, the depth of standing water in the deepest, usually central, part of the basin determines how many zones will be present. Lists of species associated with different water regimes and salinity levels are presented.Disturbances due to agricultural activities have impacted wetlands throughout the region. Drainage has eliminated many potholes, particularly in the southern and eastern parts of the region. Grazing, mowing, and burning have altered the composition of pothole vegetation. The composition of different vegetation types impacted by grazing, haying, and cultivation is presented in a series of tables. Indirect impacts of agriculture (increased sediment, nutrient, and pesticide inputs) are widespread over the region, but their impacts on the vegetation have never been studied.Because of the periodic droughts and wet periods, many palustrine and lacustrine wetlands undergo vegetation cycles associated with water-level changes produced by these wet-dry cycles. Periods of above normal precipitation can raise water levels high enough to drown out emergent vegetation or produce 'eat outs' due to increases in the size of muskrat populations that accompany periods of high water. The elimination of emergents creates a lake marsh dominated by submersed vegetation. During the next drought when the marsh bottom is exposed by receding water levels (a drawdown), seeds of emergents and mudflat annuals in the soil (the seed bank) germinate (the dry marsh stage). When the marsh refloods, ending the dry marsh stage, the emergents survive and spread vegetatively. This is the regenerating marsh. This stage continues until high water again eliminates the emergents, starting the next degenerating stage.Zonation patterns are conspicuous because each zone often is dominated by a single species that has a lifeform different from those in adjacent zones. The species composition of each zone is a function of its environment (water or moisture regime, salinity, and disturbance history). Within a zone it may take a year or more for species composition to adjust to a change of environmental conditions. These lags sometimes result in abnormal zonation patterns, particularly after a change in water level.Classification of prairie wetlands is more difficult than for most other wetland type, because of these vegetation cycles. Early attempts to classify prairie wetlands did not take the dynamic nature of their vegetation into account. Stewart and Kantrud (1971) developed a classification system for prairie potholes that recognized different phases of vegetation zones dominated by deep marsh species. It used the composition of the vegetation in the deepest part (zone) of a pothole as an indicator of its water-level regime and water chemistry. The application of the national wetland classification system of Cowardin et al. (1979) to potholes is also discussed, and lists of species that characterize the various dominance types associated with the subclasses in this system are presented.

Book chapter↗

Liquefaction-induced lateral spreading in Oceano, California, during the 2003 San Simeon Earthquake

The December 22, 2003, San Simeon, California, (M6.5) earthquake caused damage to houses, road surfaces, and underground utilities in Oceano, California. The community of Oceano is approximately 50 miles (80 km) from the earthquake epicenter. Damage at this distance from a M6.5 earthquake is unusual. To understand the causes of this damage, the U.S. Geological Survey conducted extensive subsurface exploration and monitoring of aftershocks in the months after the earthquake. The investigation included 37 seismic cone penetration tests, 5 soil borings, and aftershock monitoring from January 28 to March 7, 2004. The USGS investigation identified two earthquake hazards in Oceano that explain the San Simeon earthquake damage?site amplification and liquefaction. Site amplification is a phenomenon observed in many earthquakes where the strength of the shaking increases abnormally in areas where the seismic-wave velocity of shallow geologic layers is low. As a result, earthquake shaking is felt more strongly than in surrounding areas without similar geologic conditions. Site amplification in Oceano is indicated by the physical properties of the geologic layers beneath Oceano and was confirmed by monitoring aftershocks. Liquefaction, which is also commonly observed during earthquakes, is a phenomenon where saturated sands lose their strength during an earthquake and become fluid-like and mobile. As a result, the ground may undergo large permanent displacements that can damage underground utilities and well-built surface structures. The type of displacement of major concern associated with liquefaction is lateral spreading because it involves displacement of large blocks of ground down gentle slopes or towards stream channels. The USGS investigation indicates that the shallow geologic units beneath Oceano are very susceptible to liquefaction. They include young sand dunes and clean sandy artificial fill that was used to bury and convert marshes into developable lots. Most of the 2003 damage was caused by lateral spreading in two separate areas, one near Norswing Drive and the other near Juanita Avenue. The areas coincided with areas with the highest liquefaction potential found in Oceano. Areas with site amplification conditions similar to those in Oceano are particularly vulnerable to earthquakes. Site amplification may cause shaking from distant earthquakes, which normally would not cause damage, to increase locally to damaging levels. The vulnerability in Oceano is compounded by the widespread distribution of highly liquefiable soils that will reliquefy when ground shaking is amplified as it was during the San Simeon earthquake. The experience in Oceano can be expected to repeat because the region has many active faults capable of generating large earthquakes. In addition, liquefaction and lateral spreading will be more extensive for moderate-size earthquakes that are closer to Oceano than was the 2003 San Simeon earthquake. Site amplification and liquefaction can be mitigated. Shaking is typically mitigated in California by adopting and enforcing up-to-date building codes. Although not a guarantee of safety, application of these codes ensures that the best practice is used in construction. Building codes, however, do not always require the upgrading of older structures to new code requirements. Consequently, many older structures may not be as resistant to earthquake shaking as new ones. For older structures, retrofitting is required to bring them up to code. Seismic provisions in codes also generally do not apply to nonstructural elements such as drywall, heating systems, and shelving. Frequently, nonstructural damage dominates the earthquake loss. Mitigation of potential liquefaction in Oceano presently is voluntary for existing buildings, but required by San Luis Obispo County for new construction. Multiple mitigation procedures are available to individual property owners. These procedures typically involve either

Open-File Report↗

The use of process models to inform and improve statistical models of nitrate occurrence, Great Miami River Basin, southwestern Ohio

Statistical models of nitrate occurrence in the glacial aquifer system of the northern United States, developed by the U.S. Geological Survey, use observed relations between nitrate concentrations and sets of explanatory variables—representing well-construction, environmental, and source characteristics— to predict the probability that nitrate, as nitrogen, will exceed a threshold concentration. However, the models do not explicitly account for the processes that control the transport of nitrogen from surface sources to a pumped well and use area-weighted mean spatial variables computed from within a circular buffer around the well as a simplified source-area conceptualization. The use of models that explicitly represent physical-transport processes can inform and, potentially, improve these statistical models. Specifically, groundwater-flow models simulate advective transport—predominant in many surficial aquifers— and can contribute to the refinement of the statistical models by (1) providing for improved, physically based representations of a source area to a well, and (2) allowing for more detailed estimates of environmental variables. A source area to a well, known as a contributing recharge area, represents the area at the water table that contributes recharge to a pumped well; a well pumped at a volumetric rate equal to the amount of recharge through a circular buffer will result in a contributing recharge area that is the same size as the buffer but has a shape that is a function of the hydrologic setting. These volume-equivalent contributing recharge areas will approximate circular buffers in areas of relatively flat hydraulic gradients, such as near groundwater divides, but in areas with steep hydraulic gradients will be elongated in the upgradient direction and agree less with the corresponding circular buffers. The degree to which process-model-estimated contributing recharge areas, which simulate advective transport and therefore account for local hydrologic settings, would inform and improve the development of statistical models can be implicitly estimated by evaluating the differences between explanatory variables estimated from the contributing recharge areas and the circular buffers used to develop existing statistical models. The larger the difference in estimated variables, the more likely that statistical models would be changed, and presumably improved, if explanatory variables estimated from contributing recharge areas were used in model development. Comparing model predictions from the two sets of estimated variables would further quantify—albeit implicitly—how an improved, physically based estimate of explanatory variables would be reflected in model predictions. Differences between the two sets of estimated explanatory variables and resultant model predictions vary spatially; greater differences are associated with areas of steep hydraulic gradients. A direct comparison, however, would require the development of a separate set of statistical models using explanatory variables from contributing recharge areas. Area-weighted means of three environmental variables—silt content, alfisol content, and depth to water from the U.S. Department of Agriculture State Soil Geographic (STATSGO) data—and one nitrogen-source variable (fertilizer-application rate from county data mapped to Enhanced National Land Cover Data 1992 (NLCDe 92) agricultural land use) can vary substantially between circular buffers and volume-equivalent contributing recharge areas and among contributing recharge areas for different sets of well variables. The differences in estimated explanatory variables are a function of the same factors affecting the contributing recharge areas as well as the spatial resolution and local distribution of the underlying spatial data. As a result, differences in estimated variables between circular buffers and contributing recharge areas are complex and site specific as evidenced by differences in estimated variables for circular buffers and contributing recharge areas of existing public-supply and network wells in the Great Miami River Basin. Large differences in areaweighted mean environmental variables are observed at the basin scale, determined by using the network of uniformly spaced hypothetical wells; the differences have a spatial pattern that generally is similar to spatial patterns in the underlying STATSGO data. Generally, the largest differences were observed for area-weighted nitrogen-application rate from county and national land-use data; the basin-scale differences ranged from -1,600 (indicating a larger value from within the volume-equivalent contributing recharge area) to 1,900 kilograms per year (kg/yr); the range in the underlying spatial data was from 0 to 2,200 kg/yr. Silt content, alfisol content, and nitrogen-application rate are defined by the underlying spatial data and are external to the groundwater system; however, depth to water is an environmental variable that can be estimated in more detail and, presumably, in a more physically based manner using a groundwater-flow model than using the spatial data. Model-calculated depths to water within circular buffers in the Great Miami River Basin differed substantially from values derived from the spatial data and had a much larger range. Differences in estimates of area-weighted spatial variables result in corresponding differences in predictions of nitrate occurrence in the aquifer. In addition to the factors affecting contributing recharge areas and estimated explanatory variables, differences in predictions also are a function of the specific set of explanatory variables used and the fitted slope coefficients in a given model. For models that predicted the probability of exceeding 1 and 4 milligrams per liter as nitrogen (mg/L as N), predicted probabilities using variables estimated from circular buffers and contributing recharge areas generally were correlated but differed significantly at the local and basin scale. The scale and distribution of prediction differences can be explained by the underlying differences in the estimated variables and the relative weight of the variables in the statistical models. Differences in predictions of exceeding 1 mg/L as N, which only includes environmental variables, generally correlated with the underlying differences in STATSGO data, whereas differences in exceeding 4 mg/L as N were more spatially extensive because that model included environmental and nitrogen-source variables. Using depths to water from within circular buffers derived from the spatial data and depths to water within the circular buffers calculated from the groundwater-flow model, restricted to the same range, resulted in large differences in predicted probabilities. The differences in estimated explanatory variables between contributing recharge areas and circular buffers indicate incorporation of physically based contributing recharge area likely would result in a different set of explanatory variables and an improved set of statistical models. The use of a groundwater-flow model to improve representations of source areas or to provide more-detailed estimates of specific explanatory variables includes a number of limitations and technical considerations. An assumption in these analyses is that (1) there is a state of mass balance between recharge and pumping, and (2) transport to a pumped well is under a steady state flow field. Comparison of volumeequivalent contributing recharge areas under steady-state and transient transport conditions at a location in the southeastern part of the basin shows the steady-state contributing recharge area is a reasonable approximation of the transient contributing recharge area after between 10 and 20 years of pumping. The first assumption is a more important consideration for this analysis. A gradient effect refers to a condition where simulated pumping from a well is less than recharge through the corresponding contributing recharge area. This generally takes place in areas with steep hydraulic gradients, such as near discharge locations, and can be mitigated using a finer model discretization. A boundary effect refers to a condition where recharge through the contributing recharge area is less than pumping. This indicates other sources of water to the simulated well and could reflect a real hydrologic process. In the Great Miami River Basin, large gradient and boundary effects—defined as the balance between pumping and recharge being less than half—occurred in 5 and 14 percent of the basin, respectively. The agreement between circular buffers and volume-equivalent contributing recharge areas, differences in estimated variables, and the effect on statisticalmodel predictions between the population of wells with a balance between pumping and recharge within 10 percent and the population of all wells were similar. This indicated process-model limitations did not affect the overall findings in the Great Miami River Basin; however, this would be model specific, and prudent use of a process model needs to entail a limitations analysis and, if necessary, alterations to the model.

Ohio↗

Biomarkers of mercury exposure in two eastern Ukraine cities

This study evaluates biomarkers of mercury exposure among residents of Horlivka, a city in eastern Ukraine located in an area with geologic and industrial sources of environmental mercury, and residents of Artemivsk, a nearby comparison city outside the mercury-enriched area. Samples of urine, blood, hair, and nails were collected from study participants, and a questionnaire was administered to obtain data on age, gender, occupational history, smoking, alcohol consumption, fish consumption, tattoos, dental amalgams, home heating system, education, source of drinking water, and family employment in mines. Median biomarker mercury concentrations in Artemivsk were 0.26 μg/g-Cr (urine), 0.92 μg/L (blood), 0.42 μg/g (hair), 0.11 μg/g (toenails), and 0.09 μg/g (fingernails); median concentrations in Horlivka were 0.15 μg/g-Cr (urine), 1.01 μg/L (blood), 0.14 μg/g (hair), 0.31 μg/g (toenails), and 0.31 μg/g (fingernails). Biomarkers of mercury exposure for study participants from Horlivka and Artemivsk are low in comparison with occupationally exposed workers at a mercury recycling facility in Horlivka and in comparison with exposures known to be associated with clinical effects. Blood and urinary mercury did not suggest a higher mercury exposure among Horlivka residents as compared with Artemivsk; however, three individuals living in the immediate vicinity of the mercury mines had elevated blood and urinary mercury, relative to overall results for either city. For a limited number of residents from Horlivka (N = 7) and Artemivsk (N = 4), environmental samples (vacuum cleaner dust, dust wipes, soil) were collected from their residences. Mercury concentrations in vacuum cleaner dust and soil were good predictors of blood and urinary mercury.

Journal of Occupational and Environmental Hygiene↗

Preliminary report on geophysics of the Verde River headwaters region, Arizona

This report summarizes the acquisition, data processing, and preliminary interpretation of a high-resolution aeromagnetic and radiometric survey near the confluence of the Big and Little Chino basins in the headwaters of the Verde River, Arizona. The goal of the aeromagnetic study is to improve understanding of the geologic framework as it affects groundwater flow, particularly in relation to the occurrence of springs in the upper Verde River headwaters region. Radiometric data were also collected to map surficial rocks and soils, thus aiding geologic mapping of the basin fill. Additional gravity data were collected to enhance existing coverage. Both aeromagnetic and gravity data indicate a large gradient along the Big Chino fault, a fault with Quaternary movement. Filtered aeromagnetic data show other possible faults within the basin fill and areas where volcanic rocks are shallowly buried. Gravity lows associated with Big Chino and Williamson Valleys indicate potentially significant accumulations of low-density basin fill. The absence of a gravity low associated with Little Chino Valley indicates that high-density rocks are shallow. The radiometric maps show higher radioactivity associated with the Tertiary latites and with the sediments derived from them. The surficial materials on the eastern side of the Big Chino Valley are significantly lower in radioactivity and reflect the materials derived from the limestone and basalt east of the valley. The dividing line between the low radioactivity materials to the east and the higher radioactiviy materials to the west coincides approximately with the major drainage system of the valley, locally known as Big Chino Wash. This feature is remarkably straight and is approximately parallel to the Big Chino Fault. The uranium map shows large areas with concentrations greater than 5 ppm eU, and we expect that these areas will have a significantly higher risk potential for indoor radon.

Arizona↗

Evaluation of the precipitation-runoff modeling system, Beaver Creek basin, Kentucky

The Precipitation Runoff Modeling System (PRMS) was evaluated with data from Cane branch and Helton Branch in the Beaver Creek basin of Kentucky. Because of previous studies, 10.6 years of record were available to establish a data base for the basin including 60 storms for Cane Branch and 50 storms for Helton Branch. The model was calibrated initially using data from the 1956-58 water years. Runoff predicted by the model was 94.7% of the observed runoff at Cane Branch (mined area) and 96.9% at Helton Branch (unmined area). After the model and data base were modified, the model was refitted to the 1956-58 data for Helton Branch. It then predicted 98.6% of the runoff for the 10.6-year period. The model parameters from Helton Branch were then used to simulate the Cane Branch runoff and discharge. The model predicted 102.6% of the observed runoff at Cane Branch for the 10.6 years. The simulations produced reasonable storm volumes and peak discharges. Sensitivity analysis of model parameters indicated the parameters associated with soil moisture are the most sensitive. The model was used to predict sediment concentration and daily sediment load for selected storm periods. The sediment computations indicated the model can be used to predict sediment concentrations during storm events. (USGS)

Water-Resources Investigations Report↗

Groundwater residence times in Shenandoah National Park, Blue Ridge Mountains, Virginia, USA: A multi-tracer approach

Chemical and isotopic properties of water discharging from springs and wells in Shenandoah National Park (SNP), near the crest of the Blue Ridge Mountains, VA, USA were monitored to obtain information on groundwater residence times. Investigated time scales included seasonal (wet season, April, 1996; dry season, August–September, 1997), monthly (March through September, 1999) and hourly (30-min interval recording of specific conductance and temperature, March, 1999 through February, 2000). Multiple environmental tracers, including tritium/helium-3 ( 3 H/ 3 He), chlorofluorocarbons (CFCs), sulfur hexafluoride (SF 6 ), sulfur-35 ( 35 S), and stable isotopes ( δ 18 O and δ 2 H) of water, were used to estimate the residence times of shallow groundwater discharging from 34 springs and 15 wells. The most reliable ages of water from springs appear to be based on SF 6 and 3 H/ 3 He, with most ages in the range of 0–3 years. This range is consistent with apparent ages estimated from concentrations of CFCs; however, CFC-based ages have large uncertainties owing to the post-1995 leveling-off of the CFC atmospheric growth curves. Somewhat higher apparent ages are indicated by 35 S (>1.5 years) and seasonal variation of δ 18 O (mean residence time of 5 years) for spring discharge. The higher ages indicated by the 35 S and δ 18 O data reflect travel times through the unsaturated zone and, in the case of 35 S, possible sorption and exchange of S with soils or biomass. In springs sampled in April, 1996, apparent ages derived from the 3 H/ 3 He data (median age of 0.2 years) are lower than those obtained from SF 6 (median age of 4.3 years), and in contrast to median ages from 3 H/ 3 He (0.3 years) and SF 6 (0.7 years) obtained during the late summer dry season of 1997. Monthly samples from 1999 at four springs in SNP had SF 6 apparent ages of only 1.2 to 2.5±0.8 years, and were consistent with the 1997 SF 6 data. Water from springs has low excess air (0–1 cm 3 kg −1 ) and N 2 –Ar temperatures that vary seasonally. Concentrations of He and Ne in excess of solubility equilibrium indicate that the dissolved gases are not fractionated. The seasonal variations in N 2 –Ar temperatures suggest shallow, seasonal recharge, and the excess He and Ne data suggest waters mostly confined to gas exchange in the shallow, mountain-slope, water-table spring systems. Water from wells in the fractured rock contains up to 8 cm 3 kg −1 of excess air with ages in the range of 0–25 years. Transient responses in specific conductance and temperature were observed in spring discharge within several hours of large precipitation events in September, 1999; both parameters increased initially, then decreased to values below pre-storm base-flow values. The groundwater residence times indicate that flushing rates of mobile atmospheric constituents through groundwater to streams draining the higher elevations in SNP average less than 3 years in base-flow conditions.

Virginia↗

Carbon stocks across a chronosequence of thinned and unmanaged red pine ( Pinus resinosa ) stands

Forests function as a major global C sink, and forest management strategies that maximize C stocks offer one possible means of mitigating the impacts of increasing anthropogenic CO 2 emissions. We studied the effects of thinning, a common management technique in many forest types, on age-related trends in C stocks using a chronosequence of thinned and unmanaged red pine ( Pinus resinosa ) stands ranging from 9 to 306 years old. Live tree C stocks increased with age to a maximum near the middle of the chronosequence in unmanaged stands, and increased across the entire chronosequence in thinned stands. C in live understory vegetation and C in the mineral soil each declined rapidly with age in young stands but changed relatively little in middle-aged to older stands regardless of management. Forest floor C stocks increased with age in unmanaged stands, but forest floor C decreased with age after the onset of thinning around age 40 in thinned stands. Deadwood C was highly variable, but decreased with age in thinned stands. Total ecosystem C increased with stand age until approaching an asymptote around age 150. The increase in total ecosystem C was paralleled by an age-related increase in total aboveground C, but relatively little change in total belowground C. Thinning had surprisingly little impact on total ecosystem C stocks, but it did modestly alter age-related trends in total ecosystem C allocation between aboveground and belowground pools. In addition to characterizing the subtle differences in C dynamics between thinned and unmanaged stands, these results suggest that C accrual in red pine stands continues well beyond the 60–100 year management rotations typical for this system. Management plans that incorporate longer rotations and thinning in some stands could play an important role in maximizing C stocks in red pine forests while meeting other objectives including timber extraction, biodiversity conservation, restoration, and fuel reduction goals.

Ecological Applications↗

Assessment of rangeland ecosystem conditions in Grand Canyon-Parashant National Monument, Arizona

Sustainability of dryland ecosystems depends on the functionality of soil-vegetation feedbacks that affect ecosystem processes, such as nutrient cycling, water capture and retention, soil erosion and deposition, and plant establishment and reproduction. Useful, common indicators can provide information on soil and site stability, hydrologic function, and biotic integrity. Evaluation of rangeland health thus relies on describing the condition and sustainability of these individual, measurable, and observable indicators that are linked to important ecosystem processes. This report focuses on the ~200,000 acres of the Grand Canyon-Parashant National Monument that is administered by the National Park Service (NPS)—one of the largest NPS units where livestock grazing is a permitted land-use activity. Many ecosystems in the monument are characterized by a low degree of resilience to improper grazing because of low and variable precipitation. The monument is marked by a high degree of environmental heterogeneity, including a large elevation gradient, widely differing precipitation patterns, a diversity of geologic substrates, and unique combinations of plant species. The objective of this report is to (1) increase our understanding of the underlying landscape, soil, and climate setting factors that affect Grand Canyon-Parashant National Monument dryland ecosystem structure and function (also referred to as land potential) and (2) characterize the condition of monument ecosystems in relation to management concepts, such as rangeland health. Data were analyzed by elevation zone using both univariate and multivariate approaches. Survey results document the high level of diversity within the study area, including 15 unique soil taxa and 271 species of plants. We collected three new plant species for Grand Canyon-Parashant National Monument and 17 new species for the NPS portion of the monument. Results also document a strong association between rangeland health indicators and elevation, topographic setting, and soils. Soil factors found to explain important variation across plots include the amount of exposed bedrock, soil rockiness, soil texture (and associated hydrologic properties), and soil depth. We also found that dominant species turnover across elevation may represent species’ differences in adaptation to climates, including Larrea tridentata , Coleogyne ramosissima , and Artemisia spp. Bromus rubens is the most common invasive species of concern recorded in this study, but other common invasive species are Bromus tectorum , Erodium cicutarium , and Schismus arabicus . Correlations between an index of cattle use and indicators of rangeland health suggest that areas with high cattle use have increased bare ground, decreased ground cover, increased frequency of Schismus arabicus , decreased cover of Coleogyne ramosissima and Ephedra spp., and increased cover of Gutierrezia spp. The few strong correlations observed between indicators of vascular plant community cover or abundance and indicators of cattle activity support rangeland assessment and monitoring strategies that do not rely solely on plant-based indicators are needed. This work supports management of dryland ecosystems, including Grand Canyon-Parashant National Monument, using concepts of land potential. We conclude the report with recommendations on improving existing land-potential-based classification systems, associated interpretations, and methods for moving forward with a Grand Canyon-Parashant National Monument rangeland monitoring program.

Arizona↗

Reconnaissance of the chemical quality of surface waters of the Neches River basin, Texas

The kinds and quantities of minerals dissolved in the surface water of the Neches River basin result from such environmental factors as geology, streamflow patterns and characteristics, and industrial influences. As a result of high rainfall in the basin, much of the readily soluble material has been leached from the surface rocks and soils. Consequently, the water in the streams is usually low in concentrations of dissolved minerals and meets the U.S. Public Health Service drinking-water standards. In most streams the concentration of dissolved solids is less than 250 ppm (parts per million). The Neches River drains an area of about 10,000 square miles in eastern Texas. From its source in southeast Van Zandt County the river flows in a general southeasterly direction and empties into Sabine Lake, an arm of the Gulf of Mexico. In the basin the climate ranges from moist subhumid to humid, and the average annual rainfall ranges from 46 inches is the northwest to more than 52 inches in the southeast. Annual runoff from the basin has averaged 11 inches; however, runoff rates vary widely from year to year. The yearly mean discharge of the Neches River at Evadale has ranged from 994 to 12,720 cubic feet per second. The rocks exposed in the Neches River basin are of the Quaternary and Tertiary Systems and range in age from Eocene to Recent. Throughout most of the basin the geologic formations dip generally south and southeast toward the gulf coast. The rate of dip is greater than that of the land surface; and as a result, the older formations crop out to the north of the younger formations. Water from the outcrop areas of the Wilcox Group and from the older formations of the Claiborne Group generally has dissolved-solids concentrations ranging from 100 to 250 ppm; water from the younger formations has concentrations less than 100 ppm. The northern half of the basin has soft water, with less than 60 ppm hardness. The southern half of .the basin has very soft water, usually with less than 30 ppm hardness. The chloride concentrations are less than 20 ppm in surface water in the southern half of the basin and usually range from 20 to 100 ppm in the northern half of the basin. Concentrations greater than 100 ppm are found only where pollution is occurring. The Neches River basin has an abundance of surface water, but uneven distribution of runoff makes storage projects necessary to provide dependable water supplies. The principal existing reservoirs, with the exception of Striker Creek Reservoir, contain water of excellent quality. Chemical-quality data for the Striker Creek drainage area indicate that its streams are affected by .the disposal of brines associated with oil production. Sam Rayburn Reservoir began impounding water in 1965. The water impounded should prove of acceptable quality for most uses, but municipal and industrial wastes released into the Angelina River near Lufkin may have a degrading effect on the quality of the water, especially during extended periods of low flows. Water available for storage at the many potential reservoir sites will be of good quality; but, if the proposed salt-water barrier is to impound acceptable water, the disposal of oilfield brine into Pine Island Bayou should be discontinued.

Water Supply Paper↗

Survival, cause-specific mortality, and population growth of white-tailed deer in western Virginia

Understanding the role of recruitment in population dynamics of white-tailed deer ( Odocoileus virginianus ) is important for management. In the central Appalachian Mountains, deer are part of a largely forested ecosystem that supports 3 carnivore species thought to be capable of influencing white-tailed deer recruitment: black bears ( Urus americanus ), coyotes ( Canis latrans ), and bobcats ( Lynx rufus ). Yet little is known about predation, how other environmental factors influence recruitment, or the importance of neonate survival to white-tailed deer population performance in the region. Our objectives were to identify causes of mortality for neonates, analyze effects of landscape attributes on survival of neonates, estimate survival rates for neonates and adult female white-tailed deer, and to model population growth trends based on current vital rates and hypothetical harvest and neonate survival scenarios. During 2019–2020, we captured 57 neonate deer in Bath County, Virginia, USA, by monitoring 38 pregnant females equipped with global positioning system collars and vaginal implant transmitters and by conducting transect searches for recently born neonates. We observed 37 neonate mortalities and identified cause of death using field and genetic evidence. Mortalities included 28 predation events and 9 deaths from other causes (e.g., abandonment, malnutrition, disease). Black bears accounted for 48.6% of neonate mortalities, and 64.2% of predation events ( n = 18), followed by bobcats ( n = 5) and coyotes ( n = 3). Annual survival for adult female deer was 0.871 and neonate survival to 12 weeks old was 0.310. Elevation was a significant predictor of neonate survival; mortality risk increased 20% for every 100-m increase in elevation. Models of annual population growth using observed vital rates predicted an increasing population (λ = 1.10). A 10% increase in female harvest would still result in a potential population increase of 2% (λ = 1.02), but a 20% increase in harvest rate would result in a potential 7% decline (λ = 0.93). Neonate survival was higher near fertile valley bottoms and lower along forested ridges characterized by shallow, infertile soils and limited edge or early successional forests. While predation, largely influenced by black bears, was the leading cause of neonate mortality and contributed to low neonate survival, we observed little evidence of population decline, and suggest there is opportunity for a modest increase in harvest of female deer.

Virginia↗

Trace silicon determination in biological samples by inductively coupled plasma mass spectrometry (ICP-MS): Insight into volatility of silicon species in hydrofluoric acid digests for optimal sample preparation and introduction to ICP-MS

A method for the determination of trace levels of silicon from biological materials by inductively coupled plasma mass spectrometry (ICP-MS) has been developed. The volatility of water-soluble silicon species, hexafluorosilicic acid (H 2 SiF 6 ), and sodium metasilicate (Na 2 SiO 3 ) was investigated by evaporating respective solutions (50 µg/mL silicon) in nitric acid (HNO 3 ), nitric acid + hydrochloric acid (HNO 3 + HCl), and nitric acid + hydrochloric acid + hydrofluoric acid (HNO 3 + HCl + HF) at 120 °C on a hot-block to near dryness. The loss of silicon from H 2 SiF 6 solutions was substantial (>99%) regardless of the digestion medium. Losses were also substantial (>98%) for metasilicate solutions heated in HNO 3 + HCl + HF, while no significant loss occurred in HNO 3 or HNO 3 + HCl. These results show that H 2 SiF 6 species were highly volatile and potential losses could confound accuracy at trace level determinations by ICP-MS if digestates prepared in HF are heated to eliminate HF. Among the various matrices comprising major elements, sodium appeared to be effective in reducing silicon loss from H 2 SiF 6 solutions. Excess sodium chloride (NaCl) matrix provided better stability, improving silicon recoveries by up to about 80% in evaporated HF digestates of soil and mine waste samples, but losses could not be fully prevented. To safely remove excess acids and circumvent the adverse effects of excess HF (e.g., risk of high Si background signals), a two-step digestion scheme was adopted for the preparation of biological samples containing trace silicon levels. A closed-vessel digestion was performed either in 4 mL of concentrated HNO 3 and 1 mL of concentrated HCl or 4 mL of concentrated HNO 3 , 1 mL of concentrated HCl and 1 mL of concentrated HClO 4 on a hot plate at 140 °C. Digestates were then evaporated to incipient dryness at 120 °C to remove the acids. A second closed-vessel digestion was carried out to dissolve silicates in 0.5 mL of concentrated HNO 3 and 0.5 mL of concentrated HF at 130 °C. After digestion, digestates were diluted to 10 mL. The solution containing about 5% HNO 3 and 5% HF was directly analyzed by ICP-MS equipped with an HF-inert sample introduction system. The limit of detection was about 110 µg/L for 28 Si when using the Kinetic Energy Discrimination (KED) mode. The method was used to determine silicon in various plant and tissue certified reference materials. Data were acquired for 28 Si using KED and standard (STD) modes, and 74 Ge and 103 Rh as internal standard elements. There was not any significant difference between the accuracy and precision of the results obtained with 74 Ge and 103 Rh within the same measurement mode. Precision, calculated as relative standard deviation for four replicate analyses, varied from 5.3 (tomato leaves) to 21% (peach leaves) for plant and from 2.2 (oyster tissue) to 33% (bovine liver) for tissue SRM/CRMs. Poor precision was attributed to material heterogeneity and the large particle size distribution. An analysis of lung tissue samples from those with occupational exposure to silica dust revealed that tissues possessed substantial levels of water-soluble silicates, but the most silicon was present in the particulate matter fraction.

Minerals↗