USGS ScienceSearch

SEARCH · USGS Science

Results for “Engineering Analysis”

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 397 records · Page 22Linked to original sources

Analysis of borehole-radar reflection data from Machiasport, Maine, December 2003

In December 2003, the U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, collected borehole-radar reflection logs in two boreholes in Machiasport, Maine. These bedrock boreholes were drilled as part of a hydrogeologic investigation of the area surrounding the former Air Force Radar Tracking Station site on Howard Mountain near Bucks Harbor. The boreholes, MW09 and MW10, are located approximately 50 meters (m) from, and at the site of, respectively, the locations of former buildings where trichloroethylene was used as part of defense-site operations. These areas are thought to be potential source areas for contamination that has been detected in downgradient bedrock wells. This investigation focused on testing borehole-radar methods at this site. Single-hole radar-reflection surveys were used to identify the depth, orientation, and spatial continuity of reflectors that intersect and surround the boreholes. In addition, the methods were used to (1) identify the radial depth of penetration of the radar waves in the electrically resistive bimodal volcanic formation at the site, (2) provide information for locating additional boreholes at the site, and (3) test the potential applications of borehole-radar methods for further aquifer characterization and (or) evaluation of source-area remediation efforts. Borehole-radar reflection logging uses a pair of downhole transmitting and receiving antennas to record the reflected wave amplitude and transit time of high-frequency electromagnetic waves. For this investigation, 60- and 100-megahertz antennas were used. The electromagnetic waves emitted by the transmitter penetrate into the formation surrounding the borehole and are reflected off of a material with different electromagnetic properties, such as a fracture or change in rock type. Single-hole directional radar surveys indicate the bedrock surrounding these boreholes is highly fractured, because several reflectors were identified in the radar-reflection data. There are several steeply dipping reflectors with orientations similar to the fracture patterns observed with borehole imaging techniques and in outcrops. The radar-reflection data showed that the vitrophyre in borehole MW09 was more highly fractured than the underlying gabbroic unit. The velocities of radar waves in the bedrock surrounding the boreholes were determined using single-hole vertical radar profiling. Velocities of 114 and 125 meters per microsecond were used to determine the distance to reflectors, the radial depth of penetration, and the dip of reflectors. The bimodal volcanic units appear to be ideal for radar-wave propagation. For the radar surveys collected at this site, radar reflections were detected up to 40 m into the rock from the borehole. These results indicate that boreholes could conservatively be spaced about 15-20 m apart for hole-to-hole radar methods to be effective for imaging between the boreholes and monitoring remediation. Integrated analysis of drilling and borehole-geophysical logs indicates the vitrophyric formation is more fractured than the more mafic gabbroic units in these boreholes. There does not, however, appear to be a quantifiable difference in the radar-wave penetration in these two rock units.

Maine

Finite-difference modeling and dispersion analysis of high-frequency love waves for near-surface applications

Love-wave propagation has been a topic of interest to crustal, earthquake, and engineering seismologists for many years because it is independent of Poisson's ratio and more sensitive to shear (S)-wave velocity changes and layer thickness changes than are Rayleigh waves. It is well known that Love-wave generation requires the existence of a low S-wave velocity layer in a multilayered earth model. In order to study numerically the propagation of Love waves in a layered earth model and dispersion characteristics for near-surface applications, we simulate high-frequency (>5 Hz) Love waves by the staggered-grid finite-difference (FD) method. The air-earth boundary (the shear stress above the free surface) is treated using the stress-imaging technique. We use a two-layer model to demonstrate the accuracy of the staggered-grid modeling scheme. We also simulate four-layer models including a low-velocity layer (LVL) or a high-velocity layer (HVL) to analyze dispersive energy characteristics for near-surface applications. Results demonstrate that: (1) the staggered-grid FD code and stress-imaging technique are suitable for treating the free-surface boundary conditions for Love-wave modeling, (2) Love-wave inversion should be treated with extra care when a LVL exists because of a lack of LVL information in dispersions aggravating uncertainties in the inversion procedure, and (3) energy of high modes in a low-frequency range is very weak, so that it is difficult to estimate the cutoff frequency accurately, and "mode-crossing" occurs between the second higher and third higher modes when a HVL exists. ?? 2010 Birkh??user / Springer Basel AG.

Pure and Applied Geophysics

Cost analysis of ground-water supplies in the North Atlantic region, 1970

The cost of municipal and industrial ground water (or, more specifically, large supplies of ground water) at the wellhead in the North Atlantic Region in 1970 generally ranged from 1.5 to 5 cents per thousand gallons. Water from crystalline rocks and shale is relatively expensive. Water from sandstone is less so. Costs of water from sands and gravels in glaciated areas and from Coastal Plain sediments range from moderate to very low. In carbonate rocks costs range from low to fairly high. The cost of ground water at the wellhead is low in areas of productive aquifers, but owing to the cost of connecting pipe, costs increase significantly in multiple-well fields. In the North Atlantic Region, development of small to moderate supplies of ground water may offer favorable cost alternatives to planners, but large supplies of ground water for delivery to one point cannot generally be developed inexpensively. Well fields in the less productive aquifers may be limited by costs to 1 or 2 million gallons a day, but in the more favorable aquifers development of several tens of millions of gallons a day may be practicable and inexpensive. Cost evaluations presented cannot be applied to any one specific well or specific site because yields of wells in any one place will depend on the local geologic and hydrologic conditions; however, with such cost adjustments as may be necessary, the methodology presented should have wide applicability. Data given show the cost of water at the wellhead based on the average yield of several wells. The cost of water delivered by a well field includes costs of connecting pipe and of wells that have the yields and spacings specified. Cost of transport of water from the well field to point of consumption and possible cost of treatment are not evaluated. In the methodology employed, costs of drilling and testing, pumping equipment, engineering for the well field, amortization at 5% percent interest, maintenance, and cost of power are considered. The report includes an analysis of test drilling costs leading to a production well field. The discussion shows that test drilling is a relatively low cost item and that more than a minimum of test holes in a previously unexplored area is, above all, simple insurance in keeping down costs and may easily result in final lower costs for the system. Use of the jet drill for testing is considered short sighted and may result in higher total costs and possibly failure to discover good aquifers. Economic development of ground water supplies will depend on obtaining qualified hydrologic and engineering advice, on carrying out adequate test drilling, and on utilizing high-quality (at times, more costly) material.

Water Supply Paper

Geophysical logging For hydrogeology

Geophysical logging is the measurement and analysis of electrical, acoustic, nuclear, and other physical properties in a borehole using wireline or direct push technology. Geophysical logging is one of the primary methods of collecting subsurface information for hydrogeologic investigations. Groundwater scientists and engineers should have a basic understanding of borehole geophysics and how it is applied in the characterization of aquifer systems. In this book, the different types of geophysical logs and equipment used in hydrogeologic investigations are reviewed. Applications and analysis methods are illustrated by specific examples in both unconsolidated deposits and bedrock settings. Key references for further exploration of the methods and applications are provided. Hands-on exercises are included for students to develop and practice log-analysis skills.

Book

Possible reestablishment of sawtimber production following mining in the Osceola National Forest, an analysis of the state-of-the-art, February 9, 1978

There is a good possibility that sawtimber production could be reestablished on the Osceola National Forest. Engineering technology that is coming on-line should allow the establishment of land/water relationships that meet the requirements of desired tree species. Information available on soil-plant relationships indicate an absence of phytotoxic materials in various mining wastes and that good tree growth is possible.

Florida

Coastal Louisiana ecosystem assessment and restoration program: The role of ecosystem forecasting in evaluating restoration planning in the Mississippi River Deltaic Plain

The development of ecosystem management plans to restore and rehabilitate natural resources requires an understanding of how specific ecological mechanisms regulate the structure and function of ecosystems. To achieve restoration goals, comprehensive plans and engineering designs must effectively change environmental drivers at the regional level to reduce stress conditions at the local environment that are responsible for ecosystem degradation. This document focuses on the Coastal Louisiana Ecosystem Assessment and Restoration (CLEAR) ecosystem forecasting framework and how it can be used to support the analysis of Louisiana’s coastal restoration plans. Specifically, the framework is designed to (1) develop and incorporate conceptual ecological models that can be used to integrate ecological needs and opportunities with engineering designs, (2) utilize wetland loss rates to describe the most likely “future without” scenario for a variety of ecosystem attributes, (3) estimate broad ecosystem responses to restoration alternatives based on processes associated with succession of geomorphic and ecological systems, and (4) calculate ecological benefits for incorporation into decision support tools associated with large-scale geomorphic and hydrologic processes. This paper provides a brief overview of the spatial framework and modular design of the CLEAR ecosystem forecasting framework and describes in greater detail the evolution of the landscape change module, concepts for its refinement, and how it was utilized in evaluating a coastal restoration alternative proposed in the Coastal Protection and Restoration Authority Preliminary Draft Master Plan. Such projections by the CLEAR forecasting framework can evaluate processes and conditions that result in sustainable coastal ecosystems with habitat functions that support higher trophic levels.

Louisiana

A global ensemble of ocean wave climate statistics from contemporary wave reanalysis and hindcasts

There are numerous global ocean wave reanalysis and hindcast products currently being distributed and used across different scientific fields. However, there is not a consistent dataset that can sample across all existing products based on a standardized framework. Here, we present and describe the first coordinated multi-product ensemble of present-day global wave fields available to date. This dataset, produced through the Coordinated Ocean Wave Climate Project (COWCLIP) phase 2, includes general and extreme statistics of significant wave height ( H s ), mean wave period ( T m ) and mean wave direction ( θ m ) computed across 1980–2014, at different frequency resolutions (monthly, seasonally, and annually). This coordinated global ensemble has been derived from fourteen state-of-the-science global wave products obtained from different atmospheric reanalysis forcing and downscaling methods. This data set has been processed, under a specific framework for consistency and quality, following standard Data Reference Syntax, Directory Structures and Metadata specifications. This new comprehensive dataset provides support to future broad-scale analysis of historical wave climatology and variability as well as coastal risk and vulnerability assessments across offshore and coastal engineering applications.

Scientific Data

Discharge estimation using video recordings from small unoccupied aircraft systems

Measurement of river discharge during flooding events has especially been a challenging and dangerous task in the southwestern US, where flows can be flashy, laden with sediment, and at high velocity. Small unoccupied aircraft systems (sUAS) can be deployed to access unsafe field sites and capture imagery for measuring surface flow velocity and discharge. This paper compares flow discharge estimation at eight field sites—located at or near USGS gauging stations—using time-averaged surface velocities and the turbulence dissipation rate (TDR) derived from large-scale particle image velocimetry (LSPIV) analysis of sUAS videos with conventional measurement techniques conducted by professional USGS hydrographers. Sites characteristics include both natural and engineered channels. The conventional measured discharges were treated as the reference discharges for evaluating the accuracy of the LSPIV discharge estimates. This study evaluated four approaches to estimate the depth-averaged or cross-sectional averaged velocity: constant-velocity index, logarithmic law, power-law, and the entropy method. Results showed the discharges can be accurately calculated by using any of these methods, and that choice of method depended on width to depth ratios.

Journal of Hydraulic Engineering

Volcanogenic uranium deposits: Geology, geochemical processes, and criteria for resource assessment

Felsic volcanic rocks have long been considered a primary source of uranium for many kinds of uranium deposits, but volcanogenic uranium deposits themselves have generally not been important resources. Until the past few years, resource summaries for the United States or the world generally include volcanogenic in the broad category of "other deposits" because they comprised less than 0.5 percent of past production or estimated resources. Exploration in the United States from the 1940s through 1982 discovered hundreds of prospects in volcanic rocks, of which fewer than 20 had some recorded production. Intensive exploration in the late 1970s found some large deposits, but low grades (less than about 0.10 percent U 3 O 8 ) discouraged economic development. A few deposits in the world, drilled in the 1980s and 1990s, are now known to contain large resources (>20,000 tonnes U 3 O 8 ). However, research on ore-forming processes and exploration for volcanogenic deposits has lagged behind other kinds of uranium deposits and has not utilized advances in understanding of geology, geochemistry, and paleohydrology of ore deposits in general and epithermal deposits in particular. This review outlines new ways to explore and assess for volcanogenic deposits, using new concepts of convection, fluid mixing, and high heat flow to mobilize uranium from volcanic source rocks and form deposits that are postulated to be large. Much can also be learned from studies of epithermal metal deposits, such as the important roles of extensional tectonics, bimodal volcanism, and fracture-flow systems related to resurgent calderas. Regional resource assessment is helped by genetic concepts, but hampered by limited information on frontier areas and undiscovered districts. Diagnostic data used to define ore deposit genesis, such as stable isotopic data, are rarely available for frontier areas. A volcanic environment classification, with three classes (proximal, distal, and pre-volcanic structures), permits use of geologic features on 1:500,000 to 1:100,000 scale maps. Geochemical databases for volcanic rocks are postulated to be more effective than databases for stream sediments or surface radioactivity, both of which tend to be inconsistent because of variable leaching of uranium from soils. Based on empirical associations, spatial associations with areas of wet paleoclimate, adjacent oil and gas fields, or evaporite beds are deemed positive. Most difficult to estimate is the location of depositional traps and reduction zones, in part because they are mere points at regional scale. Grade and tonnage data are reviewed and discussed for 32 deposits in the world. Experience of mining engineers and geologists in Asia suggests that tonnages could be higher than presently known in the Western Hemisphere. Geological analysis, and new data from Asia, suggest a typical or median deposit tonnage of about 5,000 tonnes U 3 O 8 , and an optimistic forecast of discoveries in the range of 5,000 to 20,000 tonnes U 3 O 8 . The likely grade of undiscovered deposits could be about 0.15 percent U 3 O 8 , based on both western and eastern examples. Volcanic terrane is under-explored, relative to other kinds of uranium deposits, and is considered a favorable frontier area for new discoveries.

Open-File Report

A probabilistic estimate of maximum acceleration in rock in the contiguous United States

This paper presents a probabilistic estimate of the maximum ground acceleration to be expected from earthquakes occurring in the contiguous United States. It is based primarily upon the historic seismic record which ranges from very incomplete before 1930 to moderately complete after 1960. Geologic data, primarily distribution of faults, have been employed only to a minor extent, because most such data have not been interpreted yet with earthquake hazard evaluation in mind. The map provides a preliminary estimate of the relative hazard in various parts of the country. The report provides a method for evaluating the relative importance of the many parameters and assumptions in hazard analysis. The map and methods of evaluation described reflect the current state of understanding and are intended to be useful for engineering purposes in reducing the effects of earthquakes on buildings and other structures. Studies are underway on improved methods for evaluating the relativ( earthquake hazard of different regions. Comments on this paper are invited to help guide future research and revisions of the accompanying map. The earthquake hazard in the United States has been estimated in a variety of ways since the initial effort by Ulrich (see Roberts and Ulrich, 1950). In general, the earlier maps provided an estimate of the severity of ground shaking or damage but the frequency of occurrence of the shaking or damage was not given. Ulrich's map showed the distribution of expected damage in terms of no damage (zone 0), minor damage (zone 1), moderate damage (zone 2), and major damage (zone 3). The zones were not defined further and the frequency of occurrence of damage was not suggested. Richter (1959) and Algermissen (1969) estimated the ground motion in terms of maximum Modified Mercalli intensity. Richter used the terms "occasional" and "frequent" to characterize intensity IX shaking and Algermissen included recurrence curves for various parts of the country in the paper accompanying his map. The first probabilistic hazard maps covering portions of the United States were by Milne and Davenport (1969a). Recently, Wiggins, Hirshberg and Bronowicki (1974) prepared a probabilistic map of maximum particle velocity and Modified Mercalli intensity for the entire United States. The maps are based on an analysis of the historical seismicity. In general, geological data were not incorporated into the development of the maps.

Open-File Report

Regional method to assess offshore slope stability.

The slope stability of some offshore environments can be evaluated by using only conventional acoustic profiling and short-core sampling, followed by laboratory consolidation and strength testing. The test results are synthesized by using normalized-parameter techniques. The normalized data are then used to calculate the critical earthquake acceleration factors or the wave heights needed to initiate failure. These process-related parameters provide a quantitative measure of the relative stability for locations from which short cores were obtained. The method is most applicable to offshore environments of gentle relief and simple subsurface structure and is not considered a substitute for subsequent site-specific analysis. -from ASCE Publications Information

Journal of Geotechnical Engineering

Nitrous oxide emissions from cropland: a procedure for calibrating the DayCent biogeochemical model using inverse modelling

DayCent is a biogeochemical model of intermediate complexity widely used to simulate greenhouse gases (GHG), soil organic carbon and nutrients in crop, grassland, forest and savannah ecosystems. Although this model has been applied to a wide range of ecosystems, it is still typically parameterized through a traditional “trial and error” approach and has not been calibrated using statistical inverse modelling (i.e. algorithmic parameter estimation). The aim of this study is to establish and demonstrate a procedure for calibration of DayCent to improve estimation of GHG emissions. We coupled DayCent with the parameter estimation (PEST) software for inverse modelling. The PEST software can be used for calibration through regularized inversion as well as model sensitivity and uncertainty analysis. The DayCent model was analysed and calibrated using N2O flux data collected over 2 years at the Iowa State University Agronomy and Agricultural Engineering Research Farms, Boone, IA. Crop year 2003 data were used for model calibration and 2004 data were used for validation. The optimization of DayCent model parameters using PEST significantly reduced model residuals relative to the default DayCent parameter values. Parameter estimation improved the model performance by reducing the sum of weighted squared residual difference between measured and modelled outputs by up to 67 %. For the calibration period, simulation with the default model parameter values underestimated mean daily N2O flux by 98 %. After parameter estimation, the model underestimated the mean daily fluxes by 35 %. During the validation period, the calibrated model reduced sum of weighted squared residuals by 20 % relative to the default simulation. Sensitivity analysis performed provides important insights into the model structure providing guidance for model improvement.

Water, Air, & Soil Pollution

Improving tsunami resiliency: California's Tsunami Policy Working Group

California has established a Tsunami Policy Working Group to facilitate development of policy recommendations for tsunami hazard mitigation. The Tsunami Policy Working Group brings together government and industry specialists from diverse fields including tsunami, seismic, and flood hazards, local and regional planning, structural engineering, natural hazard policy, and coastal engineering. The group is acting on findings from two parallel efforts: The USGS SAFRR Tsunami Scenario project, a comprehensive impact analysis of a large credible tsunami originating from an M 9.1 earthquake in the Aleutian Islands Subduction Zone striking California’s coastline, and the State’s Tsunami Preparedness and Hazard Mitigation Program. The unique dual-track approach provides a comprehensive assessment of vulnerability and risk within which the policy group can identify gaps and issues in current tsunami hazard mitigation and risk reduction, make recommendations that will help eliminate these impediments, and provide advice that will assist development and implementation of effective tsunami hazard risk communication products to improve community resiliency.

California

Megaporosity and permeability of Thalassinoides-dominated ichnofabrics in the Cretaceous karst-carbonate Edwards-Trinity aquifer system, Texas

Current research has demonstrated that trace fossils and their related ichnofabrics can have a critical impact on the fluid-flow properties of hydrocarbon reservoirs and groundwater aquifers. Most petroleum-associated research has used ichnofabrics to support the definition of depositional environments and reservoir quality, and has concentrated on siliciclastic reservoir characterization and, to a lesser degree, carbonate reservoir characterization (for example, Gerard and Bromley, 2008; Knaust, 2009). The use of ichnology in aquifer characterization has almost entirely been overlooked by the hydrologic community because the dynamic reservoir-characterization approach has not caught on with hydrologists and so hydrology is lagging behind reservoir engineering in this area (de Marsily and others, 2005). The objective of this research is to show that (1) ichnofabric analysis can offer a productive methodology for purposes of carbonate aquifer characterization, and (2) a clear relation can exist between ichnofabrics and groundwater flow in carbonate aquifers.

Texas

Some stream waters of the Western United States, with chapters on sediment carried by the Rio Grande and the industrial application of water analyses

A systematic study of the waters likely to be utilized on the Reclamation Service projects was made in order to determine the influence of the salinity of the waters on the growth of vegetation and the effect of suspended matter in silting canals and reservoirs. The work was begun early in 1905, under the direction of Thomas H. Means, engineer, and was continued during 1906 and until May, 1907, under the direction of W. H. Heileman, engineer. The analyses were made in a laboratory established at quarters provided by the University of California at Berkeley, Cal., by C. H. Stone, P. L. McCreary, F. M. Eaton, O. J. Hawley, W. C. Riddell, F. T. Berry, H. A. Burns, J. H. Hampson, J. A. Pearce, and M. Vaygouny, the greater part of the work being that of the first five named. C. H. Stone was chemist in charge at the beginning of the investigations and is chiefly responsible for the plan of the analytical work and the methods of analysis. The results of the investigations were prepared for publication under instructions from F. H. Newell, Director of the United States Reclamation Service, by Herman Stabler, assistant engineer, who assembled and checked- the analyses, compiled the accompanying stream-flow data from records of the United States Geological Survey, and computed daily discharge of suspended matter and dissolved solids, under the supervision of D. W. Murphy, engineer in charge of Washington office engineering.

Water Supply Paper

Liquefaction evidence for at least two strong Holocene paleo-earthquakes in central and southwestern Illinois, USA

Two strong mid-Holocene earthquakes in Illinois have been documented by paleoliquefaction features such as clastic dikes, sills, and detachments of fine-grained sediment that sunk into liquefied sand. At least one paleo-earthquake occurred in central Illinois about 35 km NE of Springfield, probably sometime between 5,900 and 7,400 yr BP. Dike widths are as much as 0.4 m near the energy center of the earthquake. Outward from this center, dike widths attenuate and ultimately disappear at about 35 km. More than one paleo-earthquake is probably represented by liquefaction features near Springfield. Another paleo-earthquake that appears to have been centered about 65 km ESE of St. Louis, Missouri, occurred near 5,700 yr BP. The energy center is inferred as being in Illinois, and most likely near lowermost Shoal Creek where the meizoseismal region is defined by dikes as wide as 0.5 m and by a regional abundance of dikes. Dikes from this earthquake probably extend at least as far as 35 km from its inferred energy center. The earthquake near Shoal Creek and one earthquake near Springfield almost certainly exceeded M 6. The paleomagnitudes can be more closely bracketed by geotechnical testing and analysis, when used in conjunction with existing data.

Illinois

The Virtual Data Center Tagged-Format Tool: Introduction and executive summary

This Report introduces and summarizes the new Virtual Data Center (VDC) Tagged Format (VTF) Tool, which was developed by a diverse group of seismologists, earthquake engineers, and information technology professionals for internal use by the COSMOS VDC and other interested parties for the exchange, archiving, and analysis of earthquake strong-ground-motion data.

Open-File Report

Analysis of regional rainfall-runoff parameters for the Lake Michigan Diversion hydrological modeling

The Lake Michigan Diversion Accounting (LMDA) system has been developed by the U.S. Army Corps of Engineers, Chicago District (USACE-Chicago) and the State of Illinois as a part of the interstate Great Lakes water regulatory program. The diverted Lake Michigan watershed is a 673-square-mile watershed that is comprised of the Chicago River and Calumet River watersheds. They originally drained into Lake Michigan, but now flow to the Mississippi River watershed via three canals constructed in the Chicago area in the early twentieth century. Approximately 393 square miles of the diverted watershed is ungaged, and the runoff from the ungaged portion of the diverted watershed has been estimated by the USACE-Chicago using the Hydrological Simulation Program-FORTRAN (HSPF) program. The accuracy of simulated runoff depends on the accuracy of the parameter set used in the HSPF program. Nine parameter sets comprised of the North Branch, Little Calumet, Des Plaines, Hickory Creek, CSSC, NIPC, 1999, CTE, and 2008 have been developed at different time periods and used by the USACE-Chicago. In this study, the U.S. Geological Survey and the USACE-Chicago collaboratively analyzed the parameter sets using nine gaged watersheds in or adjacent to the diverted watershed to assess the predictive accuracies of selected parameter sets. Six of the parameter sets, comprising North Branch, Hickory Creek, NIPC, 1999, CTE, and 2008, were applied to the nine gaged watersheds for evaluating their simulation accuracy from water years 1996 to 2011. The nine gaged watersheds were modeled by using the three LMDA land-cover types (grass, forest, and hydraulically connected imperviousness) based on the 2006 National Land Cover Database, and the latest meteorological and precipitation data consistent with the current (2014) LMDA modeling framework. Results indicate that the North Branch and Hickory Creek parameter sets, which belong to the original calibration group, attained an overall “satisfactory” rating on monthly runoff volumes based on the three performance statistics selected, but the annual and over-the-period runoff volumes were generally underestimated. Parameter sets CTE and 2008 attained a similar satisfactory rating on monthly runoff volumes but the annual and over-the-period runoff volumes were overestimated in general. Although the percent bias was improved, the CTE and 2008 parameter sets also had increased residuals in monthly runoff volumes and decreased quality of the model fit to the measured streamflows relative to the North Branch and Hickory Creek parameter sets. The NIPC and 1999 parameter sets, on the other hand, had larger percent bias and residuals in monthly runoff volumes, and underestimated the annual and over-the-period runoff volumes. Recalibration of the HSPF parameters to the updated inputs and land covers was completed on two representative watershed models selected from the nine by using a manual method (HSPEXP) and an automatic method (PEST). The objective of the recalibration was to develop a regional parameter set that improves the accuracy in runoff volume prediction for the nine study watersheds. Knowledge about flow and watershed characteristics plays a vital role for validating the calibration in both manual and automatic methods. The best performing parameter set was determined by the automatic calibration method on a two-watershed model. Applying this newly determined parameter set to the nine watersheds for runoff volume simulation resulted in “very good” ratings in five watersheds, an improvement as compared to “very good” ratings achieved for three watersheds by the North Branch parameter set.

Illinois, Indiana