Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Natural Resource Modelling”

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

Assessment experimental semivariogram uncertainty in the presence of a polynomial drift

The semivariogram, which measures the spatial variability between experimental data, is generally used as a structural input in all two-point geostatistical procedures. However, in most geoscience applications, experimental semivariograms are usually computed from a limited number of sparsely spaced measurements, which results in uncertainty associated with the semivariance values estimated for a specified number of lags. More importantly, considering a spatial variable modelled by a nonstationary random field, uncertainty is not only in the experimental semivariogram of the residuals, but also in the coefficients of the drift model estimated from the available experimental data. Therefore, when assessing the reliability of an experimental semivariogram (or estimated semivariances) in the nonstationary case, both aforementioned uncertainties should be taken into account. The aim of this paper is to extend the “Generalised Bootstrap” procedure to the nonstationary model by propagating the uncertainty associated with the estimated drift coefficients into the uncertainty in the experimental semivariogram of the residuals. The proposed methodology is demonstrated in a case study using abundant geophysical measurements characterised by a nonstationary random function. Two scenarios are evaluated in the case study: (1) it is assumed that the drift coefficients can be estimated without any uncertainty, and (2) uncertainty of the drift coefficients is taken into account. We have explained the methodology that allows to assess the uncertainty of the semivariogram lag estimates in the presence of the drift in the mean. Considering the second scenario, uncertainty is obviously larger than the case where uncertainty of the drift in the mean is ignored. This evaluation should be considered in applications where the data is often rather limited, such as subsurface hydrology (i.e. porosity, transmissivity), soil science (i.e. heavy metal content, soil moisture) and mining (i.e. scoping or pre-feasibility stage of the project). In fact, modern geostatistics should provide not only the semivariogram estimates but also estimation of its uncertainty.

Natural Resources Research↗

Prefabricated Zero Ascend Omnispecies (ZAO) modular fish passage modules using advanced manufacturing techniques

The Zero Ascend Omnispecies (ZAO) fish passage attraction module “ZAO Attractor” is a novel system designed to facilitate upstream passage for a broad spectrum of fish species, with a particular focus on American shad and river herring (alewife and/or blueback herring). Combining the use of prefabricated, modular components with advanced manufacturing, the ZAO aims to provide a cost-effective and scalable solution for fish passage attractions across various hydropower sites. The system features a surface level entry, so fish do not need to climb. In addition, it is designed to be adjusted and repositioned to create the most effective attraction flows by managing the outflows from hydropower turbines. Budget Period 1 (BP1) focused on the design and modeling of the ZAO Attractor and demonstrated the engineering feasibility of a modular, prefabricated, hydrodynamically efficient structure with surface level entry and passage. During BP1, the hypothesis was that certain structural adaptations, such as the addition of long fins, could direct turbine outflows to create a desirable downstream flow. Simulations, however, were inconclusive, revealing the complex nature of fish preferences for hydraulic conditions. It was discovered that there were limited published studies on the direct correlation between quantified flow patterns and fish attraction, emphasizing the need for controlled studies to observe and measure these effects. In Budget Period 2 (BP2), controlled studies with live fish were prioritized to address these complexities. A partially scaled ZAO Attractor was evaluated in laboratory conditions using native river water to assess its capability to attract and guide fish. Despite challenges with flume flow conditions and the impact of flume dimensions on flow directionality and velocity, the studies provided insights into fish behavior and preferences, particularly the requirement for a consistent flow stream for migration cues. Adjustments to the fin structure demonstrated potential effectiveness in guiding fish toward the module. However, the hydraulic conditions within the flume environment did not accurately reflect typical in-river scenarios, therefore the results were inconclusive with respect to the optimal system configuration for in-river use. In the future, conducting comprehensive field studies can enhance understanding of the interactions between fish behavior and flow conditions created by the ZAO Attractor. Additionally, a design review is advised, particularly regarding the load scenarios throughout the seasons and positioning of the ZAO Attractor. These studies should continue to refine the design to achieve its intended benefits: economical and effective fish attraction and passage that can be adapted to specific site conditions while allowing for continued energy generation at sites up for re-licensure and adding hydropower generation to previously untapped resources in non-powered dams and new stream reaches. In conclusion, the ZAO fish passage attraction module presents a promising approach to fish passage at hydropower stations, with the potential to balance environmental concerns with renewable energy development. However, further field validation and optimization are necessary to ensure its effectiveness and practicality in real-world applications.

Final Technical Report↗

Modeling, mapping, and measuring the risk of freshwater invasive species across Alaska

Freshwater ecosystems of the Alaskan Arctic and Subarctic provide resources that are culturally, ecologically, and economically invaluable. Presently, these regions are relatively free of the impacts from invasive species compared to southern latitudes. To date, there have been relatively few verified introductions of aquatic invasive species (AIS) to freshwater ecosystems in Alaska. The expanding list and distribution of AIS has led to significant negative ecological and economic impacts (e.g., waterweed Elodea nuttalli ; E. canadensis and northern pike Esox Lucius introduced outside its native range in Alaska). Escalating human activity across Alaskan lands and waters, coupled with rapidly shifting environmental conditions, increases the potential for new species introductions and subsequent establishment. Creating a proactive framework for well-informed decision-making and action can improve the effectiveness of prevention efforts and bolster decision support tools that help resource managers direct limited resources. Prioritizing AIS that may be introduced and become established, as well as the locations at highest risk of invasion, is foundational to building a proactive invasive species management framework in Alaska. This project sought to identify and prioritize AIS known to be invasive in the contiguous United States, evaluate current and future habitat suitability for AIS in Alaska, and assess potential for AIS to be transported to habitats across Alaska, utilizing similar assessment methods as implemented for Bering Sea marine invasive species and non-native plants in Alaska. To accomplish this goal, the objectives of the project were to: 1) develop a formal ranked list of potential AIS to freshwater systems of Alaska; 2) assess the level of establishment risk for potential AIS by developing habitat suitability models for waterbodies across Alaska; and 3), identify potential pathways and specific vectors for high-risk AIS to invade Alaska and develop a framework for how vector analysis will be completed to understand transport risk. Overall, our goal is horizon scanning which is defined by Roy et al. (2019) as “a systematic examination of potential threats and opportunities, within a given context, and likely future developments, which are at the margin of current thinking and planning.” The scans include pathway analyses and risk screening of species present at pathway origin points, with a focus on identifying species at high risk of being introduced, becoming established, spreading, and causing harm. We refined a list of 28 AIS from a list of hundreds based on characterizations of species’ invasiveness and species’ proximity to Alaska (USGS 2020; GBIF 2022). Next, we evaluated the relative invasiveness of individual species to create an initial AIS ranking. We sought to characterize habitat suitability of AIS by selecting variables that were continental in scale, covering North America to include Alaska as well as the lower 48 states comparing natural discharge, sub-basin average terrain slope (degrees), average silt fraction, average organic carbon, lithological class, and human footprint in sub-basin in 2009. We estimated AIS habitat suitability across the entire state of Alaska using the physiological tolerances of the AIS (Appendix 2). We also evaluated pathways and vectors for the introduction of AIS (Appendix 2). Many pathways and vectors considered did not meet the criteria for Alaska or freshwater systems. Of the 28 ranked species that we categorized as very high, high, and moderate levels of invasiveness; all three risk groups included fish and mollusks (Appendix 2). One commonality of the very high-invasiveness-ranked species was the availability of Ecological Risk Screening Summary documents (USFWS, 2022) produced by U.S. Fish and Wildlife Service (USFWS), except for the goldfish ( Carassius auratus ) and the New Zealand mudsnail ( Potamopyrgus antipodarum ). The Ecological Risk Screening Summary is now available for New Zealand mudsnails. In general, fish species often ranked very high or high in invasiveness and included sportfish and aquarium fish, suggesting the importance of pathways such as aquarium trade, fishing industry, intentional (but illegal) introductions of sportfishes and aquarium fishes for establishment. The technique we used for habitat suitability models necessitated aquatic environmental datasets that were continental in scale, which was often interpolated from very coarse resolution source data layers, particularly in Alaska. Better spatial data representing aquatic environments would likely improve this approach. While the lack of introductions in Alaska and nearby provinces and states is encouraging, the lack of occurrence data for the focal species also created complications for habitat suitability modeling. Despite the challenges, the habitat suitability models indicated limited suitability for warmwater species while some species, such as Brook trout (S alvelinus fontinalis ), have high habitat suitability across Alaska no matter what threshold approach is taken. Some environmental predictors were more important than others. Specifically, the most important predictor variable, ‘frost free days,’ was critical for 15 out of 28 species as expected due to harsh winter conditions in Arctic and Subarctic regions. The second most important predictor was ‘subbasin land surface runoff’, a variable that indicates the amount of discharge and runoff, while the third most important predictor was ‘snow cover’ another indication of winter conditions. Overall, the ability to understand the effect of future climate scenarios on the establishment of AIS was challenging. A detailed dataset of freshwater temperatures and water chemistry (e.g., pH, calcium) would greatly improve the ability to predict invasiveness of freshwater species to Alaska’s ecosystems on a regional basis. Future studies may benefit from a more focused geographic scope examining a group of subbasins or a regional basin rather than the entire state. These drainages could be selected based upon the mostly likely locations of introduction pathways. The two most prevalent pathway risks for AIS are in-state transfer and stowaways/contaminants. Although there are examples of introductions from other pathways, the risk is somewhat mitigated by Alaska’s climate and regulations. However, variable application of protocols for inspection and cleaning of fishing gear, watercraft, and other similar items while traveling into Alaska as well as transferring from waterbody to waterbody within the state creates a substantial risk in introducing invasive species. We plot cumulative invasive vulnerability for all subbasins and for the top 10% of subbasins (Appendix 3).

Alaska↗

Water quality, bed-sediment quality, and simulation of potential contaminant transport in Foster Creek, Berkeley County, South Carolina, 1991-93

Foster Creek, a freshwater tidal creek in Berkeley County, South Carolina, is located in an area of potential contaminant sources from residential, commercial, light industrial, and military activities. The creek is used as a secondary source of drinking water for the surrounding Charleston area. Foster Creek meets most of the freshwater- quality requirements of State and Federal regulatory agencies, but often contains low concentrations of dissolved oxygen and has been characterized as eutrophic. Investigations of water- and bed-sediment quality were made between 1991 and 1993 to assess the effects of anthropogenic sources of contamination on Foster Creek. Low-flow surface-water samples were generally free of toxic compounds with the exception of laboratory artifacts and naturally occurring trace metals. Storm-runoff samples generally contained very low concentrations (near detection limits) of a small number of volatile and semivolatile organics and naturally occurring trace metals. Concentrations of toxic compounds in excess of current (1995) South Carolina Department of Health and Environmental Control and U.S. Environmental Protection Agency regulations were not detected in surface-water samples collected from Foster Creek. Chemical analyses of streambed sediments indicated minimal anthropogenic effects on sediment quality. The particle-tracking option of the U.S. Geological Survey one-dimensional unsteady-flow model (BRANCH) indicated that as the simulated volume of rainfall runoff increased in the Foster Creek Basin, simulated particles in Foster Creek were transported greater distances. Simulating flow through the Bushy Park Dam (also known as Back River Dam) had little effect on particle movement in Foster Creek. Simulating typical withdrawal rates at a water-supply intake resulted in a slight attraction of particles toward the intake during conditions of relatively low runoff. These withdrawals had a greater influence on particles downstream of the intake than on those upstream of the intake. Simulations confirmed earlier findings which suggested that the creek would not flush during baseflow conditions, with the exception of the lower 1-mile reach, where flushing results from tidal movements. According to the simulations, Foster Creek will fully flush if a 2-year, 7-day storm occurs. Flushing appears to be affected more by the total volume of storm runoff than by typical municipal withdrawals or tidal effects.

South Carolina↗

Assessing denitrification in groundwater using natural gradient tracer tests with 15N: In situ measurement of a sequential multistep reaction

Denitrification was measured within a nitrate‐contaminated aquifer on Cape Cod, Massachusetts, using natural gradient tracer tests with 15 N nitrate. The aquifer contained zones of relatively high concentrations of nitrite (up to 77 μM) and nitrous oxide (up to 143 μM) and has been the site of previous studies examining ground water denitrification using the acetylene block technique. Small‐scale (15–24 m travel distance) tracer tests were conducted by injecting 15 N nitrate and bromide as tracers into a depth interval that contained nitrate, nitrite, nitrous oxide, and excess nitrogen gas. The timing of the bromide breakthrough curves at down‐gradient wells matched peaks in 15 N abundance above background for nitrate, nitrite, nitrous oxide, and nitrogen gas after more than 40 days of travel. Results were simulated with a one‐dimensional transport model using linked reaction kinetics for the individual steps of the denitrification reaction pathway. It was necessary to include within the model spatial variations in background concentrations of all nitrogen oxide species. The model indicated that nitrite production (0.036–0.047 μmol N (L aquifer) −1 d −1 ) was faster than the subsequent denitrification steps (0.013–0.016 μmol N (L aquifer) −1 d −1 for nitrous oxide and 0.013–0.020 μmol N (L aquifer) −1 d −1 for nitrogen gas) and that the total rate of reaction was slower than indicated by both acetylene block tracer tests and laboratory incubations. The rate of nitrate removal by denitrification was much slower than the rate of transport, indicating that nitrate would migrate several kilometers down‐gradient before being completely consumed.

Water Resources Research↗

Potential for, and possible effects of, artificial recharge in Carson Valley, Douglas County, Nevada

Rapid population growth in Carson Valley, west- central Nevada, requires a dependable municipal water source. Artificial recharge of aquifers using available flow of the Carson River is one way to increase the amount of water in underground storage and maintain a dependable ground-water supply. Ground water can be artificially recharged by routing excess surface water or, after proper treatment, routing wastewater to infiltration basins or injection wells. Withdrawal wells would remove stored water when needed. As a first step, maps showing areas in Carson Valley with high, low, moderate and unknown potential for artificial recharge were developed on the basis of the distribution of geologic units, depth to water, specific yield, infiltration rate, and location of natural recharge and discharge. For recharge by means of infiltration, areas totaling 5,700 acres have high potential, 23,900 acres have moderate potential, and 6,200 acres have low potential. For recharge through injection, areas totaling 7,800 acres have high potential and 43,500 acres have moderate potential; 23,000 acres have unknown potential because data are lacking on subsurface conditions. A ground-water-flow model was used to assess the possible results of artificial recharge. Simulations with no accompanying ground-water withdrawal show that, when recharge by injection is simulated near the valley floor, heads in the semiconfined aquifer increase over much of the valley, floor; only about 20 percent of the recharged water is stored in the aquifer after 5 years and as much as 80 percent is lost to streamflow and evapotranspiration. When recharge is simulated on the eastern side of the valley, 80 percent of the recharged water remains in storage after 5 years. When recharge is simulated near the valley floor, more water is lost to discharge than when recharge is on the eastern side of the valley. When recharge is applied for long periods without accompanying withdrawal, recharged water moves downgradient to discharge areas. The recharge water that discharges to the surface-water system could in turn replenish base flow of the Carson River and benefit downstream users.

Water-Resources Investigations Report↗

Deep well injection of brine from Paradox Valley, Colorado: Potential major precipitation problems remediated by nanofiltration

Groundwater brine seepage into the Dolores River in Paradox Valley, Colorado, increases the dissolved solids load of the Colorado River annually by ∼2.0 × 10 8 kg. To abate this natural contamination, the Bureau of Reclamation plans to pump ∼3540 m 3 /d of brine from 12 shallow wells located along the Dolores River. The brine, with a salinity of 250,000 mg/L, will be piped to the deepest (4.9 km) disposal well in the world and injected mainly into the Mississippian Leadville Limestone. Geochemical modeling indicates, and water-rock experiments confirm, that a huge mass of anhydrite (∼1.0 × 10 4 kg/d) likely will precipitate from the injected brine at downhole conditions of 120°C and 500 bars. Anhydrite precipitation could increase by up to 3 times if the injected brine is allowed to mix with the highly incompatible formation water of the Leadville Limestone and if the Mg in this brine dolomitizes the calcite of the aquifer. Laboratory experiments demonstrate that nanofiltration membranes, which are selective to divalent anions, provide a new technology that remediates the precipitation problem by removing ∼98% of dissolved SO 4 from the hypersaline brine. The fluid pressure used (50 bars) is much lower than would be required for traditional reverse osmosis membranes because nanofiltration membranes have a low rejection efficiency (5–10%) for monovalent anions. Our results indicate that the proportion of treatable brine increases from ∼60% to >85% with the addition of trace concentrations of a precipitation inhibitor and by blending the raw brine with the effluent stream.

Water Resources Research↗

Flood characteristics of Oklahoma streams techniques for calculating magnitude and frequency of floods in Oklahoma, with compilations of flood data through 1971

The 2-, 5-, 10-, 25-, 50-, and 100-year recurrence interval floods are related to basin and climatic parameters for natural streams in Oklahoma by multiple regression techniques through the mathematical model, Q x =aA b S c P d , where Q x is peak discharge for recurrence interval x, A is contributing drainage area, S is main channel slope, P is mean annual precipitation, and a, b, c, and d are regression constants and coefficients. One equation for each recurrence interval applies statewide for all natural streams of less than 2,500 mil (6,500 km 2 ), except where manmade works, such as dams, flood-detention structures, levees, channelization, and urban development, appreciably affect flood runoff. The equations can be used to estimate flood frequency of a stream at an ungaged site if drainage area size, main channel slope, and mean annual precipitation are known. At or near gaged sites, a weighted average of the regression results and the gaging station data is recommended. Individual relations of flood magnitude to contributing drainage area are given for all or parts of the main stems of the Arkansas, Salt Fork Arkansas, Cimarron, North Canadian, Canadian, Washita, North Fork Red, and Red Rivers. Parts of some of these streams, and all of the Neosho and Verdigris Rivers are not included because the effects of. major regulation from large reservoirs cannot be evaluated within the scope of the report. Graphical relations of maximum floods of record for eastern and western Oklahoma provide a guide to maximum probable floods. A random sampling of the seasonal occurrence of floods indicated about two-thirds of all annual floods in Oklahoma occur during. April through July. Less than one-half of one percent of annual floods occur in December. A compilation of flood records at all gaging sites in Oklahoma and some selected sites in adjacent States is given in an appendix. Basin and climatic parameters and log-Pearson Type III frequency data and statistics are given for most station records. A second appendix gives a reprint of the U.S. Water Resources Council Bulletin 15 which describes procedures for fitting a log-Pearson Type III distribution to gaging station data.

Oklahoma↗

Maps of estimated nitrate and arsenic concentrations in basin-fill aquifers of the southwestern United States

Human-health concerns and economic considerations associated with meeting drinking-water standards motivated a study of the vulnerability of basin-fill aquifers to nitrate contamination and arsenic enrichment in the southwestern United States. Statistical models were developed by using the random forest classifier algorithm to predict concentrations of nitrate and arsenic across a model grid representing about 190,600 square miles of basin-fill aquifers in parts of Arizona, California, Colorado, Nevada, New Mexico, and Utah. The statistical models, referred to as classifiers, reflect natural and human-related factors that affect aquifer vulnerability to contamination and relate nitrate and arsenic concentrations to explanatory variables representing local- and basin-scale measures of source and aquifer susceptibility conditions. Geochemical variables were not used in concentration predictions because they were not available for the entire study area. The models were calibrated to assess model accuracy on the basis of measured values. Only 2 percent of the area underlain by basin-fill aquifers in the study area was predicted to equal or exceed the U.S. Environmental Protection Agency drinking-water standard for nitrate as N (10 milligrams per liter), whereas 43 percent of the area was predicted to equal or exceed the standard for arsenic (10 micrograms per liter). Areas predicted to equal or exceed the drinking-water standard for nitrate include basins in central Arizona near Phoenix; the San Joaquin Valley, the Santa Ana Inland, and San Jacinto Basins of California; and the San Luis Valley of Colorado. Much of the area predicted to equal or exceed the drinking-water standard for arsenic is within a belt of basins along the western portion of the Basin and Range Physiographic Province that includes almost all of Nevada and parts of California and Arizona. Predicted nitrate and arsenic concentrations are substantially lower than the drinking-water standards in much of the study area-about 93 percent of the area underlain by basin-fill aquifers was less than one-half the standard for nitrate as N (5.0 milligrams per liter), and 50 percent was less than one-half the standard for arsenic (5.0 micrograms per liter). The predicted concentrations and the improved understanding of the susceptibility and vulnerability of southwestern basin-fill aquifers to nitrate contamination and arsenic enrichment can be used by water managers as a qualitative tool to assess and protect the quality of groundwater resources in the Southwest.

Arizona, California, Colorado, Idaho, Nevada, New ↗

Ground-water resources and simulated effects of withdrawals in the East Shore area of Great Salt Lake, Utah

The ground-water resources in the East Shore area of Great Salt Lake, Utah, were studied to better define the ground-water system; to document changes in ground-water levels, quality, and storage; and to simulate effects of an increase in ground-water withdrawals. The East Shore aquifer system is in basin-fill deposits, and is primarily a confined system with unconfined parts near the mountain front. Recharge to and discharge from the East Shore aquifer system were estimated to average about 160,000 acre-feet per year during 1969-84, with minor amounts of water being removed from storage during that period. Major sources of ground-water recharge are seepage from surface water in natural channels and irrigation canals, and subsurface inflow from consolidated rock to the basin-fill deposits. Discharge of ground water is primarily to wells, water courses, springs, and as diffuse seepage to Great Salt Lake. Average annual surface-water inflow to the study area was estimated to be 860,000 acre-feet for the period 1969-84. Annual withdrawal of ground water for municipal and industrial use increased from about 10,000 acre-feet in 1960 to more than 30,000 acre-feet in 1980 to supply a population that increased from 175,000 in 1960 to 290,000 in 1980. Long-term trends of ground-water levels indicate a steady decline at most observation wells since 1952, despite near normal or increased precipitation since the late 1960's. Water levels declined as much as 50 feet near the principal pumping center in the east-central part of the study area. They declined as much as 35 feet more than five miles from the pumping center. The increase in withdrawals and subsequent water-level declines have caused about 700 wells within 30 square miles to cease flowing since 1954. A numerical model of the East Shore aquifer system in the Weber Delta area was constructed and calibrated using water-level data and changes in ground-water withdrawals for 1955-85. Predictive simulations were made based on doubling the 1980-84 rate of municipal and industrial withdrawals for 20 years, and using both average and below-average recharge rates. The simulations indicated water-level declines of an additional 35 to 50 feet near the principal pumping center; a decrease in natural discharge to drains, evapotranspiration, and Great Salt Lake; and a decrease in ground-water storage of 80,000 to 115,000 acre-feet after 20 years.

Utah↗

Intra-community implications of implementing multiple tsunami-evacuation zones in Alameda, California

Tsunami-evacuation planning in coastal communities is typically based on maximum evacuation zones for a single scenario or a composite of sources; however, this approach may over-evacuate a community and overly disrupt the local economy and strain emergency-service resources. To minimize the potential for future over-evacuations, multiple evacuation zones based on arrival time and inundation extent are being developed for California coastal communities. We use the coastal city of Alameda, California (USA), as a case study to explore population and evacuation implications associated with multiple tsunami-evacuation zones. We use geospatial analyses to estimate the number and type of people in each tsunami-evacuation zone and anisotropic pedestrian evacuation models to estimate pedestrian travel time out of each zone. Results demonstrate that there are tens of thousands of individuals in tsunami-evacuation zones on the two main islands of Alameda, but they will likely have sufficient time to evacuate before wave arrival. Quality of life could be impacted by the high number of government offices, schools, day-care centers, and medical offices in certain evacuation zones and by potentially high population density at one identified safe area after an evacuation. Multi-jurisdictional evacuation planning may be warranted, given that many at-risk individuals may need to evacuate to neighboring jurisdictions. The use of maximum evacuation zones for local tsunami sources may be warranted given the limited amount of available time to confidently recommend smaller zones which would result in fewer evacuees; however, this approach may also result in over-evacuation and the incorrect perception that successful evacuations are unlikely.

California↗

Vegetation classification and mapping, Vicksburg National Military Park, Mississippi

Executive Summary The National Park Service (NPS) Gulf Coast Inventory and Monitoring Network, with the support of the National Park Service Vegetation Inventory completed vegetation classification and mapping for Vicksburg National Military Park (VICK), in Warren County, Mississippi, from 2004 to 2010. VICK is one of more than 250 NPS units to be covered by the Vegetation Inventory. Methods and procedures follow those of the NPS Vegetation Inventory as of August, 2012 ( http://science.nature.nps.gov/im/inventory/veg/index.cfm ). Ecologists collected floristic and environmental data from 45 vegetation field plots and classified and described 13 plant community types corresponding to US National Vegetation Classification (USNVC) associations from these data. This classification was used to map 712 hectares (1,698 acres) of Vicksburg National Military Park, in ten map classes corresponding to vegetation associations and two non-vegetated land cover classes, from a digital orthomosaic created from 1:24,000 color infrared aerial photographs and digital elevation model (DEM) data. Methods used to delineate stands varied according to vegetation type and were a combination of incorporating stand boundaries from an existing map, new “heads-up” interpretation and modeling of digital environmental spatial data. A thematic accuracy assessment collected vegetation data from 118 field observations stratified across the vegetation classes. The overall accuracy of the vegetation map was estimated to be 60.5% at the thematic level of USNVC association and 92.4% at the thematic level of USNVC group. Estimates of areas of vegetation types, as adjusted for mapping error, were calculated for all types mapped and/or observed during the thematic accuracy assessment. An important finding of this study is that the composition of the forests of VICK is substantially different from that described by earlier, more qualitative, accounts. Implications of these findings and possible applications of the data for resource monitoring and management are presented. This report summarizes the methods and general results of the Vegetation Inventory for Vicksburg National Military Park. Appendices include descriptions of vegetation types, a field key to the types, a list of plant species, and accuracy assessment contingency tables. This report is supplemented by the project data, which include: spatial database containing map polygon attributes. Geospatial products are in Universal Transverse Mercator (UTM) projection, Zone 15, using the North American Datum of 1983. all field data (plot, observation point, and accuracy assessment point) stored in a Microsoft Access database metadata for all digital products The data and reports produced by this investigation reside at: https://irma.nps.gov/App/ and/or at: http://biology.usgs.gov/npsveg/

Mississippi↗

Evaluation of salinity and nutrient conditions in the Heart River Basin, North Dakota, 1970–2020

The Heart River Basin is predominantly an agricultural basin in western North Dakota and is approximately 3,350 square miles. The U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service and the Grant County Soil Conservation District, completed a study to assess spatial and temporal patterns of water quality in the Heart River Basin. The purpose of this report is to describe the methods and results of a study to evaluate salinity and nutrients in the Heart River Basin in western North Dakota. Water-quality and streamflow data used in the study were compiled from 1970 to 2020 using the National Water Quality Monitoring Council Water Quality Portal and National Water Information System. Changes in streamflow characteristics were investigated at three sites from 1970 to 2020, and changes in water quality were investigated at four sites from 1974 to 2019. Streamflow analysis indicated decreasing streamflow from 1970 until the late 1990s followed by increasing streamflow through 2020, with the largest increase in the 7-day minimum streamflow or base flow. For the historical water-quality trend period (1974–2019), total dissolved solids, sulfate, sodium, chloride, and sodium adsorption ratio concentrations have increased since the mid-1970s through 2019. Potassium concentrations during the historical period remained mostly constant with some small fluctuations. Calcium and magnesium concentrations increased since the mid-1970s at all sites, except for a decrease at one site between 1974 and 1999. During the recent trend period (1999–2019), increasing concentrations in total dissolved solids, sulfate, sodium, chloride, calcium, magnesium, and sodium adsorption ratios were observed across the Heart River Basin. The magnitude of the increases was smaller at tributary sites compared to main-stem sites. During the recent period, potassium was mostly constant, although small (−0.9 milligram per liter or less) decreases on tributaries and minor (1.3 milligrams per liter) increases on the main-stem sites were detected. Unlike dissolved ion concentrations, significant increases in nutrient concentrations were not detected from 1999 to 2019, but nitrate plus nitrite concentrations most likely decreased upstream from Lake Tschida. Inverse modeling for period 1 (1974–99) in model zone 1 (Heart River reach from site 5 to site 6) had eight reasonable models that indicated the clay mineral-water interactions and dissolution of evaporites control the geochemistry. Results of the inverse modeling for period 2 (1999–2019) in model zone 1 also had eight reasonable models that indicated that the dissolution of evaporites was the major geochemical control. Results of the geochemical modeling for period 1 (1974–99) in model zone 2 (Heart River and Sweetbriar Creek reach from sites 20 and 21 to site 22) produced seven reasonable models, and the geochemical control of the system was the dissolution of sulfate evaporite minerals. Geochemical modeling results for period 2 (1999–2019) in model zone 2 produced 11 reasonable models and was also controlled by the dissolution of sulfate evaporite minerals. Differences between the two model zones indicated that geology controls some of the water-quality changes in the Heart River Basin. Loads were estimated for total dissolved solids, sulfate, sodium, and chloride and total phosphorus. Annual loads estimated for the Heart River from 2013 through 2020 at the Heart River site upstream from Lake Tschida (site 5) and near Mandan (site 22) were generally greatest in 2014 and least in 2016 for total dissolved solids, sulfate, sodium, and chloride. Most of the annual loads of total dissolved solids, sulfate, sodium, and chloride are delivered in March through July in the Heart River at these sites and are likely from snowmelt and spring and summer rains. The mean annual yields of total dissolved solids and sodium from 2013 to 2020 generally were largest in Big Muddy Creek (site 18), whereas yields of sulfate and chloride were largest at Sweetbriar Creek (site 21) compared to the other selected sites in the Heart River Basin. Larger yields of total dissolved solids, sulfate, sodium, and chloride at sites located on Big Muddy Creek and Sweet Briar Creek in the lower Heart River Basin were likely a result of differences in geology and soils upstream from the selected sites. A mass balance of total dissolved solids, sulfate, sodium, and chloride was estimated for the lower Heart River Basin, specifically the reach below Lake Tschida to Mandan (site 7 to site 22). Intervening flow was the largest contributor to the dissolved ion loads in the lower Heart River Basin and is an important part of understanding the transport of dissolved ions in the basin. The intervening load can include groundwater discharge, irrigation return flow, local runoff, and input from smaller ephemeral tributaries. Tributaries in the lower Heart River Basin contributed portions of the total dissolved solids, sulfate, sodium, and chloride loads at the Heart River near Mandan (site 22) that generally were proportional to the streamflow contributions. Annual loads for total phosphorus between 2013 and 2020 at the Heart River site upstream from Lake Tschida (site 5) and near Mandan (site 22) generally were largest in 2019 and smallest in 2016. Most of the total phosphorus loads for main-stem sites 5 and 22 were transported in March, April, and June, likely from snowmelt and early summer rains. The mean annual yields of total phosphorus for 2013–20 were largest on the main-stem site upstream from Lake Tschida (site 5) and Sweetbriar Creek (site 21), whereas the smallest yields were in Big Muddy Creek (site 18). Much of the phosphorus that enters Lake Tschida from the upper basin does not get transported downstream to the lower basin, and much of the phosphorus in the lower basin was attributed to intervening flow.

North Dakota↗

Simulation of Fish, Mud, and Crystal Lakes and the shallow ground-water system, Dane County, Wisconsin

A new MODFLOW lake package (LAK3) that simulates ground-water/lake interaction was used in simulation of Fish, Mud and Crystal Lakes?three shallow seepage lakes located in northwestern Dane County, Wis. The simulations were done to help determine the cause of increasing lake stages and provide a tool to estimate the effect of pumping water from Fish lake on future lake stages. The ground-water-flow model was developed using a telescopic-mesh refinement of the Dane and southwestern Columbia Counties regional model previously developed by the U.S. Geological Survey and the Wisconsin Geological and Natural History Survey. The parameter estimation model, UCODE, was coupled to the steadystate ground-water model to automate and optimize the calibration procedure. The steady-state model was calibrated to measured ground-water levels, Spring Creek streamflow measured at Lodi, and Fish and Crystal Lake stages. The results of the steady-state model were used as initial conditions in a transient simulation beginning in 1966 and ending in 1998. Recharge based on annual baseflow in Black Earth Creek, runoff based on measured coefficients, and precipitation and evaporation from the lake surfaces, were varied during the transient simulation. Measured Fish Lake stage was matched to simulated stage to calibrate the transient model. Model results suggest that the increase in regional ground-water recharge resulted in increased ground-water flow to the lake, which in turn resulted in increased lake stages. Simulation results of withdrawal of water from Fish Lake at 500 gallons per minute, assuming 1990?98 climatic conditions, indicate that after 1 year of pumping the stage of Fish and Mud Lakes would be reduced more than 1 foot and the stage of Crystal Lake would be reduced by less than 0.2 foot. When pumping is stopped, the lake stages would recover to near pre-pumping levels within about 3 years. When pumping is extended to 5 years, Fish and Mud Lake stage would be reduced by a maximum of 3.8 feet and Crystal Lake stage is reduced a maximum of 0.8 feet. After 4 years of recovery, Fish and Mud Lake stages are within 0.9 foot of prepumping levels and Crystal Lake stage is within 0.7 foot.

Wisconsin↗

Geochemical and isotopic composition of ground water, with emphasis on sources of sulfate, in the upper Floridan aquifer and intermediate aquifer system in southwest Florida

In southwest Florida, sulfate concentrations in water from the Upper Floridan aquifer and overlying intermediate aquifer system are commonly above 250 milligrams per liter (the drinking water standard), particularly in coastal areas. Possible sources of sulfate include dissolution of gypsum from the deeper part of the Upper Floridan aquifer or the middle confining unit, saltwater in the aquifer, and saline waters from the middle confining unit and Lower Floridan aquifer. The sources of sulfate and geochemical processes controlling ground-water composition were evaluated for the Peace and Myakka River Basins and adjacent coastal areas of southwest Florida. Samples were collected from 63 wells and a saline spring, including wells finished at different depth intervals of the Upper Floridan aquifer and intermediate aquifer system at about 25 locations. Sampling focused along three ground-water flow paths (selected based on a predevelopment potentiometric-surface map). Ground water was analyzed for major ions, selected trace constituents, dissolved organic carbon, and stable isotopes (delta deuterium, oxygen-18, carbon-13 of inorganic carbon, and sulfur-34 of sulfate and sulfide); the ratio of strontium-87 to strontium-86 was analyzed for waters along one of the flow paths. Chemical and isotopic data indicate that dedolomitization reactions (gypsum and dolomite dissolution and calcite precipitation) control the chemical composition of water in the Upper Floridan aquifer in inland areas. This is confirmed by mass-balance modeling between wells in the shallowest interval in the aquifer along the flow paths. However, gypsum occurs deeper in the aquifer than these wells. Upwelling of sulfate-rich water that previously dissolved gypsum in deeper parts of the aquifer is a more likely source of sulfate than gypsum dissolution in shallow parts of the aquifer. This deep ground water moves to shallower zones in the aquifer discharge area. Saltwater from the Upper Floridan aquifer has not dissolved significant amounts of gypsum compared to fresher water in the aquifer. This is consistent with a shallow seawater source for the saltwater, rather than a deeper source from the underlying middle confining unit or Lower Floridan aquifer, which would have elevated sulfate concentrations. Ion exchange and dolomitization may be important reactions for saltwater in the aquifer. According to geochemical modeling, the freshwater end member for water in the saltwater mixing zone in the southwestern part of the study area is not upgradient water from the Upper Floridan aquifer that dissolved gypsum. Instead, this water appears to be isolated from the regional freshwater flow system and may be part of a more localized flow system. The chemical and isotopic composition of water in the intermediate aquifer system is controlled by differences in extent of reactions with aquifer minerals, upward leakage from the Upper Floridan aquifer, and saltwater mixing. In inland areas, water generally is characterized by relatively low sulfate concentrations (less than 250 milligrams per liter) and differences in extent of carbonate mineral dissolution. Some inland waters have elevated chloride concentrations, which may be related to evaporation prior to recharge. In coastal Sarasota County and in isolated inland areas, water from the intermediate aquifer system has high sulfate concentrations characteristic of dedolomitization waters from the Upper Floridan aquifer. The chemical and isotopic composition of these waters is controlled by upward leakage from the Upper Floridan aquifer, which naturally occurs in the discharge area but may be locally enhanced by pumping or interconnection of wells open to both aquifer systems. In western Charlotte County, the waters are dominated by sodium and chloride, and their compositions are consistent with mixing between saltwater and inland intermediate aquifer system water that has not been influenced by discharge from the

Florida↗

Evaluation of ground-water flow and solute transport in the Lompoc area, Santa Barbara County, California

Ground-water quality in the Lompoc area, especially in the Lompoc plain, is only marginally acceptable for most uses. Demand for ground water has increased for municipal use since the late 1950's and has continued to be high for irrigation on the Lompoc plain, the principal agricultural area in the Santa Ynez River basin. As use has increased, the quality of ground water has deteriorated in some areas of the Lompoc plain. The dissolved-solids concentration in the main zone of the upper aquifer beneath most of the central and western plains has increased from less than 1,000 milligrams per liter in the 1940's to greater than 2,000 milligrams per liter in the 1960's. Dissolved- solids concentration have remained relatively constant since the 1960's. A three-dimensional finite-difference model was used to simulate ground-water flow in the Lompoc area and a two-dimensional finite-element model was used to simulate solute transport to gain a better understanding of the ground-water system and to evaluate the effects of proposed management plans for the ground-water basin. The aquifer system was simulated in the flow model as four horizontal layers. In the area of the Lompoc plain, the layers represent the shallow, middle, and main zones of the upper aquifer, and the lower aquifer. For the Lompoc upland and Lompoc terrace, the four layers represent the lower aquifer. The solute transport model was used to simulate dissolved-solids transport in the main zone of the upper aquifer beneath the Lompoc plain. The flow and solute-transport models were calibrated to transient conditions for 1941-88. A steady-state simulation was made to provide initial conditions for the transient-state simulation by using long-term average (1941-88) recharge rates. Model- simulated hydraulic heads generally were within 5 feet of measured heads in the main zone for transient conditions. Model-simulated dissolved- solids concentrations for the main zone generally differed less than 200milligrams per liter from concentrations in 1988. During 1941-88 about 1,096,000 acre-feet of water was pumped from the aquifer system. Average pumpage for this period (22,830 acre-feet per year) exceeded pumpage for the steady-state simulation by 16,590 acre-feet per year. The results of the transient simulation indicate that about 60 percent of this increase in pumpage was contributed by increased recharge, 28 percent by decreased natural discharge from the system (primarily discharge to the Santa Ynez River and transpiration), and 13 percent was withdrawn from storage. Total simulated downward leakage from the middle zone to the main zone in the central plain and upward leakage from the consolidated rocks to the main zone significantly increased in response to increased pumpage, which increased from about 6,240 to 30,870 acre-feet per year from 1941 to 1988. Average dissolved-solid concentration in the middle zone in 1987-88 ranged from 2,000 to 3,000 milligrams per liter beneath the northeastern plain and the dissolved-solids concentration of two samples from the consolidated rocks beneath the western plain averaged 4,300 milligrams per liter. Because the dissolved-solids concentration for the middle zone and the consolidated rocks is higher than the simulated steady-state dissolved-solids concentration of the main zone, the increase in the leakage from these two sources resulted in increased dissolved-solids concentration in the main zone during the transient period. The model results indicate that the main source of increased dissolved- solids concentration in the northeastern and central plains was downward leakage from the middle zone; whereas, upward leakage from the consolidated rocks was the main source of the increased dissolved-solids concentrations in the northwestern and western plains. The models were used to estimate changes in hydraulic head and in dissolved-solids concentration resulting from three proposed management alternatives: (1) average recharge

Water-Resources Investigations Report↗

The ground-water system in the LaGrange Aquifer near LaGrange, southeastern Wyoming

Groundwater is being developed from the La Grange aquifer in southeastern Wyoming. It consists of saturated permeable alluvium that is hydraulically connected with most of the underlying White River Group. In the area of principal interest east of Horse Creek, Hawk Springs Reservoir and 14 adjacent wells used to supplement surface-water supply in the reservoir are in a natural discharge area. Upgradient of the reservoir there are 28 irrigation wells in about a 6-square-mile area. In this area, water levels declined between 3 and 12 feet from 1973 to 1978 causing concern about the effects of well pumpage on the hydrologic system. A digital model was developed and used to simulate the two-dimensional groundwater flow system in the unconfined La Grange aquifer. Transient simulations were made for 1973-78 using 12 time periods and for 1978-80 using 25 time periods. The calibrated digital model was used to simulate four 6-month pumping alternatives including three hypothetical alternatives for the area of principal interest east of Horse Creek. The reservoir altitude was held constant, approximating a reservoir volume 7,000 acre-feet. Pumping alternative 1 simulated historic conditions for 1973-78 including monthly recharge from precipitation which was included in the next three pumping alternatives. For pumping alternative 1, the calculated rate of discharge at the end of the 6-month simulation was 5.2 cubic feet per second from the aquifer to the reservoir. At the end of the 6-month simulations for pumping alternatives 2 and 3, the calculated rate of discharge to the reservoir was decreased to 0.4 cubic feet per second by pumpage from the 14 wells and to 3.8 cubic feet per second by pumpage from the 28 irrigation wells. For pumping alternative 4, pumpage from the total 42 wells resulted in a 1.0 cubic feet per second loss from the reservoir to the aquifer. (USGS)

Water-Resources Investigations Report↗

Ground-water hydrology of the Toppenish Creek basin, Yakima Indian Reservation, Washington

A groundwater flow model was constructed for the Toppenish Creek basin aquifer system in eastern Washington. Flow was simulated in three aquifer units: (1) the confined old valley fill and shallow basalt (unit 2); (2) the underlying primary basalt (unit 3); and (3) the deep basalt (unit 4). Water levels in the overlying unconfined alluvial aquifer (unit 1) were held fixed. Calibrated transmissivities ranged from 0.01 to 0.48 foot squared per second. Calibrated storage coefficients were 0.0004 to 0.006. The confining-bed leakance ranged from 2.0 x 10 to the minus 11th power to 2.5 x 10 to the minus 10th power feet per second per foot. Under steady-state conditions (1954) annual natural recharge was about 29,000 acre-ft underflow from adjacent basins. Annual pumpage increased from less than 500 acre-ft in 1954 to an average of 19,600 acre-ft for 1971 and 1972. Pumpage caused simulated declines in unit 3 of up to 95 ft for 1955-72. Projected annual declines from 1973-77 using 1971-72 pumpage were about 0.5-1.5 ft in unit 2 and 0.2-1.5 ft in unit 3. The corresponding declines from 1978-82 were 0.2-1.2 ft and 0.2-0.8 ft, respectively. Using 1971-72 pumpage plus 12 ,400 acre-ft per year from unit 3, the calculated annual declines from 1978-82 were 1-36 ft in unit 2 and 4-20 ft in unit 3. (USGS)

Water-Resources Investigations Report↗