Geology topics

Charles S. Melching

Publications and source records attributed to Charles S. Melching.

11 recordsLinked to original sources

USGS2000

Equations for estimating Clark unit-hydrograph parameters for small rural watersheds in Illinois

Equations for estimating the time of concentration ( T C ) and storage coefficient ( R ) of the Clark unit-hydrograph method were developed for small rural watersheds [0.02–2.3 square miles (mi 2 )] in Illinois. The equations will provide State and local engineers and planners with more accurate methods to estimate the T C and R for use in simulating discharge hydrographs on small rural watersheds when designing stormwater-management facilities and other hydraulic structures, determining flood-plain boundaries, and assessing the safety of structures in rivers. The rainfall and runoff data from gaged small rural watersheds (0.02–2.3 mi 2 ) with insignificant amounts of impervious land cover in Illinois were used to develop the equations. Equations were developed on the basis of data for 121 storms that occurred in 39 watersheds. Data for 29 storms in 18 watersheds were used to verify the equations. T C and R were determined by calibrating available rainfall and runoff data, using the U.S. Army Corps of Engineers Flood Hydrograph Package HEC-1. The mathematical relations between watershed and storm characteristics, and T C and R were determined by multiple-linear regression of the logarithms of the values. Main-channel length and slope were identified as important watershed characteristics for estimating T C and R . The estimation equations had coefficients of determination of 0.73 and 0.64 for the logarithms of T C and R , respectively. When storm characteristics were added in the regression of hydrograph parameters utilizing length and slope, only minimal increases to the coefficient of determination resulted. Thus, storm characteristics were not considered further in development of the equations. Simulation of the measured discharge hydrographs for the verification storms utilizing T C and R obtained from the estimation equations yielded good results. The error in peak discharge for 21 of the 29 verification storms was less than 25 percent, and the error in time-to-peak discharge for 18 of the 29 verification storms also was less than 25 percent. Therefore, applying the estimation equations to determine T C and R for design-storm simulation may result in reliable design hydrographs, as long as the physical characteristics of the watersheds under consideration are within the range of those characteristics for the watersheds in this study [area: 0.02–2.3 mi 2 , main-channel length: 0.17–3.4 miles, main-channel slope: 10.5–229 feet per mile, and insignificant percentage of impervious cover].

Illinois↗

USGS1998

Regional rainfall-runoff relations for simulation of streamflow for watersheds in Du Page County, Illinois

Rainfall and streamflow data collected from July 1986 through September 1993 were utilized to calibrate and verify a continuous-simulation rainfall-runoff model for three watersheds (11.8--18.0 square miles in area) in Du Page County. Classification of land cover into three categories of pervious (grassland, forest/wetland, and agricultural land) and one category of impervious subareas was sufficient to accurately simulate the rainfall-runoff relations for the three watersheds. Regional parameter sets were obtained by calibrating jointly all parameters except fraction of ground-water inflow that goes to inactive ground water (DEEPFR), interflow recession constant (IRC), and infiltration (INFILT) for runoff from all three watersheds. DEEPFR and IRC varied among the watersheds because of physical differences among the watersheds. Two values of INFILT were obtained: one representing the rainfall-runoff process on the silty and clayey soils on the uplands and lake plains that characterize Sawmill Creek, St. Joseph Creek, and eastern Du Page County; and one representing the rainfall-runoff process on the silty soils on uplands that characterize Kress Creek and parts of western Du Page County. Regional rainfall-runoff relations, defined through joint calibration of the rainfall-runoff model and verified for independent periods, presented in this report, allow estimation of runoff for watersheds in Du Page County with an error in the total water balance less than 4.0 percent; an average absolute error in the annual-flow estimates of 17.1 percent with the error rarely exceeding 25 percent for annual flows; and correlation coefficients and coefficients of model-fit efficiency for monthly flows of at least 87 and 76 percent, respectively. Close reproduction of the runoff-volume duration curves was obtained. A frequency analysis of storm-runoff volume indicates a tendency of the model to undersimulate large storms, which may result from underestimation of the amount of impervious land cover in the watershed and errors in measuring rainfall for convective storms. Overall, the results of regional calibration and verification of the rainfall-runoff model indicate the simulated rainfall-runoff relations are adequate for stormwater-management planning and design for watersheds in Du Page County.

Water-Resources Investigations Report↗

USGS1998

Intercomparison of principal hydrometric instruments; Third phase, Evaluation of ultrasonic velocity meters for flow measurement in streams, canals, and estuaries

As part of the World Meteorological Organization ( WMO ) project Intercomparison of Principal Hydrometric Instruments , Third Phase , a questionnaire was prepared by the U.S. Geological Survey (USGS) on the application of Ultrasonic Velocity Meters (UVM's) for flow measurement in streams , canals , and estuaries . In 1996, this questionnaire was distributed internationally by the WMO and USGS, and distributed within the United States by the USGS. Completed questionnaires were returned by 26 agencies in 7 countries (Canada, France, Germany, The Netherlands, Switzerland, the United Kingdom, and the United States). The completed questionnaires described geometric and streamflow conditions, system configurations, and reasons for applying UVM systems for 260 sites, thus providing information on the applicability of UVM systems throughout the world. The completed questionnaires also provided information on operational issues such as (1) methods used to determine and verify UVM ratings, (2) methods used to determine the mean flow velocity for UVM systems, (3) operational reliability of UVM systems, (4) methods to estimate missing data, (5) common problems with UVM systems and guidelines to mitigate these problems, and (6) personnel training issues. The completed questionnaires also described a few unique or novel applications of UVM systems. In addition to summarizing the completed questionnaires, this report includes a brief overview of UVM application and operation, and a short summary of current (1998) information from UVM system manufacturers regarding system cost and capabilities. On the basis of the information from the completed questionnaires and provided by the manufacturers, the general applicability of UVM systems is discussed. In the finalisation of this report the financial support provided by the US National Committee for Scientific Hydrology is gratefully acknowledged.

Report↗

USGS1997

Equations for estimating synthetic unit-hydrograph parameter values for small watersheds in Lake County, Illinois

Design hydrographs computed from design storms, simple models of abstractions (interception, depression storage, and infiltration), and synthetic unit hydrographs provide vital information for stormwater, flood-plain, and water-resources management throughout the United States. Rainfall and runoff data for small watersheds in Lake County collected between 1990 and 1995 were studied to develop equations for estimation of synthetic unit-hydrograph parameters on the basis of watershed and storm characteristics. The synthetic unit-hydrograph parameters of interest were the time of concentration (TC) and watershed-storage coefficient (R) for the Clark unit-hydrograph method, the unit-graph lag (UL) for the Soil Conservation Service (now known as the Natural Resources Conservation Service) dimensionless unit hydrograph, and the hydrograph-time lag (TL) for the linear-reservoir method for unit-hydrograph estimation. Data from 66 storms with effective-precipitation depths greater than 0.4 inches on 9 small watersheds (areas between 0.06 and 37 square miles (mi2)) were utilized to develop the estimation equations, and data from 11 storms on 8 of these watersheds were utilized to verify (test) the estimation equations. The synthetic unit-hydrograph parameters were determined by calibration using the U.S. Army Corps of Engineers Flood Hydrograph Package HEC-1 (TC, R, and UL) or by manual analysis of the rainfall and run-off data (TL). The relation between synthetic unit-hydrograph parameters, and watershed and storm characteristics was determined by multiple linear regression of the logarithms of the parameters and characteristics. Separate sets of equations were developed with watershed area and main channel length as the starting parameters. Percentage of impervious cover, main channel slope, and depth of effective precipitation also were identified as important characteristics for estimation of synthetic unit-hydrograph parameters. The estimation equations utilizing area had multiple correlation coefficients of 0.873, 0.961, 0.968, and 0.963 for TC, R, UL, and TL, respectively, and the estimation equations utilizing main channel length had multiple correlation coefficients of 0.845, 0.957, 0.961, and 0.963 for TC, R, UL, and TL, respectively. Simulation of the measured hydrographs for the verification storms utilizing TC and R obtained from the estimation equations yielded good results without calibration. The peak discharge for 8 of the 11 storms was estimated within 25 percent and the time-to-peak discharge for 10 of the 11 storms was estimated within 20 percent. Thus, application of the estimation equations to determine synthetic unit-hydrograph parameters for design-storm simulation may result in reliable design hydrographs; as long as the physical characteristics of the watersheds under consideration are within the range of those for the watersheds in this study (area: 0.06-37 mi2, main channel length: 0.33-16.6 miles, main channel slope: 3.13-55.3 feet per mile, and percentage of impervious cover: 7.32-40.6 percent). The estimation equations are most reliable when applied to watersheds with areas less than 25 mi2.

Open-File Report↗

USGS1997

Full Equations (FEQ) model for the solution of the full, dynamic equations of motion for one-dimensional unsteady flow in open channels and through control structures

The Full EQuations (FEQ) model is a computer program for solution of the full, dynamic equations of motion for one-dimensional unsteady flow in open channels and through control structures. A stream system that is simulated by application of FEQ is subdivided into stream reaches (branches), parts of the stream system for which complete information on flow and depth are not required (dummy branches), and level-pool reservoirs. These components are connected by special features; that is, hydraulic control structures, including junctions, bridges, culverts, dams, waterfalls, spillways, weirs, side weirs, and pumps. The principles of conservation of mass and conservation of momentum are used to calculate the flow and depth throughout the stream system resulting from known initial and boundary conditions by means of an implicit finite-difference approximation at fixed points (computational nodes). The hydraulic characteristics of (1) branches including top width, area, first moment of area with respect to the water surface, conveyance, and flux coefficients and (2) special features (relations between flow and headwater and (or) tail-water elevations, including the operation of variable-geometry structures) are stored in function tables calculated in the companion program, Full EQuations UTiLities (FEQUTL). Function tables containing other information used in unsteady-flow simulation (boundary conditions, tributary inflows or outflows, gate settings, correction factors, characteristics of dummy branches and level-pool reservoirs, and wind speed and direction) are prepared by the user as detailed in this report. In the iterative solution scheme for flow and depth throughout the stream system, an interpolation of the function tables corresponding to the computational nodes throughout the stream system is done in the model. FEQ can be applied in the simulation of a wide range of stream configurations (including loops), lateral-inflow conditions, and special features. The accuracy and convergence of the numerical routines in the model are demonstrated for the case of laboratory measurements of unsteady flow in a sewer pipe. Verification of the routines in the model for field data on the Fox River in northeastern Illinois also is briefly discussed. The basic principles of unsteady-flow modeling and the relation between steady flow and unsteady flow are presented. Assumptions and the limitations of the model also are presented. The schematization of the stream system and the conversion of the physical characteristics of the stream reaches and a wide range of special features into function tables for model applications are described. The modified dynamic-wave equation used in FEQ for unsteady flow in curvilinear channels with drag on minor hydraulic structures and channel constrictions determined from an equivalent energy slope is developed. The matrix equation relating flows and depths at computational nodes throughout the stream system by the continuity (conservation of mass) and modified dynamic-wave equations is illustrated for four sequential examples. The solution of the matrix equation by Newton's method is discussed. Finally, the input for FEQ and the error messages and warnings issued are presented.

Water-Resources Investigations Report↗

USGS1997

Full equations utilities (FEQUTL) model for the approximation of hydraulic characteristics of open channels and control structures during unsteady flow

The Full EQuations UTiLities (FEQUTL) model is a computer program for computation of tables that list the hydraulic characteristics of open channels and control structures as a function of upstream and downstream depths; these tables facilitate the simulation of unsteady flow in a stream system with the Full Equations (FEQ) model. Simulation of unsteady flow requires many iterations for each time period computed. Thus, computation of hydraulic characteristics during the simulations is impractical, and preparation of function tables and application of table look-up procedures facilitates simulation of unsteady flow. Three general types of function tables are computed: one-dimensional tables that relate hydraulic characteristics to upstream flow depth, two-dimensional tables that relate flow through control structures to upstream and downstream flow depth, and three-dimensional tables that relate flow through gated structures to upstream and downstream flow depth and gate setting. For open-channel reaches, six types of one-dimensional function tables contain different combinations of the top width of flow, area, first moment of area with respect to the water surface, conveyance, flux coefficients, and correction coefficients for channel curvilinearity. For hydraulic control structures, one type of one-dimensional function table contains relations between flow and upstream depth, and two types of two-dimensional function tables contain relations among flow and upstream and downstream flow depths. For hydraulic control structures with gates, a three-dimensional function table lists the system of two-dimensional tables that contain the relations among flow and upstream and downstream flow depths that correspond to different gate openings. Hydraulic control structures for which function tables containing flow relations are prepared in FEQUTL include expansions, contractions, bridges, culverts, embankments, weirs, closed conduits (circular, rectangular, and pipe-arch shapes), dam failures, floodways, and underflow gates (sluice and tainter gates). The theory for computation of the hydraulic characteristics is presented for open channels and for each hydraulic control structure. For the hydraulic control structures, the theory is developed from the results of experimental tests of flow through the structure for different upstream and downstream flow depths. These tests were done to describe flow hydraulics for a single, steady-flow design condition and, thus, do not provide complete information on flow transitions (for example, between free- and submerged-weir flow) that may result in simulation of unsteady flow. Therefore, new procedures are developed to approximate the hydraulics of flow transitions for culverts, embankments, weirs, and underflow gates.

Water-Resources Investigations Report↗

USGS1996

Simulation of water quality for Salt Creek in northeastern Illinois

Water-quality processes in the Salt Creek watershed in northeastern Illinois were simulated with a computer model. Selected waste-load scenarios for 7-day, 10-year low-flow conditions were simulated in the stream system. The model development involved the calibration of the U.S. Environmental Protection Agency QUAL2E model to water-quality constituent concentration data collected by the Illinois Environmental Protection Agency (IEPA) for a diel survey on August 29-30, 1995, and the verification of this model with water-quality constituent concentration data collected by the IEPA for a diel survey on June 27-28, 1995. In-stream measurements of sediment oxygen demand rates and carbonaceous biochemical oxygen demand (CBOD) decay rates by the IEPA and traveltime and reaeration-rate coefficients by the U.S. Geological Survey facilitated the development of a model for simulation of water quality in the Salt Creek watershed. In general, the verification of the calibrated model increased confidence in the utility of the model for water-quality planning in the Salt Creek watershed. However, the model was adjusted to better simulate constituent concentrations measured during the June 27-28, 1995, diel survey. Two versions of the QUAL2E model were utilized to simulate dissolved oxygen (DO) concentrations in the Salt Creek watershed for selected effluent discharge and concentration scenarios for water-quality planning: (1) the QUAL2E model calibrated to the August 29-30, 1995, diel survey, and (2) the QUAL2E model adjusted to the June 27-28, 1995, diel survey. The results of these simulations indicated that the QUAL2E model adjusted to the June 27-28, 1995, diel survey simulates reliable information for water-quality planning. The results of these simulations also indicated that to maintain DO concentrations greater than 5 milligrams per liter (mg/L) throughout most of Salt Creek for 7-day, 10-year low-flow conditions, the sewage-treatment plants (STP's) must discharge effluent with CBOD and total ammonia as nitrogen concentrations substantially below the permit limits. If the STP's discharge effluent with CBOD and total ammonia as nitrogen concentrations at the permit limits for 7-day, 10-year low-flow conditions, DO concentrations less than 5 mg/L are expected for all of Salt Creek downstream from Fullerton Avenue (river mile 23.1).

Open-File Report↗

USGS1995

Relation of precipitation quality to storm type, and deposition of dissolved chemical constituents from precipitation in Massachusetts, 1983-85

Precipitation samples were collected for 83 storms at a rural inland site in Princeton, Mass., and 73 storms at a rural coastal site in Truro, Mass., to examine the quality of precipitation from storms and relate quality to three storm types (oceanic cyclone, continental cyclone, and cold front). At the inland site, Princeton, ranked-means of precipitation depth, storm duration, specific conductance, and concentrations and loads of hydrogen, sulfate, aluminum, bromide, and copper ions were affected by storm type. At the coastal site, Truro, ranked means of precipitation depth, storm duration, and concentrations and loads of calcium, chloride, magnesium, potassium, and sodium ions were affected by storm type. Precipitation chemistry at the coastal site was 85 percent oceanic in orgin, whereas precipitation 72 kilometers inland was 60 percent hydrogen, nitrate, and sulfate ions, reflecting fossil-fuel combustion. Concentrations and loads for specific conductance and 9 chemical constituents on an annual and seasonal basis were determined from National Atmospheric Deposition Program data for spring 1983 through winter 1985 at Quabbin (rural, inland), Waltham (suburban, inland) and Truro (rural, coastal), Massachusetts. Concentrations of magnesium, potassium, sodium, and chloride concentrations were highest at the coast and much lower inland, with very little difference between Waltham and Quabbin. Loads of ammonium, nitrate, sulfate, and hydrogen are highest at Quabbin and are about equal at Waltham and Truro. About twice as much nitrate and hydrogen and about 35 percent more sulfate is deposited at Quabbin than at Waltham or Truro; this pattern indicates that the interior of Massachusetts receives more acidic precipitation than do the eastern or the coastal areas of Massachusetts.

Water-Resources Investigations Report↗

USGS1995

Regional rainfall-runoff relations for simulation of streamflow for watersheds in Lake County, Illinois

Rainfall and streamflow data collected in Lake County, Ill., from March 1990 through September 1993 were used to (1) calibrate a rainfall-runoff model for an area encompassing three watersheds (individual areas of 17.2, 35.7, and 37.0 mi 2 (square miles) and (2) verify the regional model parameter set obtained from the calibration by applying the parameter set to rainfall-runoff models for an additional small (6.3 mi 2 ) watershed and a large (59.6 mi 2 ) watershed. In addition, rainfall and streamflow data collected from April 1991 through September 1993 were used to calibrate the rainfall-runoff model for three single land-use watersheds (38.2-305 acres), called hydrologic response units (HRU's). Significant differences were found between the best parameters used in the HRU models and in the larger watershed models. The main channels in the HRU's are intermittent streams; thus, the parameters in the HRU models were selected such that a fluctuating water table could be simulated; runoff from the larger watersheds is not as sensitive to the effects of a fluctuating water table. Classification of land cover into two pervious subareas (forest and grass) and one impervious subarea (including parking lots, streets, and rooftops, among others) was sufficient to simulate the rainfall-runoff relations for all watersheds accurately. The model parameters presented in this report, which were refined through regional calibration and verified for watersheds not considered in the calibration, allow simulation of runoff in watersheds in Lake County, Ill., with approximately 93-percent accuracy in the total water balance, an average absolute error in the annual- flow estimates of 10.9 percent (and an error rarely exceeding 25 percent for annual flow), and monthly water balances with correlation coefficients of 93 percent and coefficients of model-fit efficiency of 86 percent. The models closely reproduced the partial-duration series of runoff and storm-runoff frequencies for the modeled watersheds.

Illinois↗

USGS1993

Comparison, analysis, and estimation of discharge data from two acoustic velocity meters on the Chicago Sanitary and Ship Canal at Romeoville, Illinois

The acoustic velocity meter (AVM) on the Chicago Sanitary and Ship Canal (the Canal) at Romeoville, Ill., provides vital information for the accounting of the diversion of water from Lake Michigan. A detailed analysis of the discharge record on the Canal at Romeoville was done by the U.S. Geological Survey to establish the most accurate estimates of discharge for water years 1986-91. The analysis involved (1) checking the consistency of the discharges estimated by two different AVM's installed at Romeoville for consecutive time periods by statistical and regression analyses, (2) adjusting the discharge record to account for corrections to the width and depth of the Canal determined by field measurements, and (3) development of equations for estimating discharge on days when the AVM was inoperative using discharge estimates made by the Metropolitan Water Reclamation District of Greater Chicago at the lock, powerhouse, and controlling works at Lockport, Ill. No signi- ficant difference in the discharge estimates made by the two AVM's could be documented. The estimation equations combined regression analysis with physical principles of the outlet-works operation. The estimation equations simulated the verification period of October 1, 1991, to May 31, 1992, within 0.22, 5.15, and 0.66 percent for the mean, standard deviation, and skewness coefficient, respectively. Discharges were recalculated for the corrected width and depth, estimated for the periods of AVM inoperation, and entered into the discharge record for the station.

Water-Resources Investigations Report↗

USGS1992

Improved first-order uncertainty method for water-quality modeling

Uncertainties are unavoidable in water-quality modeling and subsequent management decisions. Monte Carlo simulation and first-order uncertainty analysis (involving linearization at central values of the uncertain variables) have been frequently used to estimate probability distributions for water-quality model output due to their simplicity. Each method has its drawbacks: Monte Carlo simulation's is mainly computational time; and first-order analysis are mainly questions of accuracy and representativeness, especially for nonlinear systems and extreme conditions. An improved (advanced) first-order method is presented, where the linearization point varies to match the output level whose exceedance probability is sought. The advanced first-order method is tested on the Streeter-Phelps equation to estimate the probability distribution of critical dissolved-oxygen deficit and critical dissolved oxygen using two hypothetical examples from the literature. The advanced first-order method provides a close approximation of the exceedance probability for the Streeter-Phelps model output estimated by Monte Carlo simulation using less computer time - by two orders of magnitude - regardless of the probability distributions assumed for the uncertain model parameters.

Journal of Environmental Engineering↗