USGS Science⌕ Search

SEARCH · USGS Science

Results for “Nature Climate Change”

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

Preliminary photointerpretation map of landslide and other surficial deposits of the Mount Hamilton quadrangle and parts of the Mount Boardman and San Jose quadrangles, Alameda and Santa Clara Counties, California

The nine San Francisco Bay region counties lie within a geologically active, young, and dynamic part of the central and northern Coast Ranges of California. Significant movements of the earth's crust are occurring here at the present time, posing numerous problems to urbanization, including some of special concern. Geological processes such as fault movements, earthquakes, land subsidence, landsliding, slow downslope movement of bedrock and surficial materials, coastal and stream erosion, flooding, and sedimentation are all potentially hazardous. Because of these factors, an understanding of the operation of physical processes in the bay region is desirable for harmonious, efficient, and safe land-use planning, particularly now, with greatly expanded pressures for urban growth. This map presents preliminary information about one aspect of the physical environment necessary to sound land-use planning--the nature and distribution of surficial deposits. Because surficial deposits are common and well developed in much of the bay region, it is useful to know how and why they have formed, as well as what properties they possess. When maps like this are used in combination with other types of environmental information, such as data on soils, bedrock geology, slopes, vegetation, climatic variation, seismic response, and hydrology, it should be easier to arrive at sound decisions regarding the physical aspects of land use. The U.S. Geological Survey is studying many of these factors in the bay region and hopes to provide the community with much of the required information as part of its San Francisco Bay Region Study in cooperation with the Department of Housing and Urban Development. The representation of surficial deposits on this map reflects the way in which a geologist, working exclusively with aerial photographs, interpreted the origin of various elements of the present landscape. The deposits shown here have not been examined in the field. However, by viewing overlapping vertical aerial photographs through a stereoscope, the geologist sees a three-dimensional relief model of the ground surface and can study and interpret the origins of landforms with considerable ease. In fact, for mapping surficial deposits, particularly in reconnaissance-type studies, photointerpretation has advantages over both ground observations and laboratory studies of surficial materials. Of course, better information can be provided when all aspects of the study are integrated. These preliminary photointerpretation maps are only the first stage in a detailed study of surficial deposits in the bay region, but they should provide land-use planners with immediately useful information about the regional distribution of landslide and other surficial deposits. This map indicates the dominant surficial processes that have probably been operative by showing the distribution of different types of surficial deposits. Processes such as weathering, erosion, sedimentation, and the slow as well as rapid downslope movement of earth materials have constantly reshaped the land surface in the past and will continue to in the future, although at varying rates. The processes are interrelated to varying degrees. For example, crustal uplift of the Coast Ranges will lead to increased erosion and downcutting by streams that in turn generally results in increased deposition of sediments in river valleys, lakes, and shoreline areas. Older flood planes and river deposits may be eroded, leaving elevated terrace deposits. In addition, downcutting by streams may cause adjacent slopes to become unstable, thereby increasing the possibility of slope failures. Man's activities can alter natural physical processes in many ways. Simple acts such as overwatering a lawn or placing a septic tank drainfield in ground that is marginally stable may weaken the bedrock and surficial materials enough to induce landsliding. Relatively stable areas may be made unstable as a result of construction activities that involve cutting or oversteepening of natural slopes. Engineers, builders, conservationists, and others concerned with land use must evaluate the potential effects of all types of development, and maps that show the nature and distribution of surficial deposits should provide much of the basic information they need. This map, then, shows the cumulative effects of various processes that have yielded surficial deposits up to the time the photographs used for photointerpretation were taken. It does not indicate directly areas where processes will be most active, nor does it show the rate at which they will operate. However, knowledge of the history of geologic events is a key to understanding and predicting the evolution of an area, even where man's activities significantly change the character of the land. Almost all new landslides, for example, occur in areas with a history of landslide activity.

California↗

Precise interpolar phasing of abrupt climate change during the last ice age

The last glacial period exhibited abrupt Dansgaard–Oeschger climatic oscillations, evidence of which is preserved in a variety of Northern Hemisphere palaeoclimate archives 1 . Ice cores show that Antarctica cooled during the warm phases of the Greenland Dansgaard–Oeschger cycle and vice versa 2 , 3 , suggesting an interhemispheric redistribution of heat through a mechanism called the bipolar seesaw 4 , 5 , 6 . Variations in the Atlantic meridional overturning circulation (AMOC) strength are thought to have been important, but much uncertainty remains regarding the dynamics and trigger of these abrupt events 7 , 8 , 9 . Key information is contained in the relative phasing of hemispheric climate variations, yet the large, poorly constrained difference between gas age and ice age and the relatively low resolution of methane records from Antarctic ice cores have so far precluded methane-based synchronization at the required sub-centennial precision 2 , 3 , 10 . Here we use a recently drilled high-accumulation Antarctic ice core to show that, on average, abrupt Greenland warming leads the corresponding Antarctic cooling onset by 218 ± 92 years (2 σ ) for Dansgaard–Oeschger events, including the Bølling event; Greenland cooling leads the corresponding onset of Antarctic warming by 208 ± 96 years. Our results demonstrate a north-to-south directionality of the abrupt climatic signal, which is propagated to the Southern Hemisphere high latitudes by oceanic rather than atmospheric processes. The similar interpolar phasing of warming and cooling transitions suggests that the transfer time of the climatic signal is independent of the AMOC background state. Our findings confirm a central role for ocean circulation in the bipolar seesaw and provide clear criteria for assessing hypotheses and model simulations of Dansgaard–Oeschger dynamics.

Nature↗

Groundwater quality and age of secondary bedrock aquifers in the glaciated portion of eastern Nebraska, 2016–18

The Eastern Nebraska Water Resources Assessment (ENWRA) project was initiated in 2006 to assist water managers by developing a hydrogeologic framework and water budget for the glaciated portion of eastern Nebraska. Within the ENWRA area, the primary groundwater sources for municipal, domestic, and irrigation water needs are provided by withdrawals from alluvial, buried paleovalley, and the High Plains aquifer (where present). Generally, other bedrock aquifers are considered a secondary water source. However, in some areas, such as parts of Sarpy and Nemaha Counties, these secondary bedrock aquifers are the only source of water within glaciated upland areas. To improve the understanding of the quality, geochemistry, and age of groundwater from bedrock aquifers, the U.S. Geological Survey (USGS), in cooperation with the ENWRA group, which includes the Lewis and Clark, Lower Elkhorn, Lower Platte North, Lower Platte South, Nemaha, and Papio-Missouri River Natural Resources Districts, designed a study to sample 31 wells completed in the secondary bedrock aquifers and analyze samples for major ions, physical properties, nutrients, stable isotopes, and selected age tracers. Of the 31 samples collected for this report, 22 samples were collected from the Dakota aquifer contained in the Dakota Sandstone, 3 from the Niobrara aquifer contained in the Niobrara Formation of Colorado Group, and 6 from Paleozoic aquifers contained in undifferentiated Paleozoic-age units. The results of this study indicate that major ion data collected from the Dakota aquifer can be used for assessing the quality, recharge source, and age of groundwater. Calcium bicarbonate dominant samples were characterized as modern or mixed, indicating that, in these areas, groundwater is unconfined and is recharged by precipitation and (or) surface water. If groundwater extraction rates exceed recharge rates, total dissolved solid concentrations may increase as a result of upwelling of groundwater from deeper units or formations, which can adversely affect groundwater quality. Sampling results presented in this report indicate water quality is good, but that groundwater in the Dakota aquifer with calcium bicarbonate water type may be vulnerable to surface contamination. In contrast, groundwater sampled from the Dakota aquifer, having a dominant water type other than calcium bicarbonate, generally has low dissolved oxygen and nitrate concentrations, and higher concentrations of total dissolved solids and trace elements, including iron and strontium. The geochemical characteristics of noncalcium bicarbonate samples from the Dakota aquifer indicated confining conditions and limited groundwater recharge from local precipitation. Apparent groundwater ages estimated from radiocarbon (carbon-14) sampling of noncalcium bicarbonate samples from the Dakota aquifer indicated that the time of groundwater recharge to the Dakota aquifer occurred during Pleistocene time. Depleted stable isotopes results indicate recharge during a colder climate. Groundwater under confined conditions is not easily recharged from precipitation or surface water. Future groundwater-level monitoring in locations where the Dakota aquifer appears to be confined could provide information to evaluate whether groundwater supplies remain sufficient to meet future municipal, domestic, and irrigation needs. For the Niobrara aquifer and Paleozoic aquifers, the dominant water type was not a diagnostic indicator of recharge source, age, and groundwater quality as with the Dakota aquifer. Most likely this is because the host formation was dominated by calcium-carbonate-rich rocks; however, few samples were collected from these aquifers to be able to confirm this interpretation. Samples collected from wells completed in the Niobrara aquifer and Paleozoic aquifers and characterized as calcium sulfate water type have statistically significantly higher concentrations of total dissolved solids compared to other samples from the Niobrara aquifer and Paleozoic aquifers characterized as calcium bicarbonate. Given that six of the nine of samples collected from the Niobrara and Paleozoic aquifers indicated modern recharge, these secondary bedrock aquifers are reliant on precipitation to sustain groundwater levels and may be vulnerable to a multiyear drought. Well yields of the Niobrara and Paleozoic aquifers are dependent on the presence of secondary porosity and these units offer little storage. Samples collected from wells completed in Paleozoic aquifers were the most isotopically enriched and similar to modern precipitation and had the highest concentrations of nitrate, indicating that groundwater is affected by agricultural activities. Future groundwater sampling would be beneficial to characterize groundwater-quality changes within the Niobrara and Paleozoic aquifers over time.

Nebraska↗

The National Ocean Biodiversity Strategy

President Biden has been clear that the ocean is central to life on Earth. As he has proclaimed, “the ocean powers millions of jobs; feeds and sustains us; and is a rejuvenating source of inspiration, exploration, and recreation.” The Biden-Harris Administration has worked hard to fulfill the President’s goal to protect and conserve at least 30% of U.S. waters by 2030. The ocean faces increased threats from warming, overfishing, increased acidity, and loss of biodiversity. It is now more important than ever to sustain the many benefits that the ocean, coasts, and Great Lakes provide, including food, a favorable climate, recreation, physical and mental health, and for many, a sense of cultural identity. Ocean life represents an irreplaceable heritage, the foundation of a habitable planet, and a vast trove of resources. Keeping our ocean healthy requires reliable information on the changing status of these living organisms, the drivers of biodiversity change, and options for effectively addressing those drivers. Over 2 million species are estimated to live in the ocean, yet only about 240,000 species have been described by scientists. Most of the ocean’s benefits result from those diverse species interacting with one another and the environment they create. To protect and conserve the ocean, we as a nation need to make better use of existing knowledge and prioritize acquiring new biodiversity knowledge to enable better policy and management decisions. The ability to monitor ocean species and habitats has expanded dramatically over the past decade, with innovations in technology, genomics, taxonomy, big data management and sharing, artificial intelligence, and machine learning. Yet large fractions of the U.S. ocean remain almost unknown. The National Ocean Biodiversity Strategy (strategy) reflects the urgent need to leverage these advances. The goals of this strategy must be guided by the nation’s diverse voices and ways of knowing, in order to maximize effective and equitable stewardship of the ocean’s diverse life and its benefits to people. The strategy is intended as a guiding document for government to advance three overarching goals: ● Goal 1: Drive delivery of ocean biodiversity knowledge at the national scale. Objectives include developing an Implementation Plan for achieving the strategy’s three goals; establishing a coordination mechanism to manage the implementation; and documenting gaps in biodiversity knowledge and the benefits of ocean biodiversity to people and economies. ● Goal 2: Strengthen tools and institutions to deliver ocean biodiversity knowledge. Objectives include establishing a robust information pipeline to support indicators and dynamic maps of ocean biodiversity, from the coasts to the deep sea. This pipeline should include expanded observing systems and comprehensive data management; science and technology solutions to accelerate the availability of biodiversity information; and plans to leverage previous investments to rebuild and expand the nation’s human capital and infrastructure to sustain foundational taxonomy and biodiversity science. ● Goal 3: Protect, conserve, restore, and sustainably use ocean biodiversity. Objectives include expanding the collection, delivery, and use of biodiversity knowledge to inform actions that advance ocean protection, conservation, restoration, and sustainable development. Government should lead in establishing and incentivizing diverse partnerships across scales and sectors to implement those actions and should educate and involve the public to discover and value the nation’s diverse ocean life. Achieving these goals will require commitments across society: new federal and private investments, coordination across sectors to address climate and equity challenges, and engagement of Indigenous Knowledge holders and frontline communities as full partners throughout planning and implementation. The Subcommittee on Ocean Science and Technology (SOST) IWG-Biodiversity will begin developing an Implementation Plan to describe and direct specific actions to implement the strategy. Successful implementation of the strategy will harmonize and expand collection and delivery of timely knowledge on ocean life to all of society. The strategy will also enable evidence-based management and protection of the ocean. Advancing the strategy will build human and institutional capital and partnerships that support both existing mandates and new needs to rebuild and sustain biodiversity, achieve healthy ocean ecosystems, and manage living resources. Implementing the strategy will deliver knowledge for monitoring, modeling, forecasting, and assessments that support food security, public health, and cultural values, and that more effectively protect, conserve, and restore nature.

Report↗

Status of whitebarkpine in the Greater Yellowstone Ecosystem: A step-trend analysis comparing 2004-2007 to 2008-2011

Whitebark pine (Pinus albicaulis) is a foundation and keystone species in upper subalpine environments of the northern Rocky Mountains that strongly influences the biodiversity and productivity of high-elevation ecosystems (Tomback et al. 2001, Ellison et al. 2005). Throughout its historic range, whitebark pine has decreased significantly as a major component of high-elevation forests. As a result, it is critical to understand the challenges to whitebark pine—not only at the tree and stand level, but also as these factors influence the distribution of whitebark pine across the Greater Yellowstone Ecosystem (GYE). In 2003, the National Park Service (NPS) Greater Yellowstone Inventory & Monitoring Network identified whitebark pine as one of twelve significant natural resource indicators or vital signs to monitor (Jean et al. 2005, Fancy et al. 2009) and initiated a long-term, collaborative monitoring program. Partners in this effort include the U.S. Geological Survey, U.S. Forest Service, and Montana State University with representatives from each comprising the Greater Yellowstone Whitebark Pine Monitoring Working Group. The objectives of the monitoring program are to assess trends in (1) the proportion of live, whitebark pine trees (>1.4-m tall) infected with white pine blister rust (blister rust); (2) to document blister rust infection severity by the occurrence and location of persisting and new infections; (3) to determine mortality of whitebark pine trees and describe potential factors contributing to the death of trees; and (4) to assess the multiple components of the recruitment of understory whitebark pine into the reproductive population. In this report we summarize the past eight years (2004-2011) of whitebark pine status and trend monitoring in the GYE. Our study area encompasses six national forests (NF), two national parks (NP), as well as state and private lands in portions of Wyoming, Montana, and Idaho; this area is collectively described as the GYE here and in other studies. The sampling design is a probabilistic, twostage cluster design with stands of whitebark pine as the primary units and 10x50 m belt transects as the secondary units. Primary sampling units (stands) were selected randomly from a sample frame of approximately 10,770 mapped pure and mixed whitebark pine stands ≥2.0 hectares in the GYE (Dixon 1997, Landenburger 2012). From 2004 through 2007 (monitoring transect establishment or initial time-step), we established 176 permanent belt transects (secondary sampling units=176) in 150 whitebark pine stands and permanently marked approximately 4,740 individual trees >1.4 m tall to monitor long-term changes in blister rust infection and survival rates. Between 2008 and 2011 (revisit time-step), these same 176 transects were surveyed and again all previously tagged trees were observed for changes in blister rust infection and survival status. Objective 1. Using a combined ratio estimator, we estimated the proportion of live trees infected in the GYE in the initial time-step (2004-2007) to be 0.22 (0.031 SE). Following the completion of all surveys in the revisit time-step (2008-2011), we estimated the proportion of live trees infected with white pine blister rust as 0.23 (0.028 SE; Table 2). We detected no significant change in the proportion of trees infected in the GYE between the two time-steps. Objective 2. We documented blister rust canker locations as occurring in the canopy or bole. We compared changes in canker position between the initial time-step (2004-2007) and the revisit time-step (2008-2011) in order to assess changes in infection severity. This analysis included the 3,795 trees tagged during the initial time-step that were located and documented as alive at the end of the revisit time-step. At the end of the revisit time-step, we found 1,217 trees infected with blister rust. This includes the 287 newly tagged trees in the revisit time step of which 14 had documented infections. Of these 1,217 trees, 780 trees were infected with blister rust in both time steps. Trees with only canopy cankers made up approximately 43% (519 trees) of the total number of trees infected with blister rust at the end of the revisit time-step, while trees with only bole cankers comprised 20% (252 trees), and those with both canopy and bole cankers included 37% (446 trees) of the infected sample. A bole infection is considered to be more consequential than a canopy canker, as it compromises not only the overall longevity of the tree, but its functional capacity for reproductive output as well (Kendall and Arno 1990, Campbell and Antos 2000, McDonald and Hoff 2001, Schwandt and Kegley 2004). In addition to infection location, we also documented infection transition between the canopy and bole. Of the 780 live trees that were infected with blister rust in both time-steps, approximately 31% (242) maintained canopy cankers and 36% (281) retained bole infections at the end of the revisit time-step. Infection transition from canopy to bole occurred in 30% (234) of the revisit time-step trees while 3% (23) transitioned from bole to canopy infections during this period. Objective 3. To determine whitebark pine mortality, we resurveyed all belt transects to reassess the life status of permanently tagged trees >1.4 m tall. We compared the total number of live tagged trees recorded during monitoring transect establishment to the total number of resurveyed dead tagged trees recorded during the revisit time-step and identified all potential mortality-influencing conditions (blister rust, mountain pine beetle, fire and other). By the end of the revisit time-step, we observed a total of 975 dead tagged whitebark pine trees; using a ratio estimator, this represents a loss of approximately 20% (SE=4.35%) of the original live tagged tree population (GYWPMWG 2012). Objective 4. To investigate the proportion of live, reproducing tagged trees, we divided the total number of positively identified cone-bearing trees by the total number of live trees in the tagged tree sample at the end of the revisit time-step. To approximate the average density of recruitment trees per stand, trees ≤1.4 m tall were summed by stand (within the 500 m² transect area) and divided by the total number of stands. Reproducing trees made up approximately 24% (996 trees) of the total live tagged population at the end of the revisit time-step. Differentiating between whitebark pine and limber pine seedlings or saplings is problematic given the absence of cones or cone scars. Therefore, understory summaries as presented in this report may include individuals of both species when they are sympatric in a stand. The average density of small trees ≤1.4 m tall was 53 understory trees per 500 m². Raw counts of these understory individuals ranged from 0-635 small trees per belt transect. In addition, a total of 287 trees were added to the tagged tree population by the end of 2011. These newly tagged trees were individuals that upon subsequent revisits had reached a height of >1.4 m tall and subsequently added to the sample. Throughout the past decade in the GYE, monitoring has helped document shifts in whitebark pine forests; whitebark pine stands have been impacted by insect, pathogen, wildland fire, and other disturbance events. Blister rust infection is ubiquitous throughout the ecosystem and infection proportions are variable across the region. And while we have documented mortality of whitebark pine, we have also recorded considerable recruitment. We provide this first step-trend report as a quantifiable baseline for understanding the state of whitebark pine in the GYE. Many aspects of whitebark pine health are highly variable across the range of its distribution in the GYE. Through sustained implementation of the monitoring program, we will continue efforts to document and quantify whitebark pine forest dynamics as they arise under periodic upsurges in insect, pathogen, fire episodes, and climatic events in the GYE. Since its inception, this monitoring program perseveres as one of the only sustained longterm efforts conducted in the GYE with a singular purpose to track the health and status of this prominent keystone species.

Natural Resource Technical Report↗

Groundwater flow model for the Little Plover River basin in Wisconsin’s Central Sands

The Little Plover River is a groundwater-fed stream in the sand plains region of central Wisconsin. In this region, sandy sediment deposited during or soon after the last glaciation forms an important unconfined sand and gravel aquifer. This aquifer supplies water for numerous high-capacity irrigation, municipal, and industrial wells that support a thriving agricultural industry. In recent years the addition of many new wells, combined with observed diminished flows in the Little Plover and other nearby rivers, has raised concerns about the impacts of the wells on groundwater levels and on water levels and flows in nearby lakes, streams, and wetlands. Diverse stakeholder groups, including well operators, Growers, environmentalists, local land owners, and regulatory and government officials have sought a better understanding of the local groundwater-surface water system and have a shared desire to balance the water needs of the he liagricultural, industrial, and urban users with the maintenance and protection of groundwater-dependent natural resources. To help address these issues, the Wisconsin Department of Natural Resources requested that the Wisconsin Geological and Natural History Survey and U.S. Geological Survey cooperatively develop a groundwater flow model that could be used to demonstrate the relationships among groundwater, surface water, and well withdrawals and also be a tool for testing and evaluating alternative water management strategies for the central sands region. Because of an abundance of previous studies, data availability, local interest, and existing regulatory constraints the model focuses on the Little Plover River watershed, but the modeling methodology developed during this study can apply to much of the larger central sands of Wisconsin. The Little Plover River groundwater flow model simulates three-dimensional groundwater movement in and around the Little Plover River basin under steady-state and transient conditions. This model explicitly includes all high-capacity wells in the model domain and simulates seasonal variations in recharge and well pumping. The model represents the Little Plover River, and other significant streams and drainage ditches in the model domain, as fully connected to the groundwater system, computes stream base flow resulting from groundwater discharge, and routes the flow along the stream channel. A separate soil-water-balance (SWB) model was used to develop groundwater recharge arrays as input for the groundwater flow model. The SWB model uses topography, soils, land use, and climatic data to estimate recharge as deep drainage from the soil zone. The SWB model explicitly includes recharge originating as irrigation water, and computes irrigation using techniques similar to those used by local irrigation operators. The groundwater flow model uses the U.S. Geological Survey’s MODFLOW modeling code which is freely available, widely accepted, and commonly used by the groundwater community. The groundwater flow model and the SWB model use identical high-resolution numerical grids having model cells 100 feet on a side, with physical properties assigned to each grid cell. This grid allows accurate geographic placement of wells, streams, and other model features. The 3-dimensional grid has three layers; layers 1 and 2 represent the sand and gravel aquifer and layer 3 represents the underlying sandstone. The distribution of material properties in the model (hydraulic conductivity, aquifer thickness, etc.) comes from previous published geologic studies of the region, updated by calibration to recent streamflow and groundwater level data. The SWB model operates on a daily time step. The groundwater flow model was calibrated to monthly stress periods with time steps ranging from 1 to 16 days. More detailed time discretization is possible. The groundwater model was calibrated to water-level and streamflow data collected during 2013 and 2014 by adjusting model parameters (primarily hydraulic conductivity, storage, and recharge) until the model produced a conditionally optimal fit between field observations and model output, subject to consistency with previously published geologic studies. Calibration was performed under both steady and transient conditions, and used a sophisticated parameter-estimation procedure (PEST) for the calibration process and to identify important model parameters. For the Little Plover River, the two most important parameters are the global recharge multiplier and the hydraulic conductivity of the stream bed. The calibrated model produces water-level and mass-balance results that are consistent with field observations and previous studies of the area. The completed model is a powerful tool for testing and demonstrating alternative water-management scenarios. Example model applications described in this report include simulating how the cumulative impacts of pumping and land-use change have affected average baseflow in the Little Plover River. Depletion-potential mapping represents a method for predicting which wells and well locations have the greatest impact on nearby surface-water resources. The completed model is publicly available, along with a companion user’s guide to assist with its operation, at http://wgnhs.org/littleplover- river-groundwater-model.

Wisconsin↗

Hotspots of biogeochemical activity linked to aridity and plant traits across global drylands

Perennial plants create productive and biodiverse hotspots, known as fertile islands, beneath their canopies. These hotspots largely determine the structure and functioning of drylands worldwide. Despite their ubiquity, the factors controlling fertile islands under conditions of contrasting grazing by livestock, the most prevalent land use in drylands, remain virtually unknown. Here we evaluated the relative importance of grazing pressure and herbivore type, climate and plant functional traits on 24 soil physical and chemical attributes that represent proxies of key ecosystem services related to decomposition, soil fertility, and soil and water conservation. To do this, we conducted a standardized global survey of 288 plots at 88 sites in 25 countries worldwide. We show that aridity and plant traits are the major factors associated with the magnitude of plant effects on fertile islands in grazed drylands worldwide. Grazing pressure had little influence on the capacity of plants to support fertile islands. Taller and wider shrubs and grasses supported stronger island effects. Stable and functional soils tended to be linked to species-rich sites with taller plants. Together, our findings dispel the notion that grazing pressure or herbivore type are linked to the formation or intensification of fertile islands in drylands. Rather, our study suggests that changes in aridity, and processes that alter island identity and therefore plant traits, will have marked effects on how perennial plants support and maintain the functioning of drylands in a more arid and grazed world.

Nature Plants↗

Simulated effects of increased groundwater withdrawals in the Cave Springs area, Hixson, Tennessee

Concern for future water supplies in Tennessee has grown in recent years as a result of increased awareness of competing needs, the impact of droughts, and the need for more water to support growing populations. The U.S. Geological Survey conducts investigations to improve the knowledge about interactions of geology, climate, humans, and ecosystems with the water cycle, which is critical to understanding and optimizing water availability. The Hixson Utility District in Hamilton County, Tennessee, uses groundwater resources in the Cave Springs area as a water supply, withdrawing water from two well fields located at Cave Springs and Walkers Corner. Historically, Hixson Utility District has withdrawn about 5 million gallons per day (Mgal/d) at the Cave Springs well field and between 2 and 3 Mgal/d at the Walkers Corner well field. To assess the capacity of the groundwater resources in the Cave Springs area to meet future demands, four different scenarios of increased groundwater withdrawals were analyzed using computer model simulations. In the study area, groundwater is present in both regolith and bedrock. Groundwater flow in the regolith occurs as diffuse flow as recharge from precipitation moves through the regolith to discharge to streams and springs or to the underlying bedrock. Most of the bedrock in the study area has low primary porosity and permeability; however, fracturing and dissolution have produced substantial secondary porosity and permeability. Groundwater flow through the bedrock occurs as both diffuse and conduit flow. Recharge to the aquifer is from two distinct sources: direct infiltration of precipitation and losing streams. A major source of recharge to the aquifer that supplies Cave Springs is surface water that is lost from North Chickamauga Creek as it flows from the Cumberland Plateau onto the Newman Limestone. Average annual streamflow loss (groundwater recharge) from this reach of North Chickamauga Creek for the period November 2000 through June 2006 is about 18 cubic feet per second (ft 3 /s). Groundwater leaves the aquifer as either discharge to North Chickamauga Creek, Poe Branch, and Lick Branch; discharge to Chickamauga Lake; spring flow to Cave Springs or Rogers Spring; or withdrawals at the Cave Springs or Walkers Corner well fields. Using computer model simulations, four scenarios of increased groundwater withdrawals were analyzed. Each of these four scenarios are compared to a base-case simulation that uses groundwater withdrawal rates from 2012 of 5.1 Mgal/d from the Cave Springs well field and 2.7 Mgal/d from the Walkers Corner well field. Under scenarios A and B, pumpage is increased at Cave Springs by 2 Mgal/d and 5 Mgal/d, respectively, while pumpage at Walkers Corner remains unchanged. Under scenarios C and D, pumpage is increased at Walkers Corner by 2.6 Mgal/d and 4.5 Mgal/d, respectively, while pumpage at Cave Springs remains unchanged. The effects of the increased withdrawals were analyzed by comparing water budget changes of the groundwater discharges to Chickamauga Lake, North Chickamauga Creek, Cave Springs, Poe Branch, and Lick Branch/Rogers Spring for each of the scenarios and evaluating changes in groundwater levels at the well fields. Under scenarios A and B, the largest change in the water budget occurs for flow to Cave Springs with decreases of 1.9 and 4.7 ft 3 /s, respectively. Similarly, groundwater discharge to North Chickamauga Creek decreases by 1.0 ft 3 /s and 2.6 ft 3 3/s, respectively. Under scenarios C and D, the largest change in the water budget occurs for flow to Chickamauga Lake with decreases of 1.3 ft 3 /s and 2.3 ft 3 /s, respectively. Similarly, groundwater discharge to North Chickamauga Creek decreases by 1.1 ft 3 /s and 2.1 ft 3 /s, respectively. Changes in groundwater levels at the well fields were also analyzed. At the Cave Springs well field, maximum declines in groundwater levels due to additional pumpage are less than 1 foot for all scenarios. Groundwater level changes at the Cave Springs well field are small due to the highly transmissive nature of the aquifer in this location. Maximum groundwater-level declines at Walkers Corner are less than 1 foot for scenarios A and B and about 52 feet and 82 feet for scenarios C and D, respectively. Under scenarios C and D, the regional potentiometric surface shows a large cone of depression centered on the Walkers Corner well field and elongated along geologic strike.

Tennessee↗

The 2016 groundwater flow model for Dane County, Wisconsin

A new groundwater flow model for Dane County, Wisconsin, replaces an earlier model developed in the 1990s by the Wisconsin Geological and Natural History Survey (WGNHS) and the U.S. Geological Survey (USGS). This modeling study was conducted cooperatively by the WGNHS and the USGS with funding from the Capital Area Regional Planning Commission (CARPC). Although the overall conceptual model of the groundwater system remains largely unchanged, the incorporation of newly acquired high-quality datasets, recent research findings, and improved modeling and calibration techniques have led to the development of a more detailed and sophisticated model representation of the groundwater system. The new model is three-dimensional and transient, and conceptualizes the county’s hydrogeology as a 12-layer system including all major unlithified and bedrock hydrostratigraphic units and two high-conductivity horizontal fracture zones. Beginning from the surface down, the model represents the unlithified deposits as two distinct model layers (1 and 2). A single layer (3) simulates the Ordovician sandstone and dolomite of the Sinnipee, Ancell, and Prairie du Chien Groups. Sandstone of the Jordan Formation (layer 4) and silty dolostone of the St. Lawrence Formation (layer 5) each comprise separate model layers. The underlying glauconitic sandstone of the Tunnel City Group makes up three distinct layers: an upper aquifer (layer 6), a fracture feature (layer 7), and a lower aquifer (layer 8). The fracture layer represents a network of horizontal bedding-plane fractures that serve as a preferential pathway for groundwater flow. The model simulates the sandstone of the Wonewoc Formation as an upper aquifer (layer 9) with a bedding-plane fracture feature (layer 10) at its base. The Eau Claire aquitard (layer 11) includes shale beds within the upper portion of the Eau Claire Formation. This layer, along with overlying bedrock units, is mostly absent in the preglacially eroded valleys along the Yahara River valley and in northeastern Dane County. Layer 12 represents the Mount Simon sandstone as the lowermost model layer. It directly overlies the Precambrian crystalline basement rock, whose top surface forms the lower boundary of the model. The model uses the USGS MODFLOW-NWT finite-difference code, a standalone version of MODFLOW-2005 that incorporates the Newton (NWT) solver. MODFLOW-NWT improves the handling of unconfined conditions by smoothing the transition from wet to dry cells. The model explicitly simulates groundwater–surface-water interaction with streamflow routing and lake-level fluctuation. Model input included published and unpublished hydrogeologic data from recent estimates of aquifer hydraulic conductivities. A spatial groundwater recharge distribution was obtained from a recent GIS-based, soil-water-balance model for Dane County. Groundwater withdrawals from pumping were simulated for 572 wells across the entire model domain, which includes Dane County and portions of seven neighboring counties—Columbia, Dodge, Green, Iowa, Jefferson, Lafayette, and Rock. These wells withdrew an average of 60 million gallons per day (mgd) over the 5-year period from 2006 through 2010. Within Dane County, 385 wells were simulated with an average withdrawal rate of 52 mgd. Model calibration used the parameter estimation code PEST, and calibration targets included heads, stream and spring flows, lake levels, and borehole flows. Steady-state calibration focused on the period 2006 through 2010; the transient calibration focused on the 7-week drought period from late May through July 2012. This model represents a significant step forward from previous work because of its finer grid resolution, improved hydrostratigraphic discretization, transient capabilities, and more sophisticated representation of surface-water features and multi-aquifer wells. Potential applications of the model include evaluation of potential sites for and impacts of new high-capacity wells, development of wellhead protection plans, evaluating the effects of changing land use and climate on groundwater, and quantifying the relationships between groundwater and surface water.

Wisconsin↗

Temporal Differences in Flow Depth and Velocity Distributions and Hydraulic Microhabitats Near Bridges of the Lower Platte River, Nebraska, 1934-2006

As part of a collaborative study of the cumulative impacts on stream and riparian ecology of water and channel management practices in the lower Platte River, Nebraska, this report describes a study by the U.S. Geological Survey in cooperation with the Lower Platte South Natural Resources District that summarizes: (1) temporal differences in distribution of streamflow depth, velocity, and microhabitats among five discrete 11-water-year periods 1934-44, 1951-61, 1966-76, 1985-95, and 1996-2006, and (2) the effects of bridge proximity on distribution of streamflow depth, velocity, and microhabitat of the Platte River when cross sections were measured at a similar discharge. The scope of the study included the four presently (2008) active streamflow-gaging stations located near bridges over the lower Platte River at North Bend, near Leshara, near Ashland, and at Louisville, Nebraska, and the most downstream streamflow-gaging station within the central Platte River segment near Duncan, Nebraska. Generally, in cases where temporal differences in streamflow depth and velocity were evident, at least one of the water-year periods from 1934 through 1995 had deeper streamflow than the recent water-year period (1996-2006). Temporal differences in distributions of streamflow depth were not strongly associated with differences in either climatic conditions or the maximum peak flow that occurred prior to the latest discharge measurement during each period. The relative cross-sectional area of most hydraulic niches did not differ among the water-year periods. Part of this apparent uniformity likely was an artifact of the broad microhabitat classification used for this study. In cases where temporal differences in relative cross-sectional area of hydraulic niches were evidenced, the differences occurred during high- and low-flow conditions, not during median flow conditions. The temporal differences in relative cross-sectional area were found more frequently for hydraulic niches defined by moderate and fast velocities than for hydraulic niches defined by slow velocities. Generally, any significant increase or decrease in the relative cross-sectional areas of hydraulic niches during the water-year periods from 1934 through 1995 had disappeared during the most recent water-year period, 1996-2006. Deep-Swift niche was the predominant hydraulic niche for all near-bridge sites on the lower Platte River for high- and median-flow conditions. The Deep-Swift niche also was the predominant niche for the near-bridge sites near Ashland and at Louisville for low-flow conditions; for the near-bridge sites at North Bend and near Leshara, streamflow cross-sectional areas during low-flow conditions were shared among the Shallow-Moderate, Intermediate-Moderate, Intermediate-Swift, and Deep-Swift hydraulic niches. For the near-bridge site near Duncan, the site farthest downstream in the central Platte River system, the Deep-Swift hydraulic niche was predominant only during high-flow conditions; during median- and low-flow conditions the relative cross-sectional area was shared among the Shallow-Slow, Shallow-Moderate, Intermediate-Moderate, and Intermediate-Swift hydraulic niches. Significant temporal differences in the relative cross-sectional area of the Deep-Swift hydraulic niche were found for sites near the two farthest downstream bridges near Ashland and at Louisville, but only for low-flow conditions. The Deep-Swift microhabitat was of special interest because it is the preferred hydraulic habitat during the adult life of the endangered pallid sturgeon (Scaphirhynchus albus). Temporal differences in relative cross-sectional areas of the Glide low-flow geomorphic microhabitat that contained the Deep-Swift hydraulic niche also indicated that relative cross-sectional areas of the Glide during the 1951-61 and 1996-2006 water-year periods were lower than during the 1966-76 period. The temporal differences indicated that any significant temporal chang

Scientific Investigations Report↗

Handbook of acute toxicity of chemicals to fish and aquatic invertebrates : summaries of toxicity tests conducted at Columbia National Fisheries Research Laboratory, 1965-78

Acute toxicity is a major subject of research at Columbia National Fisheries Research Laboratory for evaluating the impact of toxic chemicals on fishery resources. The Laboratory has played a leading role in developing research technology for toxicity testing and data interpretation. In 1965-78, more than 400 chemicals were tested against a variety of invertebrates and fish species representative of both cold- and warm-water climates. The use of acute toxicity tests for assessing the potential hazard of chemical contaminants to aquatic organisms is well documented (Boyd 1957; Henderson et al. 1960; Sanders and Cope 1966; Macek and McAllister 1970). Static acute toxicity tests provide rapid and (within limits) reproducible concentration-response curves for estimating toxic effects of chemicals on aquatic organisms. These tests provide a database for determining relative toxicity of a large number of chemicals to a variety of species and for estimating acute effects of chemical spills on natural aquatic systems; they also assist in determining priority and design of additional toxicity studies. Acute toxicity tests usually provide estimates of the exposure concentration causing 50% mortality (LC50) to test organisms during a specified period of time. For certain invertebrates, the effective concentration is based on immobilization, or some other identifiable endpoint, rather than on lethality. The application of the LC50 has gained acceptance among toxicologists and is generally the most highly rated test for assessing potential adverse effects of chemical contaminants to aquatic life (Brungs and Mount 1978; American Institute for Biological Sciences 1978 a ). The literature contains numerous papers dealing with the acute toxicity of chemicals to freshwater organisms. However, there is a tremendous need for a concise compendium of toxicity data covering a large variety of chemicals and test species. This Handbook is a compilation of a large volume of acute toxicity data from the Columbia Laboratory and its field laboratories. It presents definitive acute toxicity data on 271 chemicals tested against a variety of freshwater invertebrates and fishes. The chemicals represent all major groups of pesticides, as well as numerous industrial chemicals. This compilation should serve as a useful database for the many agencies and organizations dealing with research and management programs concerned with the impact of chemicals on aquatic resources. The Columbia Laboratory has played a major role in developing currently used standard methodology for static acute toxicity testing. The use of standardized methodology greatly reduces variation in results. The data presented here have been carefully scrutinized to eliminate tests that failed to follow acceptable procedures. Handling of test organisms and procedures for static toxicity tests followed those described by Lennon and Walker (1964) and Macek and McAllister (1970), and conform well with those recommended by Brauhn and Schoettger (1975) and the Committee on Methods for Toxicity Tests with Aquatic Organisms (1975). The species of fish and invertebrates that were tested are listed in phylogenetic order in Tables 1 and 2. Fish were obtained from Federal and State hatcheries as either eggs or fry. Original stocks of invertebrates were collected and cultured from wild populations with no known source of contamination; these populations were replenished regularly. The invertebrates were cultured in the Laboratory by methods similar to those described by Sanders and Cope (1966). Test chemicals usually consisted of technical or analytical grade samples of known purity. Formulations of the chemicals were also tested when available. When purity of test chemicals was known, all calculated concentrations were based on percent active ingredients. Stock solutions were prepared immediately before each test, with commercial grade acetone as the carrier solvent. Occasionally, ethanol or dimethyl-formamide was substituted. Solvent concentrations did not exceed 0.5 mL/L in final dilution water. Test water (dilution water) was reconstituted from deionized water of at least 10 6 ohms resistivity by the addition of appropriate reagent grade chemicals (Marking 1969). Water was buffered to maintain a pH of 7.2 to 7.5, an alkalinity of 30 to 35 mg/L, and a hardness of 40 to 50 mg/L as CaCO 3 . Test water was mixed thoroughly and aerated before transfer into test chambers. Fish were acclimated to dilution water by gradually changing the water in acclimated tanks from 100% well water to 100% reconstituted water over a 1- to 3-day period at the desired testing temperature. Invertebrates were acclimated from well water to dilution water over a 4- to 6-h period. Toxicity tests were conducted under static conditions without aeration, and the organisms were not fed during acclimation or testing. Temperature of test solutions was maintained within ± 1°C of that required for a given test. Toxicity tests with fish were conducted in 18.9-liter (5-gal) wide-mouthed jars containing 15 liters of test solution. Fingerling fish weighing 0.2 to 1.5 g were tested at each concentration. Caution was taken not to exceed 0.8 g of test organisms per liter of solution. Duplicate test chambers were used to accommodate larger fish. Test chambers varied in size for invertebrates, depending on the species used; volume of test solution ranged from 0.25 to 4 liters. At least 10 organisms were exposed to each concentration for all definitive tests. At least six concentrations were used per toxicity test. The tests began upon initial exposure to the toxicant and continued for 96 h. Immobilization tests with invertebrates were conducted for only 48 h. The number of dead or affected organisms in each test chamber were recorded and the dead organisms were removed every 24 h; general observations on the condition of test organisms were also recorded at these times. Toxicity data were analyzed by a statistical method described by Litchfield and Wilcoxon (1949) to determine LC50 (theoretical estimate of the concentration lethal to 50% of the test animals) and 95% confidence intervals. This method is recommended by the American Public Health Association (1971) and by Sprague (1969) for determining median lethal concentrations. The procedure is easily modified for computing a single LC50 when replicate tests are performed.

Resource Publication↗

Vegetation changes over 12 years in ungrazed and grazed Conservation Reserve Program Grasslands in the central and southern plains

The Conservation Reserve Program (CRP) established under the 1985 Food Security Act has the fundamental objectives of jointly providing economic support to segments of the agricultural community and conservation of natural resources (Osborn, 1997; Heard and others, 2000). Although soil loss on highly erodable lands was the principal natural resource conservation issue addressed in the 1985 CRP, improving water quality and wildlife habitat both became important considerations as the program evolved (Farmer and others, 1988). For example, Best and others (1997) found that production of young birds on CRP fields in the Midwest was ≥15 times the production on row-crop fields because of improved habitat. The increasing importance of wildlife habitat is reflected in continuing refinement of the Environmental Benefits Index (EBI) used by the U.S. Department of Agriculture (USDA) to quantify the potential benefits of enrolling lands in CRP (Osborn, 1997; Ribaudo and others, 2001). The refinements reflect input furnished by federal, state, and non-government organizations seeking greater wildlife habitat quality on CRP lands (Roseberry and David, 1994; Hughes and others, 1995; Millenbah and others, 1996; Patterson and Best, 1996; Rodgers, 1999; Allen and others, 2001). Refinement in the EBI has changed the types of grasses planted on newly enrolled land. In early CRP signups (1 through 11), 71% of new grassland acres were planted to introduced grasses and legumes [Conservation Practice (CP) 1] while 29% of the acres were planted to native grasses (CP2) (Osborn and others, 1992). By the 27th signup in July 2004, over 34.8 million acres (14 million ha) were enrolled in the CRP. More than 73% of these lands were planted to various mixtures of introduced (CP1) or native (CP2) grasses for a minimal contract period of 10 years (USDA, 2004). Continuation of grass plantings under the 2002 Farm Bill may result in CRP lands furnishing grass dominated cover for 20 or more consecutive years. The species of grass established in seeded grasslands can have a major influence on the potential quality of wildlife habitat where vegetation is maintained over a multi-year period. Different species of grass may have comparable abilities to alleviate soil erosion but furnish dissimilar qualities of wildlife habitat (fig. 1). For example, smooth brome, an introduced cool-season grass (grass species and scientific names are presented in table 1), is highly valued for its erosion control and forage attributes (Casler and Carlson, 1995). Switchgrass, a native warm-season grass, also is valued for its soil and water conservation qualities (Moser and Vogel, 1995) but provides greater benefits for some species of wildlife (Clubine, 1995). The quality of nesting and winter cover for ring-necked pheasants (Phasianus colchicus) furnished by smooth brome on northeastern Colorado CRP lands is inferior to that provided by the taller, more robust switchgrass (Allen, 1994). Characteristics of the agricultural landscape surrounding individual CRP fields also play a role in the wildlife habitat potential of CRP plantings (Weber and others, 2002; Nusser and others, 2004). of native and seeded grasslands change in response to the presence (and absence) of physical disturbances such as fire, grazing, tillage, and haying (Hobbs and Huenneke, 1992; Millenbah and others, 1996; Allen and others, 2001; Renfrew and Ribic, 2001; Swengel and Swengel, 2001). The perpetuation of diversity in species composition and vegetation structure following disturbance sustains desirable habitat for a variety of grassland-dependent wildlife (Hall and Willig, 1994; Barnes and others, 1995; King and Savidge, 1995; Granfors and others, 1996; Herkert and others, 1996; Kurzejeski, 1996; Patterson and Best, 1996; Klute and others, 1997). Undisturbed grasslands have lower grass and forb species diversity, greater amounts of dead plant material, decreased as well as seasonally delayed productivity, and diminished structural diversity of vegetation (Peet and others, 1975; Rice and Parenti, 1978; Butler and Briske, 1988; Campa and Winterstein, 1992). Recommendations for the timing of disturbance to increase grass and forb species diversity range from 3 to 8 years following establishment of seeded grasslands in the northern Great Plains and Midwest (Duebbert and others, 1981; Higgens, 1987; Millenbah and others, 1996). The management interval, however, is affected by climatic conditions, soils, grass species, and management history of the individual stand. We quantified changes in vegetation structure and species composition across the typical 10-year contract period in undisturbed southern and central Great Plains CRP fields (fig. 2) planted to introduced and native grasses. In addition, we compared changes in vegetation in fields grazed during the emergency release of 1996 by comparing conditions prior to grazing and two and four years post grazing relative to changes in similar fields that were not grazed. Documentation of long-term changes in vegetation structure and composition for fields planted to common grass seed mixtures across a wide range of environmental conditions provides information to improve long-term wildlife habitat potential, guide program administration, and define management practices that yield economic benefits to operators while still meeting wildlife and conservation objectives. Emergency grazing provisions of the CRP are controversial. Although grazing can alter vegetative characteristics and reduce habitat quality in the short-term (Temple and others, 1999), periodic disturbance may be necessary to maintain habitat quality, and more information is needed assessing long-term effects of emergency grazing on vegetative structure and species composition.

Conference Paper↗

A three-dimensional numerical model of predevelopment conditions in the Death Valley regional ground-water flow system, Nevada and California

In the early 1990's, two numerical models of the Death Valley regional ground-water flow system were developed by the U.S. Department of Energy. In general, the two models were based on the same basic hydrogeologic data set. In 1998, the U.S. Department of Energy requested that the U.S. Geological Survey develop and maintain a ground-water flow model of the Death Valley region in support of U.S. Department of Energy programs at the Nevada Test Site. The purpose of developing this "second-generation" regional model was to enhance the knowledge an understanding of the ground-water flow system as new information and tools are developed. The U.S. Geological Survey also was encouraged by the U.S. Department of Energy to cooperate to the fullest extent with other Federal, State, and local entities in the region to take advantage of the benefits of their knowledge and expertise. The short-term objective of the Death Valley regional ground-water flow system project was to develop a steady-state representation of the predevelopment conditions of the ground-water flow system utilizing the two geologic interpretations used to develop the previous numerical models. The long-term objective of this project was to construct and calibrate a transient model that simulates the ground-water conditions of the study area over the historical record that utilizes a newly interpreted hydrogeologic conceptual model. This report describes the result of the predevelopment steady-state model construction and calibration. The Death Valley regional ground-water flow system is situated within the southern Great Basin, a subprovince of the Basin and Range physiographic province, bounded by latitudes 35 degrees north and 38 degrees 15 minutes north and by longitudes 115 and 118 degrees west. Hydrology in the region is a result of both the arid climatic conditions and the complex geology. Ground-water flow generally can be described as dominated by interbasinal flow and may be conceptualized as having two main components: a series of relatively shallow and localized flow paths that are superimposed on deeper regional flow paths. A significant component of the regional ground-water flow is through a thick Paleozoic carbonate rock sequence. Throughout the flow system, ground water flows through zones of high transmissivity that have resulted from regional faulting and fracturing. The conceptual model of the Death Valley regional ground-water flow system used for this study is adapted from the two previous ground-water modeling studies. The three-dimensional digital hydrogeologic framework model developed for the region also contains elements of both of the hydrogeologic framework models used in the previous investigations. As dictated by project scope, very little reinterpretation and refinement were made where these two framework models disagree; therefore, limitations in the hydrogeologic representation of the flow system exist. Despite limitations, the framework model provides the best representation to date of the hydrogeologic units and structures that control regional ground-water flow and serves as an important information source used to construct and calibrate the predevelopment, steady-state flow model. In addition to the hydrogeologic framework, a complex array of mechanisms accounts for flow into, through, and out of the regional ground-water flow system. Natural discharges from the regional ground-water flow system occur by evapotranspiration, springs, and subsurface outflow. In this study, evapotranspiration rates were adapted from a related investigation that developed maps of evapotranspiration areas and computed rates from micrometeorological data collected within the local area over a multiyear period. In some cases, historical spring flow records were used to derive ground-water discharge rates for isolated regional springs. For this investigation, a process-based, numerical model was developed to estimate net infiltration. This result can be used to represent recharge by assuming that water that infiltrates past the “root zone” ultimately becomes ground-water recharge. Net infiltration based on this water-balance approach, however, suggests that recharge in the region exceeds the estimated ground-water discharge for the basin by almost a factor of four. The net infiltration model contains some assumptions and was limited in regard to soil moisture, bedrock permeability, and precipitation distributions that result in these higher than expected values. Despite this limitation, the net infiltration model does provide a good starting point for estimating regional ground-water recharge in the ground-water flow model. The Death Valley regional ground-water flow system was simulated using a three-dimensional steady-state model that incorporated a nonlinear least-squares regression technique to estimate selected model parameters. The numerical modeling program MODFLOW-2000 was used in creating a finite-difference model consisting of 194 rows, 160 columns, and 15 layers. The grid cells were oriented north-south and were of uniform size, with side dimensions of 1,500 meters. The layers span thicknesses of 50 to 300 meters. The model grid encompasses about 70,000 square kilometers. The required initial model parameter values were supplied by discretization of the three-dimensional hydrogeologic framework model and digital representations of the remaining conceptual model components. The three-dimensional simulation and corresponding sensitivity analysis supported the hypothesis of interactions between a relatively shallow local and subregional flow system and a deeper dominant regional system controlled by the carbonate aquifer. During calibration of the model, techniques available in MODFLOW-2000 allowed for estimation of a series of parameters that provided a best fit to observed hydraulic heads and flows. Numerous conceptual models were evaluated to test the validity of various interpretations about the flow system. Only those conceptual model changes contributing to a significant improvement in model fit, as indicated by a reduction in the sum of squared errors, were retained in the final optimized model. The final model was evaluated to assess the likely accuracy of simulated results by comparing measured and expected quantities with simulated values. Evaluation of the model indicates that although the model is clearly an improvement on previous representations of the flow system, there is an indication of important uncertainties and model error.

California, Nevada↗

Geohydrology of the French Creek Basin and simulated effects of drought and ground-water withdrawals, Chester County, Pennsylvania

This report describes the results of a study by the U.S. Geological Survey, in cooperation with the Delaware River Basin Commission, to develop a regional ground-water-flow model of the French Creek Basin in Chester County, Pa. The model was used to assist water-resource managers by illustrating the interconnection between ground-water and surface-water systems. The 70.7-square mile French Creek Basin is in the Piedmont Physiographic Province and is underlain by crystalline and sedimentary fractured-rock aquifers. Annual water budgets were calculated for 1969-2001 for the French Creek Basin upstream of streamflow-measurement station French Creek near Phoenixville (01472157). Average annual precipitation was 46.28 in. (inches), average annual streamflow was 20.29 in., average annual base flow determined by hydrograph separation was 12.42 in., and estimated average annual ET (evapotranspiration) was 26.10 in. Estimated average annual recharge was 14.32 in. and is equal to 31 percent of the average annual precipitation. Base flow made up an average of 61 percent of streamflow. Ground-water flow in the French Creek Basin was simulated using the finite-difference MODFLOW-96 computer program. The model structure is based on a simplified two-dimensional conceptualization of the ground-water-flow system. The modeled area was extended outside the French Creek Basin to natural hydrologic boundaries; the modeled area includes 40 square miles of adjacent areas outside the basin. The hydraulic conductivity for each geologic unit was calculated from reported specific-capacity data determined from aquifer tests and was adjusted during model calibration. The model was calibrated for above-average conditions by simulating base-flow and water-level measurements made on May 1, 2001, using a recharge rate of 20 in/yr (inches per year). The model was calibrated for below-average conditions by simulating base-flow and water-level measurements made on September 11 and 17, 2001, using a recharge rate of 6.2 in/yr. Average conditions were simulated by adjusting the recharge rate until simulated streamflow at streamflow-measurement station 01472157 matched the long-term (1968-2001) average base flow of 54.1 cubic feet per second. The recharge rate used for average conditions was 15.7 in/yr. The effect of drought in the French Creek Basin was simulated using a drought year recharge rate of 8 in/yr for 3 months. After 3 months of drought, the simulated streamflow of French Creek at streamflow-measurement station 01472157 decreased 34 percent. The simulations show that after 6 months of average recharge (15.7 in/yr) following drought, streamflow and water levels recovered almost to pre-drought conditions. The effect of increased ground-water withdrawals on stream base flow in the South Branch French Creek Subbasin was simulated under average and drought conditions with pumping rates equal to 50, 75, and 100 percent of the Delaware River Basin Commission Ground Water Protected Area (GWPA) withdrawal limit (1,393 million gallons per year) with all pumped water removed from the basin. For average recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 18, 28, and 37 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. After 3 months of drought recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 27, 40, and 52 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. The effect of well location on base flow, water levels, and the sources of water to the well was simulated by locating a hypothetical well pumping 200 gallons per minute in different places in the Beaver Run Subbasin with all pumped water removed from the basin. The smallest reduction in the base flow of Beaver Run was from a well on the drainage divide between the French Creek Basin and the Marsh Creek Basin to the south; the simulated base flow of Beaver Run at the mouth was reduced 1 percent. The greatest reduction in the base flow of Beaver Creek was from a well close to Beaver Run; the simulated base flow of Beaver Run at the mouth was reduced 8 percent. The simulations showed that (1) if the contributing area of a well is in a basin, pumping will affect stream base flow and water levels in that basin whether the well is inside or outside that basin; (2) wells in different areas of a basin away from a divide produce a similar reduction in base flow; (3) a well within a basin will derive more water from diverted base flow and less water from storage than a well on or near a basin divide; and (4) the reduction in base flow at the mouth of the stream is the same for a well in the headwaters and a well downstream near the confluence. Model simulations illustrate some of the typical analyses and results that can be produced. The model was calibrated using annual values for recharge and ground-water ET and then was run using the annual values in a seasonally independent transient mode to show changes with time. The timing and relative magnitude of some of the changes simulated with the model when viewed in terms of a normal climatic year may be subject to considerable uncertainty because of the variability in seasonal recharge and ground-water ET rates. Transient model simulations for short-term periods are indicative of possible hydrologic system response and are considered an approximation.

Pennsyvlania↗

Growth form and lifespan of herbaceous species mediate the role of traits in short-term drought response

Increased climate variability is expected to intensify short-term drought events. Plants have evolved stress tolerance strategies involving trade-offs in resource conservation, mycorrhizal collaboration and plant size, yet how these strategies promote drought resistance across different herbaceous plant groups remains unknown. Leveraging 63 globally distributed grassland and shrubland sites from the International Drought Experiment, we identified plant traits linked to drought resistance in 661 populations of 421 species after 1 year of extreme drought. We assessed how traits, site precipitation and drought severity affected cover change across growth forms and lifespans, and how trait–environment interactions influenced drought resistance. Across all species, leaf N (an acquisitive trait) was associated with drought resistance, whereas in forbs, drought resistance was also associated with a conservative root trait and plant size. In addition, interactions among traits mediated drought resistance; root traits predicted performance only in concert with other traits. Environmental variables influenced trait effects on drought resistance, notably for annuals in wetter sites, suggesting that drought-escape strategies in annuals may be advantageous only under mild stress. Our study highlights variability in traits that predict drought resistance across herbaceous plant groups, emphasizing the importance of species context, environmental stress and the selection of traits in research and management.

Nature Ecology and Evolution↗

Preliminary assessment of channel stability and bed-material transport in the Tillamook Bay tributaries and Nehalem River basin, northwestern Oregon

This report summarizes a preliminary study of bed-material transport, vertical and lateral channel changes, and existing datasets for the Tillamook (drainage area 156 square kilometers [km 2 ]), Trask (451 km 2 ), Wilson (500 km 2 ), Kilchis (169 km 2 ), Miami (94 km 2 ), and Nehalem (2,207 km 2 ) Rivers along the northwestern Oregon coast. This study, conducted in coopera-tion with the U.S. Army Corps of Engineers and Oregon Department of State Lands to inform permitting decisions regarding instream gravel mining, revealed that: Study areas along the six rivers can be divided into reaches based on tidal influence and topography. The fluvial (nontidal or dominated by riverine processes) reaches vary in length (2.4-9.3 kilometer [km]), gradient (0.0011-0.0075 meter of elevation change per meter of channel length [m/m]), and bed-material composition (a mixture of alluvium and intermittent bedrock outcrops to predominately alluvium). In fluvial reaches, unit bar area (square meter of bar area per meter of channel length [m 2 /m]) as mapped from 2009 photographs ranged from 7.1 m 2 /m on the Tillamook River to 27.9 m 2 /m on the Miami River. In tidal reaches, all six rivers flow over alluvial deposits, but have varying gradients (0.0001-0.0013 m/m) and lengths affected by tide (1.3-24.6 km). The Miami River has the steepest and shortest tidal reach and the Nehalem River has the flattest and longest tidal reach. Bars in the tidal reaches are generally composed of sand and mud. Unit bar area was greatest in the Tidal Nehalem Reach, where extensive mud flats flank the lower channel. Background factors such as valley and channel confinement, basin geology, channel slope, and tidal extent control the spatial variation in the accumulation and texture of bed material. Presently, the Upper Fluvial Wilson and Miami Reaches and Fluvial Nehalem Reach have the greatest abundance of gravel bars, likely owing to local bed-material sources in combination with decreasing channel gradient and valley confinement. Natural and human-caused disturbances such as mass movements, logging, fire, channel modifications for navigation and flood control, and gravel mining also have varying effects on channel condition, bed-material transport, and distribution and area of bars throughout the study areas and over time. Existing datasets include at least 16 and 18 sets of aerial and orthophotographs that were taken of the study areas in the Tillamook Bay tributary basins and Nehalem River basin, respectively, from 1939 to 2011. These photographs are available for future assessments of long-term changes in channel condition, bar area, and vegetation establishment patterns. High resolution Light Detection And Ranging (LiDAR) surveys acquired in 2007-2009 could support future quantitative analyses of channel morphology and bed-material transport in all study areas. A review of deposited and mined gravel volumes reported for instream gravel mining sites shows that bed-material deposition tends to rebuild mined bar surfaces in most years. Mean annual deposition volumes on individual bars exceeded 3,000 cubic meters (m 3 ) on Donaldson Bar on the Wilson River, Dill Bar on the Kilchis River, and Plant and Winslow Bars on the Nehalem River. Cumulative reported volumes of bed-material deposition were greatest at Donaldson and Dill Bars, totaling over 25,000 m 3 per site from 2004 to 2011. Within this period, reported cumulative mined volumes were greatest for the Donaldson, Plant, and Winslow Bars, ranging from 24,470 to 33,940 m 3 . Analysis of historical stage-streamflow data collected by the U.S. Geological Survey on the Wilson River near Tillamook (14301500) and Nehalem River near Foss (14301000) shows that these rivers have episodically aggraded and incised, mostly following high flow events, but they do not exhibit systematic, long-term trends in bed elevation. Multiple cross sections show that channels near bridge crossings in all six study areas are dynamic with many subject to incision and aggradation as well as lateral shifts in thalweg position and bank deposition and erosion. In fluvial reaches, unit bar area declined a net 5.3-83.6 percent from 1939 to 2009. The documented reduction in bar area may be attributable to several factors, including vegetation establishment and stabilization of formerly active bar surfaces, lateral channel changes and resulting alterations in sediment deposition and erosion patterns, and streamflow and/or tide differences between photographs. Other factors that may be associated with the observed reduction in bar area but not assessed in this reconnaissance level study include changes in the sediment and hydrology regimes of these rivers over the analysis period. In tidal reaches, unit bar area increased on the Tillamook and Nehalem Rivers (98.0 and 14.7 percent, respectively), but declined a net 24.2 to 83.1 percent in the other four tidal reaches. Net increases in bar area in the Tidal Tillamook and Nehalem Reaches were possibly attributable to tidal differences between the photographs as well as sediment deposition behind log booms and pile structures on the Tillamook River between 1939 and 1967. The armoring ratio (ratio of the median grain sizes of a bar's surface and subsurface layers) was 1.6 at Lower Waldron Bar on the Miami River, tentatively indicating a relative balance between transport capacity and sediment supply at this location. Armoring ratios, however, ranged from 2.4 to 5.5 at sites on the Trask, Wilson, Kilchis, and Nehalem Rivers; these coarse armor layers probably reflect limited bed-material supply at these sites. On the basis of mapping results, measured armoring ratios, and channel cross section surveys, preliminary conclusions are that the fluvial reaches on the Tillamook, Trask, Kilchis, and Nehalem Rivers are currently sediment supply-limited in terms of bed material - that is, the transport capacity of the channel generally exceeds the supply of bed material. The relation between transport capacity and sediment is more ambiguous for the fluvial reaches on the Wilson and Miami Rivers, but transport-limited conditions are likely for at least parts of these reaches. Some of these reaches have possibly evolved from sediment supply-limited to transport-limited over the last several decades in response to changing basin and climate conditions. Because of exceedingly low gradients, all the tidal reaches are transport-limited. Bed material in these reaches, however, is primarily sand and finer grain-size material and probably transported as suspended load from upstream reaches. These reaches will be most susceptible to watershed conditions affecting the supply and transport of fine sediment. Compared to basins on the southwestern Oregon coast, such as the Chetco and Rogue River basins, these six basins likely transport overall less gravel bed material. Although tentative in the absence of actual transport measurements, this conclusion is supported by the much lower area and frequency of bars and longer tidal reaches along all the northcoast rivers examined in this study. Previous studies suggest that the expansive and largely unvegetated bars visible in the 1939 photographs are primarily associated with voluminous sedimentation starting soon after the first Tillamook Burn fire in 1933. However, USGS studies of temporal bar trends in other Oregon coastal rivers unaffected by the Tillamook Burn show similar declines in bar area over approximately the same analysis period. In the Umpqua and Chetco River basins, historical declines in bar area are associated with long-term decreases in flood magnitude. Other factors may include changes in the type and volume of large wood and riparian vegetation. Further characterization of hydrology patterns in these basins and possible linkages with climate factors related to flood peaks, such as the Pacific Decadal Oscillation, could support inferences of expected future changes in vegetation establishment and channel planform and profile. More detailed investigations of bed-material transport rates and channel morphology would support assessments of lateral and vertical channel condition and longitudinal trends in bed material. Such assessments would be most practical for the fluvial study areas on the Wilson, Kilchis, Miami, and Nehalem Rivers and relevant to several ongoing management and ecological issues pertaining to sand and gravel transport. Tidal reaches may also be logical subjects for indepth analysis where studies would be more relevant to the deposition and transport of fine sediment (and associated channel and riparian conditions and processes) rather than coarse bed material.

Oregon↗

Water-quality trends of urban streams in Independence, Missouri, 2005–18

The U.S. Geological Survey and the city of Independence, Missouri, Water Pollution Control Department has studied the water quality and ecological condition of urban streams within Independence since 2005. Selected physical properties, nutrients, chloride, fecal indicator bacteria (Escherichia coli and total coliform), total dissolved solids, and suspended-sediment concentration data for base-flow and stormflow samples were used to document temporal trends in concentrations and flow-weighted concentrations; and annual loads were computed and investigated for selected nutrients, chloride, and suspended sediment. The six study sites included in this report are located on five urban streams: Rock Creek, a tributary in the city that drains to the Missouri River; three tributaries of the Little Blue River within the city (East Fork Little Blue River, Adair Creek, and Spring Branch Creek); and two sites on the main stem of the Little Blue River (one upstream from the city and one downstream from the three tributaries). Many factors such as population, land use, and climate, and combinations of these factors contributed to the significant changes in the concentrations and transport of nutrients, chloride, fecal indicator bacteria, and suspended sediment in the urban streams within Independence. The population of Independence and the amount of developed land in the urban watersheds remained unchanged during the 2005–18 study. Differences were noted in precipitation and in streamflow during the study. Annual precipitation and streamflow were separated into two time periods within the study—period 1 (2006–10), having greater annual streamflow and precipitation, and period 2 (2011–18), having about 30 percent lower annual streamflow and less precipitation. Streamflow was an important factor in the transport of nitrogen, phosphorus, chloride, and suspended sediment from the urban watersheds. Changes in data collection methodology during the study period and improvements to the city stormwater and wastewater infrastructure also could have contributed to some of the trends. Between 2009 and 2015, more than 35 million dollars of improvements were made to stormwater and wastewater infrastructure within the city. These improvements, such as additional sewage overflow holding tanks, removal of septic tanks, and improved and expanded sanitary sewer lines and storm overflows, also could have affected the decreased nutrients and fecal indicator bacteria trends among the urban streams in the study area. Models were used for analyzing streamflow-related variability in constituent concentrations and loads to determine if the water quality changed significantly during the study period. Trends in concentration data at four sites were analyzed using a statistical package called R–QWTREND and trends in load data were analyzed at six sites using a statistical package called Weighted Regressions on Time, Discharge, and Season-Kalman filter (WRTDS–K); both developed by the U.S. Geological Survey and publicly available for use. Statistically significant trends in flow-weighted nutrient concentrations and loads generally were downward during the study period. The only nutrient compound with a statistically significant upward trend in flow-weighted concentration was dissolved orthophosphate as phosphorus at the Rock Creek site and the upstream site on the Little Blue River. A statistically significant downward trend in annual dissolved ammonia load was identified at the downstream Little Blue River site. A significant upward linear trend in annual orthophosphate as phosphorus load was identified on Adair Creek. A statistically significant upward trend in dissolved chloride concentrations was identified at the downstream Little Blue River site. Road salt application near the site during the winter could have resulted in higher concentrated runoff during wet weather conditions. Annual chloride loads significantly decreased in Adair Creek and Spring Branch Creek. The mean annual chloride load transported in the drier (2011–18) period 2 was significantly less than during the wetter (2006–10) period 1, indicating that trends in precipitation runoff are an important factor in trends in annual transport of chloride. Statistically significant downward trends in flow-weighted fecal indicator bacteria Escherichia coli (E. coli) population densities were noted for Rock Creek and the down-stream site on the Little Blue River. However, no trend was identified in E. coli population density at the upstream Little Blue River site. The downward trend in E. coli population density at the downstream site could be a result of decreased streamflow and precipitation over the study period, storage of fecal indicator bacteria in the Little Blue River streambed within the study area, die-off of fecal indicator bacteria during travel from upstream to downstream, changes in the sample collection methodology, improvements to the city’s storm-water and wastewater infrastructures, or a combination of these factors. The statistically significant downward trend in suspended-sediment concentration identified at the upstream Little Blue River site could be affected by the decreased streamflow and precipitation during the study period, by changes in sampling methods within the study period, and by the decrease in construction and urban land development upstream from the city. No statistically significant change was indicated in the annual suspended-sediment load transported from Independence to the Little Blue River during the study period. More than one-half the suspended sediment transported in the Little Blue River originated in the watershed upstream from Independence. The Little Blue River and many of its tributaries that drain Independence have been designated as recreational waters classified for whole-body contact class B and secondary contact recreation, and some have been listed as impaired for E. coli by the Missouri Department of Natural Resources from urban runoff and storm sewers. Observations were made among the available E. coli population density data for both Little Blue River sites to further understand water-quality conditions over the study period. Both Little Blue River sites had similar medians and geometric means for the recreational season (April through October) and during the full study period, both of which are greater than the regulatory population density for both recreational classes. The Little Blue River drainage area nearly doubles in size from the upstream to downstream site; therefore, the consistent geometric mean and median of E. coli population densities at the upstream and downstream Little Blue River sites could be primarily due to the larger volume of streamflow creating a dilution effect. Other possible factors could be storage of fecal indicator bacteria in stream bed sediments, die-off of fecal indicator bacteria during transport, improvements to the city’s wastewater and stormwater infrastructure, changes to sampling methodology, or a combination of these factors. Specific sources of the E. coli are currently (2019) unknown.

Missouri↗

Water-quality assessment of the Connecticut, Housatonic, and Thames River Basins study unit: Analysis of available data on nutrients, suspended sediments, and pesticides, 1972-92

This retrospective report examines available nutrient, suspended sediment, and pesticide data in surface and ground water in the Connecticut, Housatonic and Thames Rivers Study Unit of the National Water-Quality Assessment Program. The purpose of this study is to improve the understanding of natural and anthropogenic factors affecting water quality in the study unit. Waterquality data were acquired from various sources, primarily, the U.S. Geological Survey and the U.S. Environmental Protection Agency. The report examines data for water years 1972-92, focusing on 1980-92, although it also includes additional data from as early as 1905. The study unit lies within the New England Physiographic Province and altitudes range from sea level in coastal Connecticut to 6,288 feet above sea level at Mount Washington, New Hampshire. Two major aquifer types underlie the study unit unconsolidated glacial deposits and fractured bedrock. The climate generally is temperate and humid, with four distinct seasons. Average annual precipitation ranges from 34 to 65 inches. The study unit has a population of about 4.5 million, which is most highly concentrated in southwestern Connecticut and along the south-central region of the Connecticut River Valley. Surface-water-quality data were screened to provide information about sites with adequate numbers of analyses (50) over sufficiently long periods (1980-90) to enable valid statistical analyses. In order to compare effects of different types of land use on surface-water quality, examination of data required application of several statistical and graphical techniques, including mapping, histograms, boxplots, concentration-discharge plots, trend analysis, and load estimation. Spatial and temporal analysis of surface-water-quality data indicated that, with a single exception, only stations in the Connecticut water-quality network had sufficient data collected over adequately long time periods to use in detailed analyses. Ground-water nutrient and pesticide data were compiled from several Federal and State agencies, primarily the U.S. Geological Survey, U.S. Environmental Protection Agency, and Connecticut Department of Health Services. Nutrient data were available for several thousand wells; nitrite plus nitrate as nitrogen was the most commonly reported constituent. Most wells with nutrient data are in Massachusetts and Connecticut. Relative to nutrient data in ground and surface water, pesticide data are less common. Pesticide data were available for slightly more than 200 surface-water sites and less than 500 wells; about 95 percent of the wells are completed in stratified-drift or till aquifers. Data for 81 pesticide compounds were available in various data bases. 2,4-D and silvex were the most commonly detected herbicides in surface water and dieldrin and diazinon were the most commonly detected insecticides. Most surface-water pesticide samples and detections are from bed sediment, but much of the data are not recent. Ethylene dibromide (EDB), a soil fumigant used in tobacco farming was detected in 268 wells in a 50 square-mile area of north-central Connecticut; EDB contamination also was detected in wells in Massachusetts. Atrazine, an herbicide commonly used in corn farming, commonly was detected in wells installed in tilled agricultural fields. Corn herbicides were commonly detected in the northern part of the study unit, although the sampling has been less frequent than in the southern part of the study unit. Pesticides were seldom detected in public-supply wells in Connecticut. Urban sites with relatively high population densities and high concentrations of dischargers were characterized by having the highest nutrient concentrations and loads when adjusted for differences in drainage area or population. Particularly, the Pequabuck, Naugatuck, and Quinnipiac River Basins were characterized by high nutrient concentrations median total nitrogen concentrations ranged from 3.3 to 4.2 mg/L (milligrams per liter) and median total phosphorus concentrations ranged from 0.42 to 0.8 mg/L. In contrast, the predominantly forested and low density residential land-use sites, such as Saugatuck and Salmon River Basins, were characterized by low nutrient concentrations median total nitrogen ranged from 0.50 to 0.60 mg/L and median total phosphorus concentrations ranged from 0.01 to 0.02 mg/L. Estimated total nitrogen loadings in median discharge years ranged from 940 kilograms per square mile at the Salmon River near East Hampton, Conn., to 5,800 kilograms per square mile at the Naugatuck River at Beacon Falls, Conn. Water quality, in terms of nutrient concentrations and areally adjusted loadings, for sites with large drainage basins integrating a wide variety of land-use categories fell between the extremes of the urban and forested sites total nitrogen was 1,400 kilograms per square mile per year at the Connecticut River at Thompsonville, Conn. Nitrate concentrations in ground water occasionally exceeded the safe drinking-water standard of 10 mg/L as nitrogen. The greatest number of detections exceeding the standard, however, were not in public-water supplies but in shallow observation wells in agricultural settings (the most frequently sampled type of well). None of the public-supply wells in Massachusetts exceeded the standard. Although nitrate concentrations for Vermont and New Hampshire generally were low, few data were available and those were seldom reported on the basis of drainage basin, making analysis difficult. Trend analysis indicated that flow-adjusted concentrations of total and dissolved phosphorus generally decreased during the period of analysis, however, total nitrogen did not change substantially. Decreases in ammonia concentrations with time were usually accompanied by increases in nitrate, suggesting improvements in sewage treatment. The lack of adequate data from more or less exclusively agricultural areas points to the need for further study of the effects of fanning on surface-water quality in the study unit. Furthermore, additional information is needed on the rates, transformations, and movements of nutrients and other materials through and between the aquatic and terrestrial components of the study unit.

Connecticut, Massachusetts, New Hampshire, Quebec,↗