USGS Science⌕ Search

SEARCH · USGS Science

Results for “Diversity and Distributions”

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

Agriculture — A river runs through it — The connections between agriculture and water quality

Sustaining the quality of the Nation’s water resources and the health of our diverse ecosystems depends on the availability of sound water-resources data and information to develop effective, science-based policies. Effective management of water resources also brings more certainty and efficiency to important economic sectors. Taken together, these actions lead to immediate and longterm economic, social, and environmental benefits that make a difference to the lives of the almost 400 million people projected to live in the United States by 2050. In 1991, Congress established the U.S. Geological Survey National Water-Quality Assessment (NAWQA) to address where, when, why, and how the Nation’s water quality has changed, or is likely to change in the future, in response to human activities and natural factors. Since then, NAWQA has been a leading source of scientific data and knowledge used by national, regional, state, and local agencies to develop science-based policies and management strategies to improve and protect water resources used for drinking water, recreation, irrigation, energy development, and ecosystem needs. Plans for the third decade of NAWQA (2013–23) address priority water-quality issues and science needs identified by NAWQA stakeholders, such as the Advisory Committee on Water Information and the National Research Council, and are designed to meet increasing challenges related to population growth, increasing needs for clean water, and changing land-use and weather patterns. This report is one of a series of publications, The Quality of Our Nation’s Waters , which describes major findings of the NAWQA Project on water-quality issues of regional and national concern and provides science-based information for assessing and managing the quality of our groundwater resources. Other reports in this series focus on occurrence and distribution of nutrients, pesticides, and volatile organic compounds in streams and groundwater, the effects of contaminants and stream-flow alteration on the condition of aquatic communities in streams, and on the quality of groundwater from private domestic and public supply wells. Each reports builds toward a more comprehensive understanding of the quality of regional and national water resources. All NAWQA reports are available online ( https://water.usgs.gov/nawqa/bib/ ). We hope this publication will provide you with insights and information to meet your water-resource needs and will foster increased citizen awareness and involvement in the protection and restoration of our Nation’s waters. The information in this report is intended primarily for those interested or involved in resource management and protection, conservation, regulation, and policymaking at the regional and national levels.

Circular↗

The Geodetic Centroid (gCent) Catalog: Global earthquake monitoring with satellite imaging geodesy

Remote sensing geodetic observations (Interferometric Synthetic Aperture Radar [InSAR] and optical correlation [“pixel tracking”]) serve an increasingly diverse and important role in earthquake monitoring and response. This study introduces the Geodetic Centroid (gCent) catalog—an earthquake catalog derived solely from space‐based geodetic observations—and analysis of 74 earthquakes ( ⁠ M W 4.3–7.4) imaged from 1 August 2019 to 01 February 2022. For gCent, we use InSAR and optical correlation observations derived from the Sentinel‐1 satellites and various publicly available optical satellites to systematically image all global earthquakes M W 5.5 or larger and shallower than 25 km, M W 7.0 or larger at any depth, and other high‐impact earthquakes or seismic events of special interest. We invert surface displacements from successfully imaged earthquakes for the location, orientation, and dimensions of a single slipping fault patch that describes the centroid characteristics of the earthquake. These centroid models, in turn, are compiled into a catalog and used in U.S. Geological Survey/Advanced National Seismic System (ANSS) operational earthquake response products such as ShakeMaps and finite‐fault models. We provide a comparison of the gCent catalog to the ANSS Comprehensive Catalog and Global Centroid Moment Tensor (Global CMT) catalog to compare reported locations, depths, and magnitudes. We find that global earthquake catalogs not only generally provide reasonably comparable locations (within 10 km on average), but also they systematically overestimate depth that may have implications for earthquake shaking predictions based solely on earthquake origin information. Geodetic magnitudes are comparable to seismically inferred magnitudes, indicating that gCent models are unlikely to be systematically biased by the presence of postseismic deformation. We additionally highlight limitations of the gCent catalog induced by both the limitations of remote sensing imaging of earthquakes and our imposition of a simplified earthquake source description that does not include spatially distributed slip.

Bulletin of the Seismological Society of America↗

USGS response to an urban earthquake, Northridge '94

The urban centers of our Nation provide our people with seemingly unlimited employment, social, and cultural opportunities as a result of the complex interactions of a diverse population embedded in an highly-engineered environment. Catastrophic events in one or more of the natural earth systems which underlie or envelop urban environment can have radical effects on the integrity and survivability of that environment. Earthquakes have for centuries been the source of cataclysmic events on cities throughout the world. Unlike many other earth processes, the effects of major earthquakes transcend all political, social, and geomorphic boundaries and can have decided impact on cities tens to hundreds of kilometers from the epicenter. In modern cities, where buildings, transportation corridors, and lifelines are complexly interrelated, the life, economic, and social vulnerabilities in the face of a major earthquake can be particularly acute. The 1994 Northridge Earthquake was a major test for parts of what many consider the most earthquake-prepared and best-engineered metropolitan region in the United States. While the combined efforts of concerned professionals at all levels of government, academia, and the private sector have produced significant advances in our knowledge of the causes and potential effects of earthquakes, and of ways to reduce their impact, it remains unfortunately true that actual earthquakes provide important opportunities to test those advances in our knowledge. In the hours and days following the Northridge Earthquake, the four Federal member agencies of the National Earthquake Hazards Reduction Program (NEHRP), the Federal Emergency Management Agency (FEMA), the National Institute for Standards and Technology (NIST), the National Science Foundation (NSF), and the United States Geological Survey (USGS) laid out a detailed plan for collecting, analyzing, archiving, and reporting information that would benefit the Nation in future earthquake hazards reduction efforts. Congress provided a special appropriation to FEMA to carry out this plan and FEMA distributed those funds to the four NEHRP agencies. The USGS was responsible for conducting geophysical and geological investigations in support of the NEHRP Post-earthquake Study. For the past 2 years, the USGS has rigorously pursued over 40 tasks focused on the USGS Northridge Earthquake Mission. This document is a summary report of the USGS findings; additional technical reports on specific USGS tasks are appearing in various scientific journals and USGS publications.

California↗

Methods for estimating annual exceedance-probability discharges for streams in Iowa, based on data through water year 2010

A statewide study was performed to develop regional regression equations for estimating selected annual exceedance-probability statistics for ungaged stream sites in Iowa. The study area comprises streamgages located within Iowa and 50 miles beyond the State’s borders. Annual exceedance-probability estimates were computed for 518 streamgages by using the expected moments algorithm to fit a Pearson Type III distribution to the logarithms of annual peak discharges for each streamgage using annual peak-discharge data through 2010. The estimation of the selected statistics included a Bayesian weighted least-squares/generalized least-squares regression analysis to update regional skew coefficients for the 518 streamgages. Low-outlier and historic information were incorporated into the annual exceedance-probability analyses, and a generalized Grubbs-Beck test was used to detect multiple potentially influential low flows. Also, geographic information system software was used to measure 59 selected basin characteristics for each streamgage. Regional regression analysis, using generalized least-squares regression, was used to develop a set of equations for each flood region in Iowa for estimating discharges for ungaged stream sites with 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities, which are equivalent to annual flood-frequency recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, respectively. A total of 394 streamgages were included in the development of regional regression equations for three flood regions (regions 1, 2, and 3) that were defined for Iowa based on landform regions and soil regions. Average standard errors of prediction range from 31.8 to 45.2 percent for flood region 1, 19.4 to 46.8 percent for flood region 2, and 26.5 to 43.1 percent for flood region 3. The pseudo coefficients of determination for the generalized least-squares equations range from 90.8 to 96.2 percent for flood region 1, 91.5 to 97.9 percent for flood region 2, and 92.4 to 96.0 percent for flood region 3. The regression equations are applicable only to stream sites in Iowa with flows not significantly affected by regulation, diversion, channelization, backwater, or urbanization and with basin characteristics within the range of those used to develop the equations. These regression equations will be implemented within the U.S. Geological Survey StreamStats Web-based geographic information system tool. StreamStats allows users to click on any ungaged site on a river and compute estimates of the eight selected statistics; in addition, 90-percent prediction intervals and the measured basin characteristics for the ungaged sites also are provided by the Web-based tool. StreamStats also allows users to click on any streamgage in Iowa and estimates computed for these eight selected statistics are provided for the streamgage.

Iowa↗

Geologic and physiographic controls on bed-material yield, transport, and channel morphology for alluvial and bedrock rivers, western Oregon

The rivers of western Oregon have diverse forms and characteristics, with channel substrates ranging from continuous alluvial gravel to bare bedrock. Analysis of several measurable morphologic attributes of 24 valley reaches on 17 rivers provides a basis for comparing nonalluvial and alluvial channels. Key differences are that alluvial reaches have greater bar area, greater migration rates, and show systematic correlation among variables relating grain size to bed-material transport capacity. We relate these differences between channel types to bed-material transport rates as derived from a coupled regional analysis of empirical sediment yield measurements and physical experiments of clast attrition during transport. This sediment supply analysis shows that overall bed-material transport rates for western Oregon are chiefly controlled by (1) lithology and basin slope, which are the key factors for bed-material supply into the stream network, and (2) lithologic control of bed-material attrition from in-transport abrasion and disintegration. This bed-material comminution strongly affects bed-material transport in the study area, reducing transport rates by 50%–90% along the length of the larger rivers in the study area. A comparison of the bed-material transport estimates with the morphologic analyses shows that alluvial gravel-bed channels have systematic and bounding relations between bed-material transport rate and attributes such as bar area and local transport capacity. By contrast, few such relations are evident for nonalluvial rivers with bedrock or mixed-bed substrates, which are apparently more influenced by local controls on channel geometry and sediment supply. At the scale of western Oregon, the physiographic and lithologic controls on the balance between bed-material supply and transport capacity exert far-reaching influence on the distribution of alluvial and nonalluvial channels and their consequently distinctive morphologies and behaviors—differences germane for understanding river response to tectonics and environmental perturbations, as well as for implementing effective restoration and monitoring strategies.

California, Nevada, Oregon, Washington↗

Understanding the influence of nutrients on stream ecosystems in agricultural landscapes

Sustaining the quality of the Nation’s water resources and the health of our diverse ecosystems depends on the availability of sound water-resources data and information to develop effective, science-based policies. Effective management of water resources also brings more certainty and efficiency to important economic sectors. Taken together, these actions lead to immediate and long-term economic, social, and environmental benefits that make a difference to the lives of the almost 400 million people projected to live in the United States by 2050. In 1991, Congress established the U.S. Geological Survey (USGS) National Water-Quality Assessment (NAWQA) to address where, when, why, and how the Nation’s water quality has changed, or is likely to change in the future, in response to human activities and natural factors. Since then, NAWQA has been a leading source of scientific data and knowledge used by national, regional, State, and local agencies to develop science-based policies and management strategies to improve and protect water resources used for drinking water, recreation, irrigation, energy development, and ecosystem needs ( https://water.usgs.gov/nawqa/applications/ ). Plans for the third decade of NAWQA (2013–23) address priority water-quality issues and science needs identified by NAWQA stakeholders, such as the Advisory Committee on Water Information and the National Research Council, and are designed to meet increasing challenges related to population growth, increasing needs for clean water, and changing land-use and weather patterns. Excess nutrients are a pervasive problem of streams, lakes, and coastal waters. The current report, “The Quality of Our Nation’s Waters—Understanding the Effects of Nutrients on Stream Ecosystems in Agricultural Landscapes,” presents a summary of results from USGS investigations conducted from 2003 to 2011 on processes that influence nutrients and how nutrient enrichment can alter biological components of agricultural streams. This study included collecting data from 232 sites distributed among eight study areas. This report summarizes findings on processes that influence nutrients and how nutrient enrichment can alter biological communities in agricultural streams. These findings are relevant to local, State, regional, and national decision-makers involved in efforts to (1) better understand the influence of nutrients on agricultural streams, (2) develop nutrient criteria for streams and rivers, (3) reduce nutrients to streams and downstream receiving waters, and (4) develop tools for tracking nutrient and biological conditions following nutrient reduction strategies. All NAWQA reports are available online at https://water.usgs.gov/nawqa/bib/ . We hope this publication will provide you with insights and information to meet your water-resource needs and will foster increased citizen awareness and involvement in the protection and restoration of our Nation’s waters. The information in this report is intended primarily for those interested or involved in resource management and protection, conservation, regulation, and policymaking at the regional and national levels.

Central Columbia-Yakima River Basins, Great River ↗

Ground deformation and gravity for volcano monitoring

Introduction When magma accumulates or migrates, it can cause pressurization and related ground deformation. Characterization of surface deformation provides important constraints on the potential for future volcanic activity, especially in combination with seismic activity, gas emissions, and other indicators. A wide variety of techniques and instrument types have been applied to the study of ground deformation at volcanoes (sidebar, p. 2; Dzurisin, 2000, 2003, 2007). Geodetic instruments include continuously recording Global Navigation Satellite System (GNSS; of which the United States’ Global Positioning System is one example) stations (fig. D1), borehole tiltmeters, and interferometric synthetic aperture radar (InSAR) measurements (from satellites, occupied and unoccupied aircraft systems, and ground-based sensors). Additional geodetic measurements like continuous- and survey-mode gravity (fig. D2) can contribute substantially to interpreting these data. Borehole strainmeters (see chapter K , this volume, by Hurwitz and Lowenstern, 2024) also have outstanding utility for monitoring deformation, although because of cost and permitting challenges, we do not include them as part of standard volcano monitoring networks for U.S. volcanoes. Still other techniques like light detection and ranging (lidar), structure from motion, and optical satellite data can be used to derive gross topographic changes, which can be used to map volcanic deposits, infer eruption rates, and gain insights into the source processes associated with eruptive activity (see chapter G , this volume, on tracking surface changes caused by volcanic activity; Orr and others, 2024). Experience has shown that no single geodetic monitoring technique is adequate to detect and track the entire range of ground-motion patterns that occur at volcanoes, primarily because of the temporal and spatial diversity of volcano deformation (fig. D3). Similarly, the magnitude of surface deformation varies widely. Geodetic monitoring strategies should therefore include multiple techniques and instrument types to cover a wide range of spatial and temporal scales. In identifying recommendations for geodetic instrumentation for volcano monitoring networks, we attempted to maximize the diversity of instrument types to measure the full range of deformation signals and minimize their expense and number; thus, we do not include several well-known deformation-monitoring techniques in our recommendations. Extensometers, for example, measure strains over distances of a few meters and have an excellent record of success in detecting changes in preeruptive localized ground motion across existing cracks, including at Mount St. Helens, Washington (Iwatsubo and others, 1992), and Piton de la Fournaise, Réunion Island (Peltier and others, 2006). Despite being relatively inexpensive, extensometers are best used primarily when localized ground displacements (for example, ground cracks) need to be tracked, and are not necessary at all volcanoes. In considering volcano deformation monitoring strategies, two complicating factors are deserving of special attention. First, not all deformation is driven by subsurface magmatic activity—for example, at many large stratovolcanoes (for example, Mount Rainier), flank collapses and landslides are significant geologic hazards (Reid and others, 2001) that may occur even in the absence of magmatic activity. Monitoring the stability of volcanoes is thus another critical application of geodetic monitoring networks to inform hazard assessment. One of the most famous examples of edifice instability is the large flank collapse that initiated the May 18, 1980, eruption of Mount St. Helens. Deformation monitoring had detected a bulge on the north flank of the mountain in April 1980 that was expanding by several meters per day (Lipman and others, 1981). Given that flank collapses can happen at any time during a period of volcanic unrest (or even outside a period of unrest), the capability to assess edifice stability is critical. Second, although volcanoes are commonly treated as idealized structures that erupt from single points, like centralvent stratovolcanoes, many are characterized by long rift zones from which eruptions may originate, and distributed volcanic fields are characterized by broadly spaced vents. For example, linear dikes are common at Kīlauea, Mauna Loa, and between Mount Shasta and Medicine Lake in California. At Kīlauea, one of these linear dikes emerged more than 40 kilometers (km) away from the summit of the volcano during the lower East Rift Zone eruption in 2018. Other volcanic fields, like Lassen volcanic center, California, or the San Francisco Volcanic Field, Arizona, have many small vents spread over a wide area. Although the instrumentation guidelines presented in this chapter remain phrased for central-vent volcanoes, they should be modified as needed in the context of the eruptive characteristics of each individual volcanic system. Spatial analysis of geodetic network coverage could help to ensure adequate instrumentation in areas where volcanism can occur over a broad area as opposed to a central vent. As an example, consider the adjacent volcanoes Mount Shasta and Medicine Lake. If station locations are chosen based only on the distance from the centers of the volcanoes, then any geodetic anomalies between the two volcanoes—an area of potential volcanism as indicated by the presence of volcanic features—may remain undetected by ground-based instrumentation. The spatial analysis is accomplished via a grid of pressure point sources (Mogi, 1958) evenly distributed across the map area, at a depth of 5 km in this example (fig. D4). Each source is inflated until predicted deformations exceed the GNSS white noise uncertainty estimates at one site (Langbein, 2017; Murray and Svarc, 2017). This volume of detectable magma provides a measure of the quality of the coverage (fig. D4). The results indicate that, as of 2022, there is a large area between Mount Shasta and Medicine Lake volcano with existing mapped dikes in which a substantial amount of magma could intrude without being detected geodetically. Applying this style of analysis to individual volcanic systems can provide a guide for designing network geometry given the expected locations of future eruptions.

Scientific Investigations Report↗

Traces in the dark: sedimentary processes and facies gradients in the upper shale member of the Upper Devonian-Lower Mississippian Bakken Formation, Williston Basin, North Dakota, U.S.A.

Black, organic-rich rocks of the upper shale member of the Upper Devonian–Lower Mississippian Bakken Formation, a world-class petroleum source rock in the Williston Basin of the United States and Canada, contain a diverse suite of mudstone lithofacies that were deposited in distinct facies belts. The succession consists of three discrete facies associations (FAs). These comprise: 1) siliceous mudstones; 2) quartz- and carbonate-bearing, laminated mudstones; and 3) macrofossil-debris-bearing massive mudstones. These FAs were deposited in three facies belts that reflect proximal to distal relationships in this mudstone system. The macrofossil-debris-bearing massive mudstones (FA 3) occur in the proximal facies belt and contain erosion surfaces, some with overlying conodont and phosphate–lithoclast lag deposits, mudstones with abundant millimeter-scale siltstone laminae showing irregular lateral thickness changes, and shell debris. In the medial facies belt, quartz- and carbonate-bearing, laminated mudstones dominate, exhibiting sub-millimeter-thick siltstone layers with variable lateral thicknesses and localized mudstone ripples. In the distal siliceous mudstone facies belt, radiolarites, radiolarian-bearing mudstones, and quartz- and carbonate-bearing, laminated mudstones dominate. Overall, total organic carbon (TOC) contents range between about 3 and 10 wt %, with a general proximal to distal decrease in TOC content. Abundant evidence of bioturbation exists in all FAs, and the lithological and TOC variations are paralleled by changes in burrowing style and trace-fossil abundance. While two horizontal traces and two types of fecal strings are recognized in the proximal facies belt, only a single horizontal trace fossil and one type of fecal string characterize mudstones in the distal facies belt. Radiolarites intercalated into the most distal mudstones are devoid of traces and fecal strings. Bedload transport processes, likely caused by storm-induced turbidity currents, were active across all facies belts. Suspended sediment settling from near the ocean surface, however, most likely played a role in the deposition of some of the mudstones, and was probably responsible for deposition of the radiolarites. The distribution pattern of high-TOC sediments in proximal and lower-TOC deposits in some distal facies is interpreted as a function of higher accumulation rates during radiolarian depositional events leading to a decrease in suspension-derived organic carbon in radiolarite laminae. The presence of burrows in all FAs and nearly all facies in the upper Bakken shale member indicates that dysoxic conditions prevailed during its deposition. This study shows that in intracratonic high-TOC mudstone successions such as the upper Bakken shale member bed-load processes most likely dominated sedimentation, and conditions promoted a thriving infaunal benthic community. As such, deposition of the upper Bakken shale member through dynamic processes in an overall dysoxic environment represents an alternative to conventional anoxic depositional models for world-class source rocks.

North Dakota↗

Fraser River watershed, Colorado — Assessment of available water-quantity and water-quality data through water year 1997

The water-quantity and water-quality data for the Fraser River watershed through water year 1997 were compiled for ground-water and surface-water sites. In order to assess the water-quality data, the data were related to land use/land cover in the watershed. Data from 81 water-quantity and water-quality sites, which consisted of 9 ground-water sites and 72 surface-water sites, were available for analysis. However, the data were limited and frequently contained only one or two water-quality analyses per site. The Fraser River flows about 28 miles from its headwaters at the Continental Divide to the confluence with the Colorado River. Ground-water resources in the watershed are used for residential and municipal drinking-water supplies. Surface water is available for use, but water diversions in the upper parts of the watershed reduce the flow in the river. Land use/land cover in the watershed is predominantly forested land, but increasing urban development has the potential to affect the quantity and quality of the water resources. Analysis of the limited ground-water data in the watershed indicates that changes in the land use/land cover affect the shallow ground-water quality. Water-quality data from eight shallow monitoring wells in the alluvial aquifer show that iron and manganese concentrations exceeded the U.S. Environmental Protection Agency secondary maximum contaminant level. Radon concentrations from these monitoring wells exceeded the U.S. Environmental Protection Agency proposed maximum contaminant level. The proposed radon contaminant level is currently being revised. The presence of volatile organic compounds at two monitoring wells in the watershed indicates that land use affects the shallow ground water. In addition, bacteria detected in three samples are at concentrations that would be a concern for public health if the water was to be used as a drinking supply. Methylene blue active substances were detected in the ground water at some sites and are a possible indication of contamination from wastewater. Age of the alluvial ground water ranged from 10 to 30 years; therefore, results of land-management practices to improve water quality may not be apparent for many years. Surface-water-quality data for the Fraser River watershed are sparse. The surface-water-quality data show that elevated concentrations of selected constituents generally are related to specific land uses in the watershed. For one sample (about 2 percent; 1 of 53), dissolved manganese concentration exceeded the U.S. Environmental Protection Agency secondary maximum contaminant level. Two samples from two surface-water sites in the watershed exceeded the un-ionized ammonia chronic criterion. Spatial distribution of nutrient species (ammonia, nitrite, nitrate, and total phosphorus) shows that elevated concentrations occur primarily downstream from urban areas. Sites with five or more years of record were analyzed for temporal trends in concentration of nutrient species. Downward trends were identified for ammonia and nitrite for three surface-water sites. For nitrate, no trends were observed at two sites and a downward trend was observed at one site. Total phosphorus showed no trend for the site near the mouth of the Fraser River. Downward trends in the nutrient species may reflect changes in the wastewater-treatment facilities in the watershed. Bacteria sampling completed in the watershed indicates that more bacteria are present in the water near urban settings. The limited ground-water and surface-water data for the Fraser River watershed provide a general assessment of the quantity and quality of these resources. Concentrations of most water-quality constituents generally are less than ground- and surface-water-quality standards, but the presence of bacteria, some volatile organic compounds, methylene blue active substances, and increased nutrients in the water may indicate that land use is affecting the water quality.

Colorado↗

Wildfire reveals transient changes to individual traits and population responses of a native bumble bee ( Bombus vosnesenskii )

1. Fire-induced changes in the abundance and distribution of organisms, especially plants, can alter resource landscapes for mobile consumers driving bottom-up effects on their population sizes, morphologies, and reproductive potential. We expect these impacts to be most striking for obligate visitors of plants, like bees and other pollinators, but these impacts can be difficult to interpret due to the limited information provided by forager counts in the absence of survival or fitness proxies. 2. Increased bumble bee worker abundance is often coincident with the pulses of flowers that follow recent fire. However, it is unknown if observed postfire activity is due to underlying population growth or a stable pool of colonies recruiting more foragers to abundant resource patches. This distinction is necessary for determining the net impact of disturbance on bumble bees: are there population-wide responses or do just a few colonies reap the rewards? 3. We estimated colony abundance before and after fire in burned and unburned areas using a genetic mark-recapture framework. We paired colony abundance estimates with measures of body size, counts of queens, and estimates of foraging and dispersal to assess changes in worker size, reproductive output, and landscape-scale movements. 4. Higher floral abundance following fire not only increased forager abundance, but also the number of colonies from which those foragers came. Importantly, despite a larger population size we also observed increased mean worker size. Two years following fire, queen abundance was higher in both burned and unburned sites, potentially due to the dispersal of queens from burned into unburned areas. The effects of fire were transient; within two growing seasons, worker abundance was substantially reduced across the entire sampling area and body sizes were similar between burned and unburned sites. 5. Our results reveal how disturbance can temporarily release populations from resource limitation, boosting the genetic diversity, body size, and reproductive output of populations. Given that the effects of fire on bumble bees acted indirectly through pulsed resource availability, it is likely our results are generalizable to other situations, such as habitat restorations, where resource density is enhanced within the landscape.

California↗

Mapping croplands of Europe, Middle East, Russia, and Central Asia using Landsat 30-m data, machine learning algorithms and Google Earth Engine

Accurate and timely information on croplands is important for environmental, food security, and policy studies. Spatially explicit cropland datasets are also required to derive information on crop type, crop yield, cropping intensity, as well as irrigated areas. Large area defined as continental to global cropland mapping is challenging due to differential manifestation of croplands, wide range of cultivation practices and limited reference data availability. This study presents the results of a cropland extent mapping of 64 Countriescovering large parts of Europe, Middle East, Russia and Central Asia. To cover such a vast area, roughly 160,000 Landsat scenes from 3,351 footprints between 2014 and 2016 were processed within the Google Earth Engine (GEE) cloud-platform. We used the pixel-based supervised Random Forest (RF) machine learning algorithm with a set of satellite data inputs capturing diverse spectral, temporal and topographical characteristics across twelve agroecological zones (AEZs). The reference data to train the classification model were collected from very high spatial resolution imagery (VHRI) and ancillary datasets. The result is a binary map showing cultivated/non-cultivated areas ca. 2015. The map produced an overall accuracy of 94 percent with roughly 14 percent omission and commission errors for the cropland class based on a large set of independent validation samples. The map suggests the entire study area has a total 546 million hectares (Mha) of croplands occupying 18 percent of the land area. Comparison between national cropland area estimates from United Nations Food and Agricultural Organizations (FAO) and those derived from this work also showed an R-square value of 0.95. For the entire Landsat-derived 30-m product the overall accuracy was 93.8% with cropland class providing producers accuracy of 86.5% (errors of omissions = 13.5%) and users accuracy of 85.7% (errors of commissions = 14.3%). This Landsat-derived 30-m cropland product (GFSAD30) provided 10-30% greater cropland areas compared to UN FAO in the 64 Countries. Finally, the map-to-map comparison between GFSAD30 with several other cropland products revealed that the best similarity matrix was with the 30m global land cover (GLC30) product providing an overall accuracy of 88.8 percent (Kappa 0.7) with producers cropland similarity of 89.2 percent (errors of omissions = 10.8%) and users cropland similarity of 81.8 percent (errors of commissions = 8.1%). GFSAD30 captured the missing croplands in GLC30 product around significantly irrigated agricultural areas in Germany and Belgium and rainfed agriculture in Italy. This study also established that the real strength of GFSAD30 product, compared to other products, were in: 1. Identifying precise location of croplands, and 2. Capturing fragmented croplands. The cropland extent map dataset is available through NASAs Land Processes Distributed Active Archive Center (LP DAAC) at https://doi.org/10.5067/MEaSUREs/GFSAD/GFSAD30EUCEARUMECE.001, while the training and reference data as well as visualization are available at the Global Croplands website.

ISPRS Journal of Photogrammetry and Remote Sensing↗

The potential for species distribution models to distinguish source populations from sinks

1. While species distribution models (SDM) are frequently used to predict species occurrences to help inform conservation management, there is limited evidence evaluating whether habitat suitability can reliably predict intrinsic growth rates or distinguish source from sink populations. Filling this knowledge gap is critical for conservation science, as applications of SDMs for management purposes ultimately depend on these typically unobserved population or metapopulation dynamics. 2. Using regression, we associate previously published population level estimates of intrinsic growth and abundance derived from a Bayesian analysis of mark-recapture data for 17 bird species found in the contiguous United States with SDM habitat suitability estimates fitted here to opportunistic data for these same species. We then use AUC to measure how well SDMs can distinguish populations categorized as sources and sinks, depending on their intrinsic growth rates estimated from the mark-recapture data. We built SDMs using two different approaches, boosted regression trees (BRT) and Generalized Linear Models (GLM), and compared their predictive performance. Each SDM was built with presence points obtained from eBird and 10 environmental variables previously selected to model intrinsic growth rates and abundance for these species. 3. We show that SDMs built with opportunistic data are poor predictors of species demography in general; both BRT and GLM explained very little spatial variation of intrinsic growth rate and population abundance (median R2 across 17 species was close to 0.1 for both SDM methods). SDMs do, however, estimate higher suitability for source populations as compared to sinks. Out of 13 species which had both source and sink populations, both BRT and GLM had AUC values greater than 0.7 for 7 species when discriminating between sources and sinks. 4. Habitat suitability have the potential to be a useful measure to indicate a population’s ability to sustain itself as a source population, however more research on a diverse set of taxa is essential to fully explore this potential. This interpretation of habitat suitability can be particularly useful for conservation practice, and identification of explicit cases of when and how SDMs fail to match population demography can be informative for advancing ecological theory.

Journal of Animal Ecology↗

Soil moisture partitioning between under canopy and interspace environments in shrublands of the northern Chihuahuan Desert

Soil moisture is a key link between hydrologic and ecologic processes in desert shrublands. Understanding how soil moisture is spatially distributed in desert shrublands provides valuable insights into how shrubs use and impact limiting water resources, and how shrublands may respond to future meteorological and climate change. Our goals were to determine how soil moisture is partitioned between soil volumes under canopies and in the bare soil interspaces across multiple desert shrublands, and to evaluate the roles of physical soil properties, shrub-type characteristics, meteorology, and measurement resolution in influencing and observing variation in soil moisture partitioning. Utilizing two long-term soil moisture datasets (monthly resolution, 30 years, whole soil profile measurements; and 30 min resolution, 10 years, 10–30 cm measurements), we compared soil moisture partitioning across nine northern Chihuahuan Desert shrubland sites (three sites dominated by creosotebush [ Larrea tridentata ], three by honey mesquite [ Prosopis glandulosa ], and three by tarbush [ Flourensia cernua ]) in the Jornada Basin, southern New Mexico, USA. Over 30 years, monthly, whole profile data showed that soil moisture in mesquite shrublands was consistently higher in bare soil interspaces compared to under canopies, whereas soil moisture under and between shrubs was more similar in creosotebush and tarbush shrublands. Physical soil properties were linked as explanatory variables of long-term soil moisture partitioning (monthly whole profile dataset), whereas 30-minute data showed that shorter-term periods of higher precipitation promoted greater near surface soil moisture (10–30 cm) in bare soil interspaces that was not captured at monthly time steps. Thus, although the long-term average partitioning of soil moisture in these shrublands is strongly controlled by soil physical properties, soil moisture partitioning varies at shorter timescales (daily to weekly) in response to precipitation events. Moreover, shrub-type characteristics influenced soil moisture partitioning, with dense and tall mesquite shrubs having lower under canopy soil moisture than tarbush, and root architecture potentially influencing partitioning across creosotebush sites. These results illustrate diversity in soil moisture partitioning both between and within shrublands of the northern Chihuahuan Desert, and elucidate how physical soil properties, shrub-type characteristics, and meteorological variation interact to shape their soil moisture dynamics.

New Mexico↗

Critical mineral resources of the United States—Economic and environmental geology and prospects for future supply

Summary Mineral commodities are vital for economic growth, improving the quality of life, providing for national defense, and the overall functioning of modern society. Minerals are being used in larger quantities than ever before and in an increasingly diverse range of applications. With the increasing demand for a considerably more diverse suite of mineral commodities has come renewed recognition that competition and conflict over mineral resources can pose significant risks to the manufacturing industries that depend on them. In addition, production of many mineral commodities has become concentrated in relatively few countries (for example, tungsten, rare-earth elements, and antimony in China; niobium in Brazil; and platinum-group elements in South Africa and Russia), thus increasing the risk for supply disruption owing to political, social, or other factors. At the same time, an increasing awareness of and sensitivity to potential environmental and health issues caused by the mining and processing of many mineral commodities may place additional restrictions on mineral supplies. These factors have led a number of Governments, including the Government of the United States, to attempt to identify those mineral commodities that are viewed as most “critical” to the national economy and (or) security if supplies should be curtailed. This book presents resource and geologic information on the following 23 mineral commodities currently among those viewed as important to the national economy and national security of the United States: antimony (Sb), barite (barium, Ba), beryllium (Be), cobalt (Co), fluorite or fluorspar (fluorine, F), gallium (Ga), germanium (Ge), graphite (carbon, C), hafnium (Hf), indium (In), lithium (Li), manganese (Mn), niobium (Nb), platinum-group elements (PGE), rare-earth elements (REE), rhenium (Re), selenium (Se), tantalum (Ta), tellurium (Te), tin (Sn), titanium (Ti), vanadium (V), and zirconium (Zr). For a number of these commodities—for example, graphite, manganese, niobium, and tantalum—the United States is currently wholly dependent on imports to meet its needs. The first two chapters (A and B) deal with general information pertinent to the study of mineral resources. Chapters C through V describe individual mineral commodities and include an overview of current uses of the commodity, identified resources and their distribution nationally and globally, the state of current geologic knowledge, the potential for finding additional deposits nationally and globally, and geoenvironmental issues that may be related to the production and uses of the commodity. These chapters are updates of the commodity chapters published in 1973 in U.S. Geological Survey Professional Paper 820, “United States Mineral Resources.”

Professional Paper↗

Simulation of the Lower Walker River Basin hydrologic system, west-central Nevada, using PRMS and MODFLOW models

Walker Lake is a terminal lake in west-central Nevada with almost all outflow occurring through evaporation. Diversions from Walker River since the early 1900s have contributed to a substantial reduction in flow entering Walker Lake. As a result, the lake is receding, and salt concentrations have increased to a level in which Oncorhynchus clarkii henshawi (Lahontan Cutthroat trout) are no longer present, and the lake ecosystem is threatened. Consequently, there is a concerted effort to restore the Walker Lake ecosystem and fishery to a level that is more sustainable. However, Walker Lake is interlinked with the lower Walker River and adjacent groundwater system which makes it difficult to understand the full effect of upstream water-management actions on the overall hydrologic system including the lake level, volume, and dissolved-solids concentrations of Walker Lake. To understand the effects of water-management actions on the lower Walker River Basin hydrologic system, a watershed model and groundwater flow model have been developed by the U.S. Geological Survey in cooperation with the Bureau of Reclamation and the National Fish and Wildlife Foundation. The watershed model was developed using the precipitation runoff modeling system (PRMS) and the groundwater flow model was constructed using the MODular groundwater FLOW model (MODFLOW) and both were calibrated for the lower Walker River Basin. These models can be incorporated in an integrated Groundwater and Surface-water FLOW (GSFLOW) model of the lower Walker River Basin. Additionally, the MODFLOW model developed for this study is useful for efficiently simulating long-term and large-scale effects of water-management actions on groundwater hydrology, streamflow, and Walker Lake level, volume, and dissolved-solids concentrations. The lower Walker River Basin PRMS model (LWR_PRMS) was constructed using a subbasin approach to aid in development and calibration, and simulates a 30-year period from 1978 to 2007 using daily time steps. The LWR_PRMS was used to estimate the distribution of groundwater recharge specified in the MODFLOW model. The highest rates of groundwater recharge occur in the Wassuk Range beneath perennial and ephemeral stream channels, whereas lower rates of recharge occur beneath alluvial fans along mountain fronts. The total groundwater recharge estimated using PRMS was about 25,000 acre-feet per year. The lower Walker River Basin MODFLOW (LWR_MF) model simulates an 89-year period using monthly time steps. The LWR_MF was constructed with an initial steady-state simulation to represent dynamic equilibrium conditions from 1908 to 1918 and then a transient simulation representing the period 1919–2007. The model was calibrated using a combination of manual and automated methods of adjusting model parameters to minimize errors between model simulated results and weighted observations of groundwater levels, streamflows, and lake level. Hydrologic conditions simulated with the LWR_MF include the movement and change in storage of groundwater, and the water budgets for Walker River, Walker Lake, and the groundwater system. The LWR_MF computed dissolved-solids concentrations for Walker Lake using simulated lake volume and an assumed constant internal salt mass of 37.2 million tons. Effects of potential changes in water management on future conditions (scenarios) of the lower Walker River Basin hydrologic system and Walker Lake from 2011 to 2070 were evaluated. Several water-management scenarios were considered, including a baseline scenario that represents no changes in system management, improved irrigation efficiencies for the Walker River Indian Irrigation Project (WRIIP), a range of increased streamflows entering the lower Walker River Basin, and, the fallowing of fields on the WRIIP. For the baseline scenario, it was assumed that streamflow conditions from 1981 to 2010 will be repeated in the future. Results indicate that Walker Lake level and volume continue to decline but at a slower rate as the surface area of the lake becomes smaller and lake evaporation decreases. Dissolved-solids concentrations in Walker Lake continue to increase and increase much more rapidly during periods when minimal flows reach the lake due to a diminished lake volume. Alternatively, in years with high runoff, lake level increases are greater and dissolved-solids decreases are greater, compared with equivalent runoffs experienced during 1981–2010. The simulated effects of improving WRIIP efficiencies on Walker River streamflows, Walker Lake inflow, level, and dissolved-solids concentrations, and crop consumptive use, are compared with the baseline reference scenario for a range of irrigation efficiency improvements from 0 to 25 percent over 60 years. Results indicate that water is conserved through a reduction in irrigation-induced groundwater recharge and subsequent groundwater discharge through evapotranspiration. The conserved water mostly goes to increased streamflow to Walker Lake, followed by increased crop consumptive use, then increased evaporation from Weber Reservoir. The simulated effects of increased streamflows at Walker River at Wabuska streamgage (10301500) on Walker Lake inflow, level, and dissolved-solids concentrations, and crop consumptive use, are compared with the baseline scenario after 60 years under two different management methods for Weber Reservoir. Results indicate Walker Lake level and dissolved-solids concentrations stabilized with increased irrigation-season streamflow of about 40,000 acre-feet per year at the Walker River at Wabuska streamgage. Walker Lake level increased, and dissolved-solids concentration decreased, with increased flows of 50,000 acre-feet per year or more. After 60 years with additional irrigation-season streamflows of 50,000 acre-feet per year, Walker Lake level increased by about 48 feet, and lake dissolved-solids concentrations decreased by about 3,000 milligrams per liter (mg/L). With 75,000 acre-feet per year of additional streamflow, Walker Lake level increased by 70 feet, and dissolved-solids concentration decreased by 7,600 milligrams per liter. The effects of fallowing of Walker River Indian Irrigation Project fields from 2007 to 2010 on Walker Lake inflow, level, and dissolved solids were evaluated. Fallowing resulted in a near doubling of Walker River inflow to Walker Lake during this period, an increase in Walker Lake level of about 1.4 feet, and a decrease in dissolved-solids concentration of about 540 mg/L.

Nevada↗

Oligocene caldera complex and calc-alkaline tuffs and lavas of the Indian Peak volcanic field, Nevada and Utah

The Indian Peak volcanic field is representative of the more than 50,000 km 3 of ash-flow tuff and tens of calderas in the Great Basin that formed during the Oligocene-early Miocene "ignimbrite flareup" in southwestern North America. The field formed about 32 to 27 Ma in the southeastern Great Basin and consists of the centrally positioned Indian Peak caldera complex and a surrounding blanket of related ash-flow sheets distributed over an area of about 55,000 km 2 . The field has a volume on the order of 10,000 km 3 . A cluster of two obscure source areas and four calderas comprise the ∼80 x 120 km caldera complex. Only minor volumes of rhyolite and two pyroxene andesite lavas were extruded episodically throughout the lifetime of the magma system that formed the field, chiefly during its youth and old age. Six ash-flow sequences alternate between rhyolite and dacite in a volume ratio of about 1:8, and a culminating seventh is trachytic. The first, fourth, and sixth tuff units are of rhyolite that contains sparse to modest amounts of phenocrysts, chiefly plagioclase and biotite, and abundant lithic and pumice lapilli; these deposits are confined within the caldera complex and form multiple and compound cooling units that are normally zoned with respect to bulk chemical composition and crystal type, content, and size. The second, third, and fifth tuff sequences are of crystal-rich dacite that forms extensive simple cooling-unit outflow sheets and partial caldera fillings of compound cooling units. Each dacite unit contains similar amounts of plagioclase, biotite, hornblende, quartz, two pyroxenes, and Fe-Ti oxides; trace amounts of sanidine and titanite also occur in the youngest. Cognate inclusions in the dacites show only slight intra- and inter-unit differences in bulk chemical composition. The seventh eruptive sequence consists of several cooling units of trachydacite tuff containing small to modest amounts of plagioclase and two pyroxenes. These dominantly high-K calc-alkaline rocks are a record of the birth, maturation, and death of a large, open, continental magma system that was probably initiated and sustained by influx of mafic magma derived from a southward-migrating locus of magma production in the mantle. The small volumes of chemically diverse andesitic rocks were derived from separately evolving magma bodies but are modified representatives of the mantle power supply. Recurrent production of very large batches (some greater than 3,000 km 3 ) of quite uniform dacite magmas appears to have required combination of andesite magma and crustal silicic material in vigorously convecting chambers. Compositional data indicate that rhyolites are polygenetic. As the main locus of mantle magma production shifted southward, trachydacite magma could have been produced by fractionation of andesitic magma within the crust.

Nevada, Utah↗

Methods for estimating annual exceedance-probability discharges and largest recorded floods for unregulated streams in rural Missouri

Regression analysis techniques were used to develop a set of equations for rural ungaged stream sites for estimating discharges with 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities, which are equivalent to annual flood-frequency recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, respectively. Basin and climatic characteristics were computed using geographic information software and digital geospatial data. A total of 35 characteristics were computed for use in preliminary statewide and regional regression analyses. Annual exceedance-probability discharge estimates were computed for 278 streamgages by using the expected moments algorithm to fit a log-Pearson Type III distribution to the logarithms of annual peak discharges for each streamgage using annual peak-discharge data from water year 1844 to 2012. Low-outlier and historic information were incorporated into the annual exceedance-probability analyses, and a generalized multiple Grubbs-Beck test was used to detect potentially influential low floods. Annual peak flows less than a minimum recordable discharge at a streamgage were incorporated into the at-site station analyses. An updated regional skew coefficient was determined for the State of Missouri using Bayesian weighted least-squares/generalized least squares regression analyses. At-site skew estimates for 108 long-term streamgages with 30 or more years of record and the 35 basin characteristics defined for this study were used to estimate the regional variability in skew. However, a constant generalized-skew value of -0.30 and a mean square error of 0.14 were determined in this study. Previous flood studies indicated that the distinct physical features of the three physiographic provinces have a pronounced effect on the magnitude of flood peaks. Trends in the magnitudes of the residuals from preliminary statewide regression analyses from previous studies confirmed that regional analyses in this study were similar and related to three primary physiographic provinces. The final regional regression analyses resulted in three sets of equations. For Regions 1 and 2, the basin characteristics of drainage area and basin shape factor were statistically significant. For Region 3, because of the small amount of data from streamgages, only drainage area was statistically significant. Average standard errors of prediction ranged from 28.7 to 38.4 percent for flood region 1, 24.1 to 43.5 percent for flood region 2, and 25.8 to 30.5 percent for region 3. The regional regression equations are only applicable to stream sites in Missouri with flows not significantly affected by regulation, channelization, backwater, diversion, or urbanization. Basins with about 5 percent or less impervious area were considered to be rural. Applicability of the equations are limited to the basin characteristic values that range from 0.11 to 8,212.38 square miles (mi 2 ) and basin shape from 2.25 to 26.59 for Region 1, 0.17 to 4,008.92 mi 2 and basin shape 2.04 to 26.89 for Region 2, and 2.12 to 2,177.58 mi 2 for Region 3. Annual peak data from streamgages were used to qualitatively assess the largest floods recorded at streamgages in Missouri since the 1915 water year. Based on existing streamgage data, the 1983 flood event was the largest flood event on record since 1915. The next five largest flood events, in descending order, took place in 1993, 1973, 2008, 1994 and 1915. Since 1915, five of six of the largest floods on record occurred from 1973 to 2012.

Missouri↗

Discharge and sediment loads in the Boise River drainage basin, Idaho 1939-40

The Boise River project is a highly developed agricultural area comprising some 520 square miles of valley and bench lands in southwestern Idaho. Water for irrigation is obtained from the Boise River and its tributaries which are regulated by storage in Arrow Rock and Deer Flat reservoirs. Distribution of water to the farms is effected by 27 principal canals and several small farm laterals which divert directly from the river. The- New York Canal, which is the largest, not only supplies water to smaller canals and farm laterals, but also is used to fill Deer Flat Reservoir near Nampa from which water is furnished to farms in the lower valley. During the past 15 years maintenance costs in a number of those canals have increased due to deposition of sediment in them and in the river channel itself below the mouth of Moore Creek. Interest in determining the runoff and sediment loads from certain areas in the Boise River drainage basin led to an investigation by the Flood Control Coordinating Committee of the Department of Agriculture. Measurements of daily discharge and sediments loads were made by the Geological Survey at 13 stations in the drainage basin during the 18-month period ended June 30, 1940. The stations were on streams in areas having different kinds of vegetative cover and subjected to different kinds of land-use practice. Data obtained during the investigation furnish a basis for certain comparisons of runoff and sediment loads from several areas arid for several periods of time. Runoff measured at stations on the. Boise River near Twin Springs and on Moore Creek near Arrow Rock was smaller during 1939 than during 1940 and was below the average annual runoff for the period of available record. Runoff measured at the other stations on the project also was smaller during 1939 than during 1940 and probably did not exceed the average for the previous 25 years. The sediment loads measured during the spring runoff in 1939 were smaller at most stations than those measured during the spring runoff in 1940. At those stations where the flow was not affected, or only slightly affected, by upstream diversions or by placer-mining operations, the largest sadiment loads per unit of drainage area were measured in Grouse Creek during both 1939 and 1940, amounting to 3,460 and 2,490 tons per square mile, respectively, and the smallest loads per unit of drainage area were measured in Bannock Creek during 1939 and in the Boise River near Twin Springs during 1940, amounting to 14 and 83 tons per square mile, respectively. Size anaylses of a large number of samples of suspended and deposited sediments give an indication of the origin of sediments carried past some of the stations. The analyses show that most of the sediment measured at the five stations in the Moore Creek drainages basin above Idaho City consisted largely of coarse material. They show, also, that the sediment measured at the station on Moore Creek above Thorn Creek consisted almost entirely of fine material during practically the entire period of the investigation. Most of the coarse material passing the stations above Idaho City probably was retained behind the dikes or in the pools usually formed by tailings from dredging operations in the placer-mining area below Idaho City, and much of the fine material measured at the station on Moore Creek above Thorn Creek probably was contributed by placer-mining activity. During the years when the spring runoff is greater than that measured during 1939 and 1940, it is probable that the dikes and pools will be less effective in retaining coarse sediments within the placered area. Records of sediment loads measured in the New York Canal indicate that a negligible amount of sediment was deposited there during 1939, but that in 1940 from 10 to 15 percent of the total load at the gaging station consisted of coarse sediment which was later deposited on the canal bottom. Most of the fine material was doubtless carried through the canal and eventually deposited in diversion ditches and on farm land. Because the sediment carried past the station on Moore Creek above Thorn Creek consisted almost entirely of fine material, it is probable, that a considerable part of the coarse sediment carried in the New York Canal during the 1940 spring runoff period was scoured from the large bed of deposited material in the Boise River above Diversion- Dam, and that the remainder came from Grimes Creek. Arrow Rock Reservoir was not sluiced during the investigation, and it is therefore unlikely that any of the coarse sediment in the New York Canal came from the Boise River above Moore Creek during 1939 and 1940. The average dry weight of 71 samples of deposited sediments collected from several parts of the Boise River drainage basin is about 90 pounds per cubic foot. The average specific gravity of 77 samples of deposited sediments is 2.57.

Idaho↗