USGS Science⌕ Search

SEARCH · USGS Science

Results for “Information Bulletin”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65Linked to original sources

From field station to forecast: Managing data at the Alaska Volcano Observatory

The Alaska Volcano Observatory (AVO) uses multidisciplinary data to monitor and study dozens of active and potentially active volcanoes. Here, we provide an overview of internally and externally generated data types, tools and resources used in their management, and challenges faced. Data sources include the following: (1) a multiparameter (seismic, infrasound, GNSS, web cameras) ground-based monitoring network that spans 3000 km and transmits data in real time; (2) a variety of satellite-borne sensors that provide information about surface change and volcanic emissions; (3) geologic and gas field campaigns; and (4) other external data products that provide situation awareness. Each data type requires distinct acquisition, processing, storage, visualization, and archiving approaches. AVO uses a variety of externally and internally developed tools to handle individual data types as well as multidisciplinary volcanological data. A primary tool is the Geologic Database of Information on Volcanoes in Alaska (GeoDIVA), which stores detailed, searchable information on more than 140 volcanoes and over 1000 eruptions and unrest events, including images, eruption descriptions, and geologic station and sample data, metadata, and analyses. It interacts with other internal tools that store monitoring reports and other operational records. Additional data management resources used by AVO assist with alarms and alerts, state-of-health monitoring, and multiparameter visualization. Requirements for 24/7 accessibility, the ever-expanding portfolio of data, and transitioning new tools from development to operations are all challenges faced by AVO and other volcano observatories. AVO strives to meet FAIR data practices and ensure that data are available to national and international community efforts using external repositories as well as those hosted by AVO and its parent institutions.

Alaska↗

Using Landsat imagery to detect, monitor, and project net landscape change

Detailed landscape information is a necessary component to bird habitat conservation planning. The U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center has been providing information on the Earth’s surface for over 40 years via the continuous series of Landsat satellites. In addition to operating, processing, and disseminating satellite images, EROS is the home to nationwide and global landscape mapping, monitoring, and projection products, including: National Land Cover Database (NLCD) – the definitive land cover dataset for the U.S., with updates occurring at five-year intervals; Global Land Cover Monitoring – producing 30m resolution global land cover; LANDFIRE – Landscape Fire and Resource Management Planning Tools–EROS is a partner in this joint program between U.S. Department of Agriculture and Department of Interior that produces consistent, comprehensive, geospatial data and databases that describe vegetation, wildland fuel, and fire regimes across the U.S.; Land Cover Trends – a landscape monitoring and assessment effort to understand the rates, trends, causes, and consequences of contemporary U.S. land use and land cover change; and Land Use and Land Cover (LULC) Modeling – a project extending contemporary databases of landscape change forward and backward in time through moderate-resolution land cover projections.

All Bird Bulletin↗

Survival and harvest-related mortality of white-tailed deer in Massachusetts

We monitored 142 radiocollared adult (≥1.0 yr old) white-tailed deer ( Odocoileus virginianus ) in 3 study areas of Massachusetts, USA, to estimate annual survival and mortality due to legal hunting. We then applied these rates to deer harvest information to estimate deer population trends over time, and compared these to trends derived solely from harvest data estimates. Estimated adult female survival rates were similar (0.82–0.86), and uniformly high, across 3 management zones in Massachusetts that differed in landscape composition, human density, and harvest regulations. Legal hunting accounted for 16–29% of all adult female mortality. Estimated adult male survival rates varied from 0.55 to 0.79, and legal hunting accounted for 40–75% of all mortality. Use of composite hunting mortality rates produced realistic estimates for adult deer populations in 2 zones, but not for the third, where estimation was hindered by regulatory restrictions on antlerless deer harvest. In addition, the population estimates we calculated were generally higher than those derived from population reconstruction, likely due to relatively low harvest pressure. Legal harvest may not be the dominant form of deer mortality in developed landscapes; thus, estimates of populations or trends that rely solely on harvest data will likely be underestimates.

Massachusetts↗

Evaluation of a waistband for attaching external radiotransmitters to anurans

Radiotelemetry provides fine-scale temporal and spatial information about an individual's movements and habitat use; however, its use for monitoring amphibians has been restricted by transmitter mass and lack of suitable attachment techniques. We describe a novel waistband for attaching external radiotransmitters to anurans and evaluate the percentages of resulting abrasions, lacerations, and shed transmitters. We used radiotelemetry to monitor movements and habitat use of wood frogs ( Lithobates sylvaticus ) in 2006 and 2011–2013 in Maine, USA; American toads ( Anaxyrus americanus ) in 2012 in North Carolina, USA; and, wood frogs, southern leopard frogs ( L. sphenocephalus ), and green frogs ( L. clamitans ) in 2012 in South Carolina, USA. We monitored 172 anurans for 1–365 days (56.4 ± 59.4) in a single year and 1–691 days (60.5 ± 94.1) across years. Our waistband resulted in an injury percentage comparable to 7 alternative anuran waistband attachment techniques; however, 12.5% fewer anurans shed their waistband when attached with our technique. Waistband retention facilitates longer monitoring periods and, thus, provides a greater quantity of data per radiotagged individual.

Wildlife Society Bulletin↗

Structural geologic interpretations from radar imagery

Certain structural geologic features may be more readily recognized on sidelooking airborne radar (SLAR) images than on conventional aerial photographs, other remote sensor imagery, or by ground observations. SLAR systems look obliquely to one or both sides and their images resemble aerial photographs taken at low sun angle with the sun directly behind the camera. They differ from air photos in geometry, resolution, and information content. Radar operates at much lower frequencies than the human eye, camera, or infrared sensors, and thus "sees" differently. The lower frequency enables it to penetrate most clouds and some precipitation, haze, dust, and some vegetation. Radar provides its own illumination, which can be closely controlled in intensity and frequency. It is narrow band, or essentially monochromatic. Low relief and subdued features are accentuated when viewed from the proper direction. Runs over the same area in significantly different directions (more than 45° from each other), show that images taken in one direction may emphasize features that are not emphasized on those taken in the other direction; optimum direction is determined by those features which need to be emphasized for study purposes. Lineaments interpreted as faults stand out on radar imagery of central and western Nevada; folded sedimentary rocks cut by faults can be clearly seen on radar imagery of northern Alabama. In these areas, certain structural and stratigraphic features are more pronounced on radar images than on conventional photographs; thus radar imagery materially aids structural interpretation.

Alabama;Nevada↗

A revised ground-motion and intensity interpolation scheme for shakemap

We describe a weighted-average approach for incorporating various types of data (observed peak ground motions and intensities and estimates from groundmotion prediction equations) into the ShakeMap ground motion and intensity mapping framework. This approach represents a fundamental revision of our existing ShakeMap methodology. In addition, the increased availability of near-real-time macroseismic intensity data, the development of newrelationships between intensity and peak ground motions, and new relationships to directly predict intensity from earthquake source information have facilitated the inclusion of intensity measurements directly into ShakeMap computations. Our approach allows for the combination of (1) direct observations (ground-motion measurements or reported intensities), (2) observations converted from intensity to ground motion (or vice versa), and (3) estimated ground motions and intensities from prediction equations or numerical models. Critically, each of the aforementioned data types must include an estimate of its uncertainties, including those caused by scaling the influence of observations to surrounding grid points and those associated with estimates given an unknown fault geometry. The ShakeMap ground-motion and intensity estimates are an uncertainty-weighted combination of these various data and estimates. A natural by-product of this interpolation process is an estimate of total uncertainty at each point on the map, which can be vital for comprehensive inventory loss calculations. We perform a number of tests to validate this new methodology and find that it produces a substantial improvement in the accuracy of ground-motion predictions over empirical prediction equations alone.

Bulletin of the Seismological Society of America↗

Real-time earthquake detection and alerting behavior of PLUM ground-motion-based early warning in the United States

We examine the real‐time earthquake detection and alerting behavior of the Propagation of Local Undamped Motion (PLUM) earthquake early warning (EEW) algorithm and compare PLUM’s performance with the real‐time performance of the current source‐characterization‐based ShakeAlert System. In the United States (U.S.), PLUM uses a two‐station approach to detect earthquakes. Once a detection is confirmed, observed modified Mercalli intensity (MMI) distributions are forecast onto a regular grid, in which the preferred alert regions are grid cells with MMI 4.0+ forecasts. Although locations of dense station coverage allow PLUM to detect small ( M < 4.5) earthquakes typically not considered for EEW in the U.S., a PLUM detection on a small earthquake does not always generate an alert. This is because PLUM alerts are determined by current shaking distributions. If the MMI 4.0+ shaking subsides prior to detection confirmation by shaking at a second neighboring station, the prior MMI 4.0+ information will not be in the alert forecasts. Of the 432 M 3.0+ U.S. West Coast earthquakes in 2021, 33 produced ground motions large enough to be detected by PLUM. Twenty‐four generated MMI 4.0+ PLUM alerts, whereas ShakeAlert issued public EEW alerts for 13 of these earthquakes. We compare PLUM and ShakeAlert alert regions with ShakeMap and “Did You Feel It?” intensity distributions. Because PLUM alert regions surround stations observed to have strong ground motions (regardless of earthquake magnitude), PLUM alerts reliably include locations that experience significant shaking. This is not necessarily the case for ShakeAlert alert regions when there are large errors in magnitude or epicenter estimates. For two of the largest earthquakes in our real‐time dataset, the M 6.0 Antelope Valley and M 5.1 Petrolia earthquakes, the inclusion of PLUM would have improved real‐time ShakeAlert performance. Our results indicate that incorporation of PLUM into ShakeAlert will improve the robustness of the EEW system.

California, Oregon, Washington↗

Critically assessing the utility of portable lead analyzers for wildlife conservation

Lead (Pb) exposure in wildlife is a widespread management and conservation concern. Quantitative determination of Pb concentrations in wildlife tissues is the foundation for estimating exposure and risk. Development of low‐cost, portable instruments has improved access and cost‐effectiveness of determining Pb concentrations in blood samples, while also facilitating the ability for wildlife researchers to conduct near real‐time Pb testing. However, these instruments, which use anodic stripping voltammetry (ASV) methodology, may produce an analytical bias in wildlife‐blood Pb concentrations. Additionally, their simplicity invites use without appropriate quality‐assurance–quality‐control measures. Together, these factors can reduce data quality and hamper the ability to evaluate it, raising concerns about use of these instruments to inform important conservation issues. We document the extent to which this bias is addressed in the wildlife toxicology literature, develop quantitative approaches for correcting the bias, and provide recommendations to ensure robust data quality when using these instruments. Of the 25 studies we reviewed that referenced ASV use for determining Pb exposure in wildlife, only 32% acknowledged the existence of bias from the instrument. Importantly, another 20% of the studies actually reported ASV and spectroscopic‐based results together without acknowledging their lack of equivalence. Using a multispecies data set of avian blood Pb concentrations, we found that ASV‐based estimates of paired blood Pb concentrations were 30–38% lower than those from standard spectrometric‐based methods. We provide regression equations based on this analysis of 453 blood samples to allow users of ASV instruments to adjust Pb concentrations to spectrometric‐equivalent values, and propose a series of guidelines to follow when using these instruments to improve data validity.

Wildlife Society Bulletin↗

Precipitation, density, and population dynamics of desert bighorn sheep on San Andres National Wildlife Refuge, New Mexico

Understanding the determinants of population size and performance for desert bighorn sheep (Ovis canadensis mexicana) is critical to develop effective recovery and management strategies. In arid environments, plant communities and consequently herbivore populations are strongly dependent upon precipitation, which is highly variable seasonally and annually. We conducted a retrospective exploratory analysis of desert bighorn sheep population dynamics on San Andres National Wildlife Refuge (SANWR), New Mexico, 1941-1976, by modeling sheep population size as a function of previous population sizes and precipitation. Population size and trend of desert bighorn were best and well described (R 2=0.89) by a model that included only total annual precipitation as a covariate. Models incorporating density-dependence, delayed density-dependence, and combinations of density and precipitation were less informative than the model containing precipitation alone (??AlCc=8.5-22.5). Lamb:female ratios were positively related to precipitation (current year: F1,34=7.09, P=0.012; previous year: F1,33=3.37, P=0.075) but were unrelated to population size (current year. F1,34=0.04, P=0.843; previous year: F1,33 =0.14, P=0.715). Instantaneous population rate of increase (r) was related to population size (F1,33=5.55; P=0.025). Precipitation limited populations of desert bighorn sheep on SANWR primarily in a density-independent manner by affecting production or survival of lambs, likely through influences on forage quantity and quality. Habitat evaluations and recovery plans for desert bighorn sheep need to consider fundamental influences on desert bighorn populations such as precipitation and food, rather than focus solely on proximate issues such as security cover, predation, and disease. Moreover, the concept of carrying capacity for desert bighorn sheep may need re-evaluation in respect to highly variable (CV =35.6%) localized precipitation patterns. On SANWR carrying capacity for desert bighorn sheep was zero when total annual precipitation was <28.2 cm.

Wildlife Society Bulletin↗

Applications of nonergodic site response models to ShakeAlert case studies in the Los Angeles area

In this study, we explore whether the Parker and Baltay (2022) site response models for the Los Angeles (LA) basin region can improve ground‐motion forecasts in the U.S. Geological Survey ShakeAlert earthquake early warning system (hereafter ShakeAlert). We implement the peak ground acceleration and peak ground velocity site response models of Parker and Baltay (2022) in ShakeAlert via the earthquake information to ground‐motion (hereafter eqinfo2GM) module, which predicts ground motions from the estimated earthquake parameters of magnitude, rupture length, and location. The nonergodic site response models for the greater LA area were developed using ground motions from 414 M 3–7.3 earthquakes in southern California. We test nonergodic ground‐motion forecasts for five earthquakes in the LA area: the 1994 M 6.7 Northridge earthquake, the 2008 M 5.4 Chino Hills earthquake, the 2019 M 7.1 Ridgecrest earthquake, the 2020 M 4.5 South El Monte earthquake, and a synthetic M 7.8 earthquake on the southern San Andreas fault from the ShakeOut scenario, which was the basis of a statewide emergency response exercise. From the test results, we find that with the nonergodic site response applied, ShakeAlert not only alerts larger areas but can also result in longer warning times in LA region. In addition, the modified Mercalli intensity (MMI) ground‐motion predictions generated by the ShakeAlert eqinfo2GM module are improved in accuracy when compared with the corresponding ShakeMap ground‐truth MMI when the nonergodic site response model is applied.

California↗

Noise constraints on global body‐wave measurement thresholds

Intermediate sized earthquakes (≈ M 4–6.5) are often measured using the teleseismic body‐wave magnitude ( ⁠𝑚b⁠ ). 𝑚b measurements are especially critical at the lower end of this range when teleseismic waveform modeling techniques (i.e., moment tensor analysis) are difficult. The U.S. Geological Survey National Earthquake Information Center (NEIC) determines the location and magnitude of all M 5 and greater earthquakes worldwide within 20 min of the rupture time, and therefore accurate 𝑚b magnitude estimates are essential to fulfill its mission. To better understand how network geometry and noise levels affect the global response capabilities, we developed a method to spatially estimate the minimum measurable 𝑚b⁠ . To do this, we compare expected 𝑚b amplitudes at every station to the station’s background noise level. We find that using NEIC’s current network geometry and these idealized thresholds, NEIC can potentially estimate 𝑚b magnitudes down to M 4.5 globally. Low‐latitude regions in the Southern Hemisphere present the biggest opportunity to improve monitoring capabilities. However, logistically they also present the biggest hurdles for network operators. Finally, to test the resiliency of the network we removed the 20 most important stations and found the 𝑚b threshold remains 𝑚b 4.5. However, the region where only 𝑚b 4.5 and greater can be estimated increases and is again restricted to the Southern Hemisphere.

Bulletin of the Seismological Society of America↗

Thermal history of rocks in southern San Joaquin Valley, California: evidence from fission-track analysis

The theory of the fission-track method and its application to sedimentary basin analysis is illustrated by a case study in the southern San Joaquin Valley, California. Fission tracks provide a powerful tool for studying the thermal history of sedimentary basins because the two minerals most commonly used in fission-track studies, apatite and zircon, occur as detrital constituents in many sedimentary rocks, and their annealing temperatures span the main temperature range for oil generation. Fission tracks also provide information on the sedimentation record and provenance of rocks in a basin. We have used fission-track analysis to study the thermal and depositional history of the subsurface Tertiary sedimentary rocks on both sides of the active White Wolf reverse fault in the southern San Joaquin Valley. The distinctly different thermal histories of the rocks in the two structural blocks are clearly reflected in the apatite fission-track data, which suggest that rocks in the rapidly subsiding basin northwest of the fault have been near their present temperature for only about 1 m.y. compared with about 10 m.y. for rocks southeast of the fault. These estimates of heating time agree with previous estimates for these rocks. Zircon fission-track data indicate that the Tertiary sediments were derived from parent rocks of more than one age. However, from at least the Eocene to late Miocene or Pliocene, the major sediment source was rocks related to the youngest Sierra Nevada Mesozoic intrusive complexes, which are presently exposed east and south of the southern San Joaquin Valley.

California↗

Recent deformation along the offshore Malibu Coast, Dume, and related faults west of Point Dume, southern California

Offshore faults west of Point Dume, southern California, are part of an important regional fault system that extends for about 206 km, from near the city of Los Angeles westward along the south flank of the Santa Monica Mountains and through the northern Channel Islands. This boundary fault system separates the western Transverse Ranges, on the north, from the California Continental Borderland, on the south. Previous research showed that the fault system includes many active fault strands; consequently, the entire system is considered a serious potential earthquake hazard to nearby Los Angeles. We present an integrated analysis of multichannel seismic- and high-resolution seismic-reflection data and multibeam-bathymetric information to focus on the central part of the fault system that lies west of Point Dume. We show that some of the main offshore faults have cumulative displacements of 3-5 km, and many faults are currently active because they deform the seafloor or very shallow sediment layers. The main offshore fault is the Dume fault, a large north-dipping reverse fault. In the eastern part of the study area, this fault offsets the seafloor, showing Holocene displacement. Onshore, the Malibu Coast fault dips steeply north, is active, and shows left-oblique slip. The probable offshore extension of this fault is a large fault that dips steeply in its upper part but flattens at depth. High-resolution seismic data show that this fault deforms shallow sediment making up the Hueneme fan complex, indicating Holocene activity. A structure near Sycamore knoll strikes transversely to the main faults and could be important to the analysis of the regional earthquake hazard because the structure might form a boundary between earthquake-rupture segments.

Bulletin of the Seismological Society of America↗

Potash brines in the Great Salt Lake Desert, Utah

During and immediately after the war the brines of the Salduro Marsh, in the Great Salt Lake Desert, were a source of considerable potash for the domestic supply. Although no potash has been produced from these brines in the last few years, a continued interest in the area has been shown by a large number of filings, in different parts of the desert, under the potash law of October 2, 1917 (40 Stat. 297), and the regulations issued under that law by the Department of the Interior on March 21, 1918, in Circular 594 (46 L. D. 323). As the basis for appropriate action on these permit, lease, and patent filings the Department of the Interior required information concerning the mode of occurrence, distribution, and general characteristics of the potassium-bearing brines in this region. Field investigation was accordingly undertaken for the primary purpose of satisfying these requirements. The work was done by the Geological Survey and the General Land Office in cooperation with the Bureau of Mines. The Bureau of Mines furnished an autotruck and a mechanic, who also assisted in the field work, and, for a period of a month, a mine rescue car, which permitted the establishment of headquarters in parts of the desert far from sources of supplies and water. The General Land Office surveyed the boundaries of townships over the desert and located the drill holes from which samples of brine were taken. The writer directed the work for the Geological Survey, and the analyses of the samples of brine were made in the Survey's chemical laboratory. The greater part of the field work was done between April 17 and July 7, 1925, and the work was finished during the two periods August 17-26 and September 30 to October 2 . Adverse weather throughout the season hindered the work.

Utah↗

Energy resource potential of natural gas hydrates

The discovery of large gas hydrate accumulations in terrestrial permafrost regions of the Arctic and beneath the sea along the outer continental margins of the world's oceans has heightened interest in gas hydrates as a possible energy resource. However, significant to potentially insurmountable technical issues must be resolved be fore gas hydrates can be considered a viable option for affordable supplies of natural gas. The combined information from Arctic gas hydrate studies shows that, in permafrost regions, gas hydrates may exist at subsurface depths ranging from about 130 to 2000 m. The presence of gas hydrates in offshore continental margins has been inferred mainly from anomalous seismic reflectors, known as bottom-simulating reflectors, that have been mapped at depths below the sea floor ranging from about 100 to 1100 m. Current estimates of the amount of gas in the world's marine and permafrost gas hydrate accumulations are in rough accord at about 20,000 trillion m 3 . Disagreements over fundamental issues such as the volume of gas stored within delineated gas hydrate accumulations and the concentration of gas hydrates within hydrate-bearing strata have demonstrated that we know little about gas hydrates. Recently, however, several countries, including Japan, India, and the United States, have launched ambitious national projects to further examine the resource potential of gas hydrates. These projects may help answer key questions dealing with the properties of gas hydrate reservoirs, the design of production systems, and, most important, the costs and economics of gas hydrate production.

American Association of Petroleum Geologists Bulle↗

The predictive skills of elastic Coulomb rate-and-state aftershock forecasts during the 2019 Ridgecrest, California, earthquake sequence

Operational earthquake forecasting protocols commonly use statistical models for their recognized ease of implementation and robustness in describing the short-term spatiotemporal patterns of triggered seismicity. However, recent advances on physics-based aftershock forecasting reveal comparable performance to the standard statistical counterparts with significantly improved predictive skills when fault and stress field heterogeneities are considered. Here, we perform a pseudo-prospective forecasting experiment during the first month of the 2019 Ridgecrest (California) earthquake sequence. We develop seven Coulomb rate-and-state models that couple static stress change estimates with continuum mechanics expressed by the rate-and-state friction laws. Our model parametrization supports a gradually increasing complexity; we start from a preliminary model implementation with simplified slip distributions and spatially homogeneous receiver faults to reach an enhanced one featuring optimized fault constitutive parameters, finite-fault slip models, secondary triggering effects, and spatially heterogenous planes informed by pre-existing ruptures. The data-rich environment of Southern California allows us to test whether incorporating data collected in near real-time during an unfolding earthquake sequence boosts our predictive power. We assess the absolute and relative performance of the forecasts by means of statistical tests used within the Collaboratory for the Study of Earthquake Predictability (CSEP) and compare their skills against a standard benchmark ETAS model for the short (24 hours after the two Ridgecrest mainshocks) and intermediate-term (one month). Stress-based forecasts expect heightened rates along the whole near-fault region and increased expected seismicity rates in Central Garlock Fault. Our comparative model evaluation supports that faulting heterogeneities coupled with secondary triggering effects are the most critical success components behind physics-based forecasts, but also underlines the importance of model updates incorporating near real-time available aftershock data reaching better performance than ETAS models. We explore the physical basis behind our results by investigating the localized shut down of pre-existing normal faults in the Ridgecrest near-source area.

California↗

Use of relational databases to evaluate regional petroleum accumulation, groundwater flow, and CO2 sequestration in Kansas

Large-scale relational databases and geographic information system tools are used to integrate temperature, pressure, and water geo-chemistry data from numerous wells to better understand regional-scale geothermal and hydrogeological regimes of the lower Paleozoic aquifer systems in the mid-continent and to evaluate their potential for geologic CO2 sequestration. The lower Paleozoic (Cambrian to Mississippian) aquifer systems in Kansas, Missouri, and Oklahoma comprise one of the largest regional-scale saline aquifer systems in North America. Understanding hydrologic conditions and processes of these regional-scale aquifer systems provides insight to the evolution of the various sedimentary basins, migration of hydrocarbons out of the Anadarko and Arkoma basins, and the distribution of Arbuckle petroleum reservoirs across Kansas and provides a basis to evaluate CO2 sequestration potential. The Cambrian and Ordovician stratigraphic units form a saline aquifer that is in hydrologic continuity with the freshwater recharge from the Ozark plateau and along the Nemaha anticline. The hydrologic continuity with areas of freshwater recharge provides an explanation for the apparent underpressure in the Arbuckle Group. Copyright ?? 2005. The American Association of Petroleum Geologists. All rights reserved.

American Association of Petroleum Geologists Bulle↗

Hydrology of carbonate aquifers in southwestern Linn County and adjacent parts of Benton, Iowa, and Johnson Counties, Iowa

Groundwater is the major source of water in Linn County and the surrounding area. Approximately 90 percent of the groundwater production is from Silurian, Devonian, and Quaternary aquifers. The Silurian and Devonian aquifers consist of limestone and dolomite with minor shale beds, which have a regional dip to the southwest of approximately 20 feet per mile. The Silurian aquifer in east-central Iowa is confined from below by Upper Ordovician, Maquoketa Formation shales, and from above by the Kenwood Member of the Wapsipinicon Formation and the Otis and Bertram formations. The Quaternary aquifer consists of unconsolidated sand and gravel beds in the glacial drift, and in the alluvium which is associated with modern streams. The alluvium consists of lenticular beds of poorly-to well-sorted silt, sand, and gravel. The sand and gravel beds are interlayered with relatively-impermeable beds of till, silt, and clay. Water moves through the Silurian aquifer in part due to a complex distribution of porous and dense carbonate facies. Horizons containing skeletal molds in the Silurian dolomite have porosities as much as 39 percent, and are laterally equivalent to dolomites with porosities as little as less than one percent. Because of subsequent fracturing and solutional enlargement of these porous horizons, hydrologic correlation of the primary water-yielding zones is not always possible. One horizon, however, does occur approximately 70 to 105 feet above the base of the Silurian, and is the most consistently productive water-yielding unit in the area. This horizon is Informally referred to as the Farmers Creek aquifer. The potentiometric surface of the Silurian aquifer has a gradient towards the Cedar River, indicating discharge from the aquifer through the alluvium into the river. By comparison, the potentiometric surface of the overlying Devonian aquifer is equal to that of the Silurian and may range to more than 40 feet higher. Yields to individual wells completed in the Silurian and Devonian carbonate aquifers vary from less than 10 to about 500 gallons per minute. Individual wells completed in the Quaternary aquifer yield as much as 2,000 gallons per minute. Water analyses from the Devonian and Silurian aquifers indicate that they are of similar chemical quality at most locations in the study area. However, they may commonly contain concentrations of sulfate that exceed 1,000 mil grams per liter. Dissolved-solids concentrations as much as 2,350 milligrams per liter occur in the Silurian aquifer in the western and southwestern part of the study area. Water from the Quaternary aquifer generally is suitable for most uses and dissolved-solids concentrations generally are less than 750 milligrams per liter.

Iowa↗