USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrological Processes”

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 901 records · Page 50Linked to original sources

Hydrologic response of groundwater and streamflow to natural and anthropogenic drivers of change in headwaters of the upper Colorado River basin during recent wet (1982–1999) and drought (2000–2022) conditions

Study region: Headwaters of the upper Colorado River basin (UCOL), USA Study focus: Surface-water and groundwater numerical models incorporating water-use information were used to investigate changes in climate, water use, and simulated hydrologic responses of snow processes, evapotranspiration, groundwater, and streamflow during recent wet (1982–1999) and drought (2000–2022) periods in the headwater subregions of the upper Colorado River basin. New hydrologic insights for the region: Decreases in average streamflow between wet and drought periods ranged from 20 % in the Colorado River headwaters subregion to 23 % in the Gunnison River headwaters subregion. Like streamflow, average surface runoff was statistically less during the drought than the wet period, with decreases from 24–31 % in the headwaters. On a volume basis, runoff decreases were greater than streamflow decreases in both the Colorado River and Gunnison River headwaters. Although the amount of water-year groundwater discharge to streams remained nearly the same between the wet and drought periods, groundwater as a percentage of streamflow increased between the wet and drought periods, highlighting the importance of groundwater in sustaining streamflow during drought conditions. Multiple linear regression analyses revealed that snowmelt-only models were better than the best precipitation and temperature models at explaining streamflow variability from all headwater subregions for both the wet and drought periods.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

Near-surface geophysics: Environmental applications

The field of geophysics encompasses a broad and diverse compilation of methodologies that employs principles of physics to characterize properties of earth materials within the subsurface. While geophysical methods have a long history in resource exploration and studies of Earth’s interior, the subdiscipline of “near-surface geophysics” has evolved in recent decades for examination of the shallow, near-surface environment for a range of purposes ranging from archaeological or forensic investigations to assessment of geologic, hydrologic, biologic, and geochemical properties and processes. “Environmental geophysics” are near-surface geophysical studies and methods that focus on understanding natural systems (e.g., watershed hydrology, groundwater–surface water connections, biophysical processes) as well as research pertaining to anthropogenic impacts and land management, (e.g., contamination and remediation, saltwater intrusion, agricultural practices). This field can be further subdivided into subdisciplines focused on specific topics and applications, such as water resources and hydrology (hydrogeophysics) or biologic and microbial processes (biogeophysics). Studies in environmental geophysics span a range of scales, from pore-scale laboratory tests to watershed-scale or regional field experiments. Methods vary by the nature of physics employed, the specific measurement acquired, and how that data is ultimately processed and analyzed to produce interpretable results. There exists further diversity in the acquisition logistics, geometry, and timing of data collection. Geophysical data can be collected in boreholes (one-dimensional, 1-D, vertical profiles), along survey lines (two-dimensional, 2-D, cross-sections), or in dense sensor arrays or gridded profiles (three-dimensional, 3-D, models). Regarding the temporal aspect, studies can conduct one-time geophysical surveys to obtain detailed imaging of subsurface structure or use timelapse and continuous monitoring to investigate variations in subsurface properties over time. The cumulation of all possible permutations of these factors (method, acquisition geometry, survey design, and target application) results in an immense diversity among environmental geophysical studies. Nevertheless, this field remains unified in the pursuit of understanding natural and human-impacted near-surface environments through geophysical investigations. Here we highlight some key references within environmental geophysics. Resources on geophysical theory, acquisition logistics, processing and inversion workflows, and example case studies are categorized into the most common geophysical classes within Geophysical Methods. Lastly, example references for the dominant types of applications in environmental geophysical studies are catalogued in Environmental Applications.

Book chapter↗

Geomorphic and vegetative recovery processes along modified stream channels of West Tennessee

Hundreds of miles of streams in West Tennessee have been channelized or otherwise modified since the turn of century. After all or parts of a stream are straightened, dredged, or cleared, systematic hydrologic, geomorphic, and ecologic processes collectively begin to reduce energy conditions towards the premodified state. One hundred and five sites along 15 streams were studied in the Obion, Forked Deer, Hatchie, and Wolf River basins. All studied streams, except the Hatchie River, have had major channel modification along all or parts of their courses. Bank material shear-strength properties were determined through drained borehole-shear testing (168 tests) and used to interpret present critical bank conditions and factors of safety, and to estimate future channel-bank stability. Mean values of cohesive strength and angle of internal friction were 1.26 pounds per square inch and 30.1 degrees, respectively. Dendrogeomorphic analyses were made using botanical evidence of channel-bank failures to determine rates of channel widening; buried riparian stems were analyzed to determine rates of bank accretion. Channel bed-level changes through time and space were represented by a power equation. Plant ecological analyses were made to infer relative bank stability, to identify indicator species of the stage of bank recovery, and to determine patterns of vegetation development through the course of channel evolution. Quantitative data on morphologic changes were used with previously developed six-stage models of channel evolution and bank-slope development to estimate trends of geomorphic and ecologic processes and forms through time. Immediately after channel modifications, a 10- to 1%yearperiod of channel-bed degradation ensues at and upstream from the most recent modifications (area of maximum disturbance). Channel-bed lowering by &gradation was as much as 20 feet along some stream reaches. Downstream from the area of maximum disturbance, the bed was aggraded by the deposition of sediment supplied by knickpoint migration upstream; aggradation also occurred in initially degraded sites with time. Additionally, if degradation caused an increase in bank height beyond the critical limits of the bank material, a period of channel widening by mass wasting followed. Degradation knickpoints migrated upstream at rates greater than 1 mile per year; the rates attenuated with distance above the area of maximum disturbance. Channel widening rates of up to 16 feet per year were documented along some severely degraded reaches. Planar failures were generally more frequent but rotational failures dominated the most rapidly widening reaches. Total volumes of bank erosion may represent 75percent or more of the total material eroded from the channel, but this material generally exits the drainage basin. Mean factors of safety vary with the stage of channel evolution with the lowest values for planar and rotational failures occurring during the threshold stage (stage IV) 1.00 and 1.15, respectively. As channel gradients decrease, degradation ceases and then a period of “secondary aggradation ” (at lesser rates than degradation) and bank accretion begins that may fill the channel to near floodplain level. This shift@ in process represents an oscillation in channel bed-level adjustment. Streams in basins underlain by loess may require an order of magnitude more time than sand-bed streams to stabilize due to a lack of coarse-grained material (sand) for aggradation. A systematic progression of riparian species that reflects the six-stage model of channel evolution has been identified. This progression can be used to infer ambient channel stability and hydrogeomorphic conditions. Woody vegetation establishes on low- and mid-bank surfaces (the slough line, initially) at about the same time that bank accretion begins. This slough line forms at a mean temporary stability angle of 24 degrees and expands upslope with time by the accretion of sediments. Species involved in this initial revegetation are hardy, fast growing, and can tolerate moderate amounts of slope instability and sedimentation; these species include river birch, black willow, boxelder, and silver maple. Vegetation appears to enhance bank stability, and with increasing stability, species such as bald cypress, tupelo gum, and various hydric oaks, which are more characteristic of stable, premodified riparian settings, begin to establish. Detrended-correspondence analysis indicated species assemblages associated with the six stages of channel evolution and bank-slope development. Ordination of site variables based on species data such as channel widening, bank accretion, and woody vegetative cover also reflects the temporal changes identifies by the models. Long-term channel geometry was estimated from a quantitative model of bed-level change, and from documented trends in channel widening. An idealized stable channel of a major sand-bed stream may have a width/depth ratio near 10 and bank slopes of about 24 degrees. This stable channel will ultimately undergo the development of point-bars and incipient meanders, characteristic of unmodified streams.

Tennessee↗

Wetland dynamics influence mid-continent duck recruitment

Recruitment is a key factor influencing duck population dynamics. Understanding what regulates recruitment of ducks is a prerequisite to informed habitat and harvest management. Quantity of May ponds (MP) has been linked to recruitment and population size (Kaminski and Gluesing 1987, Raveling and Heitmeyer 1989). However, wetland productivity (quality) is driven by inter-annual hydrological fluctuations. Periodic drying of wetlands due to wet-dry climate cycles releases nutrients and increases invertebrate populations when wet conditions return (Euliss et al. 1999). Wetlands may also become wet or dry within a breeding season. Accordingly, inter-annual and intra-seasonal hydrologic variation potentially influence duck recruitment. Here, we examined influences of wetland quantity, quality, and intra-seasonal dynamics on recruitment of ducks. We indexed duck recruitment by vulnerability-corrected age ratios (juveniles/adult females) for mid-continent Gadwall (Anas strepera). We chose Gadwall because the majority of the continental population breeds in the Prairie Pothole Region (PPR), where annual estimates of MP exist since 1974. We indexed wetland quality by calculating change in MP (?MP) over the past two years (?MP = 0.6[MPt – MPt-1] + 0.4[MPt – MPt-2]). We indexed intra-seasonal change in number of ponds by dividing the PPR mean standardized precipitation index for July by MP (hereafter summer index). MP and ?MP were positively correlated (r = 0.65); therefore, we calculated residual ?MP (?MPr) with a simple linear regression using MP, creating orthogonal variables. Finally, we conducted a multiple regression to examine how MP, ?MPr, and summer index explained variation in recruitment of Gadwall from 1976–2010. Our model explained 67% of the variation in mid-continent Gadwall recruitment and all three hydrologic indices were positively correlated with recruitment (Figure 1). Type II semi-partial R2 estimates indicated that MP accounted for 41%, ?MPr accounted for an additional 22%, and summer index accounted for the remaining 4% of the variation in recruitment. Our results are consistent with previous findings that quantity of MP was important for explaining variation in recruitment of ducks. However, our results also indicated that considering hydrologic dynamics was important for explaining recruitment. Additionally, the index for retention of MP within breeding year also was important, despite its coarse resolution as an average of precipitation events that can vary greatly spatially and in intensity within the PPR. Our results support the idea that wetland ecosystems in the PPR are ultimately regulated through bottom-up process driven by inter- and intra-annual hydrological dynamics. However from the ducks' perspective, hydrological dynamics could influence recruitment proximately through both bottom-up and top-down processes. Specifically, hydrological fluctuations may influence predator populations, prey switching by predators, or duckling vulnerability to predators (Cox et al. 1998). We will propose a conceptual model for understanding the potential role of bottom-up and top-down regulation of duck recruitment based on different hydrological contexts. Clearly, a better understanding of ultimate and proximate factors regulating duck recruitment would improve the effectiveness and efficiency of habitat conservation for ducks. Lastly, our findings could be used to improve models that predict fall flights for the purposes of informing harvest regulations.

Conference Paper↗

Implications of hydrologic variability on the succession of plants in Great Lakes wetlands

Primary succession of plant communities directed toward a climax is not a typical occurrence in wetlands because these ecological systems are inherently dependent on hydrology, and temporal hydrologic variability often causes reversals or setbacks in succession. Wetlands of the Great Lakes provide good examples for demonstrating the implications of hydrology in driving successional processes and for illustrating potential misinterpretations of apparent successional sequences. Most Great Lakes coastal wetlands follow cyclic patterns in which emergent communities are reduced in area or eliminated by high lake levels and then regenerated from the seed bank during low lake levels. Thus, succession never proceeds for long. Wetlands also develop in ridge and swale terrains in many large embayments of the Great Lakes. These formations contain sequences of wetlands of similar origin but different age that can be several thousand years old, with older wetlands always further from the lake. Analyses of plant communities across a sequence of wetlands at the south end of Lake Michigan showed an apparent successional pattern from submersed to floating to emergent plants as water depth decreased with wetland age. However, paleoecological analyses showed that the observed vegetation changes were driven largely by disturbances associated with increased human settlement in the area. Climate-induced hydrologic changes were also shown to have greater effects on plant-community change than autogenic processes. Other terms, such as zonation, maturation, fluctuations, continuum concept, functional guilds, centrifugal organization, pulse stability, and hump-back models provide additional means of describing organization and changes in vegetation; some of them overlap with succession in describing vegetation processes in Great Lakes wetlands, but each must be used in the proper context with regard to short- and long-term hydrologic variability.

Aquatic Ecosystem Health & Management↗

Accounting for time- and space-varying changes in the gravity field to improve the network adjustment of relative-gravity data

The relative gravimeter is the primary terrestrial instrument for measuring spatially and temporally varying gravitational fields. The background noise of the instrument—that is, non-linear drift and random tares—typically requires some form of least-squares network adjustment to integrate data collected during a campaign that may take several days to weeks. Here, we present an approach to remove the change in the observed relative-gravity differences caused by hydrologic or other transient processes during a single campaign, so that the adjusted gravity values can be referenced to a single epoch. The conceptual approach is an example of coupled hydrogeophysical inversion, by which a hydrologic model is used to inform and constrain the geophysical forward model. The hydrologic model simulates the spatial variation of the rate of change of gravity as either a linear function of distance from an infiltration source, or using a 3-D numerical groundwater model. The linear function can be included in and solved for as part of the network adjustment. Alternatively, the groundwater model is used to predict the change of gravity at each station through time, from which the accumulated gravity change is calculated and removed from the data prior to the network adjustment. Data from a field experiment conducted at an artificial-recharge facility are used to verify our approach. Maximum gravity change due to hydrology (observed using a superconducting gravimeter) during the relative-gravity field campaigns was up to 2.6 μGal d −1 , each campaign was between 4 and 6 d and one month elapsed between campaigns. The maximum absolute difference in the estimated gravity change between two campaigns, two months apart, using the standard network adjustment method and the new approach, was 5.5 μGal. The maximum gravity change between the same two campaigns was 148 μGal, and spatial variation in gravity change revealed zones of preferential infiltration and areas of relatively high groundwater storage. The accommodation for spatially varying gravity change would be most important for long-duration campaigns, campaigns with very rapid changes in gravity and (or) campaigns where especially precise observed relative-gravity differences are used in the network adjustment.

Geophysical Journal International↗

Controls on surface water chemistry in the upper Merced River basin, Yosemite National Park, California

Surface water draining granitic bedrock in Yosemite National Park exhibits considerable variability in chemical composition, despite the relative homogeneity of bedrock chemistry. Other geological factors, including the jointing and distribution of glacial till, appear to exert strong controls on water composition. Chemical data from three surface water surveys in the upper Merced River basin conducted in August 1981, June 1988 and August 1991 were analysed and compared with mapped geological, hydrological and topographic features to identify the solute sources and processes that control water chemistry within the basin during baseflow. Water at most of the sampling sites was dilute, with alkalinities ranging from 26 to 77 μequiv. l −1 . Alkalinity was much higher in two subcatchments, however, ranging from 51 to 302 μequiv. l −1 . Base cations and silica were also significantly higher in these two catchments than in the rest of the watershed. Concentrations of weathering products in surface water were correlated to the fraction of each subcatchment underlain by surficial material, which is mostly glacial till. Silicate mineral weathering is the dominant control on concentrations of alkalinity, silica and base cations, and ratios of these constituents in surface water reflect the composition of local bedrock. Chloride concentrations in surface water samples varied widely, ranging from <1 to 96 μequiv. l −1 . The annual volume-weighted mean chloride concentration in the Merced River at the Happy Isles gauge from 1968 to 1990 was 26 μequiv. l −1 , which was five times higher than in atmospheric deposition (4–5 μequiv. l −1 ), suggesting that a source of chloride exists within the watershed. Saline groundwater springs, whose locations are probably controlled by vertical jointing in the bedrock, are the most likely source of the chloride. Sulphate concentrations varied much less than most other solutes, ranging from 3 to 14 μequiv. l −1 . Concentrations of sulphate in quarterly samples collected at the watershed outlet also showed relatively little variation, suggesting that sulphate may be regulated to some extent by a within-watershed process, such as sulphate adsorption.

Hydrological Processes↗

Approaches to modeling weathered regolith

Sustainable soils are a requirement for maintaining human civilizations ( Carter and Dale 1974 ; Lal 1989 ). However, as the “most complicated biomaterial on the planet” ( Young and Crawford 2004 ), soils represent one of the most difficult systems to understand and model with respect to chemical, physical, and biological coupling over time (Fig. 1 ). Despite the complexity of these interactions, certain patterns in soil properties and development are universally observed and have been used in soil science as a means for classification. Elemental, mineralogical, or isotopic concentrations in soils plotted versus depth beneath the land surface comprise such patterns. Soil depth profiles are often reported for solid soil materials, and, less frequently, for solutes in soil pore waters. These profiles cross a large range in spatial scales that traditionally have been studied by different disciplines. For example, shallow, biologically active horizons are commonly defined as the soil zone in agronomic studies whereas the mobile layer of the regolith is referred to as soil in geomorphological studies. In contrast, many geochemical studies target chemical weathering to tens or even hundreds of meters in depth, sometimes extending the definition of “soils” to include the entire regolith down to parent bedrock or alluvium. Soil profiles also exhibit a large range in temporal scales ( Amundson 2004 ; Brantley 2008b ). Solid-state profiles document chemical and mineralogical changes integrated over the time scales of evolution of regolith from protolith. This “geologic time” can vary from tens to hundreds of years for weathered material developed on moraines deposited by active glaciers ( Anderson et al. 1997 ), to millions or possibly hundreds of millions of years of regolith evolution as documented in laterites and bauxites on stable cratons ( Nahon 1986 ). In contrast, solute profiles reflect much shorter time scales corresponding to the residence time of the soil water which commonly ranges from days to decades ( Stonestrom et al. 1998 ). Factors impacting soil minerals can therefore be related to geologically old processes while those impacting pore waters are related to contemporary processes. We first discuss a geochemical frame work for modeling soil profiles, including a simple scheme that depends on the extent of enrichment or depletion. Such profiles are comprised of reaction fronts affected by chemical, hydrologic, geologic and biologic processes that control soil evolution. We then present a hierarchy of models that have been used to interpret both solid state and solute compositions in regolith. The more simple approaches to model depletion in soils, using analytical models, are first described. The most elementary of these is a linear model that calculates rate constants from the slopes of either solid or solute weathering gradients: these rate constants represent lumped parameters that describe weathering in terms of an integrated reaction rate. Two other analytical models are then presented that have been used to fit solid state elemental profiles with exponential and sigmoidal functions. All of these analytical approaches are derived for models of soils as containing a limited number of components, phases, and species. At a more complex level, numerical models are then presented to elucidate how kinetic and transport parameters as well as chemical, hydrologic, and physical soil data can be incorporated. We consider two forms of these models, first relatively simple spreadsheet calculators and then more sophisticated multi-component, multi-phase reactive-transport numerical codes. Our treatment of reactive transport modeling is relatively cursory, in recognition of the treatment in the chapter by Steefel and Maher (2009 , this volume). Because these models incorporate more phases, components, and species than the other approaches and explicitly model the more fundamental reaction mechanisms involved, they generally have a greater need for parameterization. In our conclusion section, we discuss how this hierarchy of approaches can yield generalizations about soils that are often complementary.

Reviews in Mineralogy and Geochemistry↗

Seasonal and event-scale variations in solute chemistry for four Sierra Nevada catchments

Hydrobiogeochemical processes controlling stream water chemistry were examined in four small (<5 km2) catchments having contrasting bedrock lithologies in the western Sierra Nevada foothills of California. The Mediterranean climate with its cool/wet and hot/dry cycle produces strong seasonal patterns in hydrological, biological and geochemical processes. Stream water solutes fall into three general groups according to seasonal fluctuation in concentration: Strong, rainy season minimum-dry season maximum (Cl-, SO42-, base cations); weak, rainy season minimum-dry season maximum (Si); and rainy season maximum-dry season minimum (NO3- and K+). Solute dynamics in soil solutions and stream water suggest that mixing of drainage waters from bedrock and soil sources regulate stream water solute concentrations. Patterns are further altered by the leaching of solutes accumulated in the soil over the summer period of desiccation and the temporal discoupling of nutrient cycles that occurs due to differences in the timing between vegetation growth (late spring) and leaching (early winter). Solute concentrations are remarkably similar between watersheds with varying bedrock types, with the exception of nitrate, sulfate and bicarbonate. Three watersheds have nitrogen-bearing metasedimentary bedrock that contributes to elevated nitrate concentrations in stream waters. Watersheds whose bedrock includes mineralized veins of sulfide and carbonate minerals similarly have greater sulfate and bicarbonate concentrations in stream water. Hydrobiogeochemical processes are highly dynamic at the seasonal and storm-event temporal scales and spatially complex at the watershed scale making management of stream water chemical composition, such as nitrate concentrations, very challenging. ?? 2001 Elsevier Science B.V.

Journal of Hydrology↗

Expanding the role of reactive transport models in critical zone processes

Models test our understanding of processes and can reach beyond the spatial and temporal scales of measurements. Multi-component Reactive Transport Models (RTMs), initially developed more than three decades ago, have been used extensively to explore the interactions of geothermal, hydrologic, geochemical, and geobiological processes in subsurface systems. Driven by extensive data sets now available from intensive measurement efforts, there is a pressing need to couple RTMs with other community models to explore non-linear interactions among the atmosphere, hydrosphere, biosphere, and geosphere. Here we briefly review the history of RTM development, summarize the current state of RTM approaches, and identify new research directions, opportunities, and infrastructure needs to broaden the use of RTMs. In particular, we envision the expanded use of RTMs in advancing process understanding in the Critical Zone, the veneer of the Earth that extends from the top of vegetation to the bottom of groundwater. We argue that, although parsimonious models are essential at larger scales, process-based models offer tools to explore the highly nonlinear coupling that characterizes natural systems. We present seven testable hypotheses that emphasize the unique capabilities of process-based RTMs for (1) elucidating chemical weathering and its physical and biogeochemical drivers; (2) understanding the interactions among roots, micro-organisms, carbon, water, and minerals in the rhizosphere; (3) assessing the effects of heterogeneity across spatial and temporal scales; and (4) integrating the vast quantity of novel data, including “omics” data (genomics, transcriptomics, proteomics, metabolomics), elemental concentration and speciation data, and isotope data into our understanding of complex earth surface systems. With strong support from data-driven sciences, we are now in an exciting era where integration of RTM framework into other community models will facilitate process understanding across disciplines and across scales.

Earth-Science Reviews↗

Assessing impacts of roads: application of a standard assessment protocol

Adaptive management of road networks depends on timely data that accurately reflect the impacts those systems are having on ecosystem processes and associated services. In the absence of reliable data, land managers are left with little more than observations and perceptions to support management decisions of road-associated disturbances. Roads can negatively impact the soil, hydrologic, plant, and animal processes on which virtually all ecosystem services depend. The Interpreting Indicators of Rangeland Health (IIRH) protocol is a qualitative method that has been demonstrated to be effective in characterizing impacts of roads. The goal of this study were to develop, describe, and test an approach for using IIRH to systematically evaluate road impacts across large, diverse arid and semiarid landscapes. We developed a stratified random sampling approach to plot selection based on ecological potential, road inventory data, and image interpretation of road impacts. The test application on a semiarid landscape in southern New Mexico, United States, demonstrates that the approach developed is sensitive to road impacts across a broad range of ecological sites but that not all the types of stratification were useful. Ecological site and road inventory strata accounted for significant variability in the functioning of ecological processes but stratification based on apparent impact did not. Analysis of the repeatability of IIRH applied to road plots indicates that the method is repeatable but consensus evaluations based on multiple observers should be used to minimize risk of bias. Landscape-scale analysis of impacts by roads of contrasting designs (maintained dirt or gravel roads vs. non- or infrequently maintained roads) suggests that future travel management plans for the study area should consider concentrating traffic on fewer roads that are well designed and maintained. Application of the approach by land managers will likely provide important insights into minimizing impacts of road networks on key ecosystem services.

Rangeland Ecology and Management↗

Preliminary surficial geologic map database of the Amboy 30 x 60 minute quadrangle, California

The surficial geologic map database of the Amboy 30x60 minute quadrangle presents characteristics of surficial materials for an area approximately 5,000 km2 in the eastern Mojave Desert of California. This map consists of new surficial mapping conducted between 2000 and 2005, as well as compilations of previous surficial mapping. Surficial geology units are mapped and described based on depositional process and age categories that reflect the mode of deposition, pedogenic effects occurring post-deposition, and, where appropriate, the lithologic nature of the material. The physical properties recorded in the database focus on those that drive hydrologic, biologic, and physical processes such as particle size distribution (PSD) and bulk density. This version of the database is distributed with point data representing locations of samples for both laboratory determined physical properties and semi-quantitative field-based information. Future publications will include the field and laboratory data as well as maps of distributed physical properties across the landscape tied to physical process models where appropriate. The database is distributed in three parts: documentation, spatial map-based data, and printable map graphics of the database. Documentation includes this file, which provides a discussion of the surficial geology and describes the format and content of the map data, a database 'readme' file, which describes the database contents, and FGDC metadata for the spatial map information. Spatial data are distributed as Arc/Info coverage in ESRI interchange (e00) format, or as tabular data in the form of DBF3-file (.DBF) file formats. Map graphics files are distributed as Postscript and Adobe Portable Document Format (PDF) files, and are appropriate for representing a view of the spatial database at the mapped scale.

Open-File Report↗

Long-term demographic analysis of the Cape Sable seaside sparrow (1992–2021)

The Cape Sable seaside sparrow (Ammospiza maritima mirabilis) is an endangered species that has experienced a population decline of more than 60% since 1981. Despite its critical population status, a statistically robust analysis of the species’ demographic rates utilizing all data has yet to be completed (Benscoter et al. 2021). Furthermore, long-term population processes in response to hydrologic and environmental conditions have not been evaluated for this species. To address these substantial gaps in knowledge, 30 years of demographic data were synthesized to assess population dynamics of this imperiled species using an integrated population model (IPM). Three demographic data types (range-wide counts, capture-mark-recapture, nest monitoring) were incorporated into a unified IPM to evaluate demographic processes and predict population trajectories. The following were calculated: (1) annual estimates of annual population size, survival, and fecundity; (2) coefficient estimates of hydrologic and environmental variables on survival and fecundity; (3) estimates of annual population growth and their correlation with demographic rates; and (4) a Bayesian population viability analysis (BPVA) that includes predicted population size, demographic rates, and extinction risks ten years into the future.

Report↗

Preliminary analysis of regional-precipitation periodicity

Precipitation variability plays a major role in nearly every aspect of the hydrologic cycle. Precipitation is not a random event, but it occurs after a sequence of prerequi-sites has been fulfilled. Recent investigations have shown that activity of the sun can affect atmospheric vorticity, an important factor in precipitation formation. Solar activity is known to be periodic; therefore, through a complex series of physical processes, precipitation variance is solar forced to a certain degree. A preliminary analysis of precipitation periodicity was made for eight regions scattered across the central United States. Each region contained 5 to 10 stations with long-term precipitation records that were averaged to obtain yearly regional-precipitation values. Graphic analysis shows 11-year and 22-year cycles that are nearly in phase with the solar cycles. An example of the effect of cyclic precipitation is presented for the Powder River basin in Wyoming and Montana. A cycle of 22 years exhibits fluctuations of approximately 22 to 27 percent for precipitation and 38 to 50 percent for runoff. A more detailed study that investigates solar-forced precipitation cycles and their relationship to hydrologic processes is needed.

Water-Resources Investigations Report↗

Watershed Data Management (WDM) database for Salt Creek streamflow simulation, DuPage County, Illinois, water years 2005-11

The U.S. Geological Survey (USGS), in cooperation with DuPage County Stormwater Management Division, maintains a USGS database of hourly meteorologic and hydrologic data for use in a near real-time streamflow simulation system, which assists in the management and operation of reservoirs and other flood-control structures in the Salt Creek watershed in DuPage County, Illinois. Most of the precipitation data are collected from a tipping-bucket rain-gage network located in and near DuPage County. The other meteorologic data (wind speed, solar radiation, air temperature, and dewpoint temperature) are collected at Argonne National Laboratory in Argonne, Ill. Potential evapotranspiration is computed from the meteorologic data. The hydrologic data (discharge and stage) are collected at USGS streamflow-gaging stations in DuPage County. These data are stored in a Watershed Data Management (WDM) database. An earlier report describes in detail the WDM database development including the processing of data from January 1, 1997, through September 30, 2004, in SEP04.WDM database. SEP04.WDM is updated with the appended data from October 1, 2004, through September 30, 2011, water years 2005–11 and renamed as SEP11.WDM. This report details the processing of meteorologic and hydrologic data in SEP11.WDM. This report provides a record of snow affected periods and the data used to fill missing-record periods for each precipitation site during water years 2005–11. The meteorologic data filling methods are described in detail in Over and others (2010), and an update is provided in this report.

Illinois↗

PHT3D-UZF: A reactive transport model for variably-saturated porous media

A modified version of the MODFLOW/MT3DMS-based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably-saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D-UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated-zone flow (UZF1) package. A volume-averaged approach similar to the method used in UZF-MT3DMS was adopted. The PHREEQC-based computation of chemical processes within PHT3D-UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional-scale applications, UZF1 simulates downward-only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably-saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns.

Ground Water↗

Drainage network structure and hydrologic behavior of three lake-rich watersheds on the Arctic Coastal Plain, Alaska

Watersheds draining the Arctic Coastal Plain (ACP) of Alaska are dominated by permafrost and snowmelt runoff that create abundant surface storage in the form of lakes, wetlands, and beaded streams. These surface water elements compose complex drainage networks that affect aquatic ecosystem connectivity and hydrologic behavior. The 4676 km 2 Fish Creek drainage basin is composed of three watersheds that represent a gradient of the ACP landscape with varying extents of eolian, lacustrine, and fluvial landforms. In each watershed, we analyzed 2.5-m-resolution aerial photography, a 5-m digital elevation model, and river gauging and climate records to better understand ACP watershed structure and processes. We show that connected lakes accounted for 19 to 26% of drainage density among watersheds and most all channels initiate from lake basins in the form of beaded streams. Of the > 2500 lakes in these watersheds, 33% have perennial streamflow connectivity, and these represent 66% of total lake area extent. Deeper lakes with over-wintering habitat were more abundant in the watershed with eolian sand deposits, while the watershed with marine silt deposits contained a greater extent of beaded streams and shallow thermokarst lakes that provide essential summer feeding habitat. Comparison of flow regimes among watersheds showed that higher lake extent and lower drained lake-basin extent corresponded with lower snowmelt and higher baseflow runoff. Variation in baseflow runoff among watersheds was most pronounced during drought conditions in 2007 with corresponding reduction in snowmelt peak flows the following year. Comparison with other Arctic watersheds indicates that lake area extent corresponds to slower recession of both snowmelt and baseflow runoff. These analyses help refine our understanding of how Arctic watersheds are structured and function hydrologically, emphasizing the important role of lake basins and suggesting how future lake change may impact hydrologic processes.

Alaska↗