SEARCH · Geology Reports

Results for “Global Environmental Change”

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

SourceResearch connectionFrom yearThrough year

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

At least 775 records · Page 43Linked to original sources

USGS2001

All equal-area map projections are created equal, but some are more equal than others

High-resolution regional and global raster databases are currently being generated for a variety of environmental and scientific modeling applications. The projection of these data from geographic coordinates to a plane coordinate system is subject to significant areal error. Sources of error include users selecting an inappropriate projection or incorrect parameters for a given projection, algorithmic errors in commercial geographic information system (GIS) software, and errors resulting from the projection of data in the raster format. To assess the latter type of errors, the accuracy of raster projection was analyzed by two methods. First, a set of 12 one-degree by one-degree quadrilaterals placed at various latitudes was projected at several raster resolutions and compared to the projection of a vector representation of the same quadrilaterals. Second, several different raster resolutions of land cover data for Asia were projected and the total areas of 21 land cover categories were tabulated and compared. While equal-area projections are designed to specifically preserve area, the comparison of the results of the one-degree by one-degree quadrilaterals with the common equal area projections (e.g., the Mollweide) indicates a considerable variance in the one-degree area after projection. Similarly, the empirical comparison of land cover areas for Asia among various projections shows that total areas of land cover vary with projection type, raster resolution, and latitude. No single projection is best for all resolutions and all latitudes. While any of the equal-area projections tested are reasonably accurate for most applications with resolutions of eight-kilometer pixels or smaller, significant variances in accuracies appear at larger pixel sizes.

Cartography and Geographic Information Science↗

USGS2020

Millennial-scale climate and human drivers of environmental change and fire activity in a dry, mixed-conifer forest of northwestern Montana

Warm summer temperatures and longer fire seasons are promoting larger, and in some cases, more fires that are severe in low- and mid-elevation, dry mixed-conifer forests of the Northern Rocky Mountains (NRM). Long-term historical fire conditions and human influence on past fire activity are not well understood for these topographically and biophysically heterogeneous forests. We developed reconstructions of millennial-scale fire activity, vegetation change, and human presence at Black Lake, a small closed-basin lake on the Flathead Indian Reservation in the Mission Valley, Northwestern Montana, United States. Fossil pollen, charcoal, and biomarkers associated with human presence were used to evaluate the interaction between climate variability, fire activity, vegetation change and human activity for the past 7000 years. Comparisons among multiple proxies suggest climate variability acted as the primary control on fire activity and vegetation change from the early Holocene until the late Holocene when records suggest fire activity and climate variability decoupled. Specific biomarkers (5β-stanols including coprostanol and epi-coprostanol) associated with human presence indicate humans were present within the Black Lake watershed for thousands of years, although the inferred intensity of human presence is highly variable. A strong relationship between climate variability and fire activity during the early and mid-Holocene weakens during the last few thousand years, suggesting possible increased influence of humans in mediating fire activity in recent millennia, and/or a shift in the interaction between the distribution and abundance of woody fuel and fire severity. Human-set fires during the cooler and wetter late Holocene may have been aimed at maintaining important cultural resources associated with the heterogeneous mosaic of mixed conifer forests within the Black Lake watershed. The paleoenvironmental reconstruction at Black Lake corroborates archeological records that show humans were present within the Black Lake watershed for over 7000 years. Further research is needed to evaluate the evidence for this continuous presence and the possible role that people played in shaping fire regimes and vegetation within low- to mid-elevation mixed-conifer ecosystems of the NRM.

Montana↗

USGS2022

Invasion of annual grasses following wildfire corresponds to maladaptive habitat selection by a sagebrush ecosystem indicator species

Numerous wildlife species within semi-arid shrubland ecosystems across western North America are experiencing substantial habitat loss and fragmentation. These changes in habitat are often attributed to a diverse suite of factors including prolonged and increasingly severe droughts, conifer expansion, anthropogenic development, domestic and feral livestock grazing, and invasion of exotic annual grasses, which promotes increased wildfire frequency and severity. Greater sage-grouse ( Centrocercus urophasianus ; hereafter, sage-grouse) are considered an indicator of sagebrush ecosystem health and have experienced widespread population decline associated with habitat loss and degradation, as well as changes in predator communities. Our objectives were to model and map sage-grouse habitat selection and survival during the important brood-rearing life stage in relation to landscape-scale environmental predictors. Furthermore, we sought to understand impacts of wildfire and annual grass invasion on brood habitat, as these accelerated disturbance regimes are a primary cause of habitat loss within the Great Basin region of the USA. We used a hierarchical Bayesian modeling framework to estimate resource selection functions and survival for early and late brood-rearing stages of sage-grouse in relation to a broad suite of habitat characteristics evaluated at multiple spatial scales within the Great Basin from 2009 to 2019. Sage-grouse selected for greater perennial grass cover, higher relative elevations, and areas closer to springs and wet meadows during both early and late brood-rearing. Terrain characteristics, including heat load and aspect, were important in survival models, as was variation in shrub height. We also found strong evidence for higher survival for both early and late broods within previously burned areas, but survival within burned areas decreased as annual grass cover (i.e. cheatgrass , Bromus tectorum ) increased. This interaction effect demonstrates how invasion of annual grasses into burned areas, which has become prevalent in Great Basin sagebrush ecosystems, can lead to maladaptive habitat selection by brood-rearing greater sage-grouse. Understanding these complex relationships aids wildlife conservation and habitat management as wildfire and annual grass cycles continue to accelerate across western ecosystems.

California, Nevada↗

USGS2022

Grassy–herbaceous land moderates regional climate effects on honey bee colonies in the Northcentral US

The lack of seasonally sustained floral resources (i.e. pollen and nectar) is considered a primary global threat to pollinator health. However, the ability to predict the abundance of flowering resources for pollinators based upon climate, weather, and land cover is difficult due to insufficient monitoring over adequate spatial and temporal scales. Here we use spatiotemporally distributed honey bee hive scales that continuously measure hive weights as a standardized method to assess nectar intake. We analyze late summer colony weight gain as the response variable in a random forest regression model to determine the importance of climate, weather, and land cover on honey bee colony productivity. Our random forest model predicted resource acquisition by honey bee colonies with 71% accuracy, highlighting the detrimental effects of warm, wet regions in the Northcentral United States on nectar intake, as well as the detrimental effect of years with high growing degree day accumulation. Our model also predicted that grassy–herbaceous natural land had a positive effect on the summer nectar flow and that large areas of natural grassy–herbaceous land around apiaries can moderate the detrimental effects of warm, wet climates. These patterns characterize multi-scale ecological processes that constrain the quantity and quality of pollinator nutritional resources. That is, broad climate conditions constrain regional floral communities, while land use and weather act to further modify the quantity and quality of pollinator nutritional resources. Observing such broad-scale trends demonstrates the potential for utilizing hive scales to monitor the effects of climate change on landscape-level floral resources for pollinators. The interaction of climate and land use also present an opportunity to manage for climate-resilient landscapes that support pollinators through abundant floral resources under climate change.

Michigan, Minnesota, North Dakota, Ohio, Pennsylva↗

USGS2023

Spatial asynchrony and cross-scale climate interactions in populations of a coldwater stream fish

Climate change affects populations over broad geographic ranges due to spatially autocorrelated abiotic conditions known as the Moran effect. However, populations do not always respond to broad-scale environmental changes synchronously across a landscape. We combined multiple datasets for a retrospective analysis of time-series count data (5–28 annual samples per segment) at 144 stream segments dispersed over nearly 1,000 linear kilometers of range to characterize the population structure and scale of spatial synchrony across the southern native range of a coldwater stream fish (brook trout, Salvelinus fontinalis ), which is sensitive to stream temperature and flow variations. Spatial synchrony differed by life stage and geographic region: it was stronger in the juvenile life stage than in the adult life stage and in the northern sub-region than in the southern sub-region. Spatial synchrony of trout populations extended to 100–200 km but was much weaker than that of climate variables such as temperature, precipitation, and stream flow. Early life stage abundance changed over time due to annual variation in summer temperature and winter and spring stream flow conditions. Climate effects on abundance differed between sub-regions and among local populations within sub-regions, indicating multiple cross-scale interactions where climate interacted with local habitat to generate only a modest pattern of population synchrony over space. Overall, our analysis showed higher degrees of response heterogeneity of local populations to climate variation and consequently population asynchrony than previously shown based on analysis of individual, geographically restricted datasets. This response heterogeneity indicates that certain local segments characterized by population asynchrony and resistance to climate variation could represent unique populations of this iconic native coldwater fish that warrant targeted conservation. Advancing the conservation of this species can include actions that identify such priority populations and incorporate them into landscape-level conservation planning. Our approach is applicable to other widespread aquatic species sensitive to climate change.

Global Change Biology↗

USGS2005

Mapping stream habitats with a global positioning system: Accuracy, precision, and comparison with traditional methods

We tested the precision and accuracy of the Trimble GeoXT??? global positioning system (GPS) handheld receiver on point and area features and compared estimates of stream habitat dimensions (e.g., lengths and areas of riffles and pools) that were made in three different Oklahoma streams using the GPS receiver and a tape measure. The precision of differentially corrected GPS (DGPS) points was not affected by the number of GPS position fixes (i.e., geographic location estimates) averaged per DGPS point. Horizontal error of points ranged from 0.03 to 2.77 m and did not differ with the number of position fixes per point. The error of area measurements ranged from 0.1% to 110.1% but decreased as the area increased. Again, error was independent of the number of position fixes averaged per polygon corner. The estimates of habitat lengths, widths, and areas did not differ when measured using two methods of data collection (GPS and a tape measure), nor did the differences among methods change at three stream sites with contrasting morphologies. Measuring features with a GPS receiver was up to 3.3 times faster on average than using a tape measure, although signal interference from high streambanks or overhanging vegetation occasionally limited satellite signal availability and prolonged measurements with a GPS receiver. There were also no differences in precision of habitat dimensions when mapped using a continuous versus a position fix average GPS data collection method. Despite there being some disadvantages to using the GPS in stream habitat studies, measuring stream habitats with a GPS resulted in spatially referenced data that allowed the assessment of relative habitat position and changes in habitats over time, and was often faster than using a tape measure. For most spatial scales of interest, the precision and accuracy of DGPS data are adequate and have logistical advantages when compared to traditional methods of measurement. ?? 2006 Springer Science+Business Media, Inc.

Environmental Management↗

USGS2020

Biogeography of fire regimes in western US conifer forests: A trait-based approach

Aim Functional traits are a critical link between species distributions and the ecosystem processes that structure those species’ niches. Concurrent increases in the availability of functional trait data and our ability to model species distributions present an opportunity to develop functional trait biogeography, i.e. the mapping of functional traits across space. Functional trait biogeography can improve process-based predictions about the resistance of certain species assemblages to changing environmental conditions across landscape scales. We illustrate this concept by developing the first trait-based, quantitative ranking of fire resistance (adult tree survival) in North American conifer species, and mapping that fire resistance across space. Location and Time period Western Continental United States, present-day. Major taxa studied 29 common conifer tree species. Methods We compiled six traits for each species: three relating to tree morphology and three relating to litter flammability. We combined these traits into a single fire resistance score, and used community-weighted averaging to estimate the fire resistance scores of different forest communities, using interpolated species distribution and relative abundance data. Results Species associated with historically frequent fire have high fire resistance scores (e.g., Pinus ponderosa), reflected by thick bark, tall crowns, and flammable litter. Species associated with subalpine or arid conditions have low fire resistance scores (e.g., Picea engelmannii and Pinus edulis), reflected by thin bark, short stature, poor self-pruning and low litter flammability. A map of forest community fire resistance across the western US reveals agreement with independent assessments of historical fire regimes, while also identifying areas where community-wide species traits may be mismatched with historical fire regimes. Main conclusions Quantifying the functional traits that confer resistance to tree-killing fire provides a direct link between ecosystem disturbance and community resistance. Understanding this link is critical to evaluating long-term resilience of different forest types under dynamic fire regimes. Our work represents the first known spatial representation of fire-resistance traits at a regional scale, and as such provides a link between functional traits and biogeography relevant to a critical ecosystem process.

Global Ecology and Biogeography↗

USGS2005

Mercury recycling in the United States in 2000

Reclamation and recycling of mercury from used mercury- containing products and treatment of byproduct mercury from gold mining is vital to the continued, though declining, use of this metal. Mercury is reclaimed from mercury-containing waste by treatment in multistep high-temperature retorts-the mercury is volatized and then condensed for purification and sale. Some mercury-containing waste, however, may be landfilled, and landfilled material represents loss of a recyclable resource and a threat to the environment. Related issues include mercury disposal and waste management, toxicity and human health, and regulation of mercury releases in the environment. End-users of mercury-containing products may face fines and prosecution if these products are improperly recycled or not recycled. Local and State environmental regulations require adherence to the Resource Conservation and Recovery Act and the Comprehensive Environmental Response, Compensation, and Liability Act to regulate generation, treatment, and disposal of mercury-containing products. In the United States, several large companies and a number of smaller companies collect these products from a variety of sources and then reclaim and recycle the mercury. Because mercury has not been mined as a principal product in the United States since 1992, mercury reclamation from fabricated products has become the main source of mercury. Principal product mercury and byproduct mercury from mining operations are considered to be primary materials. Mercury may also be obtained as a byproduct from domestic or foreign gold-processing operations. In the early 1990s, U.S. manufacturers used an annual average that ranged from 500 to 600 metric tons of recycled and imported mercury for fabrication of automobile convenience switches, dental amalgam, fluorescent lamps, medical uses and thermometers, and thermostats. The amount now used for fabrication is estimated to be 200 metric tons per year or less. Much of the data on mercury is estimated because it is a low-volume commodity and its production, use, and disposal is difficult to track. The prices and volumes of each category of mercury-containing material may change dramatically from year to year. For example, the average price of mercury was approximately $150 per flask from 2000 until 2003 and then rose sharply to $650 per flask in fall 2004 and approximately $850 per flask in spring 2005. Since 1927, the common unit for measuring and pricing mercury has been the flask in order to conform to the system used at Almaden, Spain (Meyers, 1951). One flask weighs 34.5 kilograms, and 29 flasks of mercury are contained in a metric ton. In the United States, the chlorine-caustic soda industry, which is the leading end-user of elemental mercury, recycles most of its mercury in-plant as home scrap. Annual purchases of replacement mercury by the chlorine-caustic soda industry indicate that some mercury may be lost through evaporation to the environment, put into a landfill as industrial waste, or trapped within pipes in the plant. Impending closure of domestic and foreign mercury-cell chlorine-caustic soda plants and the shift to nonmercury technology for chlorine-caustic soda production could ultimately result in a significant volume of elemental mercury for recycling, sale, or storage. Globally, mercury is widely used in artisanal, or small-scale, gold mining. Most of that mercury is lost to the environment and is not recycled. The recycling rate for mercury was not available owing to insufficient data in 2000, and the efficiency of mercury recycling was estimated to be 62 percent.

Circular↗

USGS2005

Mercury recycling in the United States in 2000

Open-File Report↗

USGS2021

Bayesian change point quantile regression approach to enhance the understanding of shifting phytoplankton-dimethyl sulfide relationships in aquatic ecosystems

Dimethyl sulfide (DMS) serves as an anti-greenhouse gas, plays multiple roles 7 in aquatic ecosystems, and contributes to the global sulfur cycle. The chlorophyll 8 a (CHL, an indicator of phytoplankton biomass)-DMS relationship is critical for 9 estimating DMS emissions from aquatic ecosystems. Importantly, recent research has 10 identified that the CHL-DMS relationship has a breakpoint, where the relationship 11 is positive below a CHL threshold and negative at higher CHL concentrations. 12 Conventionally, mean regression methods are employed to characterize the CHL-DMS 13 relationship. However, these approaches focus on the response of mean conditions 14 and cannot illustrate responses of other parts of the DMS distribution, which could 15 be important in order to obtain a complete view of the CHL-DMS relationship. In 16 this study, for the first time, we proposed a novel Bayesian change point quantile 17 regression (BCPQR) model that integrates and inherits advantages of Bayesian change 18 point models and Bayesian quantile regression models. Our objective was to examine 19 whether or not the BCPQR approach could enhance the understanding of shifting 20 CHL-DMS relationships in aquatic ecosystems. We fitted BCPQR models at five 21 regression quantiles for freshwater lakes and for seas. We found that BCPQR models 22 could provide a relatively complete view on the CHL-DMS relationship. In particular, 23 it quantified the upper boundary of the relationship, representing the limiting effect of 24 CHL on DMS. Based on the results of paired parameter comparisons, we revealed the 25 inequality of regression slopes in BCPQR models for seas, indicating that applying 26 the mean regression method to develop the CHL-DMS relationship in seas might not 27 be appropriate. We also confirmed relationship differences between lakes and seas at 28 multiple regression quantiles. Further, by introducing the concept of DMS emission 29 potential, we found that pH was not likely a key factor leading to the change of the 30 CHL-DMS relationship in lakes. These findings cannot be revealed using piecewise 31 linear regression. We thereby concluded that the BCPQR model does indeed enhance 32 the understanding of shifting CHL-DMS relationships in aquatic ecosystems and is 33 expected to benefit efforts aimed at estimating DMS emissions. Considering that 34 shifting (threshold) relationships are not rare and that the BCPQR model can easily 35 be adapted to different systems, the BCPQR approach is expected to have great 36 potential for generalization in other environmental and ecological studies.

Water Research↗

USGS1999

Demography of the Yellowstone grizzly bears

We undertook a demographic analysis of the Yellowstone grizzly bears ( Ursus arctos ) to identify critical environmental factors controlling grizzly bear vital rates, and thereby to help evaluate the effectiveness of past management and to identify future conservation issues. We concluded that, within the limits of uncertainty implied by the available data and our methods of data analysis, the size of the Yellowstone grizzly bear population changed little from 1975 to 1995. We found that grizzly bear mortality rates are about double in years when the whitebark pine crop fails than in mast years, and that the population probably declines when the crop fails and increases in mast years. Our model suggests that natural variation in whitebark pine crop size over the last two decades explains more of the perceived fluctuations in Yellowstone grizzly population size than do other variables. Our analysis used demographic data from 202 radio-telemetered bears followed between 1975 and 1992 and accounted for whitebark pine ( Pinus albicaulis ) crop failures during 1993–1995. We used a maximum likelihood method to estimate demographic parameters and used the Akaike Information Criteria to judge the significance of various independent variables. We identified no independent variables correlated with grizzly bear fecundity. In order of importance, we found that grizzly bear mortality rates are correlated with season, whitebark pine crop size (mast vs. nonmast year), sex, management-trapping status (never management-trapped vs. mangement-trapped once or more), and age. The mortality rate of bears that were management-trapped at least once was almost double that of bears that were never management-trapped, implying a source/sink (i.e., never management-trapped/management-trapped) structure. The rate at which bears move between the source and sink, estimated as the management-trapping rate ( h ), is critical to estimating the finite rate of increase, λ͞ . We quantified h by estimating the rate at which bears that have never been management-trapped are management-trapped for the first time. It differed across seasons, was higher in nonmast than mast years, and varied with age. We calculate that λ͞ = 1.00 from 1975 to 1983 (four mast and five nonmast years) and 1.02 from 1984 to 1995 (seven mast and five nonmast years). Overall, we find that λ͞ = 1.01 ± 0.04 (mean ± 1 se ) from 1975 to 1995. Our models suggest that future management should concentrate on the threats to whitebark pine, such as those posed by white pine blister rust, global warming, and fire suppression. As is currently widely recognized by Yellowstone land managers, our model also suggests that future management must compensate for the increased grizzly bear mortality that is likely to be caused by an increasing number of humans in Yellowstone.

Idaho, Montana, Wyoming↗

USGS2013

Wetlands serve as natural sources for improvement of stream ecosystem health in regions affected by acid deposition

For over 40 years, acid deposition has been recognized as a serious international environmental problem, but efforts to restore acidified streams and biota have had limited success. The need to better understand the effects of different sources of acidity on streams has become more pressing with the recent increases in surface water organic acids, or 'brownification' associated with climate change and decreased inorganic acid deposition. Here, we carried out a large scale multi-seasonal investigation in the Adirondacks, one of the most acid-impacted regions in the United States, to assess how acid stream producers respond to local and watershed influences and whether these influences can be used in acidification remediation. We explored the pathways of wetland control on aluminum chemistry and diatom taxonomic and functional composition. We demonstrate that streams with larger watershed wetlands have higher organic content, lower concentrations of acidic anions, and lower ratios of inorganic to organic monomeric aluminum, all beneficial for diatom biodiversity and guilds producing high biomass. Although brownification has been viewed as a form of pollution, our results indicate that it may be a stimulating force for biofilm producers with potentially positive consequences for higher trophic levels. Our research also reveals that the mechanism of watershed control of local stream diatom biodiversity through wetland export of organic matter is universal in running waters, operating not only in hard streams, as previously reported, but also in acid streams. Our findings that the negative impacts of acid deposition on Adirondack stream chemistry and biota can be mitigated by wetlands have important implications for biodiversity conservation and stream ecosystem management. Future acidification research should focus on the potential for wetlands to improve stream ecosystem health in acid-impacted regions and their direct use in stream restoration, for example, through stream rechanneling or wetland construction in appropriate hydrologic settings.

Global Change Biology↗

USGS2018

Spring temperature, migration chronology, and nutrient allocation to eggs in three species of arctic‐nesting geese: Implications for resilience to climate warming

The macronutrients that Arctic herbivores invest in their offspring are derived from endogenous reserves of fat and protein (capital) that females build prior to the period of investment or from foods they consume concurrently with investment (income). The relative contribution from each source can be influenced by temporal and environmental constraints on a female's ability to forage on Arctic breeding areas. Warming temperatures and advancing Arctic phenology may alter those constraints. From 2011–2014, we examined relationships among spring temperature, timing of migration and reproduction, and the sources of nutrients females deposited in eggs for three sympatric species of geese that nested in northern Alaska. Compared to lesser snow geese ( Anser caerulescens caerulescens ) and greater white‐fronted geese ( Anser albifrons frontalis ), black brant ( Branta bernicla nigricans ) were more likely to initiate follicle development during migration, resulting in fewer days between their arrival in the Arctic and the onset of incubation, and requiring a relatively greater capital investment in eggs. Delaying follicle development until after their arrival in the Arctic provided snow geese and white‐fronted geese an opportunity to forage near their nesting area and to deposit exogenous nutrients in eggs. With warmer spring temperatures brant invested more capital in eggs but snow geese invested less capital. Brant likely used capital to meet costs associated with earlier onset of follicle development when phenology was advanced, whereas snow geese used capital to compensate for poor foraging conditions during colder Arctic springs. Global warming is likely to reduce the quality of lower‐latitude marine habitats where brant acquire endogenous reserves, and advancing Arctic phenology may increase their reliance on those reserves during reproduction. Near‐term warming in northern Alaska may improve foraging conditions and favor the reproductive strategies of some herbivores such as snow geese and white‐fronted geese that mainly invest Arctic nutrients in their offspring.

Global Change Biology↗

USGS2025

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming↗

USGS2012

Water quality and landscape processes of four watersheds in eastern Puerto Rico

Humid tropical regions occupy about a quarter of Earth's land surface, yet they contribute a substantially higher fraction of the water, solutes, and sediment discharged to the world's oceans. Nearly half of Earth's population lives in the tropics, and development stresses can potentially harm soil resources, water quality, and water supply and in addition increase landslide and flood hazards. Owing to Puerto Rico's steep topography, low water storage capacity, and dependence on trade-wind precipitation, the island's people, ecosystems, and water supply are vulnerable to extreme weather such as hurricanes, floods, and droughts. Eastern Puerto Rico offers a natural laboratory for separating geologic and land-cover influences from regional- and global-scale influences because of its various bedrock types and the changing land cover surrounding intact, mature forest of the Luquillo Experimental Forest. Accordingly, a multiyear assessment of hydrological and biogeochemical processes was designed to develop an understanding of the effects of these differences on local climate, streamflow, water quality, and ecosystems, and to form the basis for a long-term and event-based program of climate and hydrologic monitoring. Because infrequent, large storms play a major role in this landscape, we focused on high-runoff events, sampling 263 storms, including all major hurricanes from 1991 through 2005. The largest storms have profound geomorphic consequences, such as landslides, debris flows, deep gullying on deforested lands, excavation and suspension of sediment in stream channels, and delivery of a substantial fraction of annual stream sediment load. Large storms sometimes entrain ocean foam and spray causing high concentrations of seasalt-derived constituents in stream waters during the storm. Past deforestation and agricultural activities in the Cayaguás and Canóvanas watersheds accelerated erosion and soil loss, and this material continues to be remobilized during large storms. Nearly 5,000 routine and event samples were analyzed for parameters that allow determination of denudation rates based on suspended and dissolved loads; 860 of these samples were analyzed for a comprehensive suite of chemical constituents. The rivers studied are generally similar in water-quality characteristics, and windward or leeward aspect appears to exert a stronger influence on water quality than geology or land cover. Of samples analyzed for comprehensive chemistry and for sediment, 543 were collected at runoff rates greater than 1 millimeter per hour, 256 at rates exceeding 10 millimeters per hour, and 3 at rates exceeding 90 millimeters per hour. Streams have rarely been sampled during events with such high runoff rates. Rates of physical and chemical weathering are especially high, and physical denudation rates, forested watersheds included, are considerably greater than is expected for a steady-state system. The elevated physical erosion drives an increased particulate organic carbon flux, one that is large, important to the carbon cycle, and sustainable, because soil-carbon regeneration is rapid. The 15-year Water, Energy, and Biogeochemical Budget dataset, which includes discharge, field parameters, suspended sediment, major cations and anions, and nutrients, is available from the U.S. Geological Survey's National Water Information System (http://waterdata.usgs.gov/nwis). The dataset provides a baseline for characterizing future environmental change and will improve our understanding of the interdependencies of land, water, and biological resources and their responses to changes in climate and land use. Because eastern Puerto Rico resembles many tropical regions in terms of geology and patterns of development, implications from this study are transferable to other tropical regions facing deforestation, rapid land-use change, and climate change.

Professional Paper↗

USGS2009

Quality characteristics of ground water in the Ozark aquifer of northwestern Arkansas, southeastern Kansas, southwestern Missouri and northeastern Oklahoma, 2006-07

Because of water quantity and quality concerns within the Ozark aquifer, the State of Kansas in 2004 issued a moratorium on most new appropriations from the aquifer until results were made available from a cooperative study between the U.S. Geological Survey and the Kansas Water Office. The purposes of the study were to develop a regional ground-water flow model and a water-quality assessment of the Ozark aquifer in northwestern Arkansas, southeastern Kansas, southwestern Missouri, and northeastern Oklahoma (study area). In 2006 and 2007, water-quality samples were collected from 40 water-supply wells completed in the Ozark aquifer and spatially distributed throughout the study area. Samples were analyzed for physical properties, dissolved solids and major ions, nutrients, trace elements, and selected isotopes. This report presents the results of the water-quality assessment part of the cooperative study. Water-quality characteristics were evaluated relative to U.S. Environmental Protection Agency drinking-water standards. Secondary Drinking-Water Regulations were exceeded for dissolved solids (11 wells), sulfate and chloride (2 wells each), fluoride (3 wells), iron (4 wells), and manganese (2 wells). Maximum Contaminant Levels were exceeded for turbidity (3 wells) and fluoride (1 well). The Maximum Contaminant Level Goal for lead (0 milligrams per liter) was exceeded in water from 12 wells. Analyses of isotopes in water from wells along two 60-mile long ground-water flow paths indicated that water in the Ozark aquifer was at least 60 years old but the upper age limit is uncertain. The source of recharge water for the wells along the flow paths appeared to be of meteoric origin because of isotopic similarity to the established Global Meteoric Water Line and a global precipitation relation. Additionally, analysis of hydrogen-3 (3H) and carbon-14 (14C) indicated that there was possible leakage of younger ground water into the lower part of the Ozark aquifer. This may be caused by cracks or fissures in the confining unit that separates the upper and lower parts of the aquifer, poorly constructed or abandoned wells, or historic mining activities. Analyses of major ions in water from wells along the flow paths indicated a transition from freshwater in the east to saline water in the west. Generally, ground water along flow paths evolved from a calcium magnesium bicarbonate type to a sodium calcium bicarbonate or a sodium calcium chloride bicarbonate type as water moved from recharge areas in Missouri into Kansas. Much of this evolution occurred within the last 20 to 25 miles of the flow paths along a water-quality transition zone near the Kansas-Missouri State line and west. The water quality of the Kansas part of the Ozark aquifer is degraded compared to the Missouri part. Geophysical and well-bore flow information and depth-dependent water-quality samples were collected from a large-capacity (1,900-2,300 gallons per minute) municipal-supply well to evaluate vertical ground-water flow accretion and variability in water-quality characteristics at different levels. Although the 1,050-foot deep supply well had 500 feet of borehole open to the Ozark aquifer, 77 percent of ground-water flow entering the borehole came from two 20-foot thick rock layers above the 1,000-foot level. For the most part, water-quality characteristics changed little from the deepest sample to the well-head sample, and upwelling of saline water from deeper geologic formations below the well was not evident. However, more saline water may be present below the bottom of the well.

Arkansas, Kansas, Missouri, Oklahoma↗

USGS2020

Holocene paleoclimate change in the western US: The importance of chronology in discerning patterns and drivers

Sediment in lakes and meadows forms a powerful archive that can be used to reconstruct environmental change through time. Reconstructions of lake level, of chemical, biological, and hydrological conditions, and of surrounding vegetation provide detailed information about past climate conditions, both locally and regionally. Indeed, most of our current knowledge of centennial- to millennial-scale climate variability in the arid western United States, where information about past hydroclimate is particularly important, comes from such sediment-based reconstructions. The pressing need for robust, precise predictions of future conditions is a significant motivation for paleoclimate science, and current research questions frequently require Holocene reconstructions to be resolved at sub-centennial timescales. Increasingly, regional syntheses seek to identify synoptic-scale patterns similar to those defined from modern observations (seasonal, interannual, multi-decadal, etc.) or to compare with the output of climate model simulations. However, the age control on existing records, especially those more than about 20 years old, is often sufficient only for millennial-scale interpretation. Here we assess the age control for 84 published and unpublished records from lakes and meadows in the Great Basin, California, and desert southwest, and use Bayesian modeling to evaluate the 95% uncertainty ranges for the 42 best-dated records. In the Late Holocene, about half of the 42 records have <400-year mean uncertainty ranges; however, high-precision age control is especially critical for young records, used to develop an accurate understanding of a proxy’s response to known climate variations. In the Middle Holocene, records vary from 400 to >800-year mean uncertainty and records of the Early Holocene have 600- to >1400-year mean uncertainty ranges. We find that the largest control on modeled uncertainties is dating density, with at least 2 dates/kyr being optimal and suggest obtaining “range-finder” dates at the onset of a study to better predict the total number of dates needed for an adequate age model. Such a density avoids a commonly observed phenomenon of significant peaks in uncertainty arising in gaps between age control points. Analysis of the uncertainties associated with proxy shifts reveal that more than half are >400 years. Although such large uncertainties currently prevent sub-centennial interpretations in most cases, increased dating density, strategic use of limited funds (including budgeting for a 2 date/kyr minimum at the proposal stage), construction of age-depth models with Bayesian methods, and critical evaluation of chronological uncertainty will shed light on past climate variability at finer timescales, enhancing our understanding of global and regional drivers of western U.S. climate.

California, Idaho, Nevada, Oregon, Utah↗

USGS2018

Global sea-level contribution from Arctic land ice: 1971 to 2017

The Arctic Monitoring and Assessment Program (AMAP) (AMAP, 2017) identifies the Arctic as the largest regional source of land ice to global sea-level rise in the 2003 to 2014 period. Yet, this contextualization ignores the longer perspective from in-situ records of glacier mass balance. Here, using 18 (> 55 °N latitude) glacier and ice cap mass balance series in the 1971 to 2017 period, we develop a semi-empirical estimate of annual sealevel contribution from seven Arctic regions by scaling the in-situ records to GRACE averages. We contend that our estimate represents the most accurate mass balance assessment so far available before the 1992 start of satellite altimetry. We estimate the 1971 to 2017 eustatic sea-level contribution from land ice north of ~55° N to be 23.0±12.3 mm sea-level equivalent (SLE). In all regions, the cumulative sealevel rise curves exhibit an acceleration, especially after 1988. Greenland is the source of 46% of the Arctic sea-level rise contribution (10.6±7.3 mm), followed by Alaska (5.7±2.2 mm), Arctic Canada (3.2±0.7 mm) and the Russian High Arctic (1.5±0.4 mm). Our annual results exhibit co-variability over a 43 year overlap (1971 to 2013) with the alternative dataset of Marzeion et al (2015) (M15). However, we find a 1.36x lower sea-level contribution, in agreement with satellite gravimetry. The IPCC Fifth Assessment report identified constraining the pre-satellite era sealevel budget as a topic of low scientific understanding that we address and specify sealevel contributions coinciding with IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (SROCC) “present day” (2005-2015) and “recent past” (1986-2005) reference periods. We assess an Arctic land ice loss of 8.3 mm SLE during the recent past and 12.4 mm SLE during the present day.

Environmental Research Letters↗