Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Hydrology”

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

Automated calibration of a stream solute transport model: Implications for interpretation of biogeochemical parameters

The hydrologic processes of advection, dispersion, and transient storage are the primary physical mechanisms affecting solute transport in streams. The estimation of parameters for a conservative solute transport model is an essential step to characterize transient storage and other physical features that cannot be directly measured, and often is a preliminary step in the study of reactive solutes. Our study used inverse modeling to estimate parameters of the transient storage model OTIS (One dimensional Transport with Inflow and Storage). Observations from a tracer injection experiment performed on Uvas Creek, California, USA, are used to illustrate the application of automated solute transport model calibration to conservative and nonconservative stream solute transport. A computer code for universal inverse modeling (UCODE) is used for the calibrations. Results of this procedure are compared with a previous study that used a trial-and-error parameter estimation approach. The results demonstrated 1) importance of the proper estimation of discharge and lateral inflow within the stream system; 2) that although the fit of the observations is not much better when transient storage is invoked, a more randomly distributed set of residuals resulted (suggesting non-systematic error), indicating that transient storage is occurring; 3) that inclusion of transient storage for a reactive solute (Sr2+) provided a better fit to the observations, highlighting the importance of robust model parameterization; and 4) that applying an automated calibration inverse modeling estimation approach resulted in a comprehensive understanding of the model results and the limitation of input data.

Journal of the North American Benthological Societ↗

Cumulative spring discharge and survey effort influence occupancy and detection of a threatened freshwater mussel, the Suwannee Moccasinshell

Freshwater mussels (Unionidae) are among the most imperiled groups of organisms in the world, and the lack of information regarding species distributions, life-history characteristics, and ecological and biological requirements may limit the protection of remaining mussel populations. We examined the influence of hydrologic factors on the occurrence of the Suwannee Moccasinshell Medionidus walkeri, a federally threatened freshwater mussel species, endemic to the Suwannee River Basin in Georgia and Florida. We also evaluated the influence of survey effort on detection of Suwannee Moccasinshell during field surveys. We compiled all recent (2013–2016) mussel survey records in the Suwannee River Basin. We calculated cumulative discharge contributed by upstream springs for each of 220 survey locations. We combined the spring discharge predictor variable with Suwannee Moccasinshell detection and nondetection data from each survey location to develop a suite of occupancy models. Modeling results indicated that detection of Suwannee Moccasinshell during surveys was strongly and positively related to survey effort. Modeling results also indicated that sites with cumulative spring discharge inputs exceeding ∼28 cubic meters per second were most likely (i.e., predicted occupancy probabilities >0.5) to support Suwannee Moccasinshell populations. However, occupancy declined in the lowermost reaches of the Suwannee mainstem despite high spring discharge inputs, presumably due to greater tidal influences and differences in physicochemical habitat conditions. Historical localities where Suwannee Moccasinshell has presumably been extirpated are all devoid of springs in their upstream watersheds. We hypothesize that springs may buffer extremely tannic, and at times polluted, surface waters, in addition to maintaining adequate flows during periods of drought, thereby promoting the persistence of Suwannee Moccasinshell populations. Our study suggests that springs are a critical resource for Suwannee Moccasinshell and may be more important for conservation planning than was previously recognized.

Florida↗

Abundance, stable isotopic composition, and export fluxes of DOC, POC, and DIC from the Lower Mississippi River during 2006–2008

Sources, abundance, isotopic compositions, and export fluxes of dissolved inorganic carbon (DIC), dissolved and colloidal organic carbon (DOC and COC), and particulate organic carbon (POC), and their response to hydrologic regimes were examined through monthly sampling from the Lower Mississippi River during 2006–2008. DIC was the most abundant carbon species, followed by POC and DOC. Concentration and δ 13 C of DIC decreased with increasing river discharge, while those of DOC remained fairly stable. COC comprised 61 ± 3% of the bulk DOC with similar δ 13 C abundances but higher percentages of hydrophobic organic acids than DOC, suggesting its aromatic and diagenetically younger status. POC showed peak concentrations during medium flooding events and at the rising limb of large flooding events. While δ 13 C-POC increased, δ 15 N of particulate nitrogen decreased with increasing discharge. Overall, the differences in δ 13 C between DOC or DIC and POC show an inverse correlation with river discharge. The higher input of soil organic matter and respired CO 2 during wet seasons was likely the main driver for the convergence of δ 13 C between DIC and DOC or POC, whereas enhanced in situ primary production and respiration during dry seasons might be responsible for their isotopic divergence. Carbon export fluxes from the Mississippi River were estimated to be 13.6 Tg C yr −1 for DIC, 1.88 Tg C yr −1 for DOC, and 2.30 Tg C yr −1 for POC during 2006–2008. The discharge-normalized DIC yield decreased during wet seasons, while those of POC and DOC increased and remained constant, respectively, implying variable responses in carbon export to the increasing discharge.

Louisiana↗

Flow–recruitment relationships for Shoal Chub and implications for managing environmental flows

Objective Regulation of river flow regimes by dams and diversions impacts aquatic biota and ecosystems globally. However, our understanding of the ecological consequences of flow alteration and ecological benefits of flow restoration lags behind our ability to manipulate flows, and there is a need for broader development of flow–ecology relationships. Approaches for establishing flow–ecology relationships have recently shifted away from state-based methods that analyze snapshots of ecological conditions and towards rate-based methods focused on mechanisms that link hydrology with dynamics of important ecological components and processes. Methods We used a rate-based approach to validate environmental flow standards developed for the lower Brazos River, Texas, by analyzing the relationship between flow regime components and recruitment strength of imperiled Shoal Chub Macrhybopsis hyostoma , a fluvial specialist and pelagic-broadcast-spawning fish. We collected 254 age-0 Shoal Chub (9–40 mm total length), extracted their otoliths to estimate age in days, and used a generalized additive model to regress the number of captured recruits that hatched on a calendar date against flow regime metrics, such as pulse magnitude, flow rate of change, and pulse timing in relation to environmental flow standards proposed by a science advisory committee (Brazos Basin and Bay Area Expert Science Team). Result The model revealed that flow magnitude, rate of change, and timing were all significant predictors that collectively explained 60% of variation in the recruitment strength index. Hindcasting for 1919–2020 showed a general reduction in recruitment strength following commencement of flow regulation in the lower Brazos River and revealed that high recruitment correlated with years in which most or all proposed flow tiers were attained, whereas low recruitment correlated with years when less than half of the targeted tiers were attained. Conclusion Our work represents an effective validation method for environmental flow recommendations and reveals specific flow regimes that benefit an imperiled fish species.

New Mexico, Texas↗

Representing the effects of alpine grassland vegetation cover on the simulation of soil thermal dynamics by ecosystem models applied to the Qinghai-Tibetan Plateau

Soil surface temperature is a critical boundary condition for the simulation of soil temperature by environmental models. It is influenced by atmospheric and soil conditions and by vegetation cover. In sophisticated land surface models, it is simulated iteratively by solving surface energy budget equations. In ecosystem, permafrost, and hydrology models, the consideration of soil surface temperature is generally simple. In this study, we developed a methodology for representing the effects of vegetation cover and atmospheric factors on the estimation of soil surface temperature for alpine grassland ecosystems on the Qinghai-Tibetan Plateau. Our approach integrated measurements from meteorological stations with simulations from a sophisticated land surface model to develop an equation set for estimating soil surface temperature. After implementing this equation set into an ecosystem model and evaluating the performance of the ecosystem model in simulating soil temperature at different depths in the soil profile, we applied the model to simulate interactions among vegetation cover, freeze-thaw cycles, and soil erosion to demonstrate potential applications made possible through the implementation of the methodology developed in this study. Results showed that (1) to properly estimate daily soil surface temperature, algorithms should use air temperature, downward solar radiation, and vegetation cover as independent variables; (2) the equation set developed in this study performed better than soil surface temperature algorithms used in other models; and (3) the ecosystem model performed well in simulating soil temperature throughout the soil profile using the equation set developed in this study. Our application of the model indicates that the representation in ecosystem models of the effects of vegetation cover on the simulation of soil thermal dynamics has the potential to substantially improve our understanding of the vulnerability of alpine grassland ecosystems to changes in climate and grazing regimes.

Qinghai;Tibet↗

Hydrologic effects on water level changes associated with episodic fault creep near Parkfield, California

As part of the Parkfield, California, earthquake prediction experiment, water level is monitored in a well 460 m from the main trace of the San Andreas fault on Middle Mountain, in the preparation zone of the anticipated Parkfield earthquake. The well configuration allows water level to be monitored in two fluid reservoirs at depths of 85 and 250 m below land surface. During 1987, water level changes were recorded during 12 of the 18 episodes of accelerated fault creep detected by a creep meter spanning the fault trace 750 m northwest of the well. The creep-related water level changes in the shallow reservoir have durations of less than 1 day, whereas in the deeper reservoir the changes persist for as long as 2 months. These data suggest that the transient nature of the water level changes in the shallow interval is due to vertical flow to the water table and is not evidence that creep events propagate past the well. Phase leads of earth tidal constituents in the water level data from the shallow interval relative to the same constituents in the local volume strain tide support the interpretation of significant flow to the water table at periods of 1 day or less. The form of the water level changes in the deep interval is affected by horizontal flow to the well bore. This effect can be removed from the water level records using a theoretical response curve constrained by the phases of earth tidal constituents in the deep interval relative to the local volume strain tide. For the events where the signal in the shallow interval has been large enough to measure, the sizes of the simultaneous water level changes in the two reservoirs are consistent with the same amounts of volume strain occurring at both depths.

Journal of Geophysical Research Solid Earth↗

Hydrologic characterization of desert soils with varying degrees of pedogenesis: 1. field experiments evaluating plant-relevant soil water behavior

To assess the eff ect of pedogenesis on the soil moisture dynamics infl uencing the character and quality of ecological habitat, we conducted infi ltration and redistribution experiments on three alluvial deposits in the Mojave National Preserve: (i) recently deposited active wash sediments, (ii) a soil of early Holocene age, and (iii) a highly developed soil of late Pleistocene age. At each, we ponded water in a 1-m-diameter infi ltration ring for 2.3 h and monitored soil water content and matric pressure during and atier infi ltration, using probes and electrical resistivity imaging (ERI). Infi ltration and downward fl ow rates were greater in younger material, favoring deep-rooted species. Deep-rooted species tend to colonize the margins of washes, where they are unaff ected by sediment transport that inhibits colonization. The ERI results support important generalizations, for example that shallower than 0.5 m, infi ltrated water persists longer in highly developed soil, favoring shallow-rooted species. Soil moisture data for the two youngest soils suggested that saturation overshoot, which may have signifi cant but unexplored hydroecologic and pedogenic eff ects, occurred at the horizontally advancing weting front. Spatial heterogeneity of soil properties generally increased with pedogenic development. Evidence suggested that some early-stage developmental processes may promote uniformity; the intermediate- age soil appeared to have the least heterogeneity in terms of textural variation with depth, and also the least anisotropy. Lateral heterogeneity was pronounced in older soil, having a multitude of eff ects on the distribution and retention of soil water, and may facilitate certain water-conserving strategies of plants over what would be possible in a laterally homogeneous soil. ?? Soil Science Society of America.

California↗

Residence time controls on the fate of nitrogen in flow‐through lakebed sediments

For many glacial lakes with highly permeable sediments, water exchange rates control hydrologic residence times within the sediment‐water interface (SWI) and the removal of reactive compounds such as nitrate, a common pollutant in lakes and groundwater. Here we conducted a series of focused tracer injection experiments in the upper 20 cm of the naturally downwelling SWI in a flow‐through lake on Cape Cod, MA. We systematically varied residence time and reactant controls on nitrate processing, using isotopically labeled 15 N nitrate to monitor the effect of these changes on nitrate removal via denitrification. The addition of acetate, a labile carbon compound, triggered the lake SWI to switch from net production to net removal of nitrate. When acetate was combined with increased residence time created by controlled reductions in water flux, we observed a fivefold increase in nitrate removal, a 26‐fold increase in N 2 production, and a 42‐fold increase in N 2 O production. We demonstrate that water residence time is an important control on the fate of nitrate in these lake SWIs and illustrate that seasonal conditions that alter lake exchange rates and variability in lake carbon may predict dynamic nitrate removal across the SWI. Additionally, observed N 2 O production during the oxic pore water experiments paired with geophysical characterization of the sediment porosity revealed that the lake SWI has less mobile pores occupying upward of 50% of the total porosity volume, which function as reactive microzones for nitrate processing.

Journal of Geophysical Research: Biogeosciences↗

Challenges for creating a site-specific groundwater-use record for the Ozark Plateaus aquifer system (central USA) from 1900 to 2010

Hydrologic budgets to determine groundwater availability are important tools for water-resource managers. One challenging component for developing hydrologic budgets is quantifying water use through time because historical and site-specific water-use data can be sparse or poorly documented. This research developed a groundwater-use record for the Ozark Plateaus aquifer system (central USA) from 1900 to 2010 that related county-level aggregated water-use data to site-specific well locations and aquifer units. A simple population-based linear model, constrained to 0 million liters per day in 1900, provided the best means to extrapolate groundwater-withdrawal rates pre-1950s when there was a paucity of water-use data. To disaggregate county-level data to individual wells across a regional aquifer system, a programmatic hierarchical process was developed, based on the level of confidence that a well pumped groundwater for a specific use during a specific year. Statistical models tested on a subset of the best-available site-specific water-use data provided a mechanism to bracket historic groundwater use, such that groundwater-withdrawal rates ranged, on average, plus or minus 38% from modeled values. Groundwater withdrawn for public supply and domestic use accounted for between 48 and 74% of total groundwater use since 1901, highlighting that groundwater provides an important drinking-water resource. The compilation, analysis, and spatial and temporal extrapolation of water-use data remain a challenging task for water scientists, but is of paramount importance to better quantify groundwater use and availability.

Arkansas, Kansas, Missouri, Oklahoma↗

Uncertainty in critical source area predictions from watershed-scale hydrologic models

Watershed-scale hydrologic models are frequently used to inform conservation and restoration efforts by identifying critical source areas (CSAs; alternatively 'hotspots'), defined as areas that export relatively greater quantities of nutrients and sediment. The CSAs can then be prioritized or ‘targeted’ for conservation and restoration to ensure efficient use of limited resources. However, CSA simulations from watershed-scale hydrologic models may be uncertain and it is critical that the extent and implications of this uncertainty be conveyed to stakeholders and decision makers. We used an ensemble of four independently developed Soil and Water Assessment Tool (SWAT) models and a SPAtially Referenced Regression On Watershed attributes (SPARROW) model to simulate CSA locations for flow, phosphorus, nitrogen, and sediment within the ~17,000-km 2 Maumee River watershed at the HUC-12 scale. We then assessed uncertainty in CSA simulations determined as the variation in CSA locations across the models. Our application of an ensemble of models - differing with respect to inputs, structure, and parameterization - facilitated an improved accounting of CSA prediction uncertainty. We found that the models agreed on the location of a subset of CSAs, and that these locations may be targeted with relative confidence. However, models more often disagreed on CSA locations. On average, only 16%–46% of HUC-12 subwatersheds simulated as a CSA by one model were also simulated as a CSA by a different model. Our work shows that simulated CSA locations are highly uncertain and may vary substantially across models. Hence, while models may be useful in informing conservation and restoration planning, their application to identify CSA locations would benefit from comprehensive uncertainty analyses to avoid inefficient use of limited resources.

Ohio↗

Short-term variability of 7 Be atmospheric deposition and watershed response in a Pacific coastal stream, Monterey Bay, California, USA

Beryllium-7 is a powerful and commonly used tracer for environmental processes such as watershed sediment provenance, soil erosion, fluvial and nearshore sediment cycling, and atmospheric fallout. However, few studies have quantified temporal or spatial variability of 7 Be accumulation from atmospheric fallout, and parameters that would better define the uses and limitations of this geochemical tracer. We investigated the abundance and variability of 7 Be in atmospheric deposition in both rain events and dry periods, and in stream surface-water samples collected over a ten-month interval at sites near northern Monterey Bay (37°N, 122°W) on the central California coast, a region characterized by a rainy winters, dry summers, and small mountainous streams with flashy hydrology. The range of 7 Be activity in rainwater samples from the main sampling site was 1.3–4.4 Bq L −1 , with a mean (±standard deviation) of 2.2 ± 0.9 Bq L −1 , and a volume-weighted average of 2.0 Bq L −1 . The range of wet atmospheric deposition was 18–188 Bq m −2 per rain event, with a mean of 72 ± 53 Bq m −2 . Dry deposition fluxes of 7 Be ranged from less than 0.01 up to 0.45 Bq m −2 d −1 , with an estimated dry season deposition of 7 Bq m −2 month −1 . Annualized 7 Be atmospheric deposition was approximately 1900 Bq m −2 yr −1 , with most deposition via rainwater (>95%) and little via dry deposition. Overall, these activities and deposition fluxes are similar to values found in other coastal locations with comparable latitude and Mediterranean-type climate. Particulate 7 Be values in the surface water of the San Lorenzo River in Santa Cruz, California, ranged from <0.01 Bq g −1 to 0.6 Bq g −1 , with a median activity of 0.26 Bq g −1 . A large storm event in January 2010 characterized by prolonged flooding resulted in the entrainment of 7 Be-depleted sediment, presumably from substantial erosion in the watershed. There were too few particulate 7 Be data over the storm to accurately model a 7 Be load, but the results suggest enhanced watershed export of 7 Be from small, mountainous river systems compared to other watershed types.

California↗

Early vegetational changes on a forested wetland constructed for mitigation

Changes in vegetation were studied on 15 acres of a 35 acre forested wetland created as a mitigation site in Anne Arundel County, Maryland during 1994-96. Meter-square sampling on four different hydrologic elevations determined that grasses initially dominated the area, but decreased from 59 percent in 1994 to 51 percent in 1995 and 30 percent in 1996. Herbaceous non-grass plants (forbs) increased from 19 percent to 56 percent in the three-year period. Area with no plant cover decreased from 21 percent in 1994 to 11 percent in 1995, and 10 percent in 1996. Woody plants comprised 2 percent of the cover in 1994, increased to 4 percent in 1995, and remained at 4 percent in 1996. The increase of woody plants was mainly from natural regeneration (pioneer) plants. Monitoring of the transplanted trees and shrubs indicated 35 percent mortality and little growth of surviving plants. The pioneer woody plant forming most of the cover was black willow (Salix nigra). Differences in the vegetation were observed among the four elevations, although no differences were observed for the major vegetation classes between plots that were planted and those that were not planted with woody plants. Dominant grass species was redtop (Agrostis stolonifera), which comprised 51 percent of the cover in 1994 and 42 percent cover in 1995 and 23 percent in 1996. Other species that were common were bush clover (Lespedeza cuneata), Japanese clover (Lespedeza striata) and flat pea (Lathyrus sylvestris). All four of these dominant species were part of the original seed mixtures that were seeded on the site. A total of 134 species of plants was recorded on the site indicating a fairly diverse community for a newly established habitat.

Wetland Journal↗

Reservoir analysis of the Denver earthquakes: A case of induced seismicity

Injection of fluid wastes into the fractured Precambrian crystalline bedrock beneath the Rocky Mountain Arsenal near Denver triggered earthquakes in the 1960's. An analysis, based on the assumption that fluid flow in the fractured reservoir can be approximated by flow in a porous medium, is presented. The configuration and hydrologic properties of the reservoir are determined from two lines of evidence: (1) locations of earthquake hypocenters determined by seismic arrays installed at the Arsenal and (2) observed long-term decline in fluid levels in the injection well. Together these two sets of data indicate that a long, narrow reservoir, aligned in the direction N 60°W, exists. The reservoir is 3.35 km in width, extends 30.5 km to the northwest and infinitely to the southeast, and spans a depth interval from 3.7 to 7.0 km below land surface. It has a transmissivity of 1.08×10-5 m2/s and a storage coefficient of 1.0×10-5. Computed pressure buildup along the length of the reservoir is compared with the spatial distribution of earthquake epicenters. The comparison shows that earthquakes are confined to that part of the reservoir where the pressure buildup exceeds 32 bars. This critical value is interpreted as the pressure buildup above which earthquakes occur. The migration of earthquake epicenters away from the injection well, a phenomenon noted by previous investigators, can be accounted for by the outward propagation of the critical pressure buildup. The analysis is extended to examining the effects of rapid flow in fractures opened by high injection pressure. The results show that the effect is confined to a small region within 1 km of the injection well. The existence of a critical pressure buildup above which earthquakes occur is completely consistent with the theory on the role of fluid pressure in fault movement as presented by Hubbert and Rubey.

Colorado↗

Sediment fluxes from California Coastal Rivers: the influences of climate, geology, and topography

The influences of geologic and climatic factors on erosion and sedimentation processes in rivers draining the western flank of the California Coast Range are assessed. Annual suspended, bedload, and total sediment fluxes were determined for 16 river basins that have hydrologic records covering all or most of the period from 1950 to 2006 and have been relatively unaffected by flow storage, regulation, and depletion, which alter the downstream movement of water and sediment. The occurrence of relatively large annual sediment fluxes are strongly influenced by the El Nino&ndash;Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO). The frequency of relatively large annual sediment fluxes decreases from north to south during La Nina phases and increases from north to south during El Nino phases. The influence of ENSO is modulated over a period of decades by the PDO, such that relatively large annual sediment fluxes are more frequent during a La Nina phase in conjunction with a cool PDO and during an El Nino phase in conjunction with a warm PDO. Values of mean annual sediment flux, , were regressed against basin and climatic characteristics. Basin area, bedrock erodibility, basin relief, and precipitation explain 87% of the variation in from the 16 river basins. Bedrock erodibility is the most significant characteristic influencing . Basin relief is a superior predictor of compared with basin slope. is nearly proportional to basin area and increases with increasing precipitation. For a given percentage change, basin relief has a 2.3-fold greater effect on than a similar change in precipitation. The estimated natural from all California coastal rivers for the period 1950&ndash;2006 would have been approximately 85 million tons without flow storage, regulation, and depletion; the actual has been approximately 50 million tons, because of the effects of flow storage, regulation, and depletion.

California↗

Seasonal impoundment management reduces nitrogen cycling but not resilience to surface fire in a tidal wetland

Hydrology and salinity regimes of many impounded wetlands are manipulated to provide seasonal habitats for migratory waterfowl , with little-known consequences for ecosystem structure and function. Managed hydrology can alter ecosystems by directly changing soil properties and processes and by influencing plant community dynamics. Additionally, management history may influence ecosystem response to disturbance, including fires. To better understand how wetland management regime influences ecosystem response to disturbance, we quantified elevation, soil nitrogen concentrations and process rates, and plant community structure and diversity in a natural experiment following the 2018 Branscombe Fire. We measured paired burned-unburned patches in both tidally-influenced and managed, seasonally-impounded wetlands in Suisun Marsh, California, USA. Unburned ecosystem structure and nutrient cycling differed by wetland management history; unburned impounded wetlands were ∼1 m lower in elevation and plant community composition was dominated by succulents whereas the unburned tidal wetland was dominated by graminoids. Unburned impounded wetland soil nitrogen cycling (potential nitrification and denitrification) rates were <28% of those measured in unburned tidal wetland soils and soil extractable nitrate, ammonium, and dissolved inorganic phosphorus concentrations were also substantially lower in unburned impounded than unburned tidal wetlands. Despite these differences in pre-disturbance (i.e., unburned) conditions, all soil processes recovered to baseline levels within 6 months after surface fire, and we found no evidence of plant community change 1 year after fire in either wetland management type. Overall, water management history exerted stronger control on ecosystem processes and structure than surface fire disturbance. Low extractable soil nitrate and potential denitrification rates may indicate limitation of soil nitrogen removal in impounded wetlands, with implications for downstream environmental quality and eutrophication across managed landscapes.

California↗

Hydrologic modeling of two glaciated watersheds in Northeast Pennsylvania

ABSTRACT: A hydrologic modeling study, using the Hydrologic Simulation Program - FORTRAN (HSPF), was conducted in two glaciated watersheds, Purdy Creek and Ariel Creek in northeastern Pennsylvania. Both watersheds have wetlands and poorly drained soils due to low hydraulic conductivity and presence of fragipans. The HSPF model was calibrated in the Purdy Creek watershed and verified in the Ariel Creek watershed for June 1992 to December 1993 period. In Purdy Creek, the total volume of observed stream-flow during the entire simulation period was 13.36 × 10 6 m 3 and the simulated streamflow volume was 13.82 × 10 6 m 3 (5 percent difference). For the verification simulation in Ariel Creek, the difference between the total observed and simulated flow volumes was 17 percent. Simulated peak flow discharges were within two hours of the observed for 30 of 46 peak flow events (discharge greater than 0.1 m 3 /sec) in Purdy Creek and 27 of 53 events in Ariel Creek. For 22 of the 46 events in Purdy Creek and 24 of 53 in Ariel Creek, the differences between the observed and simulated peak discharge rates were less than 30 percent. These 22 events accounted for 63 percent of total volume of streamflow observed during the selected 46 peak flow events in Purdy Creek. In Ariel Creek, these 24 peak flow events accounted for 62 percent of the total flow observed during all peak flow events. Differences in observed and simulated peak flow rates and volumes (on a percent basis) were greater during the snowmelt runoff events and summer periods than for other times.

Journal of the American Water Resources Associatio↗

Nutrient and suspended-sediment concentrations in the Maumee River and tributaries during 2019 rain-induced fallow conditions

Above average precipitation from October 2018 through July 2019 in the Maumee River (R.) Basin resulted in 29% of cropland left fallow, providing a glimpse of potential effects from decreased nutrient application. Ongoing monitoring at 15 water-quality sites on the Maumee R. upstream from Defiance enabled comparison with 2017, which was hydrologically similar to 2019 in precipitation and streamflow. In 2019, nitrate (as nitrogen; NO 3 -N) for March-July was significantly less than previous years (2015–2018), but the response for phosphorus (P) was more complicated. Relative to 2017, total P (TP) was lower at 7 of 15 sites, but higher at 7, reflecting higher suspended sediment (SS). Dissolved P (DP) was generally lower, but less different than NO 3 ; DP was higher at 3 sites. DP-P:NO 3 -N was generally higher in 2019, DP-P:TP was lower, and there was less TP relative to SS. Overall, less P was in the system in 2019. However smaller streams showed a large range of difference between 2019 and 2017 for all constituents, indicating variability in land management and physiography. In contrast, all constituents were lower in 2019 in larger (>5000 km 2 ) streams, including the Maumee R. near Defiance, where the difference in NO 3 (−37%) exceeded that for TP (−16%), DP (−10%), and SS (−20%). Differences in these relations among N, P, and SS indicate that P was available from legacy sources that are more difficult to distinguish during typical agricultural production years and that some material from 2019 was stored in the system upstream from the largest sites.

Indiana, Ohio↗

Characterizing the two-dimensional thermal conductivity distribution in a sand and gravel aquifer

Both hydrologic and thermal transport properties play a significant role in the movement of heat through permeable sedimentary material; however, the thermal conductivity is rarely characterized in detailed spatial resolution. As part of a study of the movement of thermal plumes through a sand and gravel aquifer, we have constructed a two-dimensional profile of thermal conductivity. This work consisted of: (i) measuring the thermal conductivity of the soil solids, λ s , for the main stratigraphic units using the steady-state divided-bar apparatus and estimating conductivity from mineral composition; (ii) measuring the volumetric water content and porosity using crosshole ground-penetrating radar; (iii) evaluating four models used to predict the apparent thermal conductivity, λ, of variably saturated soils and selecting the best model using the information-theoretic approach, (iv) calculating the λ field on a 0.25-m square cell grid using measured data and the selected model, and (v) simulating thermal transport within the two-dimensional domain using a finite element numerical model. The apparent thermal conductivity in the saturated aquifer ranges from 2.14 to 2.69 W m −1 K −1 with a mean of 2.42 W m −1 K −1 Numerical simulations show that the heterogeneous thermal conductivity field results in increased thermal dispersion that is most pronounced at the plume front. Our values for λ and λ s may be used for glacial soils with similar mineralogy and texture. Our methods may also be used at other sites to construct the thermal conductivity distribution.

Ontario↗