SEARCH · USGS Science

Results for “Current Climate Change Reports”

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

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26Linked to original sources

USGS2003

Ecosystem history of southern and central Biscayne Bay: Summary report on sediment core analyses

During the last century, the environs of Biscayne Bay have been greatly affected by anthropogenic alteration through urbanization of the Miami/Dade County area. The sources, timing, delivery, and quality of freshwater flow into the Bay have been changed by construction of a complex canal system that controls movement of water throughout south Florida. Changes in shoreline and sub-aquatic vegetation and marine organisms have been observed and changes in water delivery are believed to be the cause. Current restoration goals are attempting to restore natural flow of fresh water into Biscayne and Florida Bays and to restore the natural fauna and flora, but first we need to determine pre-alteration baseline conditions in order to establish targets and performance measures for restoration. This research is part of an ongoing study designed to address the needs of the Biscayne Bay Coastal Wetlands Project (BBCW) of the Comprehensive Everglades Restoration Plan (CERP). By establishing the natural patterns of temporal change in salinity, water quality, vegetation, and benthic fauna in Biscayne Bay and the nearby wetlands over the last 100- 500 years the USGS, in collaboration with our partners, will provide the data necessary to set realistic targets to achieve the BBCW Project goals. Six cores from three sites in Biscayne Bay were collected in April 2002 for multidisciplinary multi-proxy analyses. This report details the results of these analyses and compares the 2002 cores to cores collected in 1997. The following are our significant findings to date: The salinity of central Biscayne Bay has become increasingly marine and increasingly stable over the last 100 years. At No Name Bank, prior to approximately 1915, the inter-decadal and decadal salinity fluctuations appear to have been greater than after 1915 when salinities stabilized at that site. Continental shelf/open marine influence on the sites has increased during the 20th century. There is no indication of inter-decadal salinity extremes or periods of hypersalinity. Freshwater and mesohaline salinities have had a minor influence on No Name bank throughout the time of deposition. At Featherbed Bank, the influence is reduced to mesohaline salinities. Card Bank has experienced relatively large swings in salinity over multi-decadal and centennial timescales, compared to central Biscayne Bay, but marine influence at the site has increased over the last century. Indications of regional scale patterns have been found, especially in the shell chemistry data and the pollen assemblages. These regional patterns indicate that the changes are not site specific and may not be limited to Biscayne Bay. Sub-aquatic vegetation has undergone bay-wide patterns of change over the last 200-500 years, which includes expansion prior to 1900 and declines during the last century in central Biscayne Bay. Thalassia appears to have increased at all three core sites sometime between 1550 and 1750 AD. A decline in Thalassia appears to have occurred after 1950 at No Name Bank and slightly earlier in the 20th century at Featherbed Bank. Card Bank does not appear to have experienced any declines in vegetation on an inter-decadal scale during the 20th century. Molluscan faunal abundance and diversity have undergone significant changes in central Biscayne Bay. Indicators of increased organic-rich sediments at No Name occur between 1869 to 1888 and between the 1930’s to 1975. These changes may correlate to human activities (settlement, population growth). These findings represent a first step towards the project’s goal to reconstruct the history of Biscayne Bay and they provide us with a working model to be tested at other sites. It is clear from our findings that Biscayne Bay has been a dynamic environment over the last 500 years, with natural changes occurring in salinity and benthic habitats. However, several significant changes have occurred in the 20th century: 1) increased stabilization of marine salinities; 2) declines in seagrass in central Biscayne Bay; 3) dramatic changes in molluscan abundance and diversity in central Biscayne Bay. The question remains - how do we better differentiate natural cycles of change from anthropogenic change within these observed trends? The preliminary implications from our research are that changes in salinity and benthic habitats have occurred naturally in Biscayne Bay on inter-decadal to centennial scales, perhaps due to climatic changes, changes in sea level, bank migrations, or a combination of factors. However, further work needs to be done to determine which components of change in the 20th century are human-induced and which are natural. By examining the historical records preserved in the sediments of Biscayne Bay, we can provide restoration trust agencies with the information necessary to set realistic targets and performance measures for Biscayne Bay.

Florida↗

USGS2010

Relation of urbanization to stream habitat and geomorphic characteristics in nine metropolitan areas of the United States

The relation of urbanization to stream habitat and geomorphic characteristics was examined collectively and individually for nine metropolitan areas of the United States?Portland, Oregon; Salt Lake City, Utah; Denver, Colorado; Dallas?Forth Worth, Texas; Milwaukee?Green Bay, Wisconsin; Birmingham, Alabama; Atlanta, Georgia; Raleigh, North Carolina; and Boston, Massachusetts. The study was part of a larger study conducted by the U.S. Geological Survey from 1999 to 2004 to examine the effects of urbanization on the physical, chemical, and biological components of stream ecosystems. The objectives of the current study were to determine how stream habitat and geomorphic characteristics relate to different aspects of urbanization across a variety of diverse environmental settings and spatial scales. A space-for-time rural-to-urban land-cover gradient approach was used. Reach-scale habitat data and geomorphic characteristic data were collected once during low flow and included indicators of potential habitat degradation such as measures of channel geometry and hydraulics, streambed substrate, low-flow reach volume (an estimate of base-flow conditions), habitat complexity, and riparian/bank conditions. Hydrologic metrics included in the analyses were those expected to be altered by increases in impervious surfaces, such as high-flow frequency and duration, flashiness, and low-flow duration. Other natural and human features, such as reach-scale channel engineering, geologic setting, and slope, were quantified to identify their possible confounding influences on habitat relations with watershed-scale urbanization indicators. Habitat and geomorphic characteristics were compared to several watershed-scale indicators of urbanization, natural landscape characteristics, and hydrologic metrics by use of correlation analyses and stepwise linear regression. Habitat and geomorphic characteristics were related to percentages of impervious surfaces only in some metropolitan areas and environmental settings. The relations between watershed-scale indicators of urbanization and stream habitat depended on physiography and climate, hydrology, pre-urban channel alterations, reach-scale slope and presence of bedrock, and amount of bank stabilization and grade control. Channels increased in size with increasing percentages of impervious surfaces in southeastern and midwestern metropolitan areas regardless of whether the pre-existing land use was forest or agriculture. The amount of enlargement depended on annual precipitation and frequency of high-flow events. The lack of a relation between channel enlargement and increasing impervious surfaces in other metropolitan areas was thought to be confounded by pre-urbanization hydrologic and channel alterations. Direct relations of channel shape and streambed substrate to urbanization were variable or lacking, probably because the type, amount, and source of sediment are dependent on the phase of urbanization. Reach-scale slope also was important for determining variations in streambed substrate and habitat complexity (percentage of riffles and runs). Urbanization-associated changes in reach-scale riparian vegetation varied geographically, partially depending on pre-existing riparian vegetation characteristics. Bank erosion increased in Milwaukee?Green Bay and Boston urban streams, and bank erosion also increased with an increase in a streamflow flashiness index. However, potential relations likely were confounded by the frequent use of channel stabilization and bank protection in urban settings. Low-flow reach volume did not decrease with increasing urbanization, but instead was related to natural landscape characteristics and possibly other unmeasured factors. The presence of intermittent bedrock in some sampled reaches likely limited some geomorphic responses to urbanization, such as channel bed erosion. Results from this study emphasize the importance of including a wide range of landscape variables at m

Scientific Investigations Report↗

USGS2004

The Blackwater NWR inundation model. Rising sea level on a low-lying coast: land use planning for wetlands

The Blackwater National Wildlife Refuge (BNWR), on the Eastern Shore of Chesapeake Bay (figure 1), occupies an area less than 1 meter above sea level. The Refuge has been featured prominently in studies of the impact of sea level rise on coastal wetlands. Most notably, the refuge has been sited by the Intergovernmental Panel on Climate Change (IPCC) as a key example of 'wetland loss' attributable to rising sea level due to global temperature increase. Comparative studies of aerial photos taken since 1938 show an expanding area of open water in the central area of the refuge. The expanding area of open water can be shown to parallel the record of sea level rise over the past 60 years. The U.S. Fish and Wildlife Service (FWS) manages the refuge to support migratory waterfowl and to preserve endangered upland species. High marsh vegetation is critical to FWS waterfowl management strategies. A broad area once occupied by high marsh has decreased with rising sea level. The FWS needs a planning tool to help predict current and future areas of high marsh available for waterfowl. 'Wetland loss' is a relative term. It is dependant on the boundaries chosen for measurement. Wetland vegetation, zoned by elevation and salinity (figure 3), respond to rising sea level. Wetlands migrate inland and upslope and may vary in areas depending on the adjacent land slopes. Refuge managers need a geospatial tool that allows them to predict future areas that will be converted to high and intertidal marsh. Shifts in location and area of coverage must be anticipated. Viability of a current marsh area is also important. When will sea level rise make short-term management strategies to maintain an area impractical? The USGS has developed an inundation model for the BNWR centered on the refuge and surrounding areas. Such models are simple in concept, but they require a detailed topographic map upon which to superimpose future sea level positions. The new system of LIDAR mapping of land and shallow water surfaces has solved this problem. Our team has developed a detailed LIDAR map of the BNWR area at a 30 centimeter (ca. 1 ft) contour interval (figure 2). The new map allows us to identify the present marsh vegetation zones and to predict the location and area of future zones on a decade-by- decade basis over the next century at increments of sea level rise on the order of 3 cm/decade (ca. 1 inch). We have developed two scenarios for the model. The first is a steady-state model that uses the historic rate of sea level rise of 3.1 mm/yr to predict marsh areas. The second is a 'global warming' scenario utilizing a conservative IPCC model with an exponentially-increasing rate of sea level rise. Under either scenario, the BNWR is progressively inundated with an expanding core of open water. Although their positions change in the future, the areas of intertidal marsh as well as those of the critical high marsh remain fairly constant until the year 2050. Beyond that time, the low-lying land surface is overtopped by rising sea level and the area is dominated by open water. Our model suggests that wetland habitat in the Blackwater area might be maintained and sustained through a combination of public and private preservation efforts through easements in combination with judicious Federal land acquisition into the predicted areas of suitable marsh formation - but for only the next 50 years. Beyond that time much of this area will become open water.

Open-File Report↗

USGS2005

Environmental setting, water budget, and stream assessment for the Broad Run watershed, Chester County, Pennsylvania

The Broad Run watershed lies almost entirely in West Bradford Township, Chester County, Pa., and drains 7.08 square miles to the West Branch Brandywine Creek. Because of the potential effect of encroaching development and other stresses on the Broad Run watershed, West Bradford Township, the Chester County Water Resources Authority, and the Chester County Health Department entered into a cooperative study with the U.S. Geological Survey to complete an annual water budget and stream assessment of overall conditions. The annual water budget quantified the basic parameters of the hydrologic cycle for the climatic conditions present from April 1, 2003, to March 31, 2004. These water-budget data identified immediate needs and (or) deficits that were present within the hydrologic cycle during that period, if present; however, an annual water budget encompassing a single year does not identify long-term trends. The stream assessment was conducted in two parts and assessed the overall condition of the watershed, an overall assessment of the fluvial-geomorphic conditions within the watershed and an overall assessment of the stream-quality conditions. The data collected will document present (2004) conditions and identify potential vulnerabilities to future disturbances. For the annual period from April 1, 2003, to March 31, 2004, determination of an annual water budget indicated that of the 67.8 inches of precipitation that fell on the Broad Run watershed, 38.8 inches drained by way of streamflow to the West Branch Brandywine Creek. Of this 38.8 inches of streamflow, local-minimum hydrograph separation techniques determined that 7.30 inches originated from direct runoff and 31.5 inches originated from base flow. The remaining precipitation went into ground-water storage (1.71 inches) and was lost to evapotranspiration (27.3 inches). Ground-water recharge for this period-35.2 inches-was based on these values and an estimated ground-water evapotranspiration rate of 2 inches. Assessment of fluvial-geomorphic conditions included large-scale mapping of stream classes within the Broad Run watershed and in-depth study of three representative stream reaches also within the Broad Run watershed. Based on the total distance of all stream reaches classified within the Broad Run watershed, 61 percent were classified as C-class, 14 percent as E-class, 13 percent as B-class, 5 percent as F-class, 4 percent as undifferentiated B- and F-class, 2 percent as G-class, and less than 1 percent as A-class. The map of stream classes indicates that the Broad Run watershed currently has no large-scale areas of stream impairment and that, generally, the stream is not entrenched and the main branch of the Broad Run has an available, functioning flood plain. Smaller tributary streams, however, showed signs of localized entrenchment due to site-specific influences such as natural stream-channel evolution, localized channelization, localized contraction due to road and driveway crossings, and (or) increased localized runoff. For example, one small reach along a tributary channel was observed to become entrenched due to runoff originating from a new housing development. Entrenched stream reaches are merely located by large-scale mapping and require individual assessment to determine potential causes of entrenchment and (or) future restorative actions. Three in-depth geomorphic study sites showed that the Broad Run can currently be considered graded or in a state of dynamic equilibrium. The sites did, however, identify certain vulnerabilities to future changes within the watershed. These vulnerabilities included disruption of the present sediment supply, including both increases and (or) reductions in the current sediment loads within the Broad Run; increases in both magnitude and duration of storm-water runoff; encroachment of development onto present flood-plain areas, and (or) alterations to riparian zones. Assessment of stream-quality conditions includ

Scientific Investigations Report↗

USGS2012

Methods of practice and guidelines for using survey-grade global navigation satellite systems (GNSS) to establish vertical datum in the United States Geological Survey

Geodetic surveys have evolved through the years to the use of survey-grade (centimeter level) global positioning to perpetuate and post-process vertical datum. The U.S. Geological Survey (USGS) uses Global Navigation Satellite Systems (GNSS) technology to monitor natural hazards, ensure geospatial control for climate and land use change, and gather data necessary for investigative studies related to water, the environment, energy, and ecosystems. Vertical datum is fundamental to a variety of these integrated earth sciences. Essentially GNSS surveys provide a three-dimensional position x, y, and z as a function of the North American Datum of 1983 ellipsoid and the most current hybrid geoid model. A GNSS survey may be approached with post-processed positioning for static observations related to a single point or network, or involve real-time corrections to provide positioning "on-the-fly." Field equipment required to facilitate GNSS surveys range from a single receiver, with a power source for static positioning, to an additional receiver or network communicated by radio or cellular for real-time positioning. A real-time approach in its most common form may be described as a roving receiver augmented by a single-base station receiver, known as a single-base real-time (RT) survey. More efficient real-time methods involving a Real-Time Network (RTN) permit the use of only one roving receiver that is augmented to a network of fixed receivers commonly known as Continually Operating Reference Stations (CORS). A post-processed approach in its most common form involves static data collection at a single point. Data are most commonly post-processed through a universally accepted utility maintained by the National Geodetic Survey (NGS), known as the Online Position User Service (OPUS). More complex post-processed methods involve static observations among a network of additional receivers collecting static data at known benchmarks. Both classifications provide users flexibility regarding efficiency and quality of data collection. Quality assurance of survey-grade global positioning is often overlooked or not understood and perceived uncertainties can be misleading. GNSS users can benefit from a blueprint of data collection standards used to ensure consistency among USGS mission areas. A classification of GNSS survey qualities provide the user with the ability to choose from the highest quality survey used to establish objective points with low uncertainties, identified as a Level I, to a GNSS survey for general topographic control without quality assurance, identified as a Level IV. A Level I survey is strictly limited to post-processed methods, whereas Level II, Level III, and Level IV surveys integrate variations of a RT approach. Among these classifications, techniques involving blunder checks and redundancy are important, and planning that involves the assessment of the overall satellite configuration, as well as terrestrial and space weather, are necessary to ensure an efficient and quality campaign. Although quality indicators and uncertainties are identified in post-processed methods using CORS, the accuracy of a GNSS survey is most effectively expressed as a comparison to a local benchmark that has a high degree of confidence. Real-time and post-processed methods should incorporate these "trusted" benchmarks as a check during any campaign. Global positioning surveys are expected to change rapidly in the future. The expansion of continuously operating reference stations, combined with newly available satellite signals, and enhancements to the conterminous geoid, are all sufficient indicators for substantial growth in real-time positioning and quality thereof.

Techniques and Methods↗

USGS2005

Advancing migratory bird conservation and management by using radar: An interagency collaboration

Migratory birds face many changes to the landscapes they traverse and the habitats they use. Wind turbines and communications towers, which pose hazards to birds and bats in flight, are being erected or proposed across the United States and offshore. Human activities can also destroy or threaten habitats critical to birds during migratory passage, and climate change appears to be altering migratory patterns. The U.S. Fish and Wildlife Service (USFWS) and other agencies are under increasing pressure to identify and evaluate movement patterns and habitats used during migration and other times. Few tools for deciphering migratory travels exist, but radar-based studies of movements and habitat use patterns in songbirds, waterfowl, and bats hold promise. The U.S. system of over 150 Doppler weather radars provides continental coverage, similar to the scale of bird migration. Although data stored from weather radar represent perhaps the second largest biological data archive in the world, use of those data is currently limited to technically savvy biologists who can handle the obscure data formats. Complementary mobile radar units and thermal and acoustic monitoring are also used in site-specific studies. Efforts to advance bird conservation and management through the use of radar arose independently in several USFWS/USGS collaborations. Recently, this coalition of scientists and resource managers identified the need to work together more closely to foster radar-related research and software development. U.S. Geological Survey (USGS) scientists at Fort Collins Science Center, National Wetlands Research Center, Northern Rocky Mountain Science Center, and Patuxent Wildlife Research Center, as well as USFWS Migratory Bird biologists across the country, are collaborating with university partners to develop a suite of products for managers. The goals are to identify migratory pathways and stopover sites for conservation, mitigation, and landscape planning; convey the importance of functional landscapes and unobstructed airspaces for migrating wildlife; enable use of radar by the wider biological, wind power, and related communities; and simplify the analysis of radar data. The long term focus is to use radar technologies to better understand movement patterns and habitat associations of migratory birds and other wildlife. Land managers and industry may use the knowledge and tools developed to optimize the siting of energy projects, other facilities, and migratory bird habitat projects. The complementary endeavors, not all of which are funded, concentrate on four fundamentals: (1) develop software, in collaboration with National Oceanic and Atmospheric Administration (NOAA) scientists, that enables biologists to access unfiltered weather data and integrate it into standard geographic information systems; (2) develop artificial intelligence-based filters that separate bird from nonbird radar echoes; (3) determine characteristics of bird migration in terms of altitude, speed and direction, daily movements, seasonality, and associations with habitats and landforms; and (4) examine specific movement patterns in relation to towers, wind generation facilities, and tall obstructions. Many technical issues make this work difficult, including complex data structures, massive data sets, digital recognition of birds, large areas not covered by weather radar, and model validation; however, progress will only be furthered by tackling the challenge. The new coalition will meets its goals by: (1) facilitating a productive collaboration with NOAA, Department of the Interior bureaus, state wildlife agencies, universities, power companies, and other potential partners; (2) building and strengthening scientific capabilities within USGS; (3) addressing key migratory bird management issues; and (4) ensuring full funding for the collaborative effort.

Open-File Report↗

USGS2015

Aquifer geometry, lithology, and water levels in the Anza–Terwilliger area—2013, Riverside and San Diego Counties, California

The population of the Anza–Terwilliger area relies solely on groundwater pumped from the alluvial deposits and surrounding bedrock formations for water supply. The size, characteristics, and current conditions of the aquifer system in the Anza–Terwilliger area are poorly understood, however. In response to these concerns, the U.S. Geological Survey, in cooperation with the High Country Conservancy and Rancho California Water District, undertook a study to (1) improve mapping of groundwater basin geometry and lithology and (2) to resume groundwater-level monitoring last done during 2004–07 in the Anza–Terwilliger area. Inversion of gravity data, including new data collected for this study, was done to estimate the thickness of the alluvial deposits that form the Cahuilla and Terwilliger groundwater basins and to understand the geometry of the underlying basement complex. After processing of the gravity data, the thickness of the alluvial aquifer materials was modeled by using all available lithology, density, and geophysical data. The thickest alluvial deposits (greater than 500 feet) are in the northern part of the study area along the south side of the San Jacinto fault zone, in the southern part of the Cahuilla groundwater basin, and in the western part of the Terwilliger groundwater basin. Through most of the area of alluvial materials, the thickness of the alluvium estimated from gravity data is less than 400 feet. Analysis of more than 900 drillers’ logs indicated that in areas having relatively thick alluvium, particularly along the San Jacinto fault zone and in the Terwilliger Valley, the alluvium is predominantly composed of sands and gravels. Fine-textured sediments appeared to be discontinuous rather than forming laterally extensive, low-permeability layers. More than 500 drillers’ logs indicated only bedrock is present, indicating that the fractured bedrock is an important source of groundwater, primarily for domestic use, in the study area. The depths of the holes drilled into the bedrock indicated that fractures potentially supplying water to wells persist in the upper few hundred feet and that the permeable zone of the fractured bedrock extends to depths greater than weathered zones in the upper part of the basement complex. Water-level data were collected from 59 wells during fall 2013. These data indicated that hydraulic head did not vary substantially with well depth and that the measured water levels in bedrock and alluvium were similar. Large offsets in groundwater altitude across the San Jacinto fault zone indicated that the fault zone is a barrier to groundwater flow in the northeastern part of the Anza Valley. On the basis of data from 33 wells, water levels mostly declined between the fall of 2006 and the fall of 2013; the median decline was 5.1 feet during this period, for a median rate of decline of about 0.7 feet/year. Based on data from 40 wells, water-level changes between fall 2004 and fall 2013 were variable in magnitude and trend, but had a median decline of 2.4 feet and a median rate of decline of about 0.3 feet/ year. These differences in apparent rates of groundwater-level change highlight the value of ongoing water-level measurements to distinguish decadal, or longer term, trends in groundwater storage often associated with climatic variability and trends. Fifty-four long-term hydrographs indicated the sensitivity of groundwater levels to climatic conditions; they also showed a general decline in water levels across the study area since 1986 and, in some cases, dating back to the 1950s.

California↗

USGS2020

Population and habitat analyses for greater sage-grouse (Centrocercus urophasianus) in the bi-state distinct population segment—2018 update

Executive Summary The Bi-State Distinct Population Segment (Bi-State DPS) of greater sage-grouse ( Centrocercus urophasianus , hereinafter “sage-grouse”) represents a genetically distinct and geographically isolated population that straddles the border between Nevada and California. The primary threat to these sage-grouse populations is the expansion of single-leaf pinyon ( Pinus monophylla ) and Utah juniper ( Juniperus osteosperma ) into sagebrush ecosystems, which fragments and reduces population connectivity and survival. Other important threats include low water availability during brood-rearing, particularly during drought, and increased predation by common ravens ( Corvus corax ), a generalist predator often associated with anthropogenic resource subsidies. Although the Bi-State DPS occurs at high elevations relative to sage-grouse range-wide, changes in historical wildfire cycles and the conversion of native shrubs to invasive annual grasslands still threaten these populations. The Bi-State DPS has undergone multiple federal status assessments and associated litigation. For example, in October of 2013, the Bi-State DPS was proposed for listing as threatened under the Endangered Species Act of 1973 by the U.S. Fish and Wildlife Service (USFWS), then withdrawn in April 2015. The withdrawal decision was challenged, and in May 2018, a Federal district court ordered the withdrawal decision to be vacated, and USFWS was required to re-open the October 2013 listing evaluation. In response, the U.S. Geological Survey (USGS), with State and Federal collaborators, embarked on a multipronged analysis to provide current and best available science regarding population status of sage-grouse within the Bi-State DPS. Using data from a long-term monitoring program, we carried out four analytical study objectives, and here, we provide preliminary results of these analyses. First, we used integrated population modeling (IPM) to predict annual population abundance and annual finite rate of population change for the Bi-State DPS, as a whole, and for each subpopulation between 1995 and 2018. Because sage-grouse exhibit population cycles (periodic increases and decreases in abundance across approximately 6- to 10-year wavelengths), we estimated trends across three nested temporal scales that represent one (11 years), two (18 years), and three (24 years) complete population cycles. These estimates of relatively long-term averaged population change account for temporal (that is, interannual) variation. Our model predicted population abundance for the Bi-State DPS during 2018 at 3,305 individuals (2,247–4,683), with the majority occupying Bodie Hills and Long Valley. The model also predicted cyclic dynamics in abundance through time with evidence of 24-year population growth and slight trends of decline over the past 18 years. Specifically, across the Bi-State DPS as a whole, we estimated annual average at 0.99, 0.99, and 1.02 over the one, two, and three population cycles, which equated to a 10.5 percent, 16.6 percent decrease, and 60.0 percent increase in abundance over the 11-, 18-, and 24-year cycles. Estimated abundance in 2018 had not reached numbers lower than those predicted during 1995. However, we found spatial variation in population trends across the three cycles. Bodie Hills subpopulation comprised the greatest (1,521) and exhibited average annual greater than 1.0 across all periods resulting in average annual increases of 7 percent. This relatively large subpopulation has grown 5 times larger than what was predicted in 1995 while experiencing cyclical dynamics within that period. Conversely, other smaller subpopulations within the Bi-State DPS exhibited average annual equal to or less than 1.0 resulting in estimated 10-year risks of extirpation ranging from 2.0 to 76.1 percent. In general, evidence of decline among smaller subpopulations was greatest for the most recent period (2008–18) compared to a period that encompassed three full population cycles (24-year). This difference coincides with an intense period of drought that began in 2012. For comparative purposes as part of this first objective, we conducted a similar analysis for populations of sage-grouse within Nevada and California but outside the Bi-State DPS. We developed a region-wide and distance-weighted IPM using lek count from Nevada Department of Wildlife (NDOW) and California Department of Fish and Wildlife (CDFW) databases and with telemetry data collected by USGS across 12 sage-grouse subpopulations. Our models predicted similar patterns in population cycling outside the Bi-State DPS but with much stronger evidence of long-term declines across 24 years. Specifically, median averaged across each year of the 11-, 18-, and 24-year periods resulted in average annual values of 0.94, 0.97, and 0.99, respectively. These values equate to 41.0 percent, 38.5 percent, and 21.3 percent declines over the corresponding periods. Second, we used lek count data in a state-space modeling framework to compare trends in population abundance across different spatial scales (that is, leks versus Bi-State DPS). This hierarchical framework allowed us to disentangle declines associated with climate conditions as opposed to other local level factors that might signal the need for management intervention. Specifically, we identified 7 leks that were both declining and recently decoupled from larger spatial scale trends, typically governed by climatic conditions (referred to as soft or hard signals). The goal of this analysis was to provide an early warning system that might have implications for conservation actions at local scales. Third, we developed phenological (spring, summer–fall, and winter) and reproductive life stage (nesting, early brood-rearing, and late-brood rearing) based resource selection functions using various environmental covariates. We report rankings of variable importance for each season and life stage, developed maps of habitat selection indices (HSI), binned categories representing low, moderate, and high classes of quality (where any category greater than or equal to low indicated selected habitat) for each phenological season and life stage, and produced composite maps by selected phenological and reproductive stage to estimate annual habitat. Fourth, we used for each lek within the Bi-State DPS to carry out a spatial analysis that quantified substantial changes in the distribution of occupied habitat across long- (24-year) and short- (11-year) term periods. Owing to differences among available datasets, the long-term analysis primarily reflected spatial shifts among subpopulations comprising the majority of the Bi-State DPS (that is, Bodie Hills and Long Valley) while the short-term analysis also quantified changes among subpopulations along the periphery. Over long and short-term periods, the overall distribution of occupied habitat (as measured by 99 percent utilization distributions intersecting any quantified habitat) was reduced by 20,573 ha and 55,492 ha, respectively. Occupied core areas (as measured by 50 percent utilization distributions intersecting any quantified habitat) over long-term periods were solely located in Bodie Hills and Long Valley. Although nearly all subpopulations experienced contractions in occupied overall and core distribution, Bodie Hills experienced spatial expansion that occurred with concomitant spatial contraction at Long Valley over both periods. Subpopulations at the northern (Pine Nuts), central (Sagehen) and southern (White-Mountains) extents of the Bi-State DPS also experienced spatial contraction over the short-term period. These findings, coupled with those of population trends, indicate long-term patterns in redistribution of sage-grouse from Long Valley and peripheral subpopulations to Bodie Hills. That is, sage-grouse subpopulations at the periphery are declining while the largest population at the core is increasing, which could have meaningful impacts on overall metapopulation persistence. We provide evidence for loss of occupied habitat (reduced distribution) given local extirpation of subpopulations. Fifth, we calculated percentages of selected phenological, life stage, and annual habitat that each subpopulation contributed to the Bi-State DPS. We then intersected these maps with a composite estimate of occupied habitat from the fourth objective and calculated percentages of selected habitat likely occupied by sage-grouse that each subpopulation contributed to the Bi-State DPS. These values provide evidence for loss of occupied habitat and subsequent reductions in spatial distribution given reductions in abundance and, in some cases, extirpation of leks within subpopulations. Lastly, we carried out an initial in-depth analysis of selection for irrigated pastures and wet meadows during the brood-rearing stage for the Long Valley subpopulation. We chose this subpopulation because it represents a population core, representing 26.5 percent of total sage-grouse within the Bi-State DPS, and has exhibited long-term declines in abundance and distribution. This subpopulation is highly sensitive to precipitation and other factors that influence water availability. Models predicted higher use of the interior portions of irrigated pastures and wet meadows during late brood-rearing period, which represented a potentially risky use of habitat that was exacerbated during periods of low moisture (for example, drought, reduced water delivery, or both). Sage-grouse typically used edges of riparian areas and pastures, largely because the interior of these mesic areas consisted of considerably less overhead concealment cover (for example, shrubs) that likely resulted in a higher risk of mortality. We found that a lack of water delivery to pastures in the form of overwinter precipitation or diversion ditches increased the movements of sage-grouse to the interior of pastures. Although further investigation of water delivery impacts on chick survival are warrented, our initial findings regarding resource selection may explain recent declines in population growth at Long Valley.

California, Nevada↗

USGS2007

Hydrology of Polk County, Florida

Local water managers usually rely on information produced at the State and regional scale to make water-resource management decisions. Current assessments of hydrologic and water-quality conditions in Polk County, Florida, commonly end at the boundaries of two water management districts (South Florida Water Management District and the Southwest Florida Water Management District), which makes it difficult for managers to determine conditions throughout the county. The last comprehensive water-resources assessment of Polk County was published almost 40 years ago. To address the need for current countywide information, the U.S. Geological Survey began a 3?-year study in 2002 to update information about hydrologic and water-quality conditions in Polk County and identify changes that have occurred. Ground-water use in Polk County has decreased substantially since 1965. In 1965, total ground-water withdrawals in the county were about 350 million gallons per day. In 2002, withdrawals totaled about 285 million gallons per day, of which nearly 95 percent was from the Floridan aquifer system. Water-conservation practices mainly related to the phosphate-mining industry as well as the decrease in the number of mines in operation in Polk County have reduced total water use by about 65 million gallons per day since 1965. Polk County is underlain by three principal hydrogeologic units. The uppermost water-bearing unit is the surficial aquifer system, which is unconfined and composed primarily of clastic deposits. The surficial aquifer system is underlain by the intermediate confining unit, which grades into the intermediate aquifer system and consists of up to two water-bearing zones composed of interbedded clastic and carbonate rocks. The lowermost hydrogeologic unit is the Floridan aquifer system. The Floridan aquifer system, a thick sequence of permeable limestone and dolostone, consists of the Upper Floridan aquifer, a middle semiconfining unit, a middle confining unit, and the Lower Floridan aquifer. The Upper Floridan aquifer provides most of the water required to meet demand in Polk County. Data from about 300 geophysical and geologic logs were used to construct hydrogeologic maps showing the tops and thicknesses of the aquifers and confining units within Polk County. Thickness of the surficial aquifer system ranges from several feet thick or less in the extreme northwestern part of the county and along parts of the Peace River south of Bartow to more than 200 feet along the southern part of the Lake Wales Ridge in eastern Polk County. Thickness of the intermediate aquifer system/intermediate confining unit is highly variable throughout the county because of past erosional processes and sinkhole formation. Thickness of the unit ranges from less than 25 feet in the extreme northwestern part of the county to more than 300 feet in southwestern Polk County. The altitude of the top of the Upper Floridan aquifer in the county ranges from about 50 feet above National Geodetic Vertical Datum of 1929 (NGVD 29) in the northwestern part to more than 250 feet below NGVD 29 in the southern part. Water levels in the Upper Floridan aquifer fluctuate seasonally, increasing during the wet season (June through September) and decreasing during the rest of the year. Water levels in the Upper Floridan aquifer also can change from year to year, depending on such factors as pumpage and climatic variations. In the southwestern part of the county, fluctuations in water use related to phosphate mining have had a major impact on ground-water levels. Hydrographs of selected wells in southwestern Polk County show a general decline in water levels that ended in the mid-1970s. This water-level decline coincides with an increase in water use associated with phosphate mining. A substantial increase in water levels that began in the mid-1970s coincides with a period of decreasing water use in the county. Despite reductions in water use since 1970, howev

Scientific Investigations Report↗

USGS2002

Trace-element deposition in the Cariaco Basin, Venezuela Shelf, under sulfate-reducing conditions: A history of the local hydrography and global climate, 20 ka to the present

A sediment core from the Cariaco Basin on the Venezuelan continental shelf, which recovered sediment that has been dated back to 20 ka (thousand years ago), was examined for its major-element-oxide and trace-element composition. Cadmium (Cd), chromium (Cr), copper (Cu), molybdenum (Mo), nickel (Ni), vanadium (V), and zinc (Zn) can be partitioned between a siliciclastic, terrigenous-derived fraction and two seawater-derived fractions. The two marine fractions are (1) a biogenic fraction represented by nutrient trace elements taken up mostly in the photic zone by phytoplankton, and (2) a hydrogenous fraction that has been derived from bottom water via adsorption and precipitation reactions. This suite of trace elements contrasts with a second suite of trace elements—barium (Ba), cobalt (Co), gallium (Ga), lithium (Li), the rare-earth elements, thorium (Th), yttrium (Y), and several of the major-element oxides—that has had solely a terrigenous source. The partitioning scheme, coupled with bulk sediment accumulation rates measured by others, allows us to determine the accumulation rate of trace elements in each of the three sediment fractions and of the fractions themselves. The current export of organic matter from the photic zone, redox conditions and advection of bottom water, and flux of terrigenous debris into the basin can be used to calculate independently trace-element depositional rates. The calculated rates show excellent agreement with the measured rates of the surface sediment. This agreement supports a model of trace-element accumulation rates in the subsurface sediment that gives a 20-kyr history of upwelling into the photic zone (that is, primary productivity), bottom-water advection and redox, and provenance. Correspondence of extrema in the geochemical signals with global changes in sea level and climate demonstrates the high degree to which the basin hydrography and provenance have responded to the paleoceanographic and paleoclimatic regimes of the last 20 kyr. The accumulation rate of the marine fraction of Mo increased abruptly at about 14.8 ka (calendar years), from less than 0.5 µg cm -2 yr -1 to greater than 4 µg cm -2 yr -1 . Its accumulation rate remained high but variable until 8.6 ka, when it decreased sharply to 1 µg cm -2 yr -1 . It continued to decrease to 4.0 ka, to its lowest value for the past 15 kyr, before gradually increasing to the present. Between 14.8 ka and 8.6 ka, its accumulation rate exhibited strong maxima at 14.4, 13.0, and 9.9 ka. The oldest maximum corresponds to melt-water pulse IA into the Gulf of Mexico. A relative minimum, centered at about 11.1 ka, corresponds to melt-water pulse IB; a strong maximum occurs in the immediately overlying sediment. The maximum at 13.0 ka corresponds to onset of the Younger Dryas cold event. This pattern to the accumulation rate of Mo (and V) can be interpreted in terms of its deposition from bottom water of the basin, the hydrogenous fraction, under SO 4 2- -reducing conditions, during times of intense bottom-water advection 14.8 ka to 11.1 ka and significantly less intense bottom-water advection 11 ka to the present. The accumulation rate of Cd shows a pattern that is only slightly different from that of Mo, although its deposition was determined largely by the rain rate of organic matter into the bottom water, a biogenic fraction whose deposition was driven by upwelling of nutrient-enriched water into the photic zone. Its accumulation exhibits only moderately high rates, on average, during both melt-water pulses. Its highest rate, and that of upwelling, occurred during the Younger Dryas, and again following melt-water pulse IB. The marine fractions of Cu, Ni, and Zn also have a strong biogenic signal. The siliciclastic terrigenous debris, however, represents the dominant source, and host, of Cu, Ni, and Zn. All four trace elements have a consid-erably weaker hydrogenous signal than biogenic signal. Accumulation rates of the terrigenous fraction, as reflected by accumulation rates of Th and Ga, show strong maxima at 16.2 and 12.7 ka and minima at 14.1 and 11.1 ka. Co, Li, REE, and Y have a similar distribution. The minima occurred during melt-water pulses IA and IB, the maxima during the Younger Dryas and the rise in sea level following the last glacial maximum.

Cariaco Basin↗

USGS2018

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗

USGS2019

The effects of management practices on grassland birds — An introduction to North American grasslands and the practices used to manage grasslands and grassland birds

The Great Plains of North America is defined as the land mass that encompasses the entire central portion of the North American continent that, at the time of European settlement, was an unbroken expanse of primarily herbaceous vegetation. The Great Plains extend from central Saskatchewan and Alberta to central Mexico and from Indiana to the Rocky Mountains. The expanses of herbaceous vegetation are often referred to as native prairie or native grasslands. Native grasslands share the characteristics of a general uniformity in vegetation structure, dominance by grasses and forbs, a near absence of trees and shrubs, annual precipitation ranging from 25 to 100 centimeters, extreme intra-annual fluctuations in temperature and precipitation, and a flat to rolling topography over which fires can spread. To the west of the Great Plains lie the sagebrush communities of the Great Basin, which extend from British Columbia and Saskatchewan to northern Arizona and New Mexico and from the eastern slopes of the Sierra Nevada and Cascade mountain ranges to western South Dakota. Sagebrush communities share similar characteristics to native grasslands, but their location east of the Rocky Mountains creates a more moderating influence from prevailing westerly winds that affect timing of peak precipitation and growth form of dominant vegetation. Native grasslands and sagebrush communities harbor a diverse array of grassland, wetland, and woodland plant and animal communities that are uniquely adapted to the natural forces of the Great Plains and Great Basin, namely the interactive forces of climate, fire, and grazing. The arrival of European settlers to North America brought profound change to native grassland and sagebrush communities, including the establishment of permanent towns and cities, the proliferation of cropland-based agricultural systems, and the suppression of wildfires. The near extirpation of bison by the 1860s paved the way for dramatic changes in the dominant grazers and a shift in the disturbance patterns that historically influenced vegetation structure. The greatest threat to native grasslands and sagebrush communities in modern times is their loss due to conversion to rowcrop agriculture and to urbanization. Concomitant with habitat loss is a precipitous decline in populations of bird species that evolved with, and are uniquely adapted to, the native grassland and sagebrush habitats. Avian population trends are linked strongly to agricultural land use. Besides outright loss of suitable breeding habitat, agricultural practices affect birds through factors such as pesticide exposure, habitat fragmentation, shifts in predator community composition, and occurrence of brood parasites. Bird populations face other stressors, such as loss of habitat to and behavioral avoidance of urbanized areas, roads, and infrastructure associated with energy production. Despite the many anthropogenic changes to North American grassland and sagebrush communities, some bird species are adaptable and opportunistic in their habitat selection and now utilize one or more human-created habitats. Human-created habitats include pastures, hayfields, agricultural terraces, crop buffer strips, field borders, grassed waterways, fencerows, road rights-of-way, airports, reclaimed coal mines, and planted wildlife cover. Fields of seeded grasslands enrolled in Federal long-term set-aside programs, such as the Conservation Reserve Program in the United States and the Permanent Cover Program in Canada, provide important nesting habitat for grassland bird species. The array of habitats used by birds makes habitat and avian management a complex undertaking, and the scale (for example, local, regional, international) at which management actions can be implemented are such that a universal approach to managing grasslands for the conservation of the entire suite of bird species does not exist. Experienced land managers recognize that it is impossible to manage for all bird species simultaneously, and thus, prioritization is necessary towards those habitats or bird species that the manager or management agency ranks highest for a specific region or management unit. The primary tools available for management are burning, grazing, mowing, herbicide application, and idling, but before choosing a particular practice, a manager will want to consider issues of seasonality, intensity, and frequency. Despite the thousands of studies that are cited in this compendium, much remains unknown about the effects of management practices on bird species. The series of species accounts in this compendium review the current state of knowledge regarding management of grassland and sagebrush bird species and summarize information on the effects of management practices on individual species. The accounts do not give definitive statements on the effects of management practices for any particular species, primarily because there are very few replicated studies in which identical management practices have been applied in the same geographical area with consistent results, which are elements necessary to provide concrete recommendations for the management of a particular species in a particular area. Documentation of the effects of management treatments on individual species through statistically sound methods that incorporate multiple years and locations will further scientists’ and land managers’ knowledge far more than 1–2-year studies that are limited in scope as well as time, but studies of that scope and breadth are rare.

Professional Paper↗

USGS1996

Surface-water hydrology and runoff simulations for three basins in Pierce County, Washington

The surface-water hydrology in Clear, Clarks, and Clover Creek Basins in central Pierce County, Washington, is described with a conceptual model of the runoff processes and then simulated with the Hydrological Simulation Program-FORTRAN (HSPF), a continuous, deterministic hydrologic model. The study area is currently undergoing a rapid conversion of rural, undeveloped land to urban and suburban land that often changes the flow characteristics of the streams that drain these lands. The complex interactions of land cover, climate, soils, topography, channel characteristics, and ground- water flow patterns determine the surface-water hydrology of the study area and require a complex numerical model to assess the impact of urbanization on streamflows. The U.S. Geological Survey completed this investigation in cooperation with the Storm Drainage and Surface Water Management Utility within the Pierce County Department of Public Works to describe the important rainfall-runoff processes within the study area and to develop a simulation model to be used as a tool to predict changes in runoff characteristics resulting from changes in land use. The conceptual model, a qualitative representation of the study basins, links the physical characteristics to the runoff process of the study basins. The model incorporates 11 generalizations identified by the investigation, eight of which describe runoff from hillslopes, and three that account for the effects of channel characteristics and ground-water flow patterns on runoff. Stream discharge was measured at 28 sites and precipitation was measured at six sites for 3 years in two overlapping phases during the period of October 1989 through September 1992 to calibrate and validate the simulation model. Comparison of rainfall data from October 1989 through September 1992 shows the data-collection period beginning with 2 wet water years followed by the relatively dry 1992 water year. Runoff was simulated with two basin models-the Clover Creek Basin model and the Clear-Clarks Basin model-by incorporating the generalizations of the conceptual model into the construction of two HSPF numerical models. Initially, the process-related parameters for runoff from glacial-till hillslopes were calibrated with numerical models for three catchment sites and one headwater basin where streamflows were continuously measured and little or no influence from ground water, channel storage, or channel losses affected runoff. At one of the catchments soil moisture was monitored and compared with simulated soil moisture. The values for these parameters were used in the basin models. Basin models were calibrated to the first year of observed streamflow data by adjusting other parameters in the numerical model that simulated channel losses, simulated channel storage in a few of the reaches in the headwaters and in the floodplain of the main stem of Clover Creek, and simulated volume and outflow of the ground-water reservoir representing the regional ground-water aquifers. The models were run for a second year without any adjustments, and simulated results were compared with observed results as a measure of validation of the models. The investigation showed the importance of defining the ground-water flow boundaries and demonstrated a simple method of simulating the influence of the regional ground-water aquifer on streamflows. In the Clover Creek Basin model, ground-water flow boundaries were used to define subbasins containing mostly glacial outwash soils and not containing any surface drainage channels. In the Clear-Clarks Basin model, ground-water flow boundaries outlined a recharge area outside the surface-water boundaries of the basin that was incorporated into the model in order to provide sufficient water to balance simulated ground-water outflows to the creeks. A simulated ground-water reservoir used to represent regional ground-water flow processes successfully provided the proper water balance of inflows and outfl

Washington↗

USGS2011

A water-budget model and assessment of groundwater recharge for the Island of Hawai'i

Concern surrounding increasing demand for groundwater on the Island of Hawaiʻi, caused by a growing population and an increasing reliance on groundwater as a source for municipal and private water systems, has prompted a study of groundwater recharge on the island using the most current data and accepted methods. For this study, a daily water-budget model for the entire Island of Hawaiʻi was developed and used to estimate mean recharge for various land-cover and rainfall conditions, and a submodel for the Kona area was developed and used to estimate historical groundwater recharge in the Kona area during the period 1984–2008. Estimated mean annual recharge on the Island of Hawaiʻi is 6,594 million gallons per day, which is about 49 percent of mean annual rainfall. Recharge is highest on the windward slopes of Mauna Loa, below the tradewind inversion, and lowest on the leeward slopes of Kohala and Mauna Kea. Local recharge maxima also occur on (1) windward Kohala, with the exception of the northern tip, (2) windward Mauna Kea below the tradewind inversion, (3) windward Kīlauea, (4) the middle elevations of southeastern Mauna Loa, and (5) the lower-middle elevations of leeward Mauna Loa and southwestern Hualālai, in the Kona area. Local recharge minima also occur on (1) Mauna Kea and Mauna Loa, above the tradewind inversion, (2) the northern tip of Kohala, (3) leeward Kīlauea, (4) the southern tip of Mauna Loa, and (5) the northwestern slopes of Mauna Loa and Hualālai. In 18 of the 24 aquifer systems on the island, estimated mean annual recharge for baseline conditions was higher than the recharge estimates used in the 2008 State of Hawaiʻi Water Resource Protection Plan (2008 WRPP). Baseline conditions for this study were 2008 land cover and mean annual rainfall from the period 1916–1983. Estimates of recharge for the Māhukona, Waimea, and Hāwī aquifer systems, however, were between 29 and 38 percent lower than the 2008 WRPP estimates, mainly because of much higher evapotranspiration estimates in this study compared to the 2008 WRPP. For the drought simulation (1991–95 rainfall), the estimates of recharge for these three aquifer systems were only 15 to 33 percent of the sustainable yields (maximum allowable pumping rates) set by the 2008 WRPP. This may be cause for concern, as these areas are experiencing a rapid growth in development and a related growth in water demand. Recent projections of change in rainfall owing to effects of ongoing climate change generally indicate a slight increase in islandwide rainfall, and estimates of annual recharge in the late 21st century are higher than baseline estimates for every aquifer system, except ʻAnaehoʻomalu. On average, these aquifer-system recharge estimates are higher by about 8 percent compared to baseline estimates. In the Kona area, estimated groundwater recharge during the period 1984–2008 was highest during 2004–8 and lowest during 1999–2003, with the 1999–2003 recharge being about 50 percent of the 2004–8 recharge. These extremes in recharge coincided with the periods of lowest and highest mean rainfall, respectively. No seasonal pattern in recharge is discernible. Spatially, the highest recharge occurred in a belt about 4 miles wide running parallel to the coast about 2 miles inland. The sensitivity of recharge estimates to input parameters is related to the climate and land-cover conditions of the particular area of study. For the wet, forested areas characteristic of the windward side of the island, recharge was most sensitive to the ratio of runoff to rainfall. For the dry, grassland areas characteristic of the northwestern leeward side of the island, recharge was most sensitive to root depth. For the Kona area, characterized by moderate rainfall and a wide variety of land cover, recharge was most sensitive to the pan coefficient and canopy-evaporation rates in

Hawai'i↗

USGS2020

Assessment of methods for soil monitoring in the Adirondack region of New York

Repeated sampling to detect changes in forest soils was rarely used before 1990, but the value of soil monitoring in understanding environmental change is becoming well established. The growing number of resampling studies has shown that sampling designs and procedures must be adapted to the objectives of the monitoring program and the soils being monitored. In the Adirondack region, current priorities include the response of soils to large increases, and more recently, large decreases in acidic deposition, and changes driven by trending climate such as altered pools of soil organic carbon, as well as other unforeseen factors that will occur in the future. To improve methods and assess the feasibility of long-term soil monitoring in the Adirondack region, the United States Geological Survey (USGS) conducted a pilot project to evaluate a new sampling method for characterizing soils on a watershed basis. Results obtained with this new approach, referred to as the ADK sampling method, was compared to methods used in previous sampling conducted in 2004 as part of the Western Adirondack Stream Survey (WASS), and also to previous high-replication pit sampling in the North and South Tributary watersheds of Buck Creek (North Buck and South Buck). The number of sampling locations and spatial distribution of sampling points within watersheds differed among the methods, although pit excavation was used to obtain samples in all cases. In addition, this investigation evaluated the use of small diameter corers as a means to measure forest floor mass with greater accuracy and precision than commonly used methods such as pit excavation. Sufficient statistical power to detect ecologically relevant changes in upper profile horizons (Oe, Oa and upper 10 cm of the B) were achieved with the ADK sampling method that utilized 18 pit excavations per watershed. The sampling locations were organized within each watershed into three study areas (six sampling locations per study area) that represented the primary types of landscape within the watershed. Sampling at 18 locations per watershed was found to be nearly as effective at detecting changes as sampling at 28 locations per watershed. Numerous significant changes (P < 0.10) were detected with both 18 and 28 sampling locations at sampling intervals of 12 to 16 years. The relationship between soil data obtained with the ADK method and stream chemistry at the base of the watershed suggested that this approach adequately characterized soil variability within the watershed for the purpose of studying soil-stream linkages. Significant changes in upper B horizon calcium (P < 0.10) and Oa horizon aluminum (P < 0.01) were detected when the data from the four WASS watersheds were combined with the two Buck Creek watersheds, which suggested that there would be value in resampling other WASS watersheds previously sampled in 2004 to support a regional assessment. Study results support small diameter cores as a useful method to monitor changes in the organic matter mass of the forest floor. This method showed high reproducibility in repeated sampling tests and lower spatial variability in sample data than traditional approaches when compared on a watershed basis. Soil coring is also faster and requires less equipment than pit excavation methods, which makes it more conducive to sampling over large areas. However, organic matter mass of the forest floor determined by coring was consistently less than the values obtained by the ADK sampling method that used pit sampling and vertical horizon measurements, and also literature values of a previous Adirondack study that utilized pit sampling in which the entire horizon was collected over a measured area. However, a high correlation (R2 = 0.87) occurred between organic matter content (expressed as Mg ha-1) determined by coring and the ADK sampling method. Differing methods with regard to where sample could be collected, and how organic matter was collected for chemical analysis were the likely reasons for differences in quantification of forest floor organic mass. Collection of forest floor cores in conjunction with the ADK method is recommended to provide improved sensitivity in detecting changes in the forest floor in proximity of where full analyses of the soil profile are being done. This duel sampling approach represents an optimized method for measuring and understanding how Adirondack soils will change in the future.

New York↗

USGS2004

A precipitation-runoff model for the analysis of the effects of water withdrawals and land-use change on streamflow in the Usquepaug–Queen River Basin, Rhode Island

The 36.1-square-mile Usquepaug–Queen River Basin in south-central Rhode Island is an important water resource. Streamflow records indicate that withdrawals may have diminished flows enough to affect aquatic habitat. Concern over the effect of withdrawals on streamflow and aquatic habitat prompted the development of a Hydrologic Simulation Program–FORTRAN (HSPF) model to evaluate the water-management alternatives and land-use change in the basin. Climate, streamflow, and water-use data were collected to support the model development. A logistic-regression equation was developed for long-term simulations to predict the likelihood of irrigation, the primary water use in the basin, from antecedent potential evapotranspiration and precipitation for generating irrigation demands. The HSPF model represented the basin by 13 pervious-area and 2 impervious-area land-use segments and 20 stream reaches. The model was calibrated to the period January 1, 2000 to September 30, 2001, at three continuous streamflow-gaging stations that monitor flow from 10, 54, and 100 percent of the basin drainage area. Hydrographs and flow-duration curves of observed and simulated discharges, along with statistics compiled for various model-fit metrics, indicate a satisfactory model performance. The calibrated HSPF model was modified to evaluate streamflow (1) under no withdrawals to streamflow under current (2000–01) withdrawal conditions under long-term (1960–2001) climatic conditions, (2) under withdrawals by the former Ladd School water-supply wells, and (3) under fully developed land use. The effects of converting from direct-stream withdrawals to ground-water withdrawals were evaluated outside of the HSPF model by use of the STRMDEPL program, which calculates the time delayed response of ground-water withdrawals on streamflow depletion. Simulated effects of current withdrawals relative to no withdrawals indicate about a 20-percent decrease in the lowest mean daily streamflows at the basin outlet, but withdrawals have little effect on flows that are exceeded less than about 90 percent of the time. Tests of alternative model structures to evaluate model uncertainty indicate that the lowest mean daily flows ranged between 3 and 5 cubic feet per second (ft3/s) without withdrawals and 2.2 to 4 ft3/s with withdrawals. Changes in the minimum daily streamflows are more pronounced, however; at the upstream streamflow-gaging station, a minimum daily flow of 0.2 ft3/s was sustained without withdrawals, but simulations with withdrawals indicate that the reach would stop flowing part of a day about 5 percent of the time. The effect on streamflow of potential ground-water withdrawals of 0.20, 0.90, and 1.78 million gallons per day (Mgal/d) at the former Ladd School near the central part of the basin were evaluated. The lowest daily mean flows in model reach 3, the main stem of the Queen River closest to the pumped wells, decreased by about 50 percent for withdrawals of 0.20 Mgal/d (from about 0.4 to 0.2 ft3/s) in comparison to current withdrawals. Reach 3 would occasionally stop flowing during part of the day at the 0.20-Mgal/d withdrawal rate because of diurnal fluctuation in streamflow. The higher withdrawal rates (0.90 and 1.78 Mgal/d) would cause reach 3 to stop flowing about 10 to 20 percent of the time, but the effects of pumping rapidly diminished downstream because of tributary inflows. Simulation results indicate little change in the annual 1-, 7-, and 30-day low flows at the 0.20 Mgal/d pumping rate, but at the 1.78 Mgal/d pumping rate, reach 3 stopped flowing for nearly a 7-day period every year and for a 30-day period about every other year. At the 0.90 Mgal/d pumping rate, reach 3 stopped flowing about every other year for a 7-day period and about once every 5 years for a 30-day period. Land-use change was simulated by converting model hydrologic-response units (HRUs) representing undeveloped areas to HRUs representing developed areas on the basis of development suitability and town zoning. About 55 percent of the basin is suitable for development; this area would accommodate about 4,300 new low-density residential homes under current zoning. Increases in storm volume and peak flows, and decreases in base flow, typically associated with urbanization, were not evident in buildout simulations because the effective impervious area was assumed to increase by only 2 percent. Under fully developed conditions, withdrawals from self-supply wells were estimated to reach 1.2 Mgal/d. Potential increases in water withdrawals for a fully developed basin have only a minor impact on the main stem streamflow, but the effects of urbanization could be more pronounced in localized areas where development is concentrated. Streamflow-depletion rates were calculated for varying distances of a pumped irrigation well from a stream. For the irrigation rates and aquifer conditions tested, streamflow depletion, relative to the pumping rate, decreases rapidly as the pumped well was moved away from the stream. Streamflow depletion, relative to the peak withdrawal rate, decreased by about 60, 80, and 90 percent by locating the pumped well 500, 1,000, and 1,500 feet from the stream, respectively.

Rhode Island↗

USGS2022

Long-term groundwater availability in the Waihe‘e, ‘Īao, and Waikapū aquifer systems, Maui, Hawai‘i

Groundwater levels have declined since the 1940s in the Wailuku area of central Maui, Hawai‘i, on the eastern flank of West Maui volcano, mainly in response to increased groundwater withdrawals. Available data since the 1980s also indicate a thinning of the freshwater lens and an increase in chloride concentrations of pumped water from production wells. These trends, combined with projected increases in demand for groundwater in central Maui, have led to concerns over groundwater availability and have highlighted a need to improve general understanding of the hydrologic effects of proposed groundwater withdrawals in the Waihe‘e, ‘Īao, and Waikapū areas of central Maui. A numerical groundwater model was constructed to simulate the flow and salinity of groundwater in central Maui. The model simulates the effects of changes in groundwater withdrawals and recharge on water levels, freshwater-lens thicknesses, and chloride concentrations of pumped water from production wells. The model incorporates updated water-budget estimates of groundwater recharge from infiltration and direct recharge, seepage in stream channels, and inflow from inland areas. Mean annual groundwater recharge from infiltration and direct recharge was estimated using a daily water-budget model and the most current data, including the distributions of monthly rainfall and potential evapotranspiration, for the study area for nine historical periods from 1926 through 2012: 1926–69, 1970–79, 1980–84, 1985–89, 1990–94, 1995–99, 2000–04, 2005–09, and 2010–12. The water-budget model also estimated groundwater recharge based on one hypothetical scenario that used 1980–2010 rainfall and 2017 land cover. For the nine historical periods, estimated recharge from infiltration and direct recharge within the area of the groundwater model ranged from 30.4 million gallons per day (Mgal/d) during 2010–12 to 98.7 Mgal/d during 1926–69. Variability in recharge during these periods mainly reflects changes in rainfall and irrigation over time. Between 2010 and 2014, streamflow restoration in previously diverted streams resulted in an estimated increase in recharge from seepage in stream channels of about 12.5 Mgal/d. Average groundwater inflow of about 39.6 Mgal/d from inland, dike-intruded areas to the main area of interest was estimated from an existing island-wide numerical groundwater-flow model, which is at a larger scale and incorporates a greater number of simplifying assumptions. The numerical groundwater model developed for this study was calibrated to 1926–2012 transient water levels, vertical salinity profiles, and chloride concentrations of water pumped by production wells in the study area. The model was then used to evaluate one future recharge and six selected withdrawal scenarios, developed in consultation with the Maui Department of Water Supply, in terms of long-term changes in water level and 50-percent ocean-water salinity surface. The groundwater model was also used to simulate the future salinity of water withdrawn by existing and proposed production wells. The simulations were run to steady-state conditions, providing an estimate of the long-term effects of changes in withdrawal and recharge on the groundwater resource. Results of the simulated future withdrawal scenarios indicate that, relative to 2017–18 rates, the scenarios’ long-term effect of increased withdrawals ultimately leads to lower water levels and a higher 50-percent ocean-water salinity surface indicating a thinning of the freshwater lens. Results also indicate that the increased withdrawals produce some groundwater with chloride concentration below 250 milligrams per liter and some groundwater with higher chloride concentration. The amount of drawdown near production wells and the quality of water withdrawn from production wells is dependent on the rate and spatial distribution of the withdrawals. The model was also used to evaluate how groundwater availability may be affected for a drier recharge scenario based on a published study of future climate. Model results of the future recharge scenario indicate that the rate of groundwater recharge is a controlling factor for (1) water levels, (2) the 50-percent ocean-water salinity surface, and (3) the quality of water withdrawn from production wells in the Wailuku area. Coupled with reduced groundwater recharge (with all other factors remaining equal), the modeled future withdrawals in the scenario would tend to cause lower water levels, a higher 50-percent ocean-water salinity surface, and increased salinity of water withdrawn from production wells. The three-dimensional numerical groundwater model developed for this study utilizes the latest available hydrologic and geologic information and is a useful tool for understanding the long-term hydrologic effects of additional groundwater withdrawals in central Maui. The model has several limitations, including its non-uniqueness and inability to account for local-scale heterogeneities. Short-term effects of changes in recharge and withdrawals—and optimization of pumping rates to meet increased demand for water with acceptable salinity—are possible conditions for future simulation analyses.

Hawaii↗

USGS2024

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York↗