USGS Science⌕ Search

SEARCH · USGS Science

Results for “Letter Report”

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

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

410 records · Page 23Linked to original sources

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗

Western U.S. deformation models for the 2023 update to the U.S. National Seismic Hazard Model

This report describes geodetic and geologic information used to constrain deformation models of the 2023 update to the National Seismic Hazard Model (NSHM), a set of deformation models to interpret these data, and their implications for earthquake rates in the western United States. Recent updates provide a much larger data set of Global Positioning System crustal velocities than used in the 2014 NSHM, as well as hundreds of new faults considered as active sources for the 2023 NSHM. These data are interpreted by four geodetic models of deformation that estimate fault slip rates and their uncertainties together with off‐fault moment release rates. Key innovations in the 2023 NSHM relative to past practice include (1) the addition of two new (in addition to two existing) deformation models, (2) the revision and expansion of the geologic slip rate database, (3) accounting for fault creep through development of a creep‐rate model that is employed by the four deformation models, and (4) accounting for time‐dependent earthquake‐cycle effects through development of viscoelastic models of the earthquake cycle along the San Andreas fault and the Cascadia subduction zone. The effort includes development of a geologic deformation model that complements the four geodetic models. The current deformation models provide a new assessment of outstanding discrepancies between geologic and geodetic slip rates, at the same time highlighting the need for both geologic and geodetic slip rates to robustly inform the earthquake rate model.

Seismological Research Letters↗

On algorithmically determined versus traditional macroseismic intensity assignments

The utility of macroseismic data, defined as the effects of earthquakes on humans and the built environment, has been increasingly recognized following the advent of online systems that now produce unprecedented volumes of macroseismic intensity information. Contributed reports from the U.S. Geological Survey “Did You Feel It?” (DYFI) system ( Wald et al ., 1999 ) are used to generate intensity values with an algorithm based on seminal work by Dengler and Dewey (1998) . The algorithm was developed initially to reproduce intensity values assigned by expert opinion using questionnaire results collected by telephone survey. In this article, I discuss reasons why intensity values from (self‐selected) DYFI responses can differ from values that would be assigned by expert opinion given more complete data from randomly selected participants. For example, with the data used by Dengler and Dewey (1998) , intensities near 4 could be determined from the percentage of people who felt shaking in each town. With less spatially rich data from self‐selected participants, this percentage often cannot be determined reliably. Audible noises are key additional diagnostic criteria for modified Mercalli intensity (MMI) 4, but, although the DYFI system includes a question about noise, following Dengler and Dewey (1998) , the DYFI algorithm does not include a noise indicator. At the upper end of the scale, as defined the DYFI algorithm yields a maximum intensity value of 9.05, nominally corresponding to peak ground acceleration of 75% g . These and other factors can result in DYFI values that are low compared to traditional MMI values assigned using expert opinion, even absent factors that can bias traditional MMI assignments. Modern ground‐motion intensity conversion equations determined using DYFI intensities are expected to be appropriate for DYFI intensities, but the results of this study suggest that biases may be introduced if DYFI and traditional intensities are assumed to be interchangeable.

Seismological Research Letters↗

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay↗

The Hf isotopic composition of ferromanganese nodules and crusts and hydrothermal manganese deposits: Implications for seawater Hf

We present Hf and Pb isotopic data, and chemical compositions of the outermost layers of marine ferromanganese deposits of different types (hydrogenous and hydrothermal) with a worldwide distribution. The Hf isotopic compositions display a broad range and refine previously reported regional differences as follows: Atlantic Ocean ɛ Hf = −4to+2, Indian Ocean ɛ Hf = +2to+4, Pacific Ocean ɛ Hf = +3to+10. The most radiogenic Hf isotopic compositions in the Pacific samples are for hydrothermal manganese deposits that also have low 207 Pb 204 Pb, demonstrating that this signature reflects a contribution from hydrothermal venting of Hf leached from oceanic volcanic rocks rather than from riverine inputs, volcanic ash, or eolian dust. Hafnium concentrations in the deposits increase from 20 ppb to 20 ppm with decreasing ɛ Hf , The Hf and Pb isotopic compositions for ferromanganese crusts define an apparent mixing trend between literature values of average continental crust and MORB. The range in ɛ Hf for ferromanganese crusts is narrower than it is for 206 Pb 204 Pb compared to the differences in isotopic composition of the sources of Hf and Pb. This is consistent with Hf having a longer residence time than Pb. The concentration of Hf in ferromanganese crusts has been found to co-vary with growth rate, and inversely correlates with Hf isotopic compositions. Hf isotope ratios may be used to determine not only the source of Hf, but possibly the source of Fe and Mn. Measurements of ɛ Hf and Hf concentrations in nodule tops, bottoms and associated sediments show that the ɛ Hf of nodules is sensitive to sedimentary oxic and sub-oxic diagenetic processes and thus most nodules may not reliably reflect the isotopic composition of Hf in seawater.

Earth and Planetary Science Letters↗

Short-period volcanic gas precursors to phreatic eruptions: Insights from Poás Volcano, Costa Rica

Volcanic eruptions involving interaction with water are amongst the most violent and unpredictable geologic phenomena on Earth. Phreatic eruptions are exceptionally difficult to forecast by traditional geophysical techniques. Here we report on short-term precursory variations in gas emissions related to phreatic blasts at Poás volcano, Costa Rica, as measured with an in situ multiple gas analyzer that was deployed at the edge of the erupting lake. Gas emitted from this hyper-acid crater lake approaches magmatic values of SO 2 /CO 2 1–6 days prior to eruption. The SO 2 flux derived from magmatic degassing through the lake is measureable by differential optical absorption spectrometry (sporadic campaign measurements), which allows us to constrain lake gas output and input for the major gas species during eruptive and non-eruptive periods. We can further calculate power supply to the hydrothermal system using volatile mass balance and thermodynamics, which indicates that the magmatic heat flux into the shallow hydrothermal system increases from ∼27 MW during quiescence to ∼59 MW during periods of phreatic events. These transient pulses of gas and heat from the deeper magmatic system generate both phreatic eruptions and the observed short-term changes in gas composition, because at high gas flux scrubbing of sulfur by the hydrothermal system is both kinetically and thermodynamically inhibited whereas CO 2 gas is always essentially inert in hyperacid conditions. Thus, the SO 2 /CO 2 of lake emissions approaches magmatic values as gas and power supply to the sub-limnic hydrothermal system increase, vaporizing fluids and priming the hydrothermal system for eruption. Our results suggest that high-frequency real-time gas monitoring could provide useful short-term eruptive precursors at volcanoes prone to phreatic explosions.

Poás Volcano↗

Sr and Nd isotopic and trace element compositions of Quaternary volcanic centers of the Southern Andes

Isotopic compositions of samples from six Quaternary volcanoes located in the northern and southern extremities of the Southern Volcanic Zone (SVZ, 33–46°S) of the Andes and from four centers in the Austral Volcanic Zone (AVZ, 49–54°S) range for 87 Sr/ 86 Sr from 0.70280 to 0.70591 and for 143 Nd/ 144 Nd from 0.51314 to 0.51255. The ranges are significantly greater than previously reported from the southern Andes but are different from the isotopic compositions of volcanoes in the central and northern Andes. Basalts and basaltic andesites from three centers just north of the Chile Rise-Trench triple junction have 87 Sr/ 86 Sr , 143 Nd/ 144 Nd , La/Yb , Ba/La , and Hf/Lu that lie within the relatively restricted ranges of the basic magmas erupted from the volcanic centers as far north as 35°S in the SVZ of the Andes. The trace element and Sr and Nd isotopic characteristics of these magmas may be explained by source region contamination of subarc asthenosphere, with contaminants derived from subducted pelagic sediments and seawater-altered basalts by dehydration of subducted oceanic lithosphere. In the northern extremity of the SVZ between 33° and 34°S, basaltic andesites and andesites have higher 87 Sr/ 86 Sr , Rb/Cs , and Hf/Lu , and lower 143 Nd/ 144 Nd than basalts and basaltic andesites erupted farther south in the SVZ, which suggests involvement of components derived from the continental crust. In the AVZ, the most primitive sample, high-Mg andesite from the southernmost volcanic center in the Andes (54°S) has Sr and Nd isotopic compositions and KR/b and Ba/La similar to MORB. The high La/Yb of this sample suggests formation by small degrees of partial melting of subducted MORB with garnet as a residue. Samples from centers farther north in the AVZ show a regionally regular northward increase in SiO 2 , K 2 O, Rb, Ba, Ba/La , and 87 Sr/ 86 Sr/ and decrease in MgO, Sr, KR/b , Rb/Cs , and 143 Nd/ 144 Nd , suggesting increasingly greater degrees of fractional crystallization and associated intra-crustal contamination.

Earth and Planetary Science Letters↗

Multinational evaluation of genetic diversity indicators for the Kunming-Montreal Global Biodiversity Framework

Under the recently adopted Kunming-Montreal Global Biodiversity Framework, 196 Parties committed to reporting the status of genetic diversity for all species. To facilitate reporting, three genetic diversity indicators were developed, two of which focus on processes contributing to genetic diversity conservation: maintaining genetically distinct populations and ensuring populations are large enough to maintain genetic diversity. The major advantage of these indicators is that they can be estimated with or without DNA-based data. However, demonstrating their feasibility requires addressing the methodological challenges of using data gathered from diverse sources, across diverse taxonomic groups, and for countries of varying socio-economic status and biodiversity levels. Here, we assess the genetic indicators for 919 taxa, representing 5271 populations across nine countries, including megadiverse countries and developing economies. Eighty-three percent of the taxa assessed had data available to calculate at least one indicator. Our results show that although the majority of species maintain most populations, 58% of species have populations too small to maintain genetic diversity. Moreover, genetic indicator values suggest that IUCN Red List status and other initiatives fail to assess genetic status, highlighting the critical importance of genetic indicators.

Ecology Letters↗

Geochronology and paleoclimatic implications of the last deglaciation of the Mauna Kea Ice Cap, Hawaii

We present new 3He surface exposure ages on moraines and bedrock near the summit of Mauna Kea, Hawaii, which refine the age of the Mauna Kea Ice Cap during the Local Last Glacial Maximum (LLGM) and identify a subsequent fluctuation of the ice margin. The 3He ages, when combined with those reported previously, indicate that the local ice-cap margin began to retreat from its LLGM extent at 20.5??2.5ka, in agreement with the age of deglaciation determined from LLGM moraines elsewhere in the tropics. The ice-cap margin receded to a position at least 3km upslope for ~4.5-5.0kyr before readvancing nearly to its LLGM extent. The timing of this readvance at ~15.4ka corresponds to a large reduction of the Atlantic meridional overturning circulation (AMOC) following Heinrich Event 1. Subsequent ice-margin retreat began at 14.6??1.9ka, corresponding to a rapid resumption of the AMOC and onset of the B??lling warm interval, with the ice cap melting rapidly to complete deglaciation. Additional 3He ages obtained from a flood deposit date the catastrophic outburst of a moraine-dammed lake roughly coeval with the Younger Dryas cold interval, suggesting a more active hydrological cycle on Mauna Kea at this time. A coupled mass balance and ice dynamics model is used to constrain the climate required to generate ice caps of LLGM and readvance sizes. The depression of the LLGM equilibrium line altitude requires atmospheric cooling of 4.5??1??C, whereas the mass balance modeling indicates an accompanying increase in precipitation of as much as three times that of present. We hypothesize (1) that the LLGM temperature depression was associated with global cooling, (2) that the temperature depression that contributed to the readvance occurred in response to an atmospheric teleconnection to the North Atlantic, and (3) that the precipitation enhancement associated with both events occurred in response to a southward shift in the position of the inter-tropical convergence zone (ITCZ). Such a shift in the ITCZ would have allowed midlatitude cyclones to reach Mauna Kea more frequently which would have increased precipitation at high elevations and caused additional cooling. ?? 2010 Elsevier B.V.

Earth and Planetary Science Letters↗

Osmium isotope and highly siderophile element systematics of the lunar crust

Coupled 187 Os/ 188 Os and highly siderophile element (HSE: Os, Ir, Ru, Pt, Pd, and Re) abundance data are reported for pristine lunar crustal rocks 60025, 62255, 65315 (ferroan anorthosites, FAN) and 76535, 78235, 77215 and a norite clast in 15455 (magnesian-suite rocks, MGS). Osmium isotopes permit more refined discrimination than previously possible of samples that have been contaminated by meteoritic additions and the new results show that some rocks, previously identified as pristine, contain meteorite-derived HSE. Low HSE abundances in FAN and MGS rocks are consistent with derivation from a strongly HSE-depleted lunar mantle. At the time of formation, the lunar floatation crust, represented by FAN, had 1.4 ?? 0.3 pg g - 1 Os, 1.5 ?? 0.6 pg g - 1 Ir, 6.8 ?? 2.7 pg g - 1 Ru, 16 ?? 15 pg g - 1 Pt, 33 ?? 30 pg g - 1 Pd and 0.29 ?? 0.10 pg g - 1 Re (??? 0.00002 ?? CI) and Re/Os ratios that were modestly elevated ( 187 Re/ 188 Os = 0.6 to 1.7) relative to CI chondrites. MGS samples are, on average, characterised by more elevated HSE abundances (??? 0.00007 ?? CI) compared with FAN. This either reflects contrasting mantle-source HSE characteristics of FAN and MGS rocks, or different mantle-crust HSE fractionation behaviour during production of these lithologies. Previous studies of lunar impact-melt rocks have identified possible elevated Ru and Pd in lunar crustal target rocks. The new results provide no supporting evidence for such enrichments. If maximum estimates for HSE in the lunar mantle are compared with FAN and MGS averages, crust-mantle concentration ratios (D-values) must be ??? 0.3. Such D-values are broadly similar to those estimated for partitioning between the terrestrial crust and upper mantle, with the notable exception of Re. Given the presumably completely different mode of origin for the primary lunar floatation crust and tertiary terrestrial continental crust, the potential similarities in crust-mantle HSE partitioning for the Earth and Moon are somewhat surprising. Low HSE abundances in the lunar crust, coupled with estimates of HSE concentrations in the lunar mantle implies there may be a 'missing component' of late-accreted materials (as much as 95%) to the Moon if the Earth/Moon mass-flux estimates are correct and terrestrial mantle HSE abundances were established by late accretion. ?? 2009 Elsevier B.V. All rights reserved.

Earth and Planetary Science Letters↗

Conversing with Pelehonuamea: A workshop combining 1,000+ years of traditional Hawaiian knowledge with 200 years of scientific thought on Kīlauea volcanism

The events surrounding volcanic eruptions and damaging earthquakes in Hawai‘i have often been described in journals, letters, and newspapers articles in the English language; however, the Hawaiian nation was among the most literate of countries in the 19th century, and many Hawaiian-language newspapers were in circulation through all but the earliest decades of the 19th century. Any modern reconstruction of the history of Hawaiian eruptions or earthquakes should take advantage of all available sources, and so we seek to add the Hawaiian-language newspaper articles, journals, stories, and chants to the volcano and earthquake literature. These sources have been used in many recent volcanological studies. Another aspect to the mix of science and traditional Hawaiian values is that many of the volcanic summits in Hawaiʻi are considered sacred to Hawaiians. Hawaiian travelers brought the first Western missionary team to the summit of Kīlauea and advised them of the proper protocols and behaviors while in this sacred area. The missionaries dismissed this advice as native superstition and they began a campaign of aggressively stamping out customs and protocols related to the Hawaiian volcano goddess Pelehonuamea. What has survived as native knowledge of the volcanoes is a few phrases from native guides included in some of the missionaries’ journals, and a few stories. Pualani and Ku ʻ ulei Kanahele provide excellent introductions to the Pelehonuamea chants. In the 21st century, amid a reawakening of Hawaiian culture, modern Hawaiians are demanding protection of their sacred areas, and scientists must be aware of these interests. At the very least, scientists should show respect to Hawaiian values when working in these areas, and should try to minimize disruptions caused by their work. Kaeo Duarte, Peter Mills, and Scott Rowland describe taking this approach in their field work. Traditional knowledge is also contained in place names. It is important not only to preserve old place names and to recover those no longer used, but also to preserve the stories of those places. Bobby Camara talks about the joys and frustrations of getting information on and recovering Hawaiian place names. Finally, we hope that a broader interest in Hawaiian views about locations in Hawaiʻi where physical scientific work is done will be as beneficial to physical scientists as it has been to life scientists investigating Hawaiian lifeforms on land and in the ocean, and that both studies will continue to benefit the native peoples of Hawaiʻi. Note that these proceedings are transcripts of oral presentations illustrated with PowerPoint presentations or charts. Although every effort has been made to assure the accuracy of the oral presentations, there are some gaps where words are not discernible in the voice recordings and are so noted. In other places, bracketed words were added to clarify the speaker’s meaning.

Hawaii↗

Preliminary report on the 16 October 1999 M 7.1 Hector mine, California, earthquake

The M w 7.1 Hector Mine, California, earthquake occurred at 9:46 GMT on 16 October 1999. The event caused minimal damage because it was located in a remote, sparsely populated part of the Mojave Desert, approximately 47 miles east-southeast of Barstow, with epicentral coordinates 34.59°N 116.27°W and a hypocentral depth of 5 ± 3 km. Twelve foreshocks, M 1.9-3.8, preceded the mainshock during the previous twelve hours. All of these events were located close to the hypocenter of the mainshock. The Hector Mine earthquake occurred within the Eastern California Shear Zone (ECSZ). By virtue of its remote location, the societal impact of the Hector Mine earthquake was, fortunately, minimal in spite of the event's appreciable size. The ECSZ is characterized by high seismicity, a high tectonic strain rate, and a broad, distributed zone of north-northwest-trending faults (ECSZ; Figure 1 ; Dokka and Travis, 1990 ; Sauber et al., 1986 ; Sauber et al., 1994 ; Sieh et al., 1993 ). Data regarding the slip rates of faults within the ECSZ suggest that on the order of 15% of the Pacific-North American plate motion occurs along this zone ( Sauber et al., 1986 ; Wesnousky, 1986 ). Most of the faults in the ECSZ have low slip rates and long repeat times for major earthquakes, on the order of several thousands to tens of thousands of years. The occurrence of the Hector Mine earthquake within seven years and only about 30 km east of the 1992 M w 7.3 Landers earthquake suggests that the closely spaced surface faults in the ECSZ are mechanically related. The Hector Mine event involved rupture on two previously mapped fault zones—the Bullion Fault and an unnamed, more northerly-trending fault that is informally referred to in this paper as the Lavic Lake Fault (Dibblee, 1966 , 1967a , b ). Traces of the Bullion Fault exhibit evidence of Holocene displacement and were zoned as active in 1988 under California's Mquist-Priolo Earthquake Fault Zoning Act ( Hart and Bryant, 1997 ). The pattern of rupture along more than one named fault was also observed from the 1992 Landers earthquake ( Hauksson et al., 1993 ; Sieh et al., 1994). Much of the fault zone that produced the Hector Mine earthquake had been buried by relatively young stream deposits, and the fault scarps in bedrock have a subdued morphology. It appears that these faults have not experienced significant offset for perhaps 10,000 years or more ( Hart, 1987 ). Planned future investigations will refine the age of the last event on these faults. The portion of the Lavic Lake Fault that ruptured between the northern end of the Bullion Mountains and Lavic Lake had not previously been mapped. However, our field investigations have identified ancient, subdued fault scarps along portions of the 1999 rupture zone in this area. It thus appears that the entire segment of the Lavic Lake Fault that was involved in the 1999 event had ruptured in the past. As is typical for most faults within the Eastern California Shear Zone, the rate of movement along the Lavic Lake Fault may be quite slow (<1 mm/yr) and should produce earthquakes only infrequently. This event is a reminder that faults that have ruptured in late Quaternary time, but that lack evidence of Holocene displacement, can still produce earthquakes in this low-slip-rate tectonic setting. Additionally, the Hector Mine earthquake is noteworthy for a couple of other reasons. First, it clearly produced triggered seismicity over much of southern California, from the rupture zone toward the south-southwest in particular. Second, as we will discuss, the event may provide new data and insight into recently developed paradigms concerning earthquake interactions and the role of static stress changes. Questions such as these will, of course, be the subject of extensive detailed analyses in years to come. Fortunately, the Hector Mine sequence will provide one of the best data sets obtained to date for a significant earthquake in the United States. Because it occurred when major upgrades to both the regional seismic network (TriNet) and the regional geodetic network (SCIGN) were well underway, the Earth science community will have abundant high-quality data with which to explore the important and interesting questions that have been raised. In this paper, we present and discuss the basic data and preliminary results from the Hector Mine earthquake.

California↗

Results of a modeling workshop concerning preservation and protection of wetlands in North Dakota

In a recently signed letter, the Governor of North Dakota and the Assistant Secretary of the Interior for Fish and Wildlife and Parks charged a joint state-federal study group with examination of two separate questions: 1) mitigation for the Garrison Diversion Project; and 2) planning for long-range protection and preservation of fish and wildlife habitat in North Dakota. The cochair for this study group (the Secretary of the Interior's Field Representative, Denver, Colorado, and the Natural Resources Coordinator for North Dakota) further articulated the charge concerning the second of these two questions to include three steps: 1) development of a general plan for preservation and protection of migratory waterfowl and their associated wetland habitat; 2) a comprehensive analysis of alternative strategies, including opportunities and constraints, for achieving the goals articulated in Step 1; and 3) design of a coordinated state-federal public information program to assist in plan implementation. In order to obtain input from a variety of interests, the joint study group initiated step 2 activities with a five-day workshop in Bismarck, N. D.; December 8-12, 1980. The objectives of the workshop were: 1) to identify alternative strategies for preserving and enhancing waterfowl production habitat in North Dakota; 2) to identify opportunities and constraints associated with those alternatives; and 3) to promote communication and understanding of the implications of those alternatives for all affected parties. To achieve these objectives, the workshop utilized a group of concepts and techniques collectively known as Adaptive Environmental Assessment (AEA). Developed by Dr. C. S. Holling and his co-workers at the University of British Columbia, the AEA process involves planners, managers, scientists, and other interested parties in a structures atmosphere whose focus is the construction and examination of a computerized simulation model of the resource system under consideration. The modeling process is used to promote communication, identify pertinent issues, identify key data gaps and uncertainties, direct research efforts to fill those gaps, and explore the possible consequences of various management alternatives. The workshop, which was facilitated by the AEA Group of the U.S. Fish and Wildlife Service (FWS), was attended by approximately 25 invited participants representing a variety of interests concerned with the wetlands protection issue in North Dakota. During the week workshop participants conceptualized and constructed a computerized simulation model incorporating many of the hydrologic, agricultural, and wildlife aspects of the wetlands issue. During the process of constructing this model and examining its behavior, participants identified several interesting alternative strategies that may prove useful in an overall wetland protection program, along with a variety of constraints associated with each. Perhaps the most interesting of these alternatives revolve around the idea that there may be a variety of cases in which water can be retained on the land, with concommitant benefits both for wildlife and flood control, without detriment to agricultural productivity. Examples of this kind of activity include flooding of previously drained Type I wetlands in summer fallow areas, installation of smaller drains in Type I and III wetlands to reduce the rate at which water runs off in the spring, and use of strategically located gates in drainage channels associated with the state highway system to slow runoff and create wetland habitat. Several other alternatives discussed at the workshop are related to the notion of using available, uncommitted water supplied to enhance or create wetlands. Several cases were cited in which more certain water supplies would be useful in increasing waterfowl production or reducing disease problems, especially in dry years. Water for such purposes might come from a variety of current of proposed water development projects, both large and small scale. Finally, the potential for re-establishment of "unsuccessfully" drained (i.e., not consistently usable for agricultural purposes) wetlands was discussed at some length. There are apparently substantial acreages of such wetlands in North Dakota and a program to acquire and rehabilitate them might be of considerable utility. The workshop was thus successful in accomplishing its first two objectives--identification of alternative strategies, opportunities, and constraints. However, it would be naive to suppose that any of the alternatives discussed offers a complete, simple solution to the wetlands issue in North Dakota. The most important result of the workshop may therefore be that which was accomplished relative to the third objective--promotion of communication and understanding. It was gratifying and encouraging to see the spirit of communication and cooperation that developed by the end of the workshop. The fact that representatives of many of the interests concerned with the wetlands issue participated in an open exchange of ideas and information marks an important step forward. We believe that it is imperative that this cooperative attitude be maintained, and that there are a variety of ways in which this might be accomplished. Perhaps the simplest would be a small-scale pilot project and research effort to determine the effects of wetland maintenance on summer fallow areas. Such a research program would provide not only useful information, but an opportunity for many of the affected interests to begin working toward mutually acceptable solutions to the overall wetlands issue.

North Dakota↗

Characterization and occurrence of confined and unconfined aquifers in Quaternary sediments in the glaciated conterminous United States

The glacial aquifer system, which is a collection of aquifers within Quaternary sediments in the glaciated conterminous United States, is a principal aquifer that supplies groundwater that serves about 42 million people and accounts for about 5 percent of the Nation’s drinking water. This aquifer system (the area of maximum glacial advance) underlies parts of 25 States and covers 1.87×10 6 square kilometers. A hydrogeologic framework is presented that divides the glaciated United States into 17 distinct hydrogeologic terranes using a geologic approach based on previous mapping. Each hydrogeologic terrane contains Quaternary sediment that is derived from a common depositional history and can be characterized by similar texture and thickness. Characteristics of Quaternary sediments are described using attributes computed from a lithologic database of well logs compiled from 24 States (excluding Kentucky). The hydrogeologic framework presents a nationwide picture of the glacial aquifer system and provides generalizations concerning the nature of aquifers within it (for example, whether the aquifers are shallow or deep, and unconfined or confined). In this way insights can be gained from understanding the similarities and differences in distinct parts of the glacial aquifer system and how they relate to water use and quality and to aquifer vulnerability. Delineation of hydrogeologic terranes was based on an interpretation of existing geologic mapping of Quaternary sediments and the thickness of unconsolidated material. Overall thickness of Quaternary sediment was used to qualitatively rank the generalized complexity of the hydrogeologic framework in each terrane: “lower” complexity (assigned a terrane code 1), “moderate” complexity (terrane code 2), and “higher” complexity (terrane code 3). Letter designations appended to the terrane codes (for example, 1A, 1B, or 1C) differentiate terranes of similar complexity. Two unique areas, where thick, stratified, coarse-grained sediment dominates, were assigned terrane code 4. Elements of this hydrogeologic framework include a glacial environments and surficial sediments geodatabase, which includes lithologic, geomorphic, and stratigraphic characterization of Quaternary sediments based on previous mapping; a gridded database of sediment and aquifer characteristics computed from lithologic logs obtained from water-well driller records; a water-use database with information on public-water supply systems and sources of groundwater; and estimated recharge computed from a geologically based soil-water balance model. A generalized map of the bedrock geology based on previous State-level mapping is included as well. Quaternary sediment in the glaciated United States includes glacial, postglacial (Holocene) and nonglacial sediments. At land surface, 60 percent of the glacial sediment is till. Large areas of outwash and ice contact sediments are extensive throughout the Midwest but generally are confined to valleys in the Northeast and the Northwest. Lacustrine sediments were deposited in proglacial lakes adjacent to the present Great Lakes and in glacial Lake Agassiz in the eastern Dakotas and northwestern Minnesota. The median thickness of Quaternary sediment ranges from 6 to 45 meters across the 17 hydrogeologic terranes, but the maximum thickness is more than 500 meters in some areas. Quaternary sediments generally contain less than 10 percent coarse material; the median range is near zero percent under till to about 50 percent under ice contact and outwash sediments. About 80 percent of the coarse material lies within 25 to 40 meters of land surface. In most of the glaciated United States, there is a small likelihood of penetrating an aquifer-material interval containing coarse material at least 3 meters thick. A single aquifer-material interval was recorded in about 44 percent of lithologic logs, whereas about 11 percent of the logs penetrated multiple intervals. About 44 percent of water wells in the lithologic database are completed in Quaternary sediment, and many of these Quaternary water wells (42 percent) are confined by at least 7.5 meters of fine materials. About 33 percent of these Quaternary water wells are unconfined—the remainder are where only thin layers (less than 3 meters) of coarse material are present. The median depths of Quaternary water wells range from 13 to 40 meters among the 17 hydrogeologic terranes. Recharge ranges from more than 400 millimeters per year in the Northeast to 11 millimeters or less per year in the Dakotas and Montana (median value of 136 millimeters per year). Annual groundwater withdrawals compiled by county range on an areal basis from less than 1 to 370 millimeters per year, and the mean is 7.4 millimeters per year. About 36 percent of the withdrawals are for public-water supply, of which 70 percent are derived from Quaternary sediments. Groundwater withdrawals are less than 10 percent of recharge throughout most of the glaciated conterminous United States but are a larger proportion of recharge near urban areas in the Northeast and the Midwest, and in counties throughout drier parts of the Midwest. The salient characteristics of the 17 hydrogeologic terranes are presented through maps and a set of descriptive plots to facilitate visual comparisons between selected sediment and aquifer characteristics. The thickness of Quaternary sediment generally increases from the lower complexity terranes through the higher complexity terranes, consistent with their delineation. Median proportions of coarse material in Quaternary sediment and depths to aquifer-material intervals are highly variable (less than 10 to 50 percent, and 0 to 30 meters, respectively). Median thicknesses of aquifer-material intervals generally fall within a narrow range (10 to 20 meters), except in two terranes that contain thick coarse-grained sediment (30 to 35 meters). The source of water in wells varies from mostly bedrock wells in the lower complexity terranes to mostly Quaternary wells in the higher complexity terranes where the sediment is thickest. A tree diagram compiled from a hierarchical cluster analysis of a matrix composed of metrics based on sediment and aquifer characteristics, and the distribution of water wells in each terrane, indicates some groups of terranes that can be treated as comparable when analyzing groundwater flow and quality. Aquifer-material intervals indicated on maps prepared from the lithologic logs, including unconfined and confined conditions, correlate well with aquifer systems delineated on state maps for Illinois, Indiana, and North Dakota. The large scale of the study limits the resolution at which the maps can be interpreted, however, and alluvial units are not mapped correctly for some valleys in the Northeast and the Northwest. Lithologic logs used in the study are biased toward shallow depths because not all logs penetrate the entire thickness of Quaternary sediment, but this bias should not limit the utility of the sediment and aquifer descriptions because shallow depths are commonly exploited for water supply. The hydrogeologic framework will support ongoing studies of groundwater flow and quality in the U.S. Geological Survey National Water Quality Assessment program for the glaciated United States.

Scientific Investigations Report↗