USGS ScienceSearch

SEARCH · USGS Science

Results for “Modeling Report”

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 235 records · Page 13Linked to original sources

A numerical program for steady-state flow of magma-gas mixtures through vertical eruptive conduits

This report presents a model that calculates flow properties (pressure, vesicularity, and some 35 other parameters) as a function of vertical position within a volcanic conduit during a steady-state eruption. The model idealizes the magma-gas mixture as a single homogeneousfluid and calculates gas exsolution under the assumption of equilibrium conditions. These are the same assumptions on which classic conduit models (e.g. Wilson and Head, 1981) have been based. They are most appropriate when applied to eruptions of rapidly ascending magma (basaltic lava-fountain eruptions, and Plinian or sub-Plinian eruptions of intermediate or silicic magmas) that contains abundant nucleation sites (microlites, for example) for bubble growth.

Open-File Report

Trends in the quality of water in New Jersey streams, water years 1971–2011

In a study conducted by the U.S. Geological Survey in cooperation with the New Jersey Department of Environmental Protection and the Delaware River Basin Commission, trend tests were conducted on selected water-quality characteristics measured at stations on streams in New Jersey during selected periods over water years 1971‒2011. Tests were conducted on 3 nutrients (total nitrogen, filtered nitrate plus nitrite, and total phosphorus) at 28 water-quality stations. At 4 of these stations, tests were also conducted on 3 measures of major ions (specific conductance, filtered chloride, and total dissolved solids). Two methods were used to identify trends—Weighted Regressions on Time, Discharge, and Season (WRTDS) models and seasonal rank-sum tests. For this report, the use of WRTDS models included the use of the WRTDS Bootstrap Test (WBT). WRTDS models identified trends in flow-normalized annual concentrations and flow-normalized annual fluxes over water years 1980‒2011 and 2000‒11 for each nutrient, filtered chloride, and total dissolved solids. WRTDS models were developed for each nutrient at the 20 or 21 stations at which streamflow was measured or estimated. Trends in nutrient concentration were reported for these stations; trends in nutrient fluxes were reported only for 15–17 of these stations. The results of WRTDS models for water years 1980‒2011 identified more stations with downward trends in concentrations of either total nitrogen or total phosphorus than upward trends. For total nitrogen, there were downward trends at 9 stations and an upward trend at 1 station. For total phosphorus, there were downward trends at 8 stations and an upward trend at 1 station. For filtered nitrate plus nitrite, there were downward trends at 6 stations and upward trends at 6 stations. The result of the trend test in flux for a selected nutrient at a selected station (downward trend, no trend, or upward trend) usually matched the trend result in concentration. Seasonal rank-sum tests, the second method used, identified step trends in water-quality measured in different decades—1970s, 1980s, 1990s, and 2000s. Tests were conducted on all nutrients at 28 stations and on all measures of major ions at the 4 selected stations. Results of seasonal rank-sum tests between the 1980s and the 2000s identified more stations with downward trends in concentrations of total nitrogen (14) than stations with upward trends (2) and more stations with downward trends in concentrations of total phosphorus (18) than stations with upward trends (1). A combined dataset of trend results for concentrations over water years 1980‒2011 was created from the results of the two tests for the period. Results of WRTDS models were included in this combined dataset, if available. Otherwise, the results of the seasonal rank-sum tests between water-quality characteristics measured in the 1980s and 2000s were included. Trend results over water years 1980‒2011 in the combined dataset show that few of the 28 stations had upward trends in concentrations of either total nitrogen or total phosphorus. There were only 2 stations with upward trends in total nitrogen concentration and 1 station with an upward trend in total phosphorus concentration. Results for filtered nitrate plus nitrite show about the same number of stations with upward trends (9) as stations with downward trends (7). Results for all measures of major ions show upward trends at the four stations tested.

New Jersey

Coastal strategies to predict Escherichia coli concentrations for beaches along a 35 km stretch of southern Lake Michigan

To understand the fate and movement of Escherichia coli in beach water, numerous modeling studies have been undertaken including mechanistic predictions of currents and plumes and empirical modeling based on hydrometeorological variables. Most approaches are limited in scope by nearshore currents or physical obstacles and data limitations; few examine the issue from a larger spatial scale. Given the similarities between variables typically included in these models, we attempted to take a broader view of E. coli fluctuations by simultaneously examining twelve beaches along 35 km of Indiana's Lake Michigan coastline that includes five point-source outfalls. The beaches had similar E. coli fluctuations, and a best-fit empirical model included two variables: wave height and an interactive term comprised of wind direction and creek turbidity. Individual beach R2 was 0.32-0.50. Data training-set results were comparable to validation results (R2 = 0.48). Amount of variation explained by the model was similar to previous reports for individual beaches. By extending the modeling approach to include more coastline distance, broader-scale spatial and temporal changes in bacteria concentrations and the influencing factors can be characterized. ?? 2008 American Chemical Society.

Environmental Science & Technology

An empirical model for earthquake probabilities in the San Francisco Bay region, California, 2002-2031

The moment magnitude M 7.8 earthquake in 1906 profoundly changed the rate of seismic activity over much of northern California. The low rate of seismic activity in the San Francisco Bay region (SFBR) since 1906, relative to that of the preceding 55 yr, is often explained as a stress-shadow effect of the 1906 earthquake. However, existing elastic and visco-elastic models of stress change fail to fully account for the duration of the lowered rate of earthquake activity. We use variations in the rate of earthquakes as a basis for a simple empirical model for estimating the probability of M ≥6.7 earthquakes in the SFBR. The model preserves the relative magnitude distribution of sources predicted by the Working Group on California Earthquake Probabilities' ( WGCEP, 1999 ; WGCEP, 2002 ) model of characterized ruptures on SFBR faults and is consistent with the occurrence of the four M ≥6.7 earthquakes in the region since 1838. When the empirical model is extrapolated 30 yr forward from 2002, it gives a probability of 0.42 for one or more M ≥6.7 in the SFBR. This result is lower than the probability of 0.5 estimated by WGCEP ( 1988 ), lower than the 30-yr Poisson probability of 0.60 obtained by WGCEP ( 1999 ) and WGCEP ( 2002 ), and lower than the 30-yr time-dependent probabilities of 0.67, 0.70, and 0.63 obtained by WGCEP ( 1990 ), WGCEP ( 1999 ), and WGCEP ( 2002 ), respectively, for the occurrence of one or more large earthquakes. This lower probability is consistent with the lack of adequate accounting for the 1906 stress-shadow in these earlier reports. The empirical model represents one possible approach toward accounting for the stress-shadow effect of the 1906 earthquake. However, the discrepancy between our result and those obtained with other modeling methods underscores the fact that the physics controlling the timing of earthquakes is not well understood. Hence, we advise against using the empirical model alone (or any other single probability model) for estimating the earthquake hazard and endorse the use of all credible earthquake probability models for the region, including the empirical model, with appropriate weighting, as was done in WGCEP ( 2002 ).

California

User's guide to revised method-of-characteristics solute-transport model (MOC--version 31)

The U.S. Geological Survey computer model to simulate two-dimensional solute transport and dispersion in ground water (Konikow and Bredehoeft, 1978; Goode and Konikow, 1989) has been modified to improve management of input and output data and to provide progressive run-time information. All opening and closing of files are now done automatically by the program. Names of input data files are entered either interactively or using a batch-mode script file. Names of output files, created automatically by the program, are based on the name of the input file. In the interactive mode, messages are written to the screen during execution to allow the user to monitor the status and progress of the simulation and to anticipate total running time. Information reported and updated during a simulation include the current pumping period and time step, number of particle moves, and percentage completion of the current time step. The batch mode enables a user to run a series of simulations consecutively, without additional control. A report of the model's activity in the batch mode is written to a separate output file, allowing later review. The user has several options for creating separate output files for different types of data. The formats are compatible with many commercially available applications, which facilitates graphical postprocessing of model results. Geohydrology and Evaluation of Stream-Aquifer Relations in the Apalachicola-Chattahoochee-Flint River Basin, Southeastern Alabama, Northwestern Florida, and Southwestern Georgia By Lynn J. Torak, Gary S. Davis, George A. Strain, and Jennifer G. Herndon Abstract The lower Apalachieola-Chattahoochec-Flint River Basin is underlain by Coastal Plain sediments of pre-Cretaceous to Quaternary age consisting of alternating units of sand, clay, sandstone, dolomite, and limestone that gradually thicken and dip gently to the southeast. The stream-aquifer system consism of carbonate (limestone and dolomite) and elastic sediments, which define the Upper Floridan aquifer and Intermediate system, in hydraulic connection with the principal rivers of the basin and other surface-water features, natural and man made. Separate digital models of the Upper Flori-dan aquifer and Intermediate system were constructed by using the U.S. Geological Survey's MODular Finite-Element model of two dimensional ground-water flow, based on concep- tualizations of the stream-aquifer system, and calibrated to drought conditions of October 1986. Sensitivity analyses performed on the models indicated that aquifer hydraulic conductivity, lateral and vertical boundary flows, and pumpage have a strong influence on groundwater levels. Simulated pumpage increases in the Upper Floridan aquifer, primarily in the Dougherty Plain physiographic district of Georgia,. caused significant reductions in aquifer discharge to streams that eventually flow to Lake Seminole and the Apalachicola River and Bay. Simulated pumpage increases greater than 3 times the October 1986 rates caused drying ofsome stream reaches and parts of the Upper Floridan aquifer in Georgia. Water budgets prepared from simulation results indicate that ground- water discharge to streams and recharge by horizontal and vertical flow are the principal mechanisms for moving water through the flow system. The potential for changes in ground-water quality is high in areas where chemical constituents can be mobilized by these mechanisms. Less than 2 percent of ground-water discharge to streams comes from the Intermediate system; thus, it plays a minor role in the hydrodynamics of the stream- aquifer system.

Water-Resources Investigations Report

Listings of model input and selected output values for the simulation of ground-water flow near Los Alamos, North-central New Mexico

This report contains listings of model input and selected output values for the simulation of ground-water flow near Los Alamos, north- central New Mexico. This simulation was developed by the U.S. Geological Survey in cooperation with Los Alamos National Laboratory (Frenzel, P.F., 1995, Geohydrology and simulation of ground-water flow near Los Alamos, north-central New Mexico: U.S. Geological Survey Water-Resources Investigations Report 95-4091). The simulation used the U.S. Geological Survey modular flow model code (McDonald, M.G., and Harbaugh, A.W., 1988, A modular three-dimensional finite-difference ground-water flow model: Techniques of Water-Resources Investigations of the United States Geological Survey, book 6, chap. A1). The listings in this report are in compressed format on a 1.44-megabyte IBM-PC1 compatible floppy disk. Software is included for decompression to ASCII format.

Open-File Report

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences

Measured and simulated runoff to the lower Charles River, Massachusetts, October 1999–September 2000

The lower Charles River, the water body between the Watertown Dam and the New Charles River Dam, is an important recreational resource for the Boston, Massachusetts, metropolitan area, but impaired water quality has affected its use. The goal of making this resource fishable and swimmable requires a better understanding of combined-sewer-overflow discharges, non-combined-sewer-overflow stormwater runoff, and constituent loads. This report documents the modeling effort used to calculate non-combined-sewer-overflow runoff to the lower Charles River. During the 2000 water year, October 1, 1999–September 30, 2000, the U.S. Geological Survey collected precipitation data at Watertown Dam and compiled data from five other precipitation gages in or near the watershed. In addition, surface-water discharge data were collected at eight sites—three relatively homogenous land-use sites, four major tributary sites, and the Charles River at Watertown Dam, which is the divide between the upper and lower watersheds. The precipitation and discharge data were used to run and calibrate Stormwater Management Models developed for the three land-use subbasins (single-family, multi-family, and commercial), and the two tributary subbasins (Laundry and Faneuil Brooks). These calibrated models were used to develop a sixth model to simulate 54 ungaged outfalls to the lower Charles River. Models developed by the U.S. Geological Survey at gaged sites were calibrated with up to 24 storms. Each model was evaluated by comparing simulated discharge against measured discharge for all storms with appreciable precipitation and reliable discharge data. The model-fit statistics indicated that the models generally were well calibrated to peak discharge and runoff volumes. The model fit of the commercial land-use subbasin was not as well calibrated compared to the other models because the measured flows appear to be affected by variable conditions not represented in the model. A separate Stormwater Management Model of the Stony Brook Subbasin previously developed by others was evaluated with the newly collected data from this study; this model had a model fit comparable to the models developed by the U.S. Geological Survey. The total annual runoff to the lower Charles River during the 2000 water year, not including contributions from combined-sewer-overflows except from the Stony Brook Subbasin, was 16,500 million cubic feet; 92 percent of the inflow was from the Charles River above Watertown Dam, 3 percent was from the Stony Brook Subbasin, 2 percent was from the Muddy River Subbasin, and less than 1 percent was from the combined inflows of Laundry and Faneuil Brooks. The remaining ungaged drainage area contributed about 2 percent of the total annual inflow to the lower Charles River. Excluding discharge from the Charles River above Watertown Dam, total annual runoff to the lower Charles River was 1,240 million cubic feet; 39 percent was from the Stony Brook Subbasin, 27 percent was from the Muddy River, which includes runoff that drains to the Muddy River conduit, 7 percent was from the Laundry Brook Subbasin, and 4 percent was from the Faneuil Brook Subbasin. Flow from the ungaged areas composed about 23 percent of the total annual inflow to the lower Charles River, excluding discharge from the Charles River above Watertown Dam. Runoff to the lower Charles River was calculated for two design storms representing a 3-month and a 1-year event, 1.84 and 2.79 inches of total rainfall, respectively. These simulated discharges were provided to the Massachusetts Water Resources Authority for use in a receiving-water model of the lower Charles River. Total storm runoff to the lower Charles River was 111 and 257 million cubic feet for the 3-month and 1-year storms, respectively. Excluding discharge from the Charles River above Watertown Dam, total runoff to the lower Charles River was 30 and 53 million cubic feet for the 3-month and 1-year storms, respectively. Runoff from the various tributary areas for the design storms was about in the same proportion as that for the annual runoff.

Massachusetts

Uses of epistemic uncertainties in the USGS National Seismic Hazard Models

The need for US Geological Survey (USGS) National Seismic Hazard Models (NSHMs) to report estimates of epistemic uncertainties in the hazard (e.g. fractile hazard curves) in all forthcoming releases is increasing. With fractile hazard curves as potential new outputs from the USGS 2023 NSHM, a simultaneous need is to help end-users better understand these epistemic uncertainties and clarify their potential uses. In this article, we address the latter need by (1) characterizing epistemic uncertainties in two updates of the USGS NSHM (2014 for California and 2021 for Hawaii), (2) illustrating a variety of downstream applications of fractile hazard curves in both hazard and risk contexts, and (3) discussing implications from the various types of uncertainties. We found that the epistemic uncertainty in hazard is generally larger for Hawaii than for California, the epistemic uncertainty in hazard can be reasonably approximated with a lognormal distribution for most of the cases considered, and the correlation between epistemic uncertainty in hazard at two different intensity measure levels generally varies with both location and type of intensity measure. Furthermore, we developed models for readily generating approximate fractile hazard curves in California and Hawaii. Finally, given the complexities involved in the hazard modeling process, we developed an open-source interactive tool to enable a broad range of users to independently examine and potentially start using such epistemic uncertainties for their respective applications.

Earthquake Spectra

MODFLOW-2000, the U.S. Geological Survey modular ground-water model -- Documentation of the Model-Layer Variable-Direction Horizontal Anisotropy (LVDA) capability of the Hydrogeologic-Unit Flow (HUF) package

This report documents the model-layer variable-direction horizontal anisotropy (LVDA) capability of the Hydrogeologic-Unit Flow (HUF) Package of MODFLOW-2000. The LVDA capability allows the principal directions of horizontal anisotropy to be different than the model-grid row and column directions, and for the directions to vary on a cell-by-cell basis within model layers. The HUF Package calculates effective hydraulic properties for model grid cells based on hydraulic properties of hydrogeologic units with thicknesses defined independently of the model layers. These hydraulic properties include, among other characteristics, hydraulic conductivity and a horizontal anisotropy ratio. Using the LVDA capability, horizontal anisotropy direction is defined for model grid cells within which one or more hydrogeologic units may occur. For each grid cell, the HUF Package calculates the effective horizontal hydraulic conductivity along the primary direction of anisotropy using the hydrogeologic-unit hydraulic conductivities, and calculates the effective horizontal hydraulic conductivity along the orthogonal anisotropy direction using the effective primary direction hydraulic conductivities and horizontal anisotropy ratios. The direction assigned to the model layer effective primary hydraulic conductivity is specified using a new data set defined by the LVDA capability, when active, to calculate coefficients needed to solve the ground-water flow equation. Use of the LVDA capability is illustrated in four simulation examples, which also serve to verify hydraulic heads, advective-travel paths, and sensitivities calculated using the LVDA capability. This version of the LVDA capability defines variable-direction horizontal anisotropy using model layers, not the hydrogeologic units defined by the HUF Package. This difference needs to be taken into account when designing model layers and hydrogeologic units to produce simulations that accurately represent a given field problem. This might be a reason, for example, to make model layer boundaries coincide with hydrogeologic-unit boundaries in all or part of a model grid.

Open-File Report

To cross or not to cross: modeling wildlife road crossings as a binary response variable with contextual predictors

Roads are significant barriers to landscape-scale movements of individuals or populations of many wildlife taxa. The decision by an animal near a road to either cross or not cross may be influenced by characteristics of the road, environmental conditions, traits of the individual animal, and other aspects of the context within which the decision is made. We considered such factors in a mixed-effects logistic regression model describing the nightly road crossing probabilities of invasive nocturnal Brown Treesnakes ( Boiga irregularis ) through short-term radiotracking of 691 snakes within close proximity to 50 road segments across the island of Guam. All measures of road magnitude (traffic volume, gap width, surface type, etc.) were significantly negatively correlated with crossing probabilities. Snake body size was the only intrinsic factor associated with crossing rates, with larger snakes crossing roads more frequently. Humidity was the only environmental variable affecting crossing rate. The distance of the snake from the road at the start of nightly movement trials was the most significant predictor of crossings. The presence of snake traps with live mouse lures during a portion of the trials indicated that localized prey cues reduced the probability of a snake crossing the road away from the traps, suggesting that a snake's decision to cross roads is influenced by local foraging opportunities. Per capita road crossing rates of Brown Treesnakes were very low, and comparisons to historical records suggest that crossing rates have declined in the 60+ yr since introduction to Guam. We report a simplified model that will allow managers to predict road crossing rates based on snake, road, and contextual characteristics. Road crossing simulations based on actual snake size distributions demonstrate that populations with size distributions skewed toward larger snakes will result in a higher number of road crossings. Our method of modeling per capita road crossing probabilities as a binary response variable, influenced by contextual factors, may be useful for describing or predicting road crossings by individuals of other taxa provided that appropriate spatial and temporal resolution can be achieved and that potentially influential covariate data can be obtained.

Ecosphere

Simulation of groundwater and surface-water resources and evaluation of water-management alternatives for the Chamokane Creek basin, Stevens County, Washington

A three-dimensional, transient numerical model of groundwater and surface-water flow was constructed for Chamokane Creek basin to better understand the groundwater-flow system and its relation to surface-water resources. The model described in this report can be used as a tool by water-management agencies and other stakeholders to quantitatively evaluate the effects of potential increases in groundwater pumping on groundwater and surface-water resources in the basin. The Chamokane Creek model was constructed using the U.S. Geological Survey (USGS) integrated model, GSFLOW. GSFLOW was developed to simulate coupled groundwater and surface-water resources. The model uses 1,000-foot grid cells that subdivide the model domain by 102 rows and 106 columns. Six hydrogeologic units in the model are represented using eight model layers. Daily precipitation and temperature were spatially distributed and subsequent groundwater recharge was computed within GSFLOW. Streamflows in Chamokane Creek and its major tributaries are simulated in the model by routing streamflow within a stream network that is coupled to the groundwater-flow system. Groundwater pumpage and surface-water diversions and returns specified in the model were derived from monthly and annual pumpage values previously estimated from another component of this study and new data reported by study partners. The model simulation period is water years 1980-2010 (October 1, 1979, to September 30, 2010), but the model was calibrated to the transient conditions for water years 1999-2010 (October 1, 1998, to September 30, 2010). Calibration was completed by using traditional trial-and-error methods and automated parameter-estimation techniques. The model adequately reproduces the measured time-series groundwater levels and daily streamflows. At well observation points, the mean difference between simulated and measured hydraulic heads is 7 feet with a root-mean-square error divided by the total difference in water levels of 4.7 percent. Simulated streamflow was compared to measured streamflow at the USGS streamflow-gaging station-Chamokane Creek below Falls, near Long Lake (12433200). Annual differences between measured and simulated streamflow for the site ranged from -63 to 22 percent. Calibrated model output includes a 31-year estimate of monthly water budget components for the hydrologic system. Five model applications (scenarios) were completed to obtain a better understanding of the relation between groundwater pumping and surface-water resources. The calibrated transient model was used to evaluate: (1) the connection between the upper- and middle-basin groundwater systems, (2) the effect of surface-water and groundwater uses in the middle basin, (3) the cumulative impacts of claims registry use and permit-exempt wells on Chamokane Creek streamflow, (4) the frequency of regulation due to impacted streamflow, and (5) the levels of domestic and stockwater use that can be regulated. The simulation results indicated that streamflow is affected by existing groundwater pumping in the upper and middle basins. Simulated water-management scenarios show streamflow increased relative to historical conditions as groundwater and surface-water withdrawals decreased.

Washington

Simulation of Flow Regimes to Reduce Habitat for T. tubifex

Whirling disease has had a significant impact on trout fisheries of the American west by reducing the numbers and quality of rainbow trout in infected streams. A critical factor in the life cycle of the whirling disease parasite is the fine sediment that provides the optimum habitat for Tubifex tubifex, an oligochaete worm that acts as an intermediate host for the disease. This report presents a model for the simulation of flushing flows required to remove undesirable fines and sand from a pool. Undesirable fines may also need to be flushed from runs, the surface layer, and backwater areas. Well-defined links of specific particle sizes to oligochaete worm abundance is needed to justify the use of flushing flows to move sediment. An analytical method for estimating the streamflows needed to remove the fine sediment is demonstrated herein. The overall steps to follow in removing fines from a stream are: Step 1. Determine size of the sediment that is the habitat for oligochaete worms. Step 2. Determine location of the sediment that is the habitat for oligochaete worms. Step 3. Determine streamflows needed to flush (remove) the sediment that is the habitat for oligochaete worms. The case study approach is used to present the method and to demonstrate its application. The case is derived from the sediment and oligochaete worm habitat of Willow Creek, a tributary of the Upper Colorado River located in Grand County, Colo. Willow Creek Reservoir (an element of the Colorado-Big Thompson Project) controls the streamflows of the creek and is just above the study site.

Open-File Report

Simulation of daily streamflows at gaged and ungaged locations within the Cedar River Basin, Iowa, using a Precipitation-Runoff Modeling System model

The U.S. Geological Survey, in cooperation with the Iowa Department of Natural Resources, conducted a study to examine techniques for estimation of daily streamflows using hydrological models and statistical methods. This report focuses on the use of a hydrologic model, the U.S. Geological Survey's Precipitation-Runoff Modeling System, to estimate daily streamflows at gaged and ungaged locations. The Precipitation-Runoff Modeling System is a modular, physically based, distributed-parameter modeling system developed to evaluate the impacts of various combinations of precipitation, climate, and land use on surface-water runoff and general basin hydrology. The Cedar River Basin was selected to construct a Precipitation-Runoff Modeling System model that simulates the period from January 1, 2000, to December 31, 2010. The calibration period was from January 1, 2000, to December 31, 2004, and the validation periods were from January 1, 2005, to December 31, 2010 and January 1, 2000 to December 31, 2010. A Geographic Information System tool was used to delineate the Cedar River Basin and subbasins for the Precipitation-Runoff Modeling System model and to derive parameters based on the physical geographical features. Calibration of the Precipitation-Runoff Modeling System model was completed using a U.S. Geological Survey calibration software tool. The main objective of the calibration was to match the daily streamflow simulated by the Precipitation-Runoff Modeling System model with streamflow measured at U.S. Geological Survey streamflow gages. The Cedar River Basin daily streamflow model performed with a Nash-Sutcliffe efficiency ranged from 0.82 to 0.33 during the calibration period, and a Nash-Sutcliffe efficiency ranged from 0.77 to -0.04 during the validation period. The Cedar River Basin model is meeting the criteria of greater than 0.50 Nash-Sutcliffe and is a good fit for streamflow conditions for the calibration period at all but one location, Austin, Minnesota. The Precipitation-Runoff Modeling System model accurately simulated streamflow at four of six uncalibrated sites within the basin. Overall, there was good agreement between simulated and measured seasonal and annual volumes throughout the basin for calibration and validation sites. The calibration period ranged from 0.2 to 20.8 percent difference, and the validation period ranged from 0.0 to 19.5 percent difference across all seasons and total annual runoff. The Precipitation-Runoff Modeling System model tended to underestimate lower streamflows compared to the observed streamflow values. This is an indication that the Precipitation-Runoff Modeling model needs more detailed groundwater and storage information to properly model the low-flow conditions in the Cedar River Basin.

Iowa

Quantifying watershed-scale groundwater loading and in-stream fate of nitrate using high-frequency water quality data

We describe a new approach that couples hydrograph separation with high-frequency nitrate data to quantify time-variable groundwater and runoff loading of nitrate to streams, and the net in-stream fate of nitrate at the watershed-scale. The approach was applied at three sites spanning gradients in watershed size and land use in the Chesapeake Bay watershed. Results indicate that 58-73% of the annual nitrate load to the streams was groundwater-discharged nitrate. Average annual first order nitrate loss rate constants (k) were similar to those reported in both modelling and in-stream process-based studies, and were greater at the small streams (0.06 and 0.22 d -1 ) than at the large river (0.05 d -1 ), but 11% of the annual loads were retained/lost in the small streams, compared with 23% in the large river. Larger streambed area to water volume ratios in small streams result in greater loss rates, but shorter residence times in small streams result in a smaller fraction of nitrate loads being removed than in larger streams. A seasonal evaluation of k values suggests that nitrate was retained/lost at varying rates during the growing season. Consistent with previous studies, streamflow and nitrate concentration were inversely related to k. This new approach for interpreting high-frequency nitrate data and the associated findings furthers our ability to understand, predict, and mitigate nitrate impacts on streams and receiving waters by providing insights into temporal nitrate dynamics that would be difficult to obtain using traditional field-based studies.

Delaware, Pennsylvania, Maryland, New York, Virgin

The effects of fire on avian communities: Spatio-temporal attributes of the literature 19122003

We reviewed the temporal, geographic, and biogeographic distribution, as well as relevant research and publication attributes, of 512 documents addressing the effects of fire on avian communities, to provide an assessment of the scope of this literature and recommendations for future research. We summarized relevant attributes of all documents to identify patterns that were then tested against appropriate null models. Most documents reported on original research, with the literature evenly divided between studies investigating controlled fire and those reporting on uncontrolled wildfires. Conceptual reviews made up the second largest category; methodological reviews, bibliographies, and meta-analyses were rare. Although the literature examined spans nearly a century, most documents were published within the last 15 years, with new literature being added at an increasing rate. However, increases seem to be skewed towards original research at the expense of synthesis. An overwhelming majority of documents were published in peer-reviewed scientific journals and in English. Other important publication outlets included MS and PhD theses and conference proceedings. The spatial distribution of documents by continent and biogeographic domain and division differed significantly from expectations based on land area. Future research on avian community response to fire should focus on (1) continued synthesis, emphasizing methodological reviews, bibliographies, and North America; (2) increasing research efforts in areas currently underrepresented in the literature, including Africa, Asia, and South and Central America; and (3) meta-analyses. ?? 2009 IAWF.

International Journal of Wildland Fire

Thermal conditions predict intraspecific variation in senescence rate in frogs and toads

Variation in temperature is known to influence mortality patterns in ectotherms. Even though a few experimental studies on model organisms have reported a positive relationship between temperature and actuarial senescence (i.e., the increase in mortality risk with age), how variation in climate influences the senescence rate across the range of a species is still poorly understood in free-ranging animals. We filled this knowledge gap by investigating the relationships linking senescence rate, adult lifespan, and climatic conditions using long-term capture–recapture data from multiple amphibian populations. We considered two pairs of related anuran species from the Ranidae ( Rana luteiventris and Rana temporaria ) and Bufonidae ( Anaxyrus boreas and Bufo bufo ) families, which diverged more than 100 Mya and are broadly distributed in North America and Europe. Senescence rates were positively associated with mean annual temperature in all species. In addition, lifespan was negatively correlated with mean annual temperature in all species except A. boreas . In both R. luteiventris and A. boreas , mean annual precipitation and human environmental footprint both had negligible effects on senescence rates or lifespans. Overall, our findings demonstrate the critical influence of thermal conditions on mortality patterns across anuran species from temperate regions. In the current context of further global temperature increases predicted by Intergovernmental Panel on Climate Change scenarios, a widespread acceleration of aging in amphibians is expected to occur in the decades to come, which might threaten even more seriously the viability of populations and exacerbate global decline.

Proceedings of the National Academy of Sciences

Desert bird associations with broad-scale boundary length: Applications in avian conservation

1. Current understanding regarding the effects of boundaries on bird communities has originated largely from studies of forest-non-forest boundaries in mesic systems. To assess whether broad-scale boundary length can affect bird community structure in deserts, and to identify patterns and predictors of species' associations useful in avian conservation, we studied relations between birds and boundary-length variables in Chihuahuan Desert landscapes. Operationally, a boundary was the border between two adjoining land covers, and broad-scale boundary length was the total length of such borders in a large area. 2. Within 2-km radius areas, we measured six boundary-length variables. We analysed bird-boundary relations for 26 species, tested for assemblage-level patterns in species' associations with boundary-length variables, and assessed whether body size, dispersal ability and cowbird-host status were correlates of these associations. 3. The abundances or occurrences of a significant majority of species were associated with boundary-length variables, and similar numbers of species were related positively and negatively to boundary-length variables. 4. Disproportionately small numbers of species were correlated with total boundary length, land-cover boundary length and shrubland-grassland boundary length (variables responsible for large proportions of boundary length). Disproportionately large numbers of species were correlated with roadside boundary length and riparian vegetation-grassland boundary length (variables responsible for small proportions of boundary length). Roadside boundary length was associated (positively and negatively) with the most species. 5. Species' associations with boundary-length variables were not correlated with body size, dispersal ability or cowbird-host status. 6. Synthesis and applications. For the species we studied, conservationists can use the regressions we report as working models to anticipate influences of boundary-length changes on bird abundance and occurrence, and to assess avifaunal composition for areas under consideration for protection. Boundary-length variables associated with a disproportionate or large number of species can be used as foci for landscape management. Assessing the underlying causes of bird-boundary relations may improve the prediction accuracy of associated models. We therefore advocate local- and broad-scale manipulative experiments involving the boundary types with which species were correlated, as indicated by the regressions. ?? 2008 The Authors.

Journal of Applied Ecology