Geology ReportsSearch

Geology topics

Research about Illinois, Indiana

Source-linked reports with geographic coverage including Illinois, Indiana.

At least 19 recordsLinked to original sources

Dynamic estimates of geologic CO2 storage resources in the Illinois Basin constrained by reinjectivity of brine extracted for pressure management

Geologic carbon storage (GCS) is likely to be an important part of global efforts to decarbonize the energy industry. Widespread deployment of GCS relies on strategies to maximize CO 2 injection rates while minimizing reservoir pressurization that could induce seismicity and/or fluid leakage into groundwater resources. Brine extraction from CO 2 storage formations with subsurface reinjection elsewhere could mitigate pressure buildup associated with GCS. Therefore, evaluation of CO 2 storage resources should consider the injectivity of produced brine in geologic layers above or below the CO 2 storage formation. For this study, a methodology was developed to estimate brine injectivity from formation depth and thickness using flow modeling and optimization techniques. The methodology was demonstrated in the Illinois Basin, where GCS in the Mt. Simon Sandstone is ongoing. Based on pressure constraints and considering only regions of the shallower units where salinity and sealing conditions were met, maximum brine injection rates were estimated within the Mt. Simon and three overlying hydrostratigraphic layers. Results indicate that a large area exists where CO 2 injectivity could be optimized by brine extraction and reinjection.

Illinois, Indiana

Verification of multiple phosphorus analyzers for use in surface-water applications

The U.S. Geological Survey (USGS) completed a verification study of selected commercially available phosphorus analyzers for their applicability to scientific surface-water applications. In this study, the analyzers were the Hach EZ7800 TOPHO, Hach Phosphax sc, Sea-Bird Scientific HydroCycle-PO 4 , and the YSI Inc. Alyza IQ PO4. Verification tests included laboratory trials comparing analyzer results to known standards with several known concentrations of dissolved organic matter and waste production estimates. Field trials were completed at the Vermilion River near Danville, Illinois (U.S. Geological Survey station 03339000), where analyzer-measured concentrations were compared against discrete samples across a wide range of environmental conditions from November 2020 to August 2021. Data coverage was closely tracked for analyzer malfunctions and operator errors that caused missing data. Laboratory and field trials indicated that each analyzer is a viable option for scientific surface-water studies depending on environmental conditions. Because of the complexity of the analyzers, a substantial time investiture was required to get maximum data coverage including considerable site infrastructure investments and well-trained technicians. Data coverage was closely related to each analyzer’s ability to handle elevated turbidity levels.

Illinois, Indiana

Effect of uncertainty of discharge data on uncertainty of discharge simulation for the Lake Michigan Diversion, northeastern Illinois and northwestern Indiana

Simulation models of watershed hydrology (also referred to as “rainfall-runoff models”) are calibrated to the best available streamflow data, which are typically published discharge time series at the outlet of the watershed. Even after calibration, the model generally cannot replicate the published discharges because of simplifications of the physical system embedded in the model structure and uncertainties of the input data and of the estimated model parameters, which, although optimized for the given calibration data, remain uncertain. The input data errors are caused by uncertainties in the forcing data, such as precipitation and other climatological data, and in the published discharges used for calibration. In the numerical algorithms used for calibration, the published discharges are often assumed to be without error, but they are themselves uncertain, typically having been computed using ratings, which are models fitted to uncertain discharge measurements. In this study, uncertainty of published daily discharge data and how the discharge uncertainty is transmitted to the parameter values of the Hydrological Simulation Program–FORTRAN (HSPF) rainfall-runoff model and to the simulated discharge at both calibration and prediction locations were investigated for the Lake Michigan diversion in northeastern Illinois and northwestern Indiana. The HSPF model used in this study is used by the U.S. Army Corps of Engineers as part of quantifying the diversion of water from Lake Michigan by the State of Illinois. In this study, the model is calibrated jointly at two watersheds in the study area; the resulting model is considered the base model in this study. Seven other gaged watersheds in the study area are used for testing predictive simulations. A Bayesian rating curve estimation (BaRatin) approach, the BaRatin stage-period-discharge (SPD) method, was used to estimate the uncertainty of the published discharge from the calibration watersheds. To characterize the effect of the discharge uncertainty on parameter values, the HSPF model parameters were recalibrated to 17 nonrandomly selected pairs of discharge series from the BaRatin SPD analysis. To provide an indicator of the effect of parameter uncertainty to compare to the effect of discharge uncertainty, 1,000 parameter sets also were randomly generated from the estimated parameter covariance matrix of the base model. The recalibrated and random parameter sets were then used in HSPF simulations of discharge at the two calibration watersheds and at the seven prediction watersheds. Selected discharge summary statistics—the period-of-study (POS, water years 1997 to 2015) mean discharge, selected flow-duration curve (FDC) quantiles, and water year mean discharges—are used to characterize the variability between simulated and published discharge. A normalized variability index ( V N ) is used as a measure of the uncertainty of flow statistics arising from the uncertainty of the sources considered in this study. When this index is at least 1, the variability of the simulations is large enough to explain the median error between simulated and published values, although offsetting errors from other sources are also likely. When the index is appreciably less than 1, the variability of the simulations is clearly insufficient to explain the median error between simulated and published values. At the two calibration watersheds and for results of the two simulation sets considered together, the V N values ranged from 0.2 to 0.8 for POS mean discharge, from 0.3 to 0.6 in the median for a set of FDC quantiles, and from 0.1 to 0.2 in the median for water year mean discharges. These values indicate that substantial uncertainty remains unexplained. Even though two watersheds were used in calibration, that calibration was highly constrained because it was applied to the watersheds simultaneously and was subject to parameter regularization that constrained the adjustment of the parameters from their initial values. These constraints were applied to avoid overfitting to the calibration watersheds and thus to increase the likelihood that the resulting parameters would give accurate results at watersheds not used in the calibration, but they created a parameter transfer error in the calibration watershed results shown by the balancing of errors between the two watersheds. Additional remaining error sources include model structural error and meteorological forcing error to the degree that the calibration was unable to adjust the parameters to account for these errors. At the prediction watersheds, the corresponding V N values were almost always substantially lower than those values at the calibration watersheds. This result is expected because the prediction watersheds have additional uncertainty, including parameter transfer error. The work described in this report provides preliminary estimates of a limited range of sources of error in predicted discharge uncertainty. Future work would be beneficial to obtain a better statistical characterization of the effect of the uncertainty of calibration discharge series and to address additional sources of uncertainty, such as from precipitation input data used in calibration and prediction and from structural (model) errors.

Illinois, Indiana

Status of the major aquaculture carps of China in the Laurentian Great Lakes Basin

There is concern of economic and environmental damage occuring if any of the four major aquacultured carp species of China, black carp Mylopharyngodon piceus , bighead carp Hypophthalmichthys nobilis , silver carp H. molitrix , or grass carp Ctenopharyngodon idella , were to establish in the Laurentian Great Lakes. All four are reproducing in the Mississippi River Basin. We review the status of these fishes in relation to the Great Lakes and their proximity to pathways into the Great Lakes, based on captures and collections of eggs and larvae. No black carp have been captured in the Great Lakes Basin. One silver carp and one bighead carp were captured within the Chicago Area Waterway System, on the Great Lakes side of electric barriers designed to keep carp from entering the Great Lakes from the greater Mississippi River Basin. Three bighead carp were captured in Lake Erie, none later than the year 2000. By December 2019, at least 650 grass carps had been captured in the Great Lakes Basin, most in western Lake Erie, but none in Lake Superior. Grass carp reproduction has been documented in the Sandusky and Maumee rivers in Ohio, tributaries of Lake Erie. We also discuss environmental DNA (eDNA) results as an early detection and monitoring tool for bighead and silver carps. Detection of eDNA does not necessarily indicate presence of live fish, but bigheaded carp eDNA has been detected on the Great Lakes side of the barriers and in a small proportion of samples from the western basin of Lake Erie.

Illinois, Indiana

Hydrous pyrolysis of New Albany Shale: A study examining maturation changes and porosity development

The characterization of nanoscale organic structures has improved our understanding of porosity development within source-rock reservoirs, but research linking organic porosity evolution to thermal maturity has generated conflicting results. To better understand this connection, an immature (0.25% solid bitumen reflectance; BR o ) sample of the New Albany Shale was used in four isothermal hydrous pyrolysis (HP) experiment sequences at 300°, 320°, 340°, and 370°C, with residues collected periodically for a maximum of 103 days. The HP residues, along with the original immature sample and two naturally matured (1.49 and 1.56% BR o ) New Albany Shale samples were analyzed for organic petrology, total organic carbon (TOC) content, and organic porosity evaluation using correlative light and electron microscopy (CLEM). All of the HP series increased in thermal maturity with increasing duration of pyrolysis, though reflectance for each series plateaued within 25 days of maturation. Initially, TOC in the HP residues decreases (from 14.24 wt. %) with increasing thermal maturity until ∼1.0% BR o where TOC remains at ∼9–10 wt. % for all remaining residues. Qualitative CLEM observations within the 50–100 day 300° and 340°C HP sequences (0.95–1.70% BR o ), and the naturally matured samples, develop organic porosity in smaller (<5 μm in diameter), void-filling solid bitumen that occurs in spaces between clays and other fine-grained minerals. The 370°C HP residues developed significant organic porosity, relative to the other HP temperature series in all solid bitumen accumulations regardless of size. Overall, the study indicates that temperature and duration of artificial maturation play an important role in the abundance of pores in the HP residues. This work expands on our understanding of the conditions needed for the generation and development of organic porosity in the New Albany Shale and potentially in other marine source-rock petroleum systems.

Illinois, Indiana

In vitro effects-based method and water quality screening model for use in pre- and post-distribution treated waters

Recent urban public water supply contamination events emphasize the importance of screening treated drinking water quality after distribution. In vitro bioassays, when run concurrently with analytical chemistry methods, are effective tools to evaluating the efficacy of water treatment processes and water quality. We tested 49 water samples representing the Chicago Department of Water Management service areas for estrogen, (anti)androgen, glucocorticoid receptor-activating contaminants and cytotoxicity. We present a tiered screening approach suitable to samples with anticipated low-level activity and initially tested all extracts for statistically identifiable endocrine activity; performing a secondary dilution-response analysis to determine sample EC 50 and biological equivalency values (BioEq). Estrogenic activity was detected in untreated Lake Michigan intake water samples using mammalian (5/49; median: 0.21 ng E2Eq/L) and yeast cell (5/49; 1.78 ng E2Eq/L) bioassays. A highly sensitive (anti)androgenic activity bioassay was applied for the first time to water quality screening and androgenic activity was detected in untreated intake and treated pre-distribution samples (4/49; 0.93 ng DHTEq/L). No activity was identified above method detection limits in the yeast androgenic, mammalian anti-androgenic, and both glucocorticoid bioassays. Known estrogen receptor agonists were detected using HPLC/MS-MS (estrone: 0.72-1.4 ng/L; 17α-estradiol: 1.3-1.5 ng/L; 17β-estradiol: 1.4 ng/L; equol: 8.8 ng/L), however occurrence did not correlate with estrogenic bioassay results. Many studies have applied bioassays to water quality monitoring using only relatively small samples sets often collected from surface and/or wastewater effluent. However, to realistically adapt these tools to treated water quality monitoring, water quality managers must have the capacity to screen potentially hundreds of samples in short timeframes. Therefore, we provided a tiered screening model that increased sample screening speed, without sacrificing statistical stringency, and detected estrogenic and androgenic activity only in pre-distribution Chicago area samples.

Illinois, Indiana

Estimating the net costs of brine production and disposal to expand pressure-limited dynamic capacity for basin-scale CO2 storage in a saline formation

If carbon capture and storage (CCS) needs to be deployed at basin- or larger-scale, it is likely that multiple sites will be injecting carbon dioxide (CO 2 ) into the same geologic formation. This could lead to excessive pressure buildup, overlapping induced pressure fronts, and pressure interference with neighboring uses of the subsurface. Extracting the in situ brine from the storage formation could be necessary to relieve pressure constraints; control migration of the CO 2 plume, displaced brine, and the induced pressure front; and sequester more CO 2 while reducing potential risks. Such active pressure management could be very costly, and it could present a formidable economic constraint on the feasible scale of deployment of CCS. Alternatively, there may be high-injectivity zones (“storage sweet spots”) where a significant volume of CO 2 could be stored without producing brine. For simulated deployment of CO 2 storage sites across the Illinois Basin, the results of this study suggest that brine production could be required to sequester 20 % or more of the regional CO 2 emissions of major stationary sources in the Mount Simon Sandstone saline formation. In some cases, brine production could expand pressure-limited CO 2 storage capacity enough to more than compensate for the additional costs of pressure management, but only if produced brine could be cheaply reinjected onsite for disposal in an overlying geologic formation. With or without brine production, this study found that the lowest-cost deployment option was to inject CO 2 only into a potential storage sweet spot of the Mount Simon Sandstone.

Illinois, Indiana

Methane and nitrous oxide temporal and spatial variability in two midwestern USA streams containing high nitrate concentrations

Concentrations and emissions of greenhouse gases CO 2 , CH 4 , and N 2 O commonly are examined individually in aquatic environments in which each is expected to be relatively important; however, their co-occurrence and dynamic interactions in fluvial settings could provide important information about their controlling biogeochemical processes and potential contributions to global climate change. Spatial and temporal variability of CH 4 , N 2 O, and CO 2 concentrations were measured from June 1999 to September 2003 in two nitrate-rich (40–1200 μM) streams draining agricultural land in the midwestern USA that differed ~13-fold in flow. Seasonal (biweekly), diel (hourly), and transport-oriented (reach-scale) sampling approaches were compared. Dissolved gas concentrations exceeded atmospheric equilibrium values up to 700- and 16-fold, for CH 4 and N 2 O, respectively. Mean concentrations were higher in the larger stream than in the smaller stream. In both streams, CH 4 emissions were generally higher in summer-fall and negatively correlated with flow and NO 3 − concentration while N 2 O emissions were generally higher in winter/spring and positively correlated with flow and NO 3 − . In the small stream, diel variations in the concentrations, emissions, and isotopic compositions of CH 4 , N 2 O, and NO 2 − resulted from diel variations in sources, sinks, and air-water gas exchange velocities. Seasonal mean total (CH 4 + N 2 O) area-normalized emission rates, expressed as CO 2 warming potential equivalents, were similar for the two streams, but the total reach-scale emission rate for the larger stream, including CO 2 , was about 2.9 times that of the smaller stream (131.6 vs 46.0 kg CO 2 equivalents km −1 day −1 , respectively). The CH 4 contribution to this flux was 9–28%, despite the relatively high NO 3 − and O 2 concentrations in the streams, indicating contributions from upwelling groundwater or reactions in streambed sediment.

Illinois, Indiana

Concentrations of lead and other inorganic constituents in samples of raw intake and treated drinking water from the municipal water filtration plant and residential tapwater in Chicago, Illinois, and East Chicago, Indiana, July–December 2017

The U.S. Geological Survey (USGS) Environmental Health Mission Area (EHMA) is providing comprehensive science on sources, movement, and transformation of contaminants and pathogens in watershed and aquifer drinking-water supplies and in built water and wastewater infrastructure (referred to as the USGS Water and Wastewater Infrastructure project) in the Greater Chicago Area and elsewhere in the United States, to fill data gaps identified by stakeholders and collaborators in drinking water and public health. EHMA Water and Wastewater Infrastructure research specifically provides insight into natural factors in the environment as well as those water-infrastructure components and processes (such as source-water corrosivity, treatment, plumbing, and so forth) that might influence human exposure to chemical and microbial contaminants at the residential tap. This infrastructure-exposure research role is fulfilled uniquely by the USGS and not by the U.S. Environmental Protection Agency (EPA), other agencies, or municipalities that focus on regulatory and policy activities and related compliance. The USGS approach to assessing the possible links between human health and chemical contaminant and pathogen exposure in drinking water is conducted in collaboration with public health experts and includes comprehensive characterization of the presence/absence and concentrations of more than 500 organic and 27 inorganic chemical constituents at the point of use (tap). Laboratory results for lead and other inorganic contaminants in Chicago, Illinois, and East Chicago, Indiana, residential tapwater are being released to ensure the timely release of quality-assured data to participants in the study. Concentrations of lead and other inorganic constituents were assessed in drinking water at the point of use (kitchen tap or filter) in 45 residential locations and in two locations within each of the two Chicago water purification plants and the two East Chicago water filtration plants during July–December 2017. Three methods were used for analyzing lead. The most sensitive method had a reporting limit of 0.020 micrograms per liter (µg/L). When using the most sensitive analytical method, lead was detected in 39 of 45 residential tapwater samples, with concentrations ranging from less than 0.020 µg/L to 5.31 µg/L (median of the detected values = 0.481 µg/L). Concentrations of lead also were detected in Lake Michigan intake water at all water purification/filtration plant facilities at concentrations ranging from 0.083 to 0.330 µg/L, but were not detected above the reporting limit in any samples of treated, pre-distribution drinking water at any of the water purification/filtration plant facilities. Because the USGS Water and Wastewater Infrastructure project in the Greater Chicago Area is focused on the potential human exposure to a broad suite of organic and inorganic contaminants in drinking water and is not focused specifically on lead, the sampling protocol did not include “first-draw,” stagnant sampling and samples were collected with point-of-use treatment in place, if present. Thus, the lead results reported herein are not appropriate for assessment of compliance with the EPA 1991 Lead and Copper Rule. Information resources for lead mitigation and water filtration are provided.

Illinois, Indiana

Comparison of benthos and plankton for Waukegan Harbor Area of Concern, Illinois, and Burns Harbor-Port of Indiana non-Area of Concern, Indiana, in 2015

During two seasonal sampling events in spring (June) and fall (August) of 2015, the U.S. Geological Survey collected benthos (benthic invertebrates) and plankton (zooplankton and phytoplankton) at three sites each in the Waukegan Harbor Area of Concern (AOC) in Illinois and in Burns Harbor-Port of Indiana, a non-AOC comparison site in Indiana. The study was done in cooperation with the U.S. Environmental Protection Agency and the Illinois Department of Natural Resources. Samples were collected concurrently for physical and chemical parameters (specific conductance, temperature, pH, dissolved oxygen, chlorophyll- a , total and volatile suspended solids in water samples; particle size and volatile-on-ignition solids of sediment in dredge samples). The purpose of the study was to assess whether or not aquatic communities at the AOC were degraded in comparison to communities at the non-AOC, which was presumed to be less impaired than the AOC. Benthos were collected by using Hester-Dendy artificial substrate samplers and a Ponar® dredge sampler to collect composited grabs of bottom sediment; zooplankton were collected by using tows from depth to the surface with a 63-micrometer mesh plankton net; phytoplankton were collected by using whole water samples composited from set depth intervals. Aquatic communities at the AOC and the non-AOC were compared by use of univariate statistical analyses with metrics such as taxa richness (number of unique taxa), diversity, and a multimetric Index of Biotic Integrity (IBI, for artificial-substrate samples only) as well as by use of multivariate statistical analyses of taxa relative abundances. Although benthos communities at Waukegan Harbor AOC were not rated as degraded in comparison to the non-AOC, metrics for zooplankton and phytoplankton communities did show some impairment for the 2015 sampling. Across seasons, benthos richness and diversity were significantly higher and rated as less degraded at the AOC compared to the non-AOC; however, benthos IBIs were not significantly different. Multivariate comparisons revealed that the benthos communities in the AOC and non-AOC were significantly different, but these comparisons do not address current degradation in either harbor. The dominant taxa in dredge samples were oligochaete worms in both harbors, but there were differences in the relative abundances of Dreissena as well as oligochaete and midge taxa. Although zooplankton richness and diversity in the AOC were lower and rated as more degraded in spring, these metrics were rated as less degraded in fall compared to the non-AOC, effectively balancing out so that there was no difference across seasons. Multivariate comparisons also indicated that zooplankton communities in the AOC were significantly different from those in the non-AOC for spring only but not across seasons, possibly because of lower water temperatures in spring at Waukegan Harbor than at the non-AOC site. The spring zooplankton community in Waukegan Harbor was dominated in density and biomass by the rotifer Synchaeta . Across seasons, diatom richness was significantly higher and rated as less degraded in the AOC than the non-AOC because of spring values, whereas soft algae richness was significantly lower and rated as more degraded in the AOC because of fall values. Spring richness of combined phytoplankton (soft algae and diatoms) was significantly higher in the AOC than in the non-AOC. Neither diatom diversity nor soft algae diversity differed significantly between the harbors, but combined phytoplankton diversity across seasons, if replicates were included, was significantly lower and rated as more degraded in the AOC than in the non-AOC. Multivariate tests indicated that the combined phytoplankton communities in the harbors were not significantly different across seasons. Significant differences were not found between harbors for chlorophyll- a , suspended solids, algal densities, or biomass.

Illinois, Indiana

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 &ldquo;satisfactory&rdquo; 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 &ldquo;very good&rdquo; ratings in five watersheds, an improvement as compared to &ldquo;very good&rdquo; ratings achieved for three watersheds by the North Branch parameter set.

Illinois, Indiana

Estimation of regional flow-duration curves for Indiana and Illinois

Flow-duration curves (FDCs) of daily streamflow are useful for many applications in water resources planning and management but must be estimated at ungaged sites. One common technique for estimating FDCs at ungaged sites in a given region is to use equations obtained by linear regression of FDC quantiles against multiple basin characteristics that can be computed by means of a geographic information system (GIS) computer program. In this study, such regional regression equations for estimating FDC quantiles were computed at the 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 25, 30, 40, 50, 60, 70, 75, 80, 90, 95, 98, 99, 99.5, 99.8, and 99.9-percent exceedance probabilities for rural, unregulated streams in Indiana and Illinois with temporally stationary records, using data through September 30, 2007. The approach used accounts for censored values below 0.01 cubic feet per second, which are observed at exceedance probabilities as low as 70 percent (that is, occurring at least 30 percent of the time). The basin characteristics used are suitable for computation by the USGS Web-based application, StreamStats, and are available for all U.S. Environmental Protection Agency (EPA) Region V states and the larger Great Lakes area, with some specific local exceptions. Indiana and Illinois were each divided into three regions, and a different set of equations for estimating FDC quantiles was computed for each region. The error of estimation of the FDC quantiles, measured as the mean square residual in log space converted to a percentage of the quantile, varies somewhat among regions and varies strongly with exceedance probability, with a minimum error of 10 to 20 percent at an exceedance probability of 5 or 10 percent, but rises to 17 to 38 percent at the high-flow end of the FDCs (the 0.1-percent quantile) and 100 to 745 percent at the low-flow end. For comparison, errors of estimation also were computed for FDC quantiles estimated by linear regression on drainage area alone and by using the drainage-area ratio (DAR) method. Three criteria, the nearest basin centroid and two others termed “strict” and “broad”, were used to select index stations for the DAR method. The “strict” and “broad” criteria put conditions on the basin centroid distance and the range of their drainage-area ratios, and the errors were averaged for all index station pairs satisfying each criterion. The use of the simpler DAR method usually resulted in higher errors of estimation compared to the linear regression equations with multiple basin characteristics, except occasionally in the case of the DAR method with the strict index station selection criterion, a criterion that is rarely possible to satisfy in practice. An example application of the estimated equations to one gaged and a few ungaged locations in a watershed in the study area is included to illustrate the steps required. These steps are the computation of the basin characteristics and, using those characteristics together with the estimated equations, the computation of the FDC quantiles and their uncertainties.

Illinois, Indiana

Tracking change over time: River flooding

Landsat satellites have been capturing images of Earth from space since 1972. These images provide a long-term record of natural and human-induced changes on the global landscape. Comparing images from multiple years reveals slow and subtle changes as well as rapid and devastating ones. Landsat images are available from the Internet at no charge. Using the free software MultiSpec, students can track changes to the landscape over time&mdash;just like remote sensing scientists do! The objective of the Tracking Change Over Time lesson plan is to get students excited about studying the changing Earth. Intended for students in grades 5-8, the lesson plan is flexible and may be used as a student self-guided tutorial or as a teacher-led class lesson. Enhance students' learning of geography, map reading, earth science, and problem solving by seeing landscape changes from space. The objective of the Tracking Change Over Time lesson plan is to get students excited about studying the changing Earth. Intended for students in grades 5-8, the lesson plan is flexible and may be used as a student self-guided tutorial or as a teacher-led class lesson. Enhance students' learning of geography, map reading, earth science, and problem solving by seeing landscape changes from space.

Illinois, Indiana

Revised magnitude-bound relation for the Wabash Valley seismic zone of the central United States

Seismic hazard assessment in the central United States, and in particular the Wabash Valley seismic zone of Indiana-Illinois, frequently relies on empirical estimates of paleoearthquake magnitudes ( M ). In large part these estimates have been made using the magnitude-bound method. Existing region-specific magnitude-bound relations rely heavily on only a few historical earthquakes in the central United States and eastern Canada that induced reported liquefaction features. Recent seismological studies have suggested smaller magnitudes than previously presumed for some of these earthquakes, however, and other studies have reinterpreted site-to-source distances to liquefaction features associated with some of these earthquakes. In this paper, we re-examine historical earthquakes ( M > ∼5) that occurred in the central and eastern United States and eastern Canada; some of these earthquakes triggered liquefaction and others did not. Based on our findings, we reinterpret the region-specific magnitude-bound relation for the Wabash Valley. Using this revised magnitude-bound relation, we present magnitude estimates for four prehistoric earthquakes that occurred in the Wabash Valley seismic zone during Holocene time.

Illinois, Indiana

The New Albany Shale petroleum system, Illinois basin - Data and map image archive from the material-balance assessment

The data files and explanations presented in this report were used to generate published material-balance approach estimates of amounts of petroleum 1) expelled from a source rock, and the sum of 2) petroleum discovered in-place plus that lost due to 3) secondary migration within, or leakage or erosion from a petroleum system. This study includes assessment of cumulative production, known petroleum volume, and original oil in place for hydrocarbons that were generated from the New Albany Shale source rocks. More than 4.00 billion barrels of oil (BBO) have been produced from Pennsylvanian-, Mississippian-, Devonian-, and Silurian-age reservoirs in the New Albany Shale petroleum system. Known petroleum volume is 4.16 BBO; the average recovery factor is 103.9% of the current cumulative production. Known petroleum volume of oil is 36.22% of the total original oil in place of 11.45 BBO. More than 140.4 BBO have been generated from the Upper Devonian and Lower Mississippian New Albany Shale in the Illinois Basin. Approximately 86.29 billion barrels of oil that was trapped south of the Cottage Grove fault system were lost by erosion of reservoir intervals. The remaining 54.15 BBO are 21% of the hydrocarbons that were generated in the basin and are accounted for using production data. Included in this publication are 2D maps that show the distribution of production for different formations versus the Rock-Eval pyrolysis hydrogen-indices (HI) contours, and 3D images that show the close association between burial depth and HI values.The primary vertical migration pathway of oil and gas was through faults and fractures into overlying reservoir strata. About 66% of the produced oil is located within the generative basin, which is outlined by an HI contour of 400. The remaining production is concentrated within 30 miles (50 km) outside the 400 HI contour. The generative basin is subdivided by contours of progressively lower hydrogen indices that represent increased levels of thermal maturity and generative capacity of New Albany Shale source rocks. The generative basin was also divided into seven oil-migration catchments. The catchments were determined using a surface-flow hydrologic model with contoured HI values as input to the model.

Illinois, Indiana

Use of isotopes to identify sources of ground water, estimate ground-water-flow rates, and assess aquifer vulnerability in the Calumet region of northwestern Indiana and northeastern Illinois

Isotope data collected in the Calumet Region of northwestern Indiana and northeastern Illinois, one of the most heavily industrialized regions of the United States, indicated that water in the surficial Calumet aquifer is well mixed. The Calumet aquifer is recharged areally by precipitation and locally may be recharged by surface water. The residence time of ground water in the Calumet aquifer is approximately 5 to 15 years. Ground-water-flow rates through the Calumet aquifer are estimated to be 400 to 2,300 feet per year. The permeable deposits, shallow water table, lack of an overlying confining unit, and proximity to numerous contaminant sources indicate that the Calumet aquifer is vulnerable to contamination. Isotopic data indicate that ground water in the confining unit underlying the Calumet aquifer is derived from a variety of sources that include Lake Michigan, modern precipitation, and discharge from the basal Silurian-Devonian bedrock aquifer. The source and apparent age of the water are variable and appear to be affected locally by various geologic and hydraulic factors. The vertical ground-water-flow rate through the unweathered part of the confining unit is about 0.20 feet per year and is about 6.3 feet per year through the weathered part. The data indicate the weathered part of the confining unit may be more vulnerable to contamination than the unweathered confining unit. Ground water in the basal Silurian-Devonian aquifer is derived from Lake Michigan, glacial-age water, and modern precipitation. Post-1952 recharge has occurred in the vicinity of Stony Island, Ill. Ground-water recharge of the Silurian-Devonian aquifer may be occurring near Calumet Harbor. The Silurian-Devonian aquifer is vulnerable to contamination where the confining unit is thin or absent or where the integrity of the confining unit has been compromised.

Illinois, Indiana

Concentrations of Escherichia coli in streams in the Kankakee and lower Wabash River watersheds in Indiana, June-September 1999

Water samples collected from 58 surface- water sites in the Kankakee and Lower Wabash River Watersheds from June through September 1999 were analyzed for concentrations of Escherichia coli bacteria. Each sitewas sampled five times in a 30-day period. Twentynine sites were sampled during June and July, and 29 different sites were sampled during August and September. A five-sample geometric mean of concentrations was computed for each site. Concentrations of Escherichia coli (E. coli) in 126 of the 289 samples exceeded the State of Indiana single-sample standard of 235 colonies per 100 milliliters for waters used for recreation. Concentrations in samples from 38 of the 58 sites exceeded the State of Indiana standard for a five-sample geometric mean of 125 colonies per 100 milliliters for waters used for recreation. Ten of the 58 sites were at or near U.S. Geological Survey streamflow-gaging stations. Based on records from the streamflowgaging stations, 18 percent of the samples collected at these sites were collected at streamflows above the median daily discharge for each station. E. coli concentrations and turbidity measurements collected during 1999 were analyzed in concert with similar concentration and turbidity data collected in 1998 at streams within the Upper Wabash River Watershed in Indiana to investigate the relation between concentrations of bacteria and turbidity. The analysis indicated a statistically significant correlation between concentrations of E. coli and turbidity. If the turbidity was greater than 83 nephelometric turbidity units, the E. coli concentration always exceeded the singlesample standard. If, however, the turbidity was less than 83 nephelometric turbidity units, concentrations of E. coli were not always below the single-sample standard.

Illinois, Indiana

Hydrology and geochemistry of a slag-affected aquifer and chemical characteristics of slag-affected ground water, northwestern Indiana and northeastern Illinois

Slag is a by-product of steel manufacturing and a ubiquitous fill material in northwestern Indiana. Ground water associated with slag deposits generally is characterized by high pH and elevated concentrations of many inorganic water-quality constituents. The U.S. Geological Survey, in cooperation with the Indiana Department of Environmental Management, conducted a study in northwestern Indiana from June 1995 to September 1996 to improve understanding of the effects of slag deposits on the water quality of a glacial-outwash aquifer. The Bairstow Landfill, a slag-fill deposit overlying the Calumet aquifer near Hammond, Indiana, was studied to represent conditions in slag-deposit settings that are common in northwestern Indiana. Ground water from 10 observation wells, located in four nests at the site, and surface water from the adjacent Lake George were analyzed for values of field-measured parameters and concentrations of major ions, nutrients, trace elements, and bulk properties. Solid-phase samples of slag and aquifer sediment collected during drilling were examined with X-ray diffraction and geochemical digestion and analysis. Concentrations of calcium, potassium, sodium, and sulfate were highest in wells screened partly or fully in slag. Potassium concentrations in ground water ranged from 2.9 to 120 milligrams per liter (mg/L), were highest in water from slag deposits, and decreased with depth. The highest concentrations for aluminum, barium, molybdenum, nickel, and selenium were in water from the slag. Silica concentrations were highest in wells screened directly beneath the slag-aquifer interface, and magnesium concentrations were highest in intermediate and deep aquifer wells. Silica concentrations in shallow and intermediate aquifer wells ranged from 27 to 41 mg/L and were about 10 times greater than those in water from slag deposits. The highest concentrations for chromium, lead, and zinc were in ground water from immediately below the slag-aquifer interface. The solid-phase analyses indicated that calcite, dolomite, and quartz generally were present throughout the slag-aquifer system; barian celestite, cristobalite, manganese-bearing calcite, and minrecordite were present in fewer samples. Trace elements that are liberated from the slag may be incorporated as impurities during precipitation of major minerals, sorbed onto clays and other grainsize fractions not analyzed as part of this study, or present in low-abundance minerals that were not detected by the X-ray analysis. Mass-balance and speciation programs were used to identify geochemical processes that may be occurring as water infiltrates through the slag, flows into the aquifer, and discharges into Lake George. The geochemical models indicate that precipitation of calcite may be occurring where slag-affected water enters the aquifer. Models also indicate that dolomite precipitation and clay-mineral dissolution may be occurring at the slag-aquifer interface; however, dolomite precipitation is generally believed to require geologically long time periods. Silica may be dissolving where slag-affected ground water enters the aquifer and may be precipitating where slag-affected ground water discharges to the lakebed of Lake George. In addition to the site-specific study, a statistical analysis of regional water quality was done to compare ground water in wells affected and unaffected by slag. When com-pared to wells in background locations in the Calumet aquifer, wells screened in slag across northwestern Indiana and northeastern Illinois generally had relatively higher pH and specific-conductance values and relatively higher concentrations of alkalinity, dissolved solids, suspended solids, total organic carbon, calcium, potassium, sodium, chloride, aluminum, barium, and possibly magnesium, sulfate, chromium, cobalt, copper, cyanide, manganese, mercury, nickel, and vanadium. When compared to wells in slag and wells in background locations, ground water from immediately beneath or immediately downgradient from slag had relatively high concentrations of arsenic and silica. Water-quality characteristics in ground water at the Bairstow Landfill were similar to water-quality characteristics in slag-contact and slag-affected wells throughout northwestern Indiana.

Illinois, Indiana