Search USGSSearch

SEARCH · Search USGS

Results for “Scientific Investigations Map”

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,135 records · Page 63Linked to original sources

Flood-inundation maps of the Blue River and major tributaries in and near Kansas City, Missouri, 2023–25

Digital flood-inundation maps for 35.5 miles of the Blue River, in and near Kansas City, Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the City of Kansas City. Maps were also developed for a combined 7.3 miles of Brush Creek (4.4 miles) and Indian Creek (2.9 miles), two primary tributaries to the Blue River in the study reach. The flood-inundation maps, available through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding, depth of water above first-floor elevations in selected buildings, water velocity, water-surface elevations, and time of travel corresponding to selected water levels or gage heights (hereafter referred to as “stages”) at nine reference USGS streamgages—five on the Blue River (USGS station numbers 06893590, 06893578, 06893500, 06893150, and 06893100), two on Brush Creek (USGS station numbers 06893562 and 06893557), one on Indian Creek (USGS station number 06893390), and one on the Missouri River (USGS station number 06893000). Near-real-time stages at these streamgages may be obtained from USGS Water Data for the Nation at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at seven of these sites (USGS station numbers 06893000, 06893590, 06893578, 06893500, 06893150, 06893557, and 06893390). Flood profiles were computed for eight map reaches by means of two-dimensional hydraulic models. The models were calibrated using the USGS-developed stage–streamflow relations at each reference streamgage. Two map reaches on the lower Blue River also include the effects of backwater from the Missouri River on flood inundation extent and hydraulic characteristics. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “Action Stage” or near bankfull to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root-mean-square error) to delineate the area flooded at each water level and the associated hydraulic characteristics. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that could be critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Kansas, Missouri

Summary of Hydrologic Data for the Tuscarawas River Basin, Ohio, with an Annotated Bibliography

The Tuscarawas River Basin drains approximately 2,600 square miles in eastern Ohio and is home to 600,000 residents that rely on the water resources of the basin. This report summarizes the hydrologic conditions in the basin, describes over 400 publications related to the many factors that affect the groundwater and surface-water resources, and presents new water-quality information and a new water-level map designed to provide decisionmakers with information to assist in future data-collection efforts and land-use decisions. The Tuscarawas River is 130 miles long, and the drainage basin includes four major tributary basins and seven man-made reservoirs designed primarily for flood control. The basin lies within two physiographic provinces-the Glaciated Appalachian Plateaus to the north and the unglaciated Allegheny Plateaus to the south. Topography, soil types, surficial geology, and the overall hydrology of the basin were strongly affected by glaciation, which covered the northern one-third of the basin over 10,000 years ago. Within the glaciated region, unconsolidated glacial deposits, which are predominantly clay-rich till, overlie gently sloping Pennsylvanian-age sandstone, limestone, coal, and shale bedrock. Stream valleys throughout the basin are filled with sands and gravels derived from glacial outwash and alluvial processes. The southern two-thirds of the basin is characterized by similar bedrock units; however, till is absent and topographic relief is greater. The primary aquifers are sand- and gravel-filled valleys and sandstone bedrock. These sands and gravels are part of a complex system of aquifers that may exceed 400 feet in thickness and fill glacially incised valleys. Sand and gravel aquifers in this basin are capable of supporting sustained well yields exceeding 1,000 gallons per minute. Underlying sandstones within 300 feet of the surface also provide substantial quantities of water, with typical well yields of up to 100 gallons per minute. Although hydraulic connection between the sandstone bedrock and the sands and gravels in valleys is likely, it has not been assessed in the Tuscarawas River Basin. In 2001, the major land uses in the basin were approximately 40 percent forested, 39 percent agricultural, and 17 percent urban/residential. Between 1992 and 2001, forested land use decreased by 2 percent with correspondingly small increases in agricultural and urban land uses, but from 1980 to 2005, the 13-county area that encompasses the basin experienced a 7.1-percent increase in population. Higher population density and percentages of urban land use were typical of the northern, headwaters parts of the basin in and around the cities of Akron, Canton, and New Philadelphia; the southern area was rural. The basin receives approximately 38 inches of precipitation per year that exits the basin through evapotranspiration, streamflow, and groundwater withdrawals. Recharge to groundwater is estimated to range from 6 to 10 inches per year across the basin. In 2000, approximately 89 percent of the 116 million gallons per day of water used in the basin came from groundwater sources, whereas 11 percent came from surface-water sources. To examine directions of groundwater flow in the basin, a new dataset of water-level contours was developed by the Ohio Department of Natural Resources. The contours were compiled on a map that shows that groundwater flows from the uplands towards the valleys and that the water-level surface mimics surface topography; however, there are areas where data were too sparse to adequately map the water-level surface. Additionally, little is known about deep groundwater that may be flowing into the basin from outside the basin and groundwater interactions with surface-water bodies. Many previous reports as well as new data collected as part of this study show that water quality in the streams and aquifers in the Tuscarawas River Basin has been degraded by urban, suburban, and rural

Scientific Investigations Report

Monitoring changes in Landsat thermal features in urban and non-urban interfaces from 1986 to 2023 in two international urban centers: Implications for climate and global issues

Rapid urbanization is reshaping thermal environments worldwide, with the strongest impacts occurring at the interface between urban and non-urban areas. Impervious surfaces, as key indicators of urban expansion, are critical for monitoring urban growth and assessing surface urban heat island (SUHI) effects. Land use and land cover change (LULCC) provides an essential link between urban dynamics and their environmental and societal consequences. Here, we integrated the U.S. Geological Survey (USGS) Climate Global Issues (CGI) Land Cover Product with Landsat thermal time-series to investigate SUHI evolution in two contrasting metropolitan regions: Wuhan, China, and Brasília, Brazil. Using data spanning 1986–2023, we analyzed the relationships between land cover, Landsat-based land surface temperature (LST), and SUHI intensity, and identified persistent thermal hotspots. Results demonstrate that the land cover data utilized increases the accuracy of impervious surface mapping along urban–rural gradients. Average SUHI intensities were 3.4 °C in Wuhan and 3.3 °C in Brasília, with statistically significant warming trends of 0.04 °C/year and 0.01 °C/year, respectively. Maximum temperature proved to be a robust indicator of SUHI intensification, capturing long-term upward trends. Our findings highlight the important role of urban land cover dynamics in shaping temporal SUHI variability and hotspot emergence. This prototype framework demonstrates the scientific and policy value of combining long-term land cover monitoring information with satellite thermal monitoring to quantify and track SUHI at city scale, supporting sustainable urban planning and climate adaptation strategies.

Remote Sensing

Coastal-change and glaciological maps of Antarctica

Changes in the area and volume of polar ice sheets are intricately linked to changes in global climate, and the resulting changes in sea level may severely impact the densely populated coastal regions on Earth. Melting of the West Antarctic part alone of the Antarctic ice sheet could cause a sea-level rise of approximately 6 meters (m). The potential sea-level rise after melting of the entire Antarctic ice sheet is estimated to be 65 m (Lythe and others, 2001) to 73 m (Williams and Hall, 1993). In spite of its importance, the mass balance (the net volumetric gain or loss) of the Antarctic ice sheet is poorly known; it is not known for certain whether the ice sheet is growing or shrinking. In a review paper, Rignot and Thomas (2002) concluded that the West Antarctic part of the Antarctic ice sheet is probably becoming thinner overall; although the western part is thickening, the northern part is thinning. Joughin and Tulaczyk (2002), based on analysis of ice-flow velocities derived from synthetic aperture radar, concluded that most of the Ross ice streams (ice streams on the east side of the Ross Ice Shelf) have a positive mass balance. The mass balance of the East Antarctic is unknown, but thought to be in near equilibrium. Measurement of changes in area and mass balance of the Antarctic ice sheet was given a very high priority in recommendations by the Polar Research Board of the National Research Council (1986), in subsequent recommendations by the Scientific Committee on Antarctic Research (SCAR) (1989, 1993), and by the National Science Foundation's (1990) Division of Polar Programs. On the basis of these recommendations, the U.S. Geological Survey (USGS) decided that the archive of early 1970s Landsat 1, 2, and 3 Multispectral Scanner (MSS) images of Antarctica and the subsequent repeat coverage made possible with Landsat and other satellite images provided an excellent means of documenting changes in the coastline of Antarctica (Ferrigno and Gould, 1987). The availability of this information provided the impetus for carrying out a comprehensive analysis of the glaciological features of the coastal regions and changes in ice fronts of Antarctica (Swithinbank, 1988; Williams and Ferrigno, 1988). The project was later modified to include Landsat 4 and 5 MSS and Thematic Mapper (TM) (and in some areas Landsat 7 Enhanced Thematic Mapper Plus (ETM+)), RADARSAT images, and other data where available, to compare changes over a 20- to 25- or 30-year time interval (or longer where data were available, as in the Antarctic Peninsula). The results of the analysis are being used to produce a digital database and a series of USGS Geologic Investigations Series Maps consisting of 24 maps at 1:1,000,000 scale and 1 map at 1:5,000,000 scale, in both paper and digital format (Williams and others, 1995; Williams and Ferrigno, 1998; and Ferrigno and others, 2002).

IMAP

Forecasting for a fractured land: A case study of the communication and use of aftershock forecasts from the Mw 7.8 2016 Kaikōura earthquake in Aotearoa New Zealand

Operational earthquake forecasts (OEFs) are represented as time‐dependent probabilities of future earthquake hazard and risk. These probabilities can be presented in a variety of formats, including tables, maps, and text‐based scenarios. In countries such as Aotearoa New Zealand, the U.S., and Japan, OEFs have been released by scientific organizations to agencies and the public, with the intent of providing information about future earthquake hazard and risk, so that people can use this information to inform their decisions and activities. Despite questions being raised about the utility of OEF for decision‐making, past earthquake events have shown that agencies and the public have indeed made use of such forecasts. Responses have included making decisions about safe access into buildings, cordoning, demolition safety, timing of infrastructure repair and rebuild, insurance, postearthquake building standards, postevent land‐use planning, and public communication about aftershocks. To add to this body of knowledge, we undertook a survey to investigate how agencies and GNS Science staff used OEFs that were communicated following the M w 7.8 2016 Kaikōura earthquake in Aotearoa New Zealand. We found that agencies utilized OEFs in many of the ways listed previously, and we document individual employee’s actions taken in their home‐life context. Challenges remain, however, regarding the interpretation of probabilistic information and applying this to practical decision‐making. We suggest that science agencies cannot expect nontechnical users to understand and utilize forecasts without additional support. This might include developing a diversity of audience‐relevant OEF information for communication purposes, alongside advice on how such information could be utilized.

northern part of South Island

Vegetation map of the watersheds between Kawela and Kamalō Gulches, Island of Molokaʻi, Hawaiʻi

In this document we describe the methods and results of a project to produce a large-scale map of the dominant plant communities for an area of 5,118.5 hectares encompassing the Kawela and Kamalō watersheds on the island of Molokaʻi, Hawaiʻi, using digital image analysis of multi-spectral satellite imagery. Besides providing a base map of the area for land managers to use, this vegetation map serves as spatial background for the U.S. Geological Survey’s (USGS) Molokaʻi Ridge-to-Reef project, which is an interdisciplinary study of erosion and sediment transport within these watersheds. A total of 14 mapping units were identified for the Kawela-Kamalō project area. The most widespread units were the ʻŌhiʻa montane wet or mesic forest and No vegetation or very sparse grasses/shrubs communities, each present in more than 800 hectares, or 16 percent of the mapping area. Next largest were the Kiawe woodland with alien grass understory and ʻAʻaliʻi dry shrubland units, each of which covered more than 500 hectares, or more than 12 percent of the area; followed by the Mixed native mesic shrubland, ʻIlima and mixed grass dry shrubland, Mixed alien grass with ʻilima shrubs, and the Mixed alien forest with alien shrub/grass understory communities, which ranged in size from approximately 391 to 491 hectares, or 7.6 to 9.6 percent of the project site. The other six mapped units covered less than 170 hectares of the landscape. Six of the map units were dominated by native vegetation, covering a total of 2,535.2 hectares combined, or approximately 50 percent of the project area. The remaining map units were dominated by nonnative species and represent vegetation types that have resulted from invasion and establishment of plant species that had been either purposely or accidently introduced into Hawaiʻi since humans arrived in these islands more than 1,500 years ago. The preponderance of mapping units that are dominated by alien species of plants is a strong indication of how much anthropogenic disturbance has occurred in this area. The native-dominated ʻŌhiʻa forest and uluhe fern communities are probably most similar to the vegetation that was originally found in the upper part of the project area this area. Portions of the mixed mesic native shrub community still persist in the lowland mesic zone, but below that area, the vegetation is either dominated by alien species, or artificially opened by animal grazing and erosion, even in the few units that are still dominated by native species. The map produced for the Kawela to Kamalō watersheds can be used as a baseline for assessing the distribution and abundance of the various plant communities found across this landscape at the time of the imagery (2004). It can also be used to help understand the dynamics of the vegetation and other attributes of this watershed—such as erosion and surface transport of sediment, relative to current and future habitat conditions.

Hawai'i

Groundwater hydrology and chemistry of Jamestown Island, Virginia—Potential effects of tides, storm surges, and sea-level rise on archaeological, cultural, and ecological resources

As the site of the first permanent English settlement in North America in 1607, Jamestown Island, Colonial National Historical Park (COLO), Virginia, contains a rich archaeological record that extends from the Paleoindian period (15,000 to 8,000 years ago) through the 20th century. The island is located on the lower James River near the mouth of Chesapeake Bay. Jamestown Island vegetation is dominated by upland forests surrounded by tidal, freshwater-to-oligohaline marshes. Along the Virginia coast, relative sea-level rise was more than 2.5 times the global average during the 20th century. Consequently, the National Park Service (NPS) has identified COLO as one of the 25 national parks most threatened by climate change. Surface waters across the island are hydraulically connected to the laterally continuous Surficial aquifer. The land-surface altitude of the island is low, with two-thirds of the island less than 5 feet (ft) above the North American Vertical Datum of 1988 (NAVD 88). Consequently, sea-level rise, combined with tides and storm surges, threatens the island and its resources as surface-water and groundwater levels rise, saltwater enters the Surficial aquifer, and groundwater chemistry changes. The impact of sea-level rise on the island’s surface-water resources has been well studied, but groundwater effects have been largely ignored. Quantifying the effects of tides, storm surges, and sea-level rise on groundwater levels and chemistry is essential to developing an effective strategy for managing climate-induced changes. The first step in developing a response strategy includes a parkwide general risk assessment for archaeological sites on the island, so that sites can be prioritized for management actions. The U.S. Geological Survey and the NPS began a study in 2015 to develop a long-term groundwater-monitoring program to evaluate this risk and to develop an updated management strategy. The groundwater-monitoring program consists of 45 wells and piezometers in two individual clusters and three transects across the island in different hydrologic and chemical settings. Samples for water quality were collected from the wells and piezometers from October 2015 through September 2018 at variable time intervals. Results of the monitoring identified disparate hydrologic and chemical responses to saltwater intrusion across the island. Specific conductance (an indicator of salinity) of groundwater beneath several marshes responded differently to changes in James River salinity. Groundwater response to changes in James River specific conductance appeared to be controlled by land-surface altitude and slope, differences in lateral and vertical sediment characteristics, distance from surface waters, and the degree of surface water/groundwater connectivity between channels and the aquifer. Groundwater chemistry data from monitoring wells at Black Point, a low-altitude, upland setting, are in contrast with conditions observed in Island House observation wells, a high-altitude, upland setting. Specific conductance (less than 200 microsiemens per centimeter [μS/cm]) and pH (greater than 5.0) of groundwater beneath much of the uplands that characterize the Island House observation wells are typical of groundwater in noncarbonate sedimentary aquifers recharged by precipitation. At Black Point, specific conductance ranged from 2,490 to 15,200 μS/cm, and pH ranged from 3.1 to 6.6 standard units. At the Black Point observation wells, the most saline and dense water was at the water table rather than deeper in the aquifer, causing a density inversion that persisted throughout the study. The density inversion likely resulted from differences in permeability between the shallow clay and fine-grained sands and the deeper coarse-grained sand and gravel. Groundwater with the lowest pH was at the water table. As saline groundwater flows through organic sediment beneath the marshes, bacterial biodegradation of organic matter creates anoxic conditions. Continued biodegradation concomitantly reduces iron-oxide minerals in the sediment and sulfate in saline water. When oxygen is reintroduced into groundwater, iron and sulfur can reoxidize to form sulfuric acid, locally lowering the pH of the water. This report describes the groundwater monitoring network design, rationale for site selection, monitoring approach, and results of monitoring from October 2015 through September 2018. Maps of inundation at selected water-level altitudes are included to identify the risk to archaeological, cultural, and ecological resources. The monitoring results of the hydrology and chemistry data are interpreted, and the different hydrologic and chemical settings are described. The implications of the study results for management decisions are presented, and suggestions for improving the monitoring network are included.

Virginia

Coastal-change and glaciological map of the Eights Coast area, Antarctica, 1972-2001

Changes in the area and volume of polar ice sheets are intricately linked to changes in global climate, and the resulting changes in sea level may severely impact the densely populated coastal regions on Earth. Melting of the West Antarctic part alone of the Antarctic ice sheet could cause a sea-level rise of approximately 6 meters (m). The potential sea-level rise after melting of the entire Antarctic ice sheet is estimated to be 65 m (Lythe and others, 2001) to 73 m (Williams and Hall, 1993). In spite of its importance, the mass balance (the net volumetric gain or loss) of the Antarctic ice sheet is poorly known; it is not known for certain whether the ice sheet is growing or shrinking. In a review paper, Rignot and Thomas (2002) concluded that the West Antarctic part of the Antarctic ice sheet is probably becoming thinner overall; although the western part is thickening, the northern part is thinning. Joughin and Tulaczyk (2002), based on analysis of ice-flow velocities derived from synthetic aperture radar, concluded that most of the Ross ice streams (ice streams on the east side of the Ross Ice Shelf) have a positive mass balance. The mass balance of the East Antarctic is unknown, but thought to be in near equilibrium. Measurement of changes in area and mass balance of the Antarctic ice sheet was given a very high priority in recommendations by the Polar Research Board of the National Research Council (1986), in subsequent recommendations by the Scientific Committee on Antarctic Research (SCAR) (1989, 1993), and by the National Science Foundation's (1990) Division of Polar Programs. On the basis of these recommendations, the U.S. Geological Survey (USGS) decided that the archive of early 1970s Landsat 1, 2, and 3 Multispectral Scanner (MSS) images of Antarctica and the subsequent repeat coverage made possible with Landsat and other satellite images provided an excellent means of documenting changes in the coastline of Antarctica (Ferrigno and Gould, 1987). The availability of this information provided the impetus for carrying out a comprehensive analysis of the glaciological features of the coastal regions and changes in ice fronts of Antarctica (Swithinbank, 1988; Williams and Ferrigno, 1988). The project was later modified to include Landsat 4 and 5 MSS and Thematic Mapper (TM) (and in some areas Landsat 7 Enhanced Thematic Mapper Plus (ETM+)), RADARSAT images, and other data where available, to compare changes over a 20- to 25- or 30-year time interval (or longer where data were available, as in the Antarctic Peninsula). The results of the analysis are being used to produce a digital database and a series of USGS Geologic Investigations Series Maps consisting of 24 maps at 1:1,000,000 scale and 1 map at 1:5,000,000 scale, in both paper and digital format (Williams and others, 1995; Williams and Ferrigno, 1998; and Ferrigno and others, 2002).

IMAP

Coastal-change and glaciological map of the Saunders Coast area, Antarctica: 1972-1997

Changes in the area and volume of polar ice sheets are intricately linked to changes in global climate, and the resulting changes in sea level may severely impact the densely populated coastal regions on Earth. Melting of the West Antarctic part alone of the Antarctic ice sheet could cause a sea-level rise of approximately 6 meters (m). The potential sea-level rise after melting of the entire Antarctic ice sheet is estimated to be 65 m (Lythe and others, 2001) to 73 m (Williams and Hall, 1993). In spite of its importance, the mass balance (the net volumetric gain or loss) of the Antarctic ice sheet is poorly known; it is not known for certain whether the ice sheet is growing or shrinking. In a review paper, Rignot and Thomas (2002) concluded that the West Antarctic part of the Antarctic ice sheet is probably becoming thinner overall; although the western part is thickening, the northern part is thinning. Joughin and Tulaczyk (2002), based on analysis of ice-flow velocities derived from synthetic aperture radar, concluded that most of the Ross ice streams (ice streams on the east side of the Ross Ice Shelf) have a positive mass balance. The mass balance of the East Antarctic is unknown, but thought to be in near equilibrium. Measurement of changes in area and mass balance of the Antarctic ice sheet was given a very high priority in recommendations by the Polar Research Board of the National Research Council (1986), in subsequent recommendations by the Scientific Committee on Antarctic Research (SCAR) (1989, 1993), and by the National Science Foundation?s (1990) Division of Polar Programs. On the basis of these recommendations, the U.S. Geological Survey (USGS) decided that the archive of early 1970s Landsat 1, 2, and 3 Multispectral Scanner (MSS) images of Antarctica and the subsequent repeat coverage made possible with Landsat and other satellite images provided an excellent means of documenting changes in the coastline of Antarctica (Ferrigno and Gould, 1987). The availability of this information provided the impetus for carrying out a comprehensive analysis of the glaciological features of the coastal regions and changes in ice fronts of Antarctica (Swithinbank, 1988; Williams and Ferrigno, 1988). The project was later modified to include Landsat 4 and 5 MSS and Thematic Mapper (TM) (and in some areas Landsat 7 Enhanced Thematic Mapper Plus (ETM+)), RADARSAT images, and other data where available, to compare changes over a 20- to 25- or 30-year time interval (or longer where data were available, as in the Antarctic Peninsula). The results of the analysis are being used to produce a digital database and a series of USGS Geologic Investigations Series Maps consisting of 24 maps at 1:1,000,000 scale and 1 map at 1:5,000,000 scale, in both paper and digital format (Williams and others, 1995; Williams and Ferrigno, 1998; and Ferrigno and others, 2002).

IMAP

Water-Level Surface Maps of the Carbonate-Rock and Basin-Fill Aquifers in the Basin and Range Carbonate-Rock Aquifer System, White Pine County, Nevada, and Adjacent Areas in Nevada and Utah

Water-level measurements in 418 wells were used to develop a potentiometric-surface map of the carbonate-rock aquifer and a water-table map of the basin-fill aquifer in the Basin and Range carbonate-rock aquifer system (BARCAS) study area. The BARCAS study area encompasses about 13,500 square miles and includes most of White Pine County, Nevada, and smaller areas of adjacent counties in Nevada and Utah. Current and historical data from the U.S. Geological Survey National Water Information System, previous publications, and field reconnaissance were used to define water-level surfaces.

Scientific Investigations Report

Simulating groundwater flow in karst aquifers with distributed parameter models—Comparison of porous-equivalent media and hybrid flow approaches

Understanding karst aquifers, for purposes of their management and protection, poses unique challenges. Karst aquifers are characterized by groundwater flow through conduits (tertiary porosity), and (or) layers with interconnected pores (secondary porosity) and through intergranular porosity (primary or matrix porosity). Since the late 1960s, advances have been made in the development of numerical computer codes and the use of mathematical model applications towards the understanding of dual (primary [matrix] and secondary [fractures and conduits]) porosity groundwater flow processes, as well as characterization and management of karst aquifers. The Floridan aquifer system (FAS) in Florida and parts of Alabama, Georgia, and South Carolina is composed of a thick sequence of predominantly carbonate rocks. Karst features are present over much of its area, especially in Florida where more than 30 first-magnitude springs occur, numerous sinkholes and submerged conduits have been mapped, and numerous circular lakes within sinkhole depressions are present. Different types of mathematical models have been applied for simulation of the FAS. Most of these models are distributed parameter models based on the assumption that, like a sponge, water flows through connected pores within the aquifer system and can be simulated with the same mathematical methods applied to flow through sand and gravel aquifers; these models are usually referred to as porous-equivalent media models. The partial differential equation solved for groundwater flow is the potential flow equation of fluid mechanics, which is used when flow is dominated by potential energy and has been applied for many fluid problems in which kinetic energy terms are dropped from the differential equation solved. In many groundwater model codes (basic MODFLOW), it is assumed that the water has a constant temperature and density and that flow is laminar, such that kinetic energy has minimal impact on flow. Some models have been developed that incorporate the submerged conduits as a one-dimensional pipe network within the aquifer rather than as discrete, extremely transmissive features in a porous-equivalent medium; these submerged conduit models are usually referred to as hybrid models and may include the capability to simulate both laminar and turbulent flow in the one-dimensional pipe network. Comparisons of the application of a porous-equivalent media model with and without turbulence (MODFLOW-Conduit Flow Process mode 2 and basic MODFLOW, respectively) and a hybrid (MODFLOW-Conduit Flow Process mode 1) model to the Woodville Karst Plain near Tallahassee, Florida, indicated that for annual, monthly, or seasonal average hydrologic conditions, all methods met calibration criteria (matched observed groundwater levels and average flows). Thus, the increased effort required, such as the collection of data on conduit location, to develop a hybrid model and its increased computational burden, is not necessary for simulation of average hydrologic conditions (non-laminar flow effects on simulated head and spring discharge were minimal). However, simulation of a large storm event in the Woodville Karst Plain with daily stress periods indicated that turbulence is important for matching daily springflow hydrographs. Thus, if matching streamflow hydrographs over a storm event is required, the simulation of non-laminar flow and the location of conduits are required. The main challenge in application of the methods and approaches for developing hybrid models relates to the difficulty of mapping conduit networks or having high-quality datasets to calibrate these models. Additionally, hybrid models have long simulation times, which can preclude the use of parameter estimation for calibration. Simulation of contaminant transport that does not account for preferential flow through conduits or extremely permeable zones in any approach is ill-advised. Simulation results in other karst aquifers or other parts of the FAS may differ from the comparison demonstrated herein.

Florida

Flood-inundation maps for Fourmile Creek at Silver Grove, Kentucky

Digital flood-inundation maps for a 3.4-mile reach of Fourmile Creek at Silver Grove, Kentucky, were created by the U.S. Geological Survey (USGS) in cooperation with the City of Silver Grove and the U.S. Army Corps of Engineers Louisville District. Because the City of Silver Grove is subject to flooding from Fourmile Creek and the Ohio River (backwater flooding up Fourmile Creek), a set of flood-inundation maps was created, including maps for each flooding source considered independently and for possible scenarios involving flooding from both sources combined. The flood-inundation maps depict estimates of the areal extent and depth of flooding corresponding to a range of gage heights (gage height is commonly referred to as “stage,” or the water-surface elevation at a streamgage) at the USGS streamgage on Fourmile Creek at Grays Crossing at Silver Grove, Ky. (station number 03238785), and the USGS streamgage on Fourmile Creek at Highway 8 at Silver Grove, Ky. (station number 03238798). Near-real-time stages at these streamgages can be obtained from the USGS National Water Information System at https://waterdata.usgs.gov/ . The USGS streamgage on the Ohio River at Cincinnati, Ohio (station number 03255000), is also important in this study because the National Weather Service (NWS) Advanced Hydrologic Prediction Service (AHPS; https://water.weather.gov/ahps/ ) forecasts flood hydrographs for this site (NWS AHPS site CCNO1). The peak-stage information forecast by the NWS AHPS can be used in conjunction with the maps developed in this study to show predicted areas of flood inundation. Flood profiles were computed for the Fourmile Creek study reach by means of a one-dimensional, step-backwater hydraulic model (HEC-RAS) developed by the U.S. Army Corps of Engineers. The hydraulic model was calibrated by using the current stage-discharge relation (USGS rating number 1.1) at USGS streamgage 03238785, Fourmile Creek at Grays Crossing at Silver Grove, Ky. The model was then used to compute water-surface profiles for 83 combinations of flood stages on the Ohio River and Fourmile Creek ranging from approximately base flow to greater than a 2-percent annual exceedance probability flood in the model reach. An additional 50 water-surface profiles were computed for backwater-only flooding (from the Ohio River) for flood elevations (referenced to the North American Vertical Datum of 1988 [NAVD 88]) at 1-foot intervals referenced to USGS streamgage 03238798, Fourmile Creek at Highway 8 at Silver Grove, Ky.; these elevations ranged from approximately normal pool (460 ft, NAVD 88) to approximately a 0.2-percent annual exceedance probability flood (509 ft, NAVD 88) on the Ohio River. The computed water-surface profile information was then combined with a digital elevation model derived from light detection and ranging (lidar) data to delineate the approximate flooded areas. The digital flood-inundation maps are available through the USGS Flood Inundation Mapper application ( https://fim.wim.usgs.gov/fim/ ), which presents map libraries and provides detailed information on flood extent and depths for selected sites. The flood-inundation maps developed in this study, in conjunction with the real-time stage data from the USGS streamgages on Fourmile Creek at Silver Grove, Ky., and forecasted stream stages from the NWS AHPS, are intended to provide information that can help inform the public about potential flooding and provide emergency management personnel with a tool to efficiently manage emergency flood operations, such as evacuations and road closures, and assist in postflood recovery efforts.

Kentucky

Coastal-Change and Glaciological Map of the Larsen Ice Shelf Area, Antarctica, 1940-2005

Changes in the area and volume of polar ice sheets are intricately linked to changes in global climate, and the resulting changes in sea level could severely impact the densely populated coastal regions on Earth. Antarctica is Earth's largest reservoir of glacial ice. Melting of the West Antarctic part alone of the Antarctic ice sheet would cause a sea-level rise of approximately 6 meters (m), and the potential sea-level rise after melting of the entire Antarctic ice sheet is estimated to be 65 m (Lythe and others, 2001) to 73 m (Williams and Hall, 1993). The mass balance (the net volumetric gain or loss) of the Antarctic ice sheet is highly complex, responding differently to different climatic and other conditions in each region (Vaughan, 2005). In a review paper, Rignot and Thomas (2002) concluded that the West Antarctic ice sheet is probably becoming thinner overall; although it is known to be thickening in the west, it is thinning in the north. The mass balance of the East Antarctic ice sheet is thought by Davis and others (2005) to be positive on the basis of the change in satellite-altimetry measurements made between 1992 and 2003. Measurement of changes in area and mass balance of the Antarctic ice sheet was given a very high priority in recommendations by the Polar Research Board of the National Research Council (1986), in subsequent recommendations by the Scientific Committee on Antarctic Research (SCAR) (1989, 1993), and by the National Science Foundation's (1990) Division of Polar Programs. On the basis of these recommendations, the U.S. Geological Survey (USGS) decided that the archive of early 1970s Landsat 1, 2, and 3 Multispectral Scanner (MSS) images of Antarctica and the subsequent repeat coverage made possible with Landsat and other satellite images provided an excellent means of documenting changes in the cryospheric coastline of Antarctica (Ferrigno and Gould, 1987). The availability of this information provided the impetus for carrying out a comprehensive analysis of the glaciological features of the coastal regions and changes in ice fronts of Antarctica (Swithinbank, 1988; Williams and Ferrigno, 1988). The project was later modified to include Landsat 4 and 5 MSS and Thematic Mapper (TM) images [and in some areas Landsat 7 Enhanced Thematic Mapper Plus (ETM+) images], RADARSAT images, aerial photography, and other data where available, to compare changes that occurred during a 20- to 25- or 30-year time interval (or longer where data were available, as in the Antarctic Peninsula). The results of the analysis are being used to produce a digital database and a series of USGS Geologic Investigations Series Maps (I-2600) (Williams and others, 1995; Ferrigno and others, 2002; and Williams and Ferrigno, 2005) (available online at http://www.glaciers.er.usgs.gov).

IMAP

Coastal-Change and Glaciological Map of the Palmer Land Area, Antarctica: 1947-2009

Reduction in the area and volume of the two polar ice sheets is intricately linked to changes in global climate, and the resulting rise in sea level could severely impact the densely populated coastal regions on Earth. Antarctica is Earth's largest reservoir of glacial ice. Melting of the West Antarctic part alone of the Antarctic ice sheet would cause a sea-level rise of approximately 6 meters (m), and the potential sea-level rise after melting of the entire Antarctic ice sheet is estimated to be 65 m (Lythe and others, 2001) to 73 m (Williams and Hall, 1993). The mass balance (the net volumetric gain or loss) of the Antarctic ice sheet is highly complex, responding differently to different climatic and other conditions in each region (Vaughan, 2005). In a review paper, Rignot and Thomas (2002) concluded that the West Antarctic ice sheet is probably becoming thinner overall; although it is known to be thickening in the west, it is thinning in the north. The mass balance of the East Antarctic ice sheet is thought by Davis and others (2005) to be positive on the basis of the change in satellite-altimetry measurements made between 1992 and 2003. Measurement of changes in area and mass balance of the Antarctic ice sheet was given a very high priority in recommendations by the Polar Research Board of the National Research Council (1986), in subsequent recommendations by the Scientific Committee on Antarctic Research (SCAR) (1989, 1993), and by the National Science Foundation's (1990) Division of Polar Programs. On the basis of these recommendations, the U.S. Geological Survey (USGS) decided that the archive of early 1970s Landsat 1, 2, and 3 Multispectral Scanner (MSS) images of Antarctica and the subsequent repeat coverage made possible with Landsat and other satellite images provided an excellent means of documenting changes in the cryospheric coastline of Antarctica (Ferrigno and Gould, 1987). The availability of this information provided the impetus for carrying out a comprehensive analysis of the glaciological features of the coastal regions and changes in ice fronts of Antarctica (Swithinbank, 1988; Williams and Ferrigno, 1988). The project was later modified to include Landsat 4 and 5 MSS and Thematic Mapper (TM) images (and in some areas Landsat 7 Enhanced Thematic Mapper Plus (ETM+) images), RADARSAT images, aerial photography, and other data where available, to compare changes that occurred during a 20- to 25- or 30-year time interval (or longer where data were available, as in the Antarctic Peninsula). The results of the analysis are being used to produce a digital database and a series of USGS Geologic Investigations Series Maps (I-2600) (Williams and others, 1995; Swithinbank and others, 2003a,b, 2004; Ferrigno and others, 2002, 2005, 2006, 2007, 2008, and in press; and Williams and Ferrigno, 2005) (available online at http://www.glaciers.er.usgs.gov).

IMAP

Preliminary report on ground water in the Michaud Flats Project, Power County, Idaho

The Michaud Flats Project area, as here described, includes about 65 square miles in central Power County, south of the Snake River in the southeastern Snake River Plain of Idaho. The principal town and commercial center of the area is American Falls. The immediate purpose of work in the area by the U.S. Geological Survey was to investigate the possibility of developing substantial quantities of ground water for irrigating high and outlying lands in the proposed Michaud Flats Project area of the U.S. Bureau of Reclamation. Initial findings are sufficiently favorable to warrant comprehensive further investigation. Advanced study would assist proper utilization of ground-water resources and would aid ultimate evaluation of total water resources available in the area. About 10,000 acres of low-lying lands in the Michaud Flats project could be irrigated with water from the Snake River under a low-line distribution system involving a maximum pumping lift of about 200 feet above the river. An additional larger area of high and outlying lands is suitable for irrigation with water pumped from wells. If sufficient ground water is economically available, the expense of constructing and operating a costly highline distribution system for surface water could be saved. Reconnaissance of the ground-water geology of the area disclosed surface outcrops of late Cenozoic sedimentary, pyroclastic, and volcanic rocks. Well logs and test borings show that similar materials are present beneath the land surface in the zone of saturation. Ground water occurs under perched, unconfined, and confined (artesian) conditions, but the aquifers have not been adequately explored. Existing irrigation wells, 300 feet or less in depth, yield several hundred to 1,400 gallons of water a minute, with pumping drawdowns of 6 to 50 feet, and perhaps more. A few wells have been pumped out at rates of less than 800 gallons a minute. Scientific well-construction and development methods would lead to more efficient well performance. A generalized water-table contour map of the area shows that the principal general direction of ground-water movement is toward the west and northwest. The southwestern part of the American Falls Reservoir, and a segment of the Snake River below the dam, may be perched above the water table. Ground water appears to move beneath this segment of the river to the Snake River Plain on the northwest side. So far as is known, recharge to the ground-water reservoir is chiefly from local sources and from the runoff from the mountain area southeast of the project. Seepage losses from surface water spread for irrigation would contribute a substantial amount of new recharge to the ground water, but the amount of such recharge might be less than the depletion of ground water by pumping. Therefore, with ground-water irrigation a part of the project, return flow to the American Falls Reservoir might be less than it is in the existing regimen. Ground-water pumping where the ground water is not tributary to the reservoir might not deplete the reservoir appreciably, but would reduce the net supply of water available west of Neeley. Evidence indicates that at least moderate supplies of ground water can be obtained in low-lying areas southwest and northeast of American Falls, but the safe perennial yields of the aquifers cannot now be estimated. The ground-water potential in high and outlying lands is not known. It is unlikely that this potential is sufficient to supply all high and outlying lands, but the supply may be adequate for a substantial part of these lands. Thorough investigation appears to be warranted.

Open-File Report

Streams, springs, and volcanic lakes for volcano monitoring

Introduction Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume, on boreholes; Hurwitz and Lowenstern, 2024). Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO 2 /SO 2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data. Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php ). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024). Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021). Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO 2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998). If magma erupts or flows into a lake, the interaction between hot magma and cold water can be explosive (Mastin and others, 2004; Zimanowski and others, 2015) and substantially expand the area affected by the eruption. Another hazard is the breaching of crater rims by landslides triggered by volcanic and (or) seismic activity. Under some circumstances, substantial volumes of water can be displaced, leading to large floods and lahars. Late Holocene lake flooding from Aniakchak Crater in the Alaska Peninsula (Waythomas, 2022) and from Paulina Lake in Newberry Crater, Oregon (Chitwood and Jensen, 2000), caused by the failure of outlet sills, testify to the substantial hazards at lake-filled calderas. Several volcanic systems in the United States host lakes known to receive heat and gas from underlying magma. These lakes vary widely in area, depth, and chemical composition. Lakes are present at level 4 volcanoes, including Crater Lake and Newberry Volcano in Oregon; Yellowstone Caldera in Wyoming; Long Valley Caldera, Clear Lake volcanic field, Medicine Lake, and Salton Buttes in California; and Aniakchak Crater, Mount Katmai, Fisher Caldera, Mount Okmok, and Kaguyak Crater, among others, in Alaska. A water lake was present in Halemaʻumaʻu, the crater of Kīlauea, Hawai‘i (fig. F1), from October 2019 to December 2020. Level 3 volcanoes with lakes include Mono Lake volcanic field (Calif.), Mount Bachelor (Ore.), Ukinrek Maars and Mount Chiginagak (Alaska), and Soda Lake (Nevada). In addition, there are lakes at many levels 1 and 2 volcanoes. In the United States, there are no strongly acidic lakes that receive abundant input of magmatic gas, such as those found at Mount Ruapehu (New Zealand), Ijen and Kelud (Indonesia), and Poás (Costa Rica). Nevertheless, many contain fluids that provide clues to magmatic processes below. Since publication of a previous report on recommended instrumentation for volcano monitoring (Moran and others, 2008), continuous hydrologic monitoring has become increasingly feasible. However, changes in water pressure, temperature, and chemistry remain, in general, poorly studied phenomena at volcanoes (Sparks, 2003; National Academies of Sciences, Engineering, and Medicine, 2017). Recent efforts by the USGS have included the temporary study of Cascade Range volcanoes, which included frequent (15 minute to hourly) temporal sampling of temperature, depth, and conductivity (Crankshaw and others, 2018; Ingebritsen and Evans, 2019). At Yellowstone Caldera, many streamgages have now added thermistors and specific conductance sensors, allowing estimation of time-dependent chloride flux as a proxy for variations in subsurface heat flux (McCleskey and others, 2012, 2016). Efforts to better understand lakes have also accelerated, with bathymetric mapping and sampling carried out at several locations in the United States. Especially thorough work was done at Yellowstone Lake thanks to the Hydrothermal Dynamics of Yellowstone Lake (HD-YLAKE, https://hdylake.org ) project, funded primarily by the National Science Foundation. In addition to geophysical surveys and recovery of cores and other samples, HD-YLAKE investigations included remotely operated vehicle (ROV) investigations of hydrothermal vents on the lake floor (fig. F2). Data collected by the ROV provided a better understanding of the thermal and chemical influx from lake-bottom hydrothermal systems (Sohn and others, 2017). In this chapter, we focus on detecting changes in the chemistry, temperature, discharge, or water levels of streams, springs, and lakes that can be caused by seismicity, volumetric strains, or increases in gas flux associated with ascending magma. There is unavoidable overlap with other chapters of this report. Samples of water and gas can also be obtained in boreholes ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume; Hurwitz and Lowenstern, 2024), both shallow and deep. Gas monitoring ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–E, 11 p., https://doi.org/10.3133/sir20245062e.">chapter E , this volume; Lewicki and others, 2024) relies in part on samples from springs and wells, particularly where measurable gas plumes are absent. Water acts as a trigger and lubricant for landslides and sediment-rich floods, and so hydrology has obvious relevance for lahar monitoring, as discussed in of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–H, 6 p., https://doi.org/10.3133/sir20245062h. ">chapter H (this volume; Thelen and others, 2024). Shared situational awareness among scientists engaged in geophysical, gas, and hydrologic monitoring will improve overall understanding of the volcanic hazard.

Scientific Investigations Report

Reservoir evolution, downstream sediment transport, downstream channel change, and synthesis of geomorphic responses of Fall Creek and Middle Fork Willamette River to water years 2012–18 streambed drawdowns at Fall Creek Lake, Oregon

Executive Summary Chapter A. Introduction Fall Creek Dam impounds Fall Creek Lake, a 10-kilometer-long reservoir in western Oregon and is operated by the U.S. Army Corps of Engineers (USACE) primarily for flood-risk management (or flood control) in late autumn through early spring months, as well as for water quality, irrigation, recreation, and habitat in late spring through early autumn. Since 2011 (water year [WY] 2012), Fall Creek Lake has been temporarily drawn down each year to facilitate downstream passage of juvenile spring Chinook salmon ( Oncorhynchus tshawytscha ) through the 55-meter (m) high dam. This annual dam operation is temporary, typically lasting about 1–2 weeks from WY 2012 through 2020 (drawdown operations in WY 2022–24 have increased to more than 6 weeks). Drawdown of the reservoir results in lake levels being lowered to the elevation near the historical, pre-dam streambed. The annual streambed drawdowns of WY 2012–18 have improved fish passage and led the USACE to formally adopt streambed drawdowns as part of annual operations at Fall Creek Dam. However, temporarily lowering the lake to streambed creates free-flowing conditions in the reservoir that result in the erosion and episodic export of predominantly sand and finer-grained sediments (less than 2 millimeters [mm]) to the lower gravel-bed reaches of Fall Creek and the Middle Fork Willamette River. The introduction of large volumes of sand and finer-grain sediment into the dam-regulated reaches downstream from Fall Creek Dam prompted questions about the geomorphic responses to annual streambed drawdowns within Fall Creek Lake and downstream reaches along Fall Creek and the Middle Fork Willamette River. The U.S. Geological Survey (USGS) in partnership with USACE initiated a comprehensive geomorphic and sediment transport investigation to assess the coupled processes of reservoir erosion, sediment evacuation from Fall Creek Lake, and patterns of sediment transport and deposition in reaches downstream from the Fall Creek Dam that have resulted from annual streambed drawdowns. The purpose of this report is to systematically describe the processes of sediment erosion, transport, and deposition at Fall Creek Lake and geomorphic interactions between reaches upstream and downstream from Fall Creek Dam that relate to dam operations. Specifically, this report focuses on evaluating geomorphic responses to streambed drawdowns from WY 2012 through 2018 and placing drawdown-induced geomorphic responses within the broader context of physiographic and historical conditions and dam operations of Fall Creek and Middle Fork Willamette Rivers. Key objectives for this study were to characterize changes in reservoir morphology and substrate at Fall Creek Lake, describe the character and temporal pattern of sediment transport downstream from Fall Creek Dam, characterize geomorphic changes in channel reaches downstream from the Fall Creek Dam, and relate these data to the annual streambed drawdowns of WY 2012–18. This study uses multiple independent monitoring and measurement approaches to assess site, reach, and river-scale geomorphic responses to drawdowns to inform dam and reservoir management. Patterns and processes of reservoir evolution were assessed with geomorphic mapping and volumetric analyses of topography through comparison of multiple digital surface models (DSMs). Just downstream from Fall Creek Dam, analyses of sediment export from the reservoir focused on suspended sediment but also incorporated bedload analyses to assess sediment sizes. Geomorphic assessments downstream from the dam used reach-scale and site-scale approaches to document changes in channel morphology and substrate, including site measurements of sand and finer-grained sediment deposition and in-channel bed-material, volumetric change analyses from comparison of digital elevation models (DEMs), and repeat geomorphic mapping. Findings from this study inform river management and dam operations by providing an understanding of (1) coupled upstream-downstream geomorphic responses to the Fall Creek Lake streambed drawdowns, (2) geomorphic responses of Fall Creek Lake streambed drawdowns in comparison to drawdowns at other large dams, (3) controls on reservoir erosion and downstream geomorphic responses, and (4) implications for future hydrogeomorphic changes that may result from continued drawdowns and monitoring activities to assess those changes. Chapter B. Reservoir Morphology and Evolution Related to Dam Operations at Fall Creek Lake To understand the volume and distribution of sediment accumulation in Fall Creek Lake since dam closure in 1965, decadal-scale sedimentation patterns (spanning approximately 1965–2016) are evaluated using a combination of storage curve analyses and geomorphic mapping. Short-term (drawdown event-scale) patterns of erosion, sedimentation, and sediment export downstream are evaluated using a combination of geomorphic mapping and change detection analyses that quantify the distribution and total volume of sediment erosion and deposition within Fall Creek Lake. Geomorphic mapping of reservoir topography and analyses of historical datasets reveals four categories of landforms and sediment processes within Fall Creek Lake related to lake level operations: lacustrine sedimentation expressed in the reservoir floor, fluvial erosion and deposition within historical stream channels during streambed drawdowns, channel-like features created by erosion within the reservoir floor during streambed drawdowns, and erosion on reservoir hillslopes. Where the reservoir floor is mapped for this study as pelagic (deep water), deposition up to 3 meters (m) thick by lacustrine processes and burial of pre-dam topography with deposits thinning toward the edges of the valley floor and upstream areas of reservoir are observed. Despite over 50 years of sediment accumulation since dam construction, the main stream channels of Fall and Winberry Creeks (or reservoir thalwegs) through the reservoir are well defined, though their distinct morphology is likely influenced by a long history of recurring historical drawdowns to or near streambed since dam construction. Unregulated streamflow and sediment transport through the reservoir primarily are confined to these channels during the streambed drawdown periods. Erosional channel-like features created by drawdowns are carved through underlying, unconsolidated reservoir floor sediments and are most prominent in the lower reservoir below minimum conservation pool (the low pool elevation during winter flood season); sediment generated from the formation of these drawdown channels is more likely to be transported through and out of the reservoir than sediment deposits along the reservoir hillslopes at the valley margins that are separated from main channels by areas of low-gradient reservoir floor. Morphologic changes in the lower reservoir topography between January 2012 and November 2016 indicate overall net erosion of about 129,500 cubic meters (m 3 ). The most prominent geomorphic changes occurred along the main channels of Fall and Winberry Creeks near the Fall Creek Dam where incision, lateral migration, and slumping banks resulted in vertical and lateral adjustments to channel position, whereas most changes fell below the detectable limit on higher-elevation reservoir floor surfaces except where erosion occurred along features mapped as drawdown channels. Chapter C. Sediment Delivery from Fall Creek Lake and Transport through Downstream Reaches USGS implemented a sediment monitoring program in WY 2013–18 to evaluate the quantity and character of reservoir sediment exported from Fall Creek Lake during streambed drawdowns. Turbidity and suspended sediments were monitored annually autumn through spring to span the WY 2013–18 streambed drawdowns; however, unequal monitoring timeframes each year reduced the ability to compare results and factors affecting sediment export from the reservoir difficult between years. These data were originally measured to develop regressions and compute suspended-sediment loads (SSL). Bedload sediment monitoring from a cableway at the Fall Creek streamgage was completed in the autumn-winter of WY 2013 and 2017. The limited number of samples and presumed variability in sediment supply from the reservoir precluded construction of streamflow and bedload discharge relations to compute more than instantaneous bedload. Sand and finer-grained silts and clays were transported from the reservoir in suspension, though some coarser grains (up to 32 mm) were also mobilized and transported downstream from the dam as bedload. Observations of increased sediment transport downstream from Fall Creek Dam coincided with lake levels approaching about 3 m (10 feet [ft] or elevation 690 ft) above the streambed regulating outlets. Suspended-sediment loads computed for the full monitoring periods WY 2013–18 at the Fall Creek streamgage, located 1.4 kilometers (km) downstream from Fall Creek Dam, range from 54,700 metric tons (t) in WY 2013 to 13,900 t in WY 2018. Although the total annual SSL varied from year to year, the overall seasonal patterns of suspended sediment transport throughout each year were similar during monitoring in WY 2013-18. Suspended-sediment loads were low prior to the drawdown, then increased rapidly as lake levels lowered and approached the streambed. In the weeks following the drawdown period, as pool levels were increased, SSL remained slightly elevated above pre-drawdown levels but generally declined through the following winter and spring except during streamflow-driven pulses of suspended-sediment transport. WY 2013 had the greatest total computed SSL for each streambed drawdown and partial-year monitoring period. SSL computed for the partial-year period have generally decreased since WY 2013 and have varied by about 6,800 t with the exception of WY 2014. WY 2014 SSL reflects anomalously low sediment export due to low streamflows and freezing conditions that stabilized reservoir floor deposits. Bedload measurements in the short 1.4-km reach between Fall Creek Dam and the Fall Creek streamgage showed an inverse correlation between bedload transport rates and discharge, which probably reflects diminishing supply of coarse-sized sediment. Sand was more abundant (60–100 percent) than gravel in bedload samples confirming sand and finer-grained sediment dominated sediment evacuated from the reservoir during streambed drawdowns at Fall Creek Lake. Chapter D. Geomorphic Responses to Fall Creek Lake Streambed Drawdowns Downstream from Fall Creek Dam In the days, weeks, and months following streambed drawdown operations at Fall Creek Dam through WY 2018, sites downstream from the dam displayed a variety of geomorphic responses to reservoir sediment delivery within the main channel and overbank areas. Evaluation of streambed elevations at two streamgages located 1.4 km downstream from the dam on Fall Creek and 16.3 km downstream from the dam on the Middle Fork Willamette River indicated the effects of drawdown sediment on bed elevations were modest and transient. Repeat particle size measurements (October 2015 and September 2016) at five sites along Fall Creek and the Middle Fork Willamette River showed similar grain-sized distributions that do not reveal substantial deposition of fine-grained sediment related to the WY 2016 streambed drawdown. Altogether, these findings indicate that transport capacity in the main, low-flow channels of Fall Creek and Middle Fork Willamette River during WY 2012–18 was sufficient to mobilize and evacuate reservoir sediments from streambed drawdowns or other bank material and tributary sources. However, other monitoring for this study indicate low-velocity zones in off-channel areas are prime locations for sand and finer-grain sediment deposition. Patterns of overbank sediment accumulation indicate that the magnitude and timing of overbank deposition on bars and low-elevation floodplain varies with proximity to the dam, geomorphic setting, streamflows, and other factors. Sand and finer-grained reservoir sediments carried as suspended-sediment load in the reaches downstream from Fall Creek Dam were deposited in overbank areas as observed with clay-horizon markers during WY 2016–17. Overbank deposition quantified with Geomorphic Change Detection (GCD) software evaluated landform-scale patterns of erosion and deposition using repeat light detection and ranging (lidar) surveys at two sites in the Upper Fall Creek reach and one site in the Jasper reach for 3 years (2012–15) and one site in the Clearwater reach for 6 years (2009–15). Deposition thickness and spatial patterns from the GCD analysis were variable; some sites had dispersed but measurable deposition while at others, deposition was highly localized and exceeded 1 m in depth. Patterns of overbank deposition illustrate interactions among bar morphology, local hydraulics, and suspended-sediment transport dynamics that can create patches of highly localized deposition. The measured deposition at the two Fall Creek GCD sites likely resulted from reservoir sediments released from Fall Creek Lake during streambed drawdowns in WY 2016 and 2017 because the limited sediment inputs from bank material (geomorphically laterally stable reach) or tributaries (no significant tributaries) provided few other sediment sources. On the Middle Fork Willamette River, observed patterns of overbank deposition could reflect sediment sourced from upstream tributaries, bank erosion, or Fall Creek Lake streambed drawdown operations. Despite the introduction of several thousand tons of reservoir sediment delivered from the Fall Creek Lake streambed drawdowns to below-dam river corridors, reach-scale mapping of channel features downstream from Fall Creek Dam shows minimal evidence of changes in channel planform or landforms that can be attributed to a drawdowns in WY 2012–16. On Upper Fall Creek reach, widespread increases in gravel bars or other in-channel sediment did not result from the five streambed drawdowns. The main changes attributable to sediment releases from Fall Creek Lake were localized increases in vegetated bar area, particularly on channel margin areas where sand and finer-grain sediment was deposited and rapidly colonized by vegetation. The area of mapped secondary water features decreased between 2005 and 2016, but that may be due to lower discharges depicted in the 2016 aerial photographs and less mapped area of inundation. Primary changes along the Lower Fall Creek reach include a 6.4 percent decrease in area of secondary water features between 2011 and 2016, and a nearly twofold increase in the area of unvegetated bars. Immediately downstream from the Fall Creek confluence, there were negligible changes in the location and areas of vegetated bars and the main wetted channel between 2005 and 2016, and local increases in bar area cannot be attributed solely to deposition of reservoir sediments from Fall Creek Lake because (1) areas along the Middle Fork Willamette River just upstream from the Fall Creek confluence display similar type and magnitude of changes and (2) some of the increases at the confluence area pre-date the drawdowns. The cumulative effect of sediment releases from Fall Creek Lake streambed drawdowns from WY 2012 to 2016 on downstream channel planform and landforms are modest compared to the river-scale transformations and planform changes that occurred in the decades following dam construction. Chapter E. Discussion of Geomorphic Responses of Fall Creek and Middle Fork Willamette River to Streambed Drawdowns at Fall Creek Lake Multiple aspects of Fall Creek Dam infrastructure and operations exert first-order controls on the magnitudes of reservoir erosion that occur during the streambed drawdowns and ultimately determine the sediment delivery to downstream reaches. Key aspects of the dam and its operations that are most relevant to assessing geomorphic responses to streambed drawdowns include the (1) dam infrastructure, including configuration and size of regulating outlets and their proximity to the streambed which dictates the capacity and competence of the river to deliver sediment to downstream reaches and mode of sediment transport as suspended-sediment load or bedload; (2) frequency of historical drawdowns and long-term, year-round dam operations and lake level management, which partly dictate reservoir morphology and locations and magnitudes of readily erodible materials; (3) dam operations and hydroclimatic conditions during the streambed drawdown (including length of the drawdown and streamflows entering the reservoir), which directly control the timing, duration and magnitude of reservoir erosion and sediment evacuation; and (4) dam operations following the streambed drawdown operation that regulate streamflows (and thereby sediment transport conditions) downstream of Fall Creek Dam which primarily reflect interactions between hydroclimatic conditions and flood control operations. Patterns of sediment erosion and evacuation observed in this study at Fall Creek Lake from WY 2012–18 suggest that reservoir erosion during annual streambed drawdowns may remain similar or decrease in future years assuming (1) annual streambed drawdown operations are implemented in similar manner as the WY 2012–18 drawdowns (in terms of duration, late autumn or early winter implementation, rate of pool-level lowering to reach streambed, and other factors), (2) streambed drawdowns coincide with similar conditions as were observed WY 2012–18 (similar sediment yield into reservoir, low reservoir inflows, limited precipitation, moderate air temperature), and (3) no major geomorphic changes in the main reservoir channels of Fall and Winberry Creeks occur (for example, channel avulsion). Under such conditions, it is hypothesized that the stream channel within the reservoir would achieve a quasi-equilibrium state with respect to annual influx and export of sediment and aided by the substantial amount of in-channel bedrock, will remain laterally stable without erosion across reservoir deposits. Patterns of sediment transport measured at the Fall Creek streamgage downstream from Fall Creek Dam provide insight into the potential effects of future streambed drawdowns at Fall Creek Lake. Analyses of suspended sediment measured in WY 2013–18 show a major reduction in suspended-sediment loads between WY 2013 and later years, indicating streamflows transporting sediment through the reservoir to downstream reaches during streambed drawdowns have become supply limited. The 6-year suspended-sediment monitoring and sampling program is insufficient to make predictions about future sediment transport conditions because of uneven monitoring periods and varying controls on reservoir sediment erosion. It is likely that future suspended-sediment loads will be variable but similar to those observed in WY 2015–18 if operational, climatic, and geomorphological factors remain similar to those monitored WY 2015–18. Suspended-sediment loads downstream from Fall Creek Lake will likely remain highest when regulating outlets are fully open and Fall Creek is free flowing with the reservoir fully drained with little to no residual pool. Over time, it is possible that the suspended-sediment loads would reflect mobilization of reservoir sediment deposited in the previous year rather than erosion of sediment deposited years or decades earlier. Bedload is likely to remain a small fraction of the total sediment load evacuated from the reservoir and is relatively modest compared with pre-dam bedload transport rates because most coarse sediment remains trapped by the dam. If sediment releases from Fall Creek Lake and ensuing streamflow conditions follow a similar pattern in the future as was assessed in this study spanning WY 2012–18, near-term geomorphic adjustments downstream of the dam are expected to be modest. Barring major operational, climatic, and geomorphological changes, local site-scale deposition on bars, overbank areas, or off-channel features that persists several months after the streambed drawdown will likely continue to be highly variable, ranging from negligible to several centimeters of deposition. At the landform-scale, low velocity areas nearest to Fall Creek Dam will likely continue to undergo rapid deposition immediately during and after a streambed drawdown event, similar to patterns observed for WY 2012–18. Some of the sediment entering these off-channel features and margin areas may be temporarily stored, then later remobilized and dispersed farther downstream. But if newly deposited sediment persists through the following spring, there is a greater likelihood that local vegetation will establish, reinforce deposited material, and trap sediment during later drawdowns. The reach-scale geomorphic changes may become more apparent if (1) streambed drawdowns continued for several decades, and geomorphic changes were measured at decadal scales or (2) the amount of sediment introduced to downstream reaches substantially increased and (or) sediment transport capacity decreased. The continued streamflow regulation of Fall Creek Dam after sediment releases provides an opportunity to strategically manage streamflows during and after the streambed drawdowns to minimize downstream sediment impacts and ensure other operational thresholds are satisfied. This study provides a comprehensive foundation of datasets and geomorphic analyses to inform dam operations at Fall Creek Lake, monitor sediment transport downstream, and consider operational schemes for future drawdowns. The datasets from this study also provide baselines of sediment transport and geomorphic conditions to assess future changes in reservoir and downstream environments. Future monitoring could be tailored to address specific questions regarding the long-term geomorphic effects of streambed drawdowns on fluvial habitats, flood hazards, cultural resources, or downstream water quality. Future monitoring activities could focus on the relevant geomorphic processes and spatial domains within the three categories used for this study: (1) reservoir erosion and net sediment evacuation, (2) sediment delivery to downstream reaches, including magnitude and temporal pattern of sediment transport, and (3) geomorphic responses of downstream reaches to sediment delivery. Specifically, high priority future monitoring activities could include: Repeat topographic or photographic surveys in the reservoir to characterize changes occurring within individual drawdowns, to quantify sediment export, to determine temporal changes in reservoir storage, and to identify locations of erosion and deposition. Continuous, year-round turbidity monitoring supplemented with suspended-sediment measurements at a streamflow-gaging station immediately downstream from the dam to quantify sediment export. Repeat geomorphic monitoring, mapping, or modeling in downstream reaches to track changes in channel and over bank features using a combination of site- and reach-scale monitoring approaches. This could support assessments of sediment deposition and ensuing vegetation encroachment on flood hazards and habitats and examine how sediment transport and depositional processes may be affected by different sediment supply, streamflow, or dam management scenarios.

Oregon

Coastal-change and glaciological map of the Trinity Peninsula area and south Shetland Islands, Antarctica: 1843-2001: Chapter A in Coastal-change and glaciological maps of Antarctica

Changes in the area and volume of polar ice sheets are intricately linked to changes in global climate, and the resulting changes in sea level could severely impact the densely populated coastal regions on Earth. Melting of the West Antarctic part alone of the Antarctic ice sheet would cause a sea-level rise of approximately 6 meters (m). The potential sea-level rise after melting of the entire Antarctic ice sheet is estimated to be 65 m (Lythe and others, 2001) to 73 m (Williams and Hall, 1993). In addition to its importance, the mass balance (the net volumetric gain or loss) of the Antarctic ice sheet is highly complex, responding differently to different conditions in each region (Vaughan, 2005). In a review paper, Rignot and Thomas (2002) concluded that the West Antarctic ice sheet is probably becoming thinner overall; although it is thickening in the west, it is thinning in the north. Thomas and others (2004), on the basis of aircraft and satellite laser altimetry surveys, believe the thinning may be accelerating. Joughin and Tulaczyk (2002), on the basis of analysis of ice-flow velocities derived from synthetic aperture radar, concluded that most of the Ross ice streams (ice streams on the east side of the Ross Ice Shelf) have a positive mass balance, whereas Rignot and others (2004) infer even larger negative mass balance for glaciers flowing northward into the Amundsen Sea, a trend suggested by Swithinbank and others (2003a,b, 2004). The mass balance of the East Antarctic ice sheet is thought by Davis and others (2005) to be strongly positive on the basis of the change in satellite altimetry measurements made between 1992 and 2003. Measurement of changes in area and mass balance of the Antarctic ice sheet was given a very high priority in recommendations by the Polar Research Board of the National Research Council (1986), in subsequent recommendations by the Scientific Committee on Antarctic Research (SCAR) (1989, 1993), and by the National Science Foundation's (1990) Division of Polar Programs. On the basis of these recommendations, the U.S. Geological Survey (USGS) decided that the archive of early 1970s Landsat 1, 2, and 3 Multispectral Scanner (MSS) images of Antarctica and the subsequent repeat coverage made possible with Landsat and other satellite images provided an excellent means of documenting changes in the coastline of Antarctica (Ferrigno and Gould, 1987). The availability of this information provided the impetus for carrying out a comprehensive analysis of the glaciological features of the coastal regions and changes in ice fronts of Antarctica (Swithinbank, 1988; Williams and Ferrigno, 1988). The project was later modified to include Landsat 4 and 5 MSS and Thematic Mapper (TM) [and in some areas Landsat 7 Enhanced Thematic Mapper Plus (ETM+)], RADARSAT images, and other data where available, to compare changes that occurred during a 20- to 25- or 30-year time interval (or longer where data were available, as in the Antarctic Peninsula). The results of the analysis are being used to produce a digital database and a series of USGS Geologic Investigations Series Maps (I–2600) consisting of 23 maps at 1:1,000,000 scale and 1 map at 1:5,000,000 scale, in both paper and digital format (Williams and others, 1995; Williams and Ferrigno, 1998; Ferrigno and others, 2002).

IMAP