Search USGSSearch

SEARCH · Search USGS

Results for “ICES Scientific Reports”

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 199 records · Page 11Linked to original sources

Streamflow, water quality, and aquatic macroinvertebrates of selected streams in Fairfax County, Virginia, 2007-12

Efforts to mitigate the effects of urbanization on streams rely on best management practices (BMPs) that are implemented with the intent of reducing and retaining stormwater runoff. A cooperative monitoring effort between the U.S. Geological Survey and Fairfax County, Virginia, was initiated in 2007 to assess the condition of county streams and document watershed-scale responses to the implementation of BMPs. Assessment of the data collected during the first 5 years of this monitoring program focused on characterizing the hydrologic and ecological condition of 14 monitored streams. Hydrologic, chemical, and macroinvertebrate community conditions in the streams monitored were found to be consistent, overall, with conditions commonly observed in urban streams. Hydrologically, the monitored streams were found to be flashy, with flashiness positively related to road cover in the watershed. Typical pH values of streams throughout the network centered around neutrality (pH = 7) with strong daily fluctuations apparent in the continuous data. Patterns in specific conductance were largely representative of anthropogenic disturbances—watersheds having the greatest percentage of open space and estate residential land-use had the lowest typical specific conductance values, and specific conductance variability was less than what is observed in watersheds that are more intensively developed. In watersheds having greater road coverage, and more development in general, increases in specific conductance over several orders of magnitude were observed during winter months as a result of the application of de-icing salts on impervious surfaces. Dissolved oxygen conditions were typically within the range required to support healthy biological communities, although occasional departures during summer months at some sites fell below the impairment threshold for streams in Virginia. Nitrogen (N) and phosphorus (P), concentration patterns were largely consistent across the network, with few exceptions. Nitrogen concentrations in monthly samples were generally low and dominated by nitrate. Exceptions to the generally low N concentrations occurred at Captain Hickory Run, which had a median total N concentration of approximately 4.9 milligrams per liter (mg/L), compared to the network-wide median of approximately 1.7 mg/L, and at Popes Head Creek Tributary, where total N concentrations spiked to 6–8 mg/L during low-flow periods in August or September of each year. Phosphorus concentrations in monthly samples were generally low and dominated by the dissolved fraction. Two monitoring stations in the network, Flatlick Branch and Frog Branch, are notable for having median total P concentrations that were, on average, approximately three times greater than the median total P concentration of 0.02 mg/L observed at the other 12 stations in the network. Suspended-sediment and nutrient loads and yields were similar to those of urbanized watersheds in other studies, although the yields from these urbanized basins were greater than, or within, the upper quartile of yields observed throughout the Chesapeake Bay watershed. Annual suspended-sediment loads ranged from 289–10,275 tons, with a median of 1,311 tons, and corresponding yields ranged from 107–2,827 tons per square mile (ton/mi 2 ), with a median of 277 ton/mi 2 . Annual total N loads ranged from 8,014–36,413 pounds, with a median of 21,314 pounds, and corresponding yields ranged from 3,361–8,360 pounds per square mile (lb/mi 2 ), with a median of 6,200 lb/mi 2 . Annual total P loads ranged from 380–6,558 pounds, with a median of 1,874 pounds, and corresponding yields ranged from 140–1,562 lb/mi 2 , with a median of 543 lb/mi 2 . Benthic macroinvertebrate community metrics indicated that streams throughout Fairfax County are generally of poor health. One station, Castle Creek, was an exception with results indicating relatively high quality aquatic health. Six additional monitoring stations were added to the network in 2012 to improve spatial coverage throughout Fairfax County. Monitoring activities are expected to continue at all 20 stations for the foreseeable future as BMP implementation is conducted.

Virginia

A New Perspective on Mount St. Helens - Dramatic Landform Change and Associated Hazards at the Most Active Volcano in the Cascade Range

Mount St. Helens has erupted more frequently than any other volcano in the Cascade Range during the past 4,000 years. The volcano has exhibited a variety of eruption styles?explosive eruptions of pumice and ash, slow but continuous extrusions of viscous lava, and eruptions of fluid lava. Evidence of the volcano?s older eruptions is recorded in the rocks that build and the deposits that flank the mountain. Eruptions at Mount St. Helens over the past three decades serve as reminders of the powerful geologic forces that are reshaping the landscape of the Pacific Northwest. On May 18, 1980, a massive landslide and catastrophic explosive eruption tore away 2.7 cubic kilometers of the mountain and opened a gaping, north-facing crater. Lahars flowed more than 120 kilometers downstream, destroying bridges, roads, and buildings. Ash from the eruption fell as far away as western South Dakota. Reconstruction of the volcano began almost immediately. Between 1980 and 1986, 80 million cubic meters of viscous lava extruded episodically onto the crater floor, sometimes accompanied by minor explosions and small lahars. A lava dome grew to a height of 267 meters, taller than the highest buildings in the nearby city of Portland, Oregon. Crater Glacier formed in the deeply shaded niche between the 1980-86 lava dome and the south crater wall. Its tongues of ice flowed around the east and west sides of the dome. Between 1989 and 1991, multiple explosions of steam and ash rocked the volcano, possibly a result of infiltrating rainfall being heated in the still-hot interior of the dome and underlying crater floor. In September 2004, rising magma caused earthquake swarms and deformation of the crater floor and glacier, which indicated that Mount St. Helens might erupt again soon. On October 1, 2004, a steam and ash explosion signaled the beginning of a new phase of eruptive activity at the volcano. On October 11, hot rock reached the surface and began building a new lava dome immediately south of the 1980-86 lava dome. The erupting lava cleaved Crater Glacier in half and bulldozed it aside, causing thickening, crevassing, and rapid northward advance of the glacier?s east and west arms. Intermittent steam and ash explosions, some generating plumes that rose up to 11 kilometers, preceded and accompanied extrusion of the new lava dome, but ceased by early 2005. As the new dome grew, a series of large fins or spines of hot lava rose, some more than 100 meters high, and then crumbled producing sometimes spectacular rock falls. The largest of these rock falls generated dust or steam plumes that rose high above the crater rim. By February 2006, the new dome had grown to a volume similar to that of the 1980-86 lava dome; and by July 2007, the new dome had grown to a volume of 93 million cubic meters, exceeding the volume of the 1980-86 lava dome. The height of the new dome also exceeded that of the 1980-86 lava dome, and at its highest point (before collapse in 2005) reached to within 2 meters of the lowest point on the south crater rim. At this height, the new dome was taller than the Empire State Building in New York City. The new lava dome initially grew very quickly, at rates of 2 to 3 cubic meters (one small dump truck load) per second. If it had continued to grow at these rates for about 100 years, it would have replaced the volume of rock removed from the volcano during the May 18, 1980, eruption. However, the lava extrusion rate slowed throughout the eruption, and, by July 2007, it was oozing at a rate of 0.1 cubic meters per second. At that rate, it would take over 700 years to replace the volume of rock lost in 1980. Lava dome extrusion has continued into early 2008.

Scientific Investigations Map

Characterization of water quality and biological communities, Fish Creek, Teton County, Wyoming, 2007-2011

Fish Creek, an approximately 25-kilometer-long tributary to Snake River, is located in Teton County in western Wyoming near the town of Wilson. Fish Creek is an important water body because it is used for irrigation, fishing, and recreation and adds scenic value to the Jackson Hole properties it runs through. Public concern about nuisance growths of aquatic plants in Fish Creek has been increasing since the early 2000s. To address these concerns, the U.S. Geological Survey conducted a study in cooperation with the Teton Conservation District to characterize the hydrology, water quality, and biologic communities of Fish Creek during 2007–11. The hydrology of Fish Creek is strongly affected by groundwater contributions from the area known as the Snake River west bank, which lies east of Fish Creek and west of Snake River. Because of this continuous groundwater discharge to the creek, land-use activities in the west bank area can affect the groundwater quality. Evaluation of nitrate isotopes and dissolved-nitrate concentrations in groundwater during the study indicated that nitrate was entering Fish Creek from groundwater, and that the source of nitrate was commonly a septic/sewage effluent or manure source, or multiple sources, potentially including artificial nitrogen fertilizers, natural soil organic matter, and mixtures of sources. Concentrations of dissolved nitrate and orthophosphate, which are key nutrients for growth of aquatic plants, generally were low in Fish Creek and occasionally were less than reporting levels (not detected). One potential reason for the low nutrient concentrations is that nutrients were being consumed by aquatic plant life that increases during the summer growing season, as a result of the seasonal increase in temperature and larger number of daylight hours. Several aspects of Fish Creek’s hydrology contribute to higher productivity and biovolume of aquatic plants in Fish Creek than typically observed in streams of its size in Wyoming. Especially in the winter, the proportionately large, continuous gain of groundwater into Fish Creek in the perennial section keeps most of the creek free of ice. Because sunlight can still reach the streambed in Fish Creek and the water is still flowing, aquatic plants continue to photosynthesize in the winter, albeit at a lower level of productivity. Additionally, the cobble and large gravel substrate in Fish Creek provides excellent attachment points for aquatic plants, and when combined with Fish Creek’s channel stability allows rapid growth of aquatic plants once conditions allow during the spring. The aquatic plant community of Fish Creek was different than most streams in Wyoming in that it contains many different macrophytes—including macroalgae such as long streamers of Cladophora , aquatic vascular plants, and moss; most other streams in the state contain predominantly algae. From the banks of Fish Creek, the bottom of the stream sometimes appeared to be a solid green carpet. A shift was observed from higher amounts of microalgae in April/May to higher amounts macrophytes in August and October, and differences in the relative abundance of microalgae and macrophytes were statistically significant between seasons. Differences in dissolved-nitrate concentrations and in the nitrogen-to-phosphorus ratio were significantly different between seasons, as concentrations of dissolved nitrate decreased from April/May to August and October. It is likely that dissolved-nitrate concentrations in Fish Creek were lower in August and October because macrophytes were quickly utilizing the nutrient, and a negative correlation between macro-phytes and nitrate was found. Macroinvertebrates also were sampled because of their role as indicators of water quality and their documented responses to perturbation such as degradation of water quality and habitat. Statistically significant seasonal differences were noted in the macroinvertebrate community. Taxa richness and relative abundance of Ephemeroptera, Plecoptera, and Trichoptera, which tend to be intolerant of water-quality degradation, decreased from April/May to August; the same time period saw a corresponding increase in Diptera and noninsects, particularly Oligochaeta (worms) that are more tolerant. Seasonal changes in macroinvertebrate functional feeding groups were significantly different. The relative abundance of gatherer-collector and scraper feeding groups decreased from April/May to August, accompanied by an increase in filterer-collector and shredders feeding groups. Seasonal changes in feeding groups might be due to the seasonal shift in aquatic plant communities, as indicated by comparison with other streams in the area that had fewer aquatic macrophytes than Fish Creek. Statistical tests of macroinvertebrate metrics indicated few differences between years or biological sampling sites on Fish Creek, although the site farthest upstream sometimes was different not only in terms of macroinvertebrates but also in streamflow, water quality, and aquatic plants. Potential effects of contributions of additional nutrients to the Fish Creek ecosystem beyond the conditions sampled during the study period are not known. However, because virtually all of the detectable dissolved nitrate commonly was consumed by aquatic plants in August (leaving dissolved nitrate less than the reporting level in water samples), it is possible that increased nutrient contributions could cause increased growth of aquatic plants. Additional long-term monitoring of the stream, with concurrent data analysis and interpretation would be needed to determine the effects of additional nutrients on the aquatic plant community and on higher levels of the food chain.

Wyoming

Halite brine in the Onondaga Trough near Syracuse, New York: Characterization and simulation of variable-density flow

Halite brine (saturation ranging from 45 to 80 percent) lies within glacial-drift deposits that fill the Onondaga Trough, a 40-km long bedrock valley deepened by Pleistocene ice near Syracuse, N.Y. The most concentrated brine occupies the northern end of the trough, more than 15 kilometers (km) beyond the northern limit of halite beds in the Silurian Salina Group, the assumed source of salt. The chemical composition of the brine and its radiocarbon age estimated from geochemical modeling with NETPATH suggest that the brine formed through dissolution of halite by glacial melt water, and later mixed with saline bedrock water about 16,500 years ago. Transient variable-density flow simulations were conducted with SEAWAT to assess current (2005) ground-water flow conditions within the glacial drift. A transient three-dimensional (3D) model using a grid spacing of 100 meters (m) and maximum layer spacing of 30 m was used to simulate a 215-year period from 1790 to 2005. The model was calibrated to observations of water levels, chloride concentrations, and discharges of water and chloride. The model produced an acceptable match to the measured data and provided a reasonable representation of the density distribution within the brine pool. The simulated mass of chloride in storage declined steadily during the 215-year period; however, the decline was mainly due to dispersion, which is probably overestimated because of the large layer spacing. Model results suggest that saline water from waste-disposal operations associated with a chemical plant has migrated beneath the western shore of Onondaga Lake. Two-dimensional (2D) cross-sectional models of the aquifer system within the Onondaga Trough were prepared to test the plausibility of a hypothesis that the brine was derived from a relict source of halite that was dissolved by glacial melt water. The 2D models used parameter estimates obtained with the calibrated 3D model. Model results indicated the brine could have migrated from the bedded-halite subcrop area and remained in the glacial sediments at the northern end of trough for over 16,000 years, as suggested by radiocarbon dating. The 2D models also indicated that slow dissipation of brine occurs through a mixing zone formed by upward flow of freshwater over the southern end of the brine pool. The simulated depletion rate is controlled by the rate of mixing, which is limited by the specified grid resolution and the accuracy of the numerical method used to solve the advection-dispersion equation. A numerical solution obtained by using an implicit finite-difference method with upstream weighting and a 2D grid containing a column and layer spacing of 76 m and 3 m, respectively, provided an acceptable match to chloride concentration profiles measured at three locations within the Onondaga Trough.

New York

Geologic Map of MTM 35337, 40337, and 45337 Quadrangles, Deuteronilus Mensae Region of Mars

Deuteronilus Mensae, first defined as an albedo feature at lat 35.0 deg N., long 5.0 deg E., by U.S. Geological Survey (USGS) and International Astronomical Union (IAU) nomenclature, is a gradational zone along the dichotomy boundary in the northern mid-latitudes of Mars. The boundary in this location includes the transition from the rugged cratered highlands of Arabia Terra to the northern lowland plains of Acidalia Planitia. Within Deuteronilus Mensae, polygonal mesas are prominent along with features diagnostic of Martian fretted terrain, including lobate debris aprons, lineated valley fill, and concentric crater fill. Lobate debris aprons, as well as the valley and crater fill deposits, are geomorphic indicators of ground ice, and their concentration in Deuteronilus Mensae is of great interest because of their potential association with Martian climate change. The paucity of impact craters on the surfaces of debris aprons and the presence of ice-cemented mantle material imply young (for example, Amazonian) surface ages that are consistent with recent climate change in this region of Mars. North of Deuteronilus Mensae are the northern lowlands, a potential depositional sink that may have had large standing bodies of water or an ocean in the past. The northern lowlands have elevations that are several kilometers below the ancient cratered highlands with significantly younger surface ages. The morphologic and topographic characteristics of the Deuteronilus Mensae region record a diverse geologic history, including significant modification of the ancient highland plateau and resurfacing of low-lying regions. Previous studies of this region have interpreted a complex array of geologic processes, including eolian, fluvial and glacial activity, coastal erosion, marine deposition, mass wasting, tectonic faulting, effusive volcanism, and hydrovolcanism. The origin and age of the Martian crustal dichotomy boundary are fundamental questions that remain unresolved at the present time. Several scenarios for its formation, including single and multiple large impact events, have been proposed and debated in the literature. Endogenic processes whereby crust is thinned by internal mantle convection and tectonic processes have also been proposed. Planetary accretion models and isotopic data from Martian meteorites suggest that the crust formed very early in Martian history. Using populations of quasi-circular depressions extracted from the topography of Mars, other studies suggest that the age difference between the highlands and lowlands could be ~100 m.y.. Furthermore, understanding the origin and age of the dichotomy boundary has been made more complicated due to significant erosion and deposition that have modified the boundary and its adjacent regions. The resulting diversity of terrains and features is likely a combined result of ancient and recent events. Detailed geologic analyses of dichotomy boundary zones are important for understanding the spatial and temporal variations in highland evolution. This information, and comparisons to other highland regions, can help elucidate the scale of potential environmental changes. Previous geomorphic and geologic mapping investigations of the Deuteronilus Mensae region have been completed at local to global scales. The regional geology was first mapped by Lucchitta (1978) at 1:5,000,000 scale using Mariner 9 data. This study concluded that high crater flux early in Martian history formed overlapping craters and basins that were later filled by voluminous lava flows that buried the impacted surface, creating the highlands. After this period of heavy bombardment, fluvial erosion of the highlands formed the canyons and valleys, followed by dissection that created the small mesas and buttes, and later, formation of the steep escarpment marking the present-day northern highland margin. After valley dissection, mass wasting and eolian processes caused lateral retreat of mesas and buttes

Scientific Investigations Map

Chloride concentrations, loads, and yields in four watersheds along Interstate 95, southeastern Connecticut, 2008-11: factors that affect peak chloride concentrations during winter storms

Chloride (Cl - ) concentrations and loads and other water chemistry characteristics were assessed to evaluate potential effects of road-deicer applications on streamwater quality in four watersheds along Interstate 95 (I–95) in southeastern Connecticut from November 1, 2008, through September 30, 2011. Streamflow and water quality were studied in the Four Mile River, Oil Mill Brook, Stony Brook, and Jordan Brook watersheds, where developed land ranged from 9 to 32 percent. Water-quality samples were collected and specific conductance was measured continuously at paired water-quality monitoring sites, upstream and downstream from I–95. Specific conductance values were related to Cl - concentrations to assist in determining the effects of road-deicing operations on the levels of Cl - in the streams. Streamflow and water-quality data were compared with weather data and with the timing, amount, and composition of deicers applied to State highways. Grab samples were collected during winter stormwater-runoff events, such as winter storms or periods of rain or warm temperatures in which melting takes place. Grab samples were also collected periodically during the spring and summer and during base-flow conditions. The estimated Cl - concentrations at the eight water-quality monitoring sites during winter storms peaked as high as 270 milligrams per liter (mg/L) and were well below the U.S. Environmental Protection Agency (EPA) recommended acute chloride toxicity criterion of 860 mg/L and the chronic 4-day average toxicity criterion of 230 mg/L. Cl - concentrations in streams, particularly at sites downstream from I–95, peaked during increased streamflow in the winter and early spring as a result of deicers applied to roads and washed off by stormwater or meltwater. Cl - concentrations during most of the nonwinter seasons decreased during increases in streamflow because storm runoff was more dilute than base flow. However, peaks in specific conductance and estimated chloride concentrations at streams in more urbanized watersheds corresponded to peaks in streamflow well after winter snow or ice events; these delayed peaks in Cl - concentration likely resulted from deicer residue that remained in melting snow piles and on roadsides and (or) that were flushed from soils and shallow groundwater, then discharged downstream. Estimated peak Cl - concentrations varied with the type of winter storm event and were highest during or after winter storms of frozen precipitation and rain, in which the rain or meltwater effectively washed off the deicers. Estimated peak Cl - concentrations correlated positively with the duration of deicer application but generally not with streamflow. Estimated peak Cl - concentrations during the winter season were highest during low streamflow at most sites. Chloride concentrations varied considerably in shallow groundwater as a result of land-use differences. Cl - concentrations were very high (as high as 800 mg/L) in shallow groundwater downstream from I–95 at the Four Mile River site. Chloride/bromide mass concentration ratios and the proximity of a former landfill and sewage lagoon upstream indicate a likely source of Cl - is landfill leachate and possibly sewage leachate. Cl - loads in streams generally were highest in the winter and early spring. The estimated daily Cl - yield for the four monitoring sites downstream from I–95 ranged from 0.0004 ton per day per square mile for one of the least developed watersheds to 0.052 ton per day per square mile for the watershed with the highest percentage of urban development and impervious area. The estimated median contribution of Cl - load from atmospheric deposition was small and ranged from 0.07 percent of Cl - load at the Jordan Brook watershed to 0.57 percent at the Oil Mill Brook watershed. The Cl - loads in streams (outputs) were compared with Cl - load inputs, which include atmospheric deposition and deicer applications; Cl - load inputs were slightly larger than the Cl - load outputs at most of the sites during most years but do not account for the Cl - load in groundwater leaving the watersheds. A multiple linear regression model was developed to describe the variability of the natural log of peak specific conductance, as well as estimated Cl - concentrations. Five significant variables best explained the variability in the natural log of the peak specific conductance after deicing events: (1) snow on ground before deicing event; (2) winter precipitation with rain; (3) specific conductance in base flow; (4) State-operated road lane miles divided by watershed area; and (5) amount of Cl - from deicers applied to State-operated roads per lane mile. In this report, winter precipitation is defined as any type of precipitation, including frozen precipitation and rain, that occurs during the active deicing season, typically November through March. Frozen precipitation is defined here as snow, sleet, freezing rain, or any winter precipitation except rain. The addition of a lane mile in both directions on I–95 would result in an estimate of approximately 2 to 11 percent increase in Cl - input from deicers applied to I–95 and other roads maintained by Connecticut Department of Transportation. The largest estimated increase in Cl - load was in the watersheds with the greatest number miles of I–95 corridor relative to the total lane miles maintained by Connecticut Department of Transportation. On the basis of these estimates and the estimated peak Cl - concentrations during the study period, it is unlikely that the increased use of deicers on the additional lanes would lead to Cl - concentrations that exceed the aquatic habitat criteria.

Connecticut

Modeled and measured glacier change and related glaciological, hydrological, and meteorological conditions at South Cascade Glacier, Washington, balance and water years 2006 and 2007

Winter snow accumulation and summer snow and ice ablation were measured at South Cascade Glacier, Washington, to estimate glacier mass balance quantities for balance years 2006 and 2007. Mass balances were computed with assistance from a new model that was based on the works of other glacier researchers. The model, which was developed for mass balance practitioners, coupled selected meteorological and glaciological data to systematically estimate daily mass balance at selected glacier sites. The North Cascade Range in the vicinity of South Cascade Glacier accumulated approximately average to above average winter snow packs during 2006 and 2007. Correspondingly, the balance years 2006 and 2007 maximum winter snow mass balances of South Cascade Glacier, 2.61 and 3.41 meters water equivalent, respectively, were approximately equal to or more positive (larger) than the average of such balances since 1959. The 2006 glacier summer balance, -4.20 meters water equivalent, was among the four most negative since 1959. The 2007 glacier summer balance, -3.63 meters water equivalent, was among the 14 most negative since 1959. The glacier continued to lose mass during 2006 and 2007, as it commonly has since 1953, but the loss was much smaller during 2007 than during 2006. The 2006 glacier net balance, -1.59 meters water equivalent, was 1.02 meters water equivalent more negative (smaller) than the average during 1953-2005. The 2007 glacier net balance, -0.22 meters water equivalent, was 0.37 meters water equivalent less negative (larger) than the average during 1953-2006. The 2006 accumulation area ratio was less than 0.10, owing to isolated patches of accumulated snow that endured the 2006 summer season. The 2006 equilibrium line altitude was higher than the glacier. The 2007 accumulation area ratio and equilibrium line altitude were 0.60 and 1,880 meters, respectively. Accompanying the glacier mass losses were retreat of the terminus and reduction of total glacier area. The terminus retreated at a rate of about 13 meters per year during balance year 2006 and at a rate of about 8 meters per year during balance year 2007. Glacier area near the end of balance years 2006 and 2007 was 1.74 and 1.73 square kilometers, respectively. Runoff from the basin containing the glacier and from an adjacent nonglacierized basin was gaged during all or parts of water years 2006 and 2007. Air temperature, wind speed, precipitation, and incoming solar radiation were measured at selected locations on and near the glacier. Air-temperature over the glacier at a height of 2 meters generally was less than at the same altitude in the air mass away from the glacier. Cooling of the air by the glacier increased systematically with increasing ambient air temperature. Empirically based equations were developed to estimate 2-meter-height air temperature over the glacier at five sites from site altitude and temperature at a non-glacier reference site.

Scientific Investigations Report

Hydrogeologic Framework and Ground-Water Budget of the Spokane Valley-Rathdrum Prairie Aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources and Washington State Department of Ecology, investigated the hydrogeologic framework and ground-water budget of the Spokane Valley-Rathdrum Prairie (SVRP) aquifer located in northern Idaho and northeastern Washington. Descriptions of the hydrogeologic framework, water-budget components, and further data needs are provided. The SVRP aquifer, which covers about 370 square miles including the Rathdrum Prairie, Idaho, and the Spokane Valley and Hillyard Trough, Washington, is the sole source of drinking water for more than 500,000 residents. Continued growth, water-management issues, and potential effects on water availability and water quality in the aquifer and in the Spokane and Little Spokane Rivers have illustrated the need to better understand and manage the region's water resources. The SVRP aquifer consists mostly of gravels, cobbles, and boulders - deposited during a series of outburst floods resulting from repeated collapse of the ice dam that impounded ancient Glacial Lake Missoula. In most places, the SVRP aquifer is bounded by bedrock of pre-Tertiary granite or metasedimentary rocks, or Miocene basalt and associated sedimentary deposits. Discontinuous fine-grained layers are scattered throughout the SVRP aquifer at considerably different altitudes and with considerably different thicknesses. In the Hillyard Trough and the Little Spokane River Arm of the aquifer, a massive fine-grained layer with a top altitude ranging from about 1,500 to 1,700 feet and thickness ranging from about 100 to 200 feet separates the aquifer into upper and lower units. Most of the Spokane Valley part of the aquifer is devoid of fine-grained layers except near the margins of the valley and near the mouths of lakes. In the Rathdrum Prairie, multiple fine-grained layers are scattered throughout the aquifer with top altitudes ranging from about 1,700 to 2,400 feet with thicknesses ranging from 1 to more than 135 feet. The altitude of the base of the aquifer ranges from less than 1,800 feet near Lake Pend Oreille to less than 1,200 feet near the aquifer's outlet near Long Lake. The thickness of the aquifer is more than 800 feet in the northwestern part of the northern Rathdrum Prairie, through the West Channel area, and through the west-central part of the Rathdrum Prairie. In Washington, the areas of greatest thickness, more than 600 feet, are mapped in the central parts of the Spokane Valley, Spokane, and the Hillyard Trough. Recharge or inflow to the SVRP aquifer occurs from six main sources: the Spokane River, lakes, infiltration from precipitation over the aquifer, tributaries, infiltration from landscape irrigation and septic systems, and subsurface inflow. Discharge or outflow from the SVRP aquifer occurs from five main sources: the Spokane River, the Little Spokane River, pumpage, subsurface discharge to Long Lake, and infiltration of ground water to sewers. Total estimated mean annual inflow to and outflow from the SVRP aquifer is about 1,470 cubic feet per second. Several data needs were identified during this investigation that would improve the definition of the hydrogeologic framework and ground-water budget components for the SVRP aquifer study area. Deep drilling along the axis of the aquifer could determine the depth to the bottom of the aquifer where data are currently unavailable as well as identify the presence of fine-grained layers and their thickness. A more detailed analysis of the geologic and hydrologic setting near the southern ends of Spirit and Hoodoo Valleys could help determine the location of the ground-water divide between the two valleys and the Rathdrum Prairie. Better estimates of seepage into the aquifer from Coeur d'Alene Lake and Lake Pend Oreille and underflow from the aquifer to Long Lake would strengthen the recharge and discharge estimates of the aquifer. A hydrochemical study incorporating analyses of envi

Scientific Investigations Report

Geohydrology and water quality of the unconsolidated aquifers in the Enfield Creek Valley, town of Enfield, Tompkins County, New York

From 2013 to 2018, the U.S. Geological Survey, in cooperation with the Town of Enfield and the Tompkins County Planning Department, studied the unconsolidated aquifer in the Enfield Creek Valley in the town of Enfield, Tompkins County, New York. The valley will likely undergo future development as the population of Tompkins County increases and spreads out from the metropolitan areas. The Town of Enfield, Tompkins County, and the New York State Departments of Health and Environmental Conservation need geohydrologic information to help planners develop a more comprehensive approach to water-resource management in Tompkins County. The Enfield Creek Valley is underlain by an unconfined aquifer that consists of saturated alluvium, alluvial-fan deposits, and ice-contact (kame) sand and gravel. A confined aquifer of discontinuous ice-contact sand and gravel overlies bedrock. Depth to bedrock in the valley ranges from about 50 feet below land surface from just north of the Enfield Creek divide in the northern part of the aquifer to the confluence of Fivemile Creek to at least 140 feet below land surface from Fivemile Creek to where the valley orientation changes from north-south to northwest-southeast. Depth to bedrock is much shallower from the valley orientation change to the southeastern part of the aquifer because Enfield Creek has carved through overlying sediments into bedrock as the creek drops 450 feet into the Cayuga Inlet Valley. A small buried valley running south to north was identified within the Fivemile Creek drainage along the western edge of the town. However, the valley fill consists of glacial till, and no sand-and-gravel aquifer is present. The unconfined aquifers are recharged by direct infiltration of precipitation, surface runoff, and shallow subsurface flow from hillsides, and by seepage loss from streams overlying the aquifer. The confined aquifers are recharged mostly by precipitation that enters the adjacent valley walls, by groundwater flowing from bordering till or bedrock, and by flow from the bottom of the valley. Also, some recharge may be occurring where confining units are absent or from confining units with sediments of moderate permeability. Groundwater discharges to Enfield Creek, its tributaries, and wetlands and is lost through evapotranspiration from the water table or is withdrawn from domestic, commercial, and agricultural wells. About 700 individual well owners depend on the unconsolidated aquifers for their water supply. An estimated 28,300,000 gallons per year are withdrawn. Groundwater samples were collected from eight test wells drilled for this study, and six surface-water samples were collected from five locations on Enfield Creek. Of the eight wells sampled, two were finished in unconfined sand-and-gravel aquifers, two were finished in confined sand-and-gravel aquifers, and four were finished at or near the shale bedrock surface. Water quality in the study area generally met State and Federal drinking-water standards. However, some samples exceeded maximum contaminant levels for barium (25 percent of samples) and secondary maximum containment levels for chloride (25 percent), dissolved solids (25 percent of samples), iron (70 percent of samples), and manganese (75 percent of samples). Groundwater from 75 percent of the wells sampled for methane had concentrations greater than the Office of Surface Mining Reclamation and Enforcement recommended action level of 10 to 28 milligrams per liter. The two deepest wells sampled, TM1075 and TM1077, had the highest specific conductance, chloride, and sodium concentrations of all wells sampled. The chloride/bromide ratios of these samples suggest the source may represent a mixture of saline formation waters with shallow dilute groundwater and may receive recharge contribution from two tributaries overlying bedrock to the west and southwest of the aquifer. In general, the highest yields are from wells completed within about 50 feet below land surface, which may tap either type of aquifer.

New York

Modeling hydrodynamics, water temperature, and water quality in Klamath Straits Drain, Oregon and California, 2012–15

Executive Summary Located southwest of Klamath Falls, Oregon, Klamath Straits Drain is a 10.1-mile-long canal that conveys water uphill and northward through the use of pumps before discharging to the Klamath River. Klamath Straits Drain traverses an area that historically encompassed Lower Klamath Lake. Currently, the Drain receives water from farmland and from parts of the Lower Klamath Lake National Wildlife Refuge. To support water-quality improvement in Klamath Straits Drain, a hydrodynamic and water-temperature model was constructed and calibrated for calendar years 2012–15 with the two-dimensional model CE-QUAL-W2 (version 4.0). Water quality was calibrated for a subset of that time, from April 1, 2012 to March 31, 2015. Flows in calendar year 2012 were within the normal range, while calendar years 2013–15 were dry years. Significant findings from this study include: In the years studied, only limited flow entered Klamath Straits Drain at the upstream Headworks (KSDH) site. Most flow entered the Drain between KSDH and the E-EE pumps near Township Road through several irrigation channels and ditches. Few data were available to describe the quality of this water for the period of study. The E-EE and F-FF pumps along Klamath Straits Drain mainly operated automatically to keep water levels relatively steady. Ten-minute flow data at streamgage 11509340, downstream of the F-FF pumps, showed high-frequency on/off switching of the F-FF pumps. Combined with daily mean flow data from the F-FF pumps, the downstream 10-minute flow data allowed estimation of 10-minute pumping rates for the F-FF pumps. Paper pump charts showed the existence of short-term variability at the E-EE pumps; however, daily pump data were used at the E-EE pump location in the model. Water temperature in Klamath Straits Drain varied from less than 5 degrees Celsius (°C) (with occasional ice cover in December–January) to greater than 20 °C in May–September. In the years studied, specific conductance was typically 250–850 microsiemens per centimeter, higher than Klamath River specific conductance (typically 100–200 microsiemens per centimeter). Increased chlorophyll a in autumn and winter, along with supersaturated oxygen concentrations, indicated algal blooms in the Drain at that time of year. The blooms were most likely diatoms, based on the timing of blooms sampled elsewhere. Total nitrogen concentration was as much as 5.5 mg/L, with most in dissolved organic and particulate forms, and lower amounts in ammonia and nitrate+nitrite. Total phosphorus concentrations were distributed between orthophophorus (at a median concentration of 0.15 mg/L) and organic and particulate forms (at a median concentration of 0.13 mg/L). Most of the organic carbon in the Klamath Straits Drain was in dissolved rather than particulate form. Newly collected water-quality data for April 1, 2012–March 31, 2015 helped provide the impetus for this modeling study. However, a lack of some data still hindered the construction and calibration of this model. The model would benefit from additional data to describe water-quality boundary conditions, water-quality calibration data upstream of the F-FF pumps, short-term E-EE pump operations, and channel bathymetry in the reach between Highway 97 and the confluence with the Klamath River. Klamath River water mixed upstream into the Klamath Straits Drain, up to the Klamath Straits Drain F-FF pumps at Highway 97, when the F-FF pumps were not operating for periods of hours to days. The F-FF pumps were off for many days during this study, especially during dry years. The boundary between Klamath Straits Drain and the Klamath River was best modeled with an external head condition, which allows exchange of water between the river and the drain in both directions, upstream and downstream. Currently there is a flow gage, water-quality monitor, and a water-quality sampling site located downstream of the F-FF pumps, in the reach where Klamath Straits Drain water can mix with Klamath River water. To sample solely Klamath Straits Drain water, water samples would need to be collected only when the F-FF pumps are actively pumping. Alternately, the sampling location could be moved upstream of the pumps. Interpretation and use of historical water-quality data at the Klamath Straits Drain at Highway 97 site should be done in conjunction with information on pump activity to help inform whether mixing with Klamath River water may have occurred. Total 2014 (a dry year) phosphorus loads from the Drain to the Klamath River were lower and closer to total maximum daily load (TMDL) allocations, as compared to 2013, a year with greater flow and pumping. Modeled travel time through the Klamath Straits Drain, from Headworks to its confluence with the Klamath River, ranged from approximately 24 hours at high flow to 16 days or more, depending on how many days the pumps were turned off. The longer travel times are sufficient for important water-quality transformations, such as algal growth and organic-matter decomposition. This newly constructed model of the Klamath Straits Drain simulates flow, water levels, water temperature, and water quality with acceptable accuracy but with certain data limitations. This model should prove useful in evaluating potential strategies for flow and water-quality management and restoration.

California, Oregon

Modeling hydrodynamics, water temperature, and water quality in the Klamath River upstream of Keno Dam, Oregon, 2006-09

A hydrodynamic, water temperature, and water-quality model was constructed for a 20-mile reach of the Klamath River downstream of Upper Klamath Lake, from Link River to Keno Dam, for calendar years 2006-09. The two-dimensional, laterally averaged model CE-QUAL-W2 was used to simulate water velocity, ice cover, water temperature, specific conductance, dissolved and suspended solids, dissolved oxygen, total nitrogen, ammonia, nitrate, total phosphorus, orthophosphate, dissolved and particulate organic matter, and three algal groups. The Link-Keno model successfully simulated the most important spatial and temporal patterns in the measured data for this 4-year time period. The model calibration process provided critical insights into water-quality processes and the nature of those inputs and processes that drive water quality in this reach. The model was used not only to reproduce and better understand water-quality conditions that occurred in 2006-09, but also to test several load-reduction scenarios that have implications for future water-resources management in the river basin. The model construction and calibration process provided results concerning water quality and transport in the Link-Keno reach of the Klamath River, ranging from interesting circulation patterns in the Lake Ewauna area to the nature and importance of organic matter and algae. These insights and results include: * Modeled segment-average water velocities ranged from near 0.0 to 3.0 ft/s in 2006 through 2009. Travel time through the model reach was about 4 days at 2,000 ft 3 /s and 12 days at 700 ft3/s flow. Flow direction was aligned with the upstream-downstream channel axis for most of the Link-Keno reach, except for Lake Ewauna. Wind effects were pronounced at Lake Ewauna during low-flow conditions, often with circulation in the form of a gyre that rotated in a clockwise direction when winds were towards the southeast and in a counterclockwise direction when winds were towards the northwest. * Water temperatures ranged from near freezing in winter to near 30 degrees C at some locations and periods in summer; seasonal water temperature patterns were similar at the inflow and outflow. Although vertical temperature stratification was not present at most times and locations, weak stratification could persist for periods up to 1-2 weeks, especially in the downstream parts of the reach. Thermal stratification was important in controlling vertical variations in water quality. * The specific conductance, and thus density, of tributaries within the reach usually was higher than that of the river itself, so that inflows tended to sink below the river surface. This was especially notable for inflows from the Klamath Straits Drain, which tended to sink to the bottom of the Klamath River at its confluence and not mix vertically for several miles downstream. * The model was able to capture most of the seasonal changes in the algal population by modeling that population with three algal groups: blue-green algae, diatoms, and other algae. The blooms of blue-green algae, consisting mostly of Aphanizomenon flos aquae that entered from Upper Klamath Lake, were dominant, dwarfing the populations of the other two algae groups in summer. A large part of the blue-green algae population that entered this reach from upstream tended to settle out, die, and decompose, especially in the upper part of the Link-Keno reach. Diatoms reached a maximum in spring and other algae in midsummer. * Organic matter, occurring in both dissolved and particulate forms, was critical to the water quality of this reach of the Klamath River, and was strongly tied to nutrient and dissolved-oxygen dynamics. Dissolved and particulate organic matter were subdivided into labile (quickly decaying) and refractory (slowing decaying) groups for modeling purposes. The particulate matter in summer, consisting largely of dead blue-green algae, decayed quickly. Consequently, this particulate matt

Oregon

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

Science mission requirements for a globally ranging, riserless drilling vessel for U.S. Scientific Ocean Drilling

Through the collection and analysis of shallow and deep subseafloor sediments, rocks, fluids, and life, scientific ocean drilling has enriched our understanding of the complex Earth system. Among other achievements, scientific ocean drilling has documented the history of Earth’s climate, the waxing and waning of polar ice sheets, the past changes in ocean and atmospheric circulation, the existence and function of microbial life in the subseafloor, the compositional variations in Earth’s crust and underlying mantle, and the physical and chemical processes acting at subduction zones, including those associated with tsunamigenic earthquakes. Over the decades, more than 12,000 articles that depend on analyses of scientific ocean drilling samples and geophysical data have been published, many detailing breakthrough contributions to global knowledge about the Earth system. Approximately 45% of these publications were led by U.S.-affiliated authors (International Ocean Discovery Program Publication Services, 2021). Since the mid-1980s, the workhorse of this multidisciplinary, international research effort has been the riserless D/V JOIDES Resolution, operated by Texas A&M University with funding from the U.S. National Science Foundation (NSF). D/V JOIDES Resolution has conducted the vast majority of scientific ocean drilling expeditions and collected most of the scientific cores over that period, including 82% of the expeditions and 93% of the cores in the last decade alone, despite being one of three platforms that is operated within the International Ocean Discovery Program. However, D/V JOIDES Resolution is approaching the end of its useful life. With a strong commitment to continue scientific ocean drilling beyond the end of the current phase, the community developed a document outlining the research frontiers that should be pursued. Exploring Earth by Scientific Ocean Drilling: 2050 Science Framework (Koppers and Coggon, 2020) describes seven scientific strategic objectives that focus on understanding interconnections within the Earth system and five flagship initiatives that integrate these objectives into long-term research efforts that address issues facing society. Additional elements in the 2050 Science Framework, including STEM education, workforce development, technology development, and innovative applications of data analytics, will advance the goals of scientific ocean drilling. Addressing the 2050 Science Framework also requires building partnerships with allied U.S. and international science programs and strengthening existing ones. To implement a significant portion of the 2050 Science Framework, the U.S. scientific community seeks to lease or acquire a newly built, globally ranging, state-of-the art, riserless drilling vessel. The many and varied technical and human resources requirements for successful accomplishment of scientific and educational goals summarized in this document and described in detail in the 2050 Science Framework require broad community input and careful consideration. Following receipt of NSF’s formal Request for Assistance to the United States Science Support Program (USSSP), the U.S. scientific ocean drilling community conducted a one-year exercise to identify its national scientific needs and priorities in order to determine the Science Mission Requirements (SMRs) presented here. This community effort included: (1) a U.S. community-wide survey to identify the specific operational and technical capabilities critical to addressing science in the 2050 Science Framework; (2) a series of online workshops focusing on critical capabilities identified by the survey; and (3) a large in-person workshop to synthesize the results of the survey and the virtual workshops (Appendix 1). The approach was designed to reach as many participants as possible. Overall, 278 survey responses were received from U.S. community members, representing 104 unique institutions from 39 states and the District of Columbia, and 137 unique individuals participated in the workshops (Appendix 2). The results of this effort comprise two classes of SMRs: Foundational Science Mission Requirements and Primary Science Mission Requirements. Foundational SMRs define minimum criteria for a new riserless drilling vessel that can address significant portions of the 2050 Science Framework. Primary SMRs build upon the Foundational SMRs and will create more robust science opportunities and data collection capabilities, will increase progress in addressing the 2050 Science Framework objectives, and will provide more real-time ship-to-shore interaction to improve science productivity, engagement, and outreach. Modern safety and environmental standards, including meeting standards to access protected waters such as exclusive economic zones, extended continental shelves, or high latitudes, while being cognizant of the vessel’s environmental footprint. Safe and efficient operations in global locations and in water depths from 70 m to 6000+ m, with total drill string length of at least 7000 m. High-quality core and data collection from a range of key subseafloor environments. Advanced heave compensation, dynamic positioning, and drill pipe stability. Modern mud and cement/casing systems. Critical onboard measurements for safety, operational decision-making, documentation of ephemeral properties, mission-specific science, and long-term science goals that extend beyond a single expedition. Designated and appropriate space for sample and data preservation. Highly skilled onboard personnel, including technical staff for curation and core handling; support for safety, time-sensitive, and critical shipboard measurements; computer support; equipment and instrument repair; application support; and data assurance. Primary Science Mission Requirements include: Flexible shipboard space for laboratories and on deck to ensure safe, successful implementation of diverse science objectives and operations. Minimizing contamination of recovered samples. Over-the-side capabilities for science-supporting technology (e.g., remotely operated vehicles, water- column sampling, sediment-water interface sampling). Downhole logging and measurements. Expanded borehole observatory capabilities. Reliable and consistent ship-to-shore communications. NSF’s investment in a new globally ranging, riserless drilling vessel will have a powerful economic multiplier effect, including the infusion of additional science support funds in the United States for training and research, the development of new technologies and tools, and the associated scientific and technical workforce development. The skills and knowledge gained through scientific ocean drilling are translatable to careers in fields such as sustainable energy development (e.g., geothermal and offshore wind), carbon sequestration, data management and cyberinfrastructure, biotechnology, communications, science education, policy, hazard mitigation, and environmental management. The United States is a leader in a well-established and internationally collaborative scientific ocean drilling community. A modern, globally ranging, riserless drilling vessel will allow the United States to expand its leadership position, address broad scientific questions that current capabilities preclude, and cultivate equitable international, multidisciplinary collaborations that will ensure scientific ocean drilling’s future success.

Report

Monitoring lahars

Introduction Lahars, or debris flows that originate from a volcano (Pierson and Scott, 1985; Pierson, 1995), are among the most destructive, far-reaching, and persistent hazards on stratovolcanoes. Lahars may be triggered by syneruptive rapid melting of snow and ice, lake breakouts, or heavy rains in conjunction with large eruptive columns. Alternatively, lahars can follow eruptions, when clastic deposits are mobilized by heavy rainfall or lake breakouts, occurring sporadically for years to decades after large eruptions. Some lahars can travel many tens of kilometers in river drainages stemming from volcanoes, as during the 1980 eruption of Mount St. Helens (Washington) (for example, Janda and others, 1981), recent eruptions of Redoubt Volcano (Alaska) (fig. H1; Dorava and Meyer, 1994; Waythomas and others, 2013), and the 1991 eruption of Mount Pinatubo (Philippines) (Major and others, 1996; Pierson and others, 1996). Large lahars are less likely in the absence of eruptive activity, but still possible. The Electron Mudflow at Mount Rainier (approximately A.D. 1500), Wash., is an example of a potential noneruptive lahar, likely initiated by a spontaneous collapse of weak rock, that reached the Puget Lowland after it flowed dozens of kilometers without a recognized eruptive trigger (Sisson and Vallance, 2009). The extreme hazard posed by lahars was demonstrated tragically by the 1985 Nevado del Ruiz (Colombia) catastrophe that claimed the lives of more than 20,000 people (Naranjo and others, 1986). The potential to provide warnings of minutes to hours in advance of lahar arrival in a populated area (for example, Voight, 1990) is a strong reason to provide special monitoring attention to the hazard. Populated river valleys are located downstream from many very high threat and high threat volcanoes, and these areas could be affected by lahars (for example, Hoblitt and others, 1998). The volume and mobility of lahars are two characteristics that can influence the extent of downstream effects (for example, George and others, 2022). The flows that reach the farthest downstream are mobile and voluminous. Additionally, entrainment of material as a lahar travels downstream may increase the volume, and a lahar that starts small may grow to a destructive size under certain conditions. Increasingly, stratovolcanoes host recreational enthusiasts who could be affected by relatively localized geologic hazards, such as rainfall-induced debris flows, glacial outburst floods, rockfalls, and avalanches. These types of events can be common on many volcanoes, occurring seasonally in the case of debris flows and several times per year in the case of avalanches and rockfalls (for example, Allstadt and others, 2018). Many very high threat stratovolcanoes, especially within the contiguous United States, have low eruption frequencies (less than once per century), such that monitoring networks could be used more often for detection and characterization of small surface flows than for identification of volcanic unrest. Such information can be used to validate avalanche forecasts, inform rescue efforts, or notify other agencies of potentially damaged infrastructure (for example, roads, powerlines, or trails). Note that although many of these smaller surface flows create seismic and infrasound waves, the signals are typically highly distorted by the complex volcanic topography and geology. In general, the smaller the flow, the weaker the geophysical signals that it generates, and thus a denser geophysical network is required to study smaller flows (for example, Allstadt and others, 2018). Lahar detection may not be an appropriate or necessary monitoring capability for all volcanoes. Some very high threat volcanoes, like Kīlauea and Mauna Loa, have no lahar hazards currently, and thus no detection, tracking, and characterization capabilities for lahars are needed. At other very high threat volcanoes, such as Pavlof Volcano, Alaska, lahars might be common but pose minimal threat because the volcano is so remote. Ideally, the local observatory would understand the combination of hazard and risk associated with surface flows and assign monitoring and detection capabilities appropriately. Several volcano monitoring techniques (for example, Real-Time Seismic Amplitude Measurement [RSAM], amplitude-based locations, and infrasound array processing) can be adapted to also detect, characterize, and track debris flows, lahars, and other surface flows, so instrumentation installed for detecting volcanic unrest and eruptions can have multiple purposes. The utility of instrumentation for the purpose of monitoring unrest and lahars further justifies the importance and utility of a dense network of monitoring stations, even if the volcano remains quiescent.

Scientific Investigations Report

Hydrology and cycling of nitrogen and phosphorus in Little Bean Marsh: A remnant riparian wetland along the Missouri River in Platte County, Missouri, 1996–97

The lack of concurrent water-quality and hydrologic data on riparian wetlands in the Midwestern United States has resulted in a lack of knowledge about the water-quality functions that these wetlands provide. Therefore, Little Bean Marsh, a remnant riparian wetland along the Missouri River, was investigated in 1996 and 1997 primarily to determine the magnitude and character of selected water-quality benefits that can be produced in such a wetland and to identify critical processes that can be managed in remnant or restored riparian wetlands for amelioration of water quality. Little Bean Marsh averages 69 hectares in size, has a maximum depth of about 1 meter, and the majority of the marsh is covered by macrophytes. In 1997, 41 percent of the water received by Little Bean Marsh was from direct precipitation, 14 percent was from ground-water seepage, 30 percent from watershed runoff, and 15 percent was backflow from Bean Lake. Although, Little Bean Marsh was both a ground-water recharge and discharge area, discharge to the marsh was three times the recharge to ground water. Ground-water levels closely tracked marsh water levels indicating a strong hydraulic connection between ground water and the marsh. Reduced surface runoff and ground-water availability are stabilizing influences on marsh hydrology and probably contribute to the persistence of emergent vegetation. The rapid hydraulic connection between Little Bean Marsh and ground water indicates that the hydrologic regime of most wetlands along the lower Missouri River is largely a function of the altitude of the marsh bottom relative to the altitude of the water table. More water was lost from the marsh through evapotranspiration (59 percent) than all other pathways combined. This is partially because the transpiration process of abundant macrophytes can greatly contribute to the evapotranspiration above that lost from open water surfaces. Surface outflow accounted for 36 percent and ground-water seepage accounted for only 5 percent of the losses. Large residence times allows the marsh to greatly affect water quality before water escapes as ground-water recharge or surface outflow. The shallowness of Little Bean Marsh and ion exclusion during ice formation caused the highest specific conductances of 1,100 to 1,300 microsiemens per centimeter at 25 degrees Celsius to occur during the winter. This concentration of dissolved solutes under ice can make wetlands more vulnerable to toxic contaminants than deeper surface-water bodies. Dissolved oxygen was less than 5 mg/L (milligrams per liter) for 3 to 4 months and near 0 mg/L for about 1 month in summer. Despite depths of less than 1 meter, temperature stratification persisted more than 3 months during the summers of 1996 and 1997, preventing mixing and contributing to periods of anoxia. Shallow depths and extended periods of anoxia in the marsh limit the ability of some organisms to escape high-temperature stress. Turbidity in Little Bean Marsh usually was low for several reasons: sediment loadings from the largely flood-plain drainage were low, emergent vegetation shade out algae and shield the water from wind, and high concentrations of bivalent cations increase flocculation rates of inorganic suspended material. The high concentrations of bivalent cations was largely because of a substantial amount of ground-water seepage into the marsh. Dissolved organic nitrogen was the dominant nitrogen species in Little Bean Marsh. Denitrification and biotic uptake kept more than 62 percent of nitrate (NO3) and 43 percent of ammonium (NH4) concentrations in marsh samples less than a detection limit of 0.005 mg/L. This contrasts with the Missouri River where inorganic NO3 dominates. Consequently, artificial flood-plain drainage that bypasses riparian wetlands likely deliver substantially more biotically available inorganic nitrogen to receiving waters than surface water that has been routed through a remnant wetland. Average total nitrogen concentrations in Little Bean Marsh were substantially less than those at other Missouri River wetlands, roughly one-half the mean concentrations in the Missouri River, but roughly twice the average nitrogen values in reservoirs of the glaciated plains of Missouri. The largest concentrations of nearly all species of nitrogen and phosphorus and the most intense period of hypereutrophy coincided with a phytoplankton bloom and senescence of River Bulrush (Scirpus fluviatilis) and common cattail (Typha latifolia) in September 1997. The rapid leaching of nitrogen that occurs soon after macrophyte senescence combined with a recent destratification of the marsh probably provided nitrogen to the nitrogen-limited open-water areas and triggered a phytoplankton bloom. Despite the rarity of runoff events, surface runoff from the watershed, combined with atmospheric deposition, contributed more than seven times the 530 kg (kilograms) of nitrogen that escaped Little Bean Marsh in surface outflow during 1997. Atmospheric deposition alone was more than 530 kg. Seepage to ground water contained less than 1.5 percent of the nitrogen leaving the marsh in surface outflow. The slow decay rate of Scirpus fluviatilis and reducing conditions in bottom sediments make burial of organic nitrogen a substantial sink of nitrogen. Denitrification experiments indicate that denitrification rates were limited by NO3 in the water column. Consequently, decomposition and nitrification of NH4 and organic nitrogen are the rate limiting steps of nitrogen removal in Little Bean Marsh. The NO3-limited rates of denitrification also indicate that Little Bean Marsh has a large unused capacity for nitrogen removal. These data indicate that the vast extent of riparian marshes along the Missouri and Mississippi Rivers may have had a substantial role in limiting NO3 loads to the Gulf of Mexico before agricultural development of flood plains. Drainage and removal of riparian marshes may be a major cause of the increased NO3 loads to the Gulf of Mexico. Periods of anoxia had much larger effects on phosphorus release than the other variables. The largest concentrations of phosphorus occurred in late summer and corresponded with senescing macrophytes, periods of anoxia, and a large algal bloom in Little Bean Marsh. Low water levels prevented the escape of phosphorus in surface outflow during these periods of highest phosphorus concentrations. Dry weather in late summer is typical and probably makes the correspondence of low water levels, anoxia, and consequent low phosphorus release a common occurrence in marshes along the Missouri River. Little Bean Marsh retained more than 95 percent of the phosphorus it received. The amount of phosphorus in surface inflows to the marsh were more than one order of magnitude greater than that escaping in surface outflows. The long hydraulic residence time of the marsh and large contributions of iron from ground water (that provide many sorption sites for phosphorus) make the marsh an effective sediment and phosphorus trap.

Missouri

Southern Salish Sea Habitat Map Series: Admiralty Inlet

In 2010 the Environmental Protection Agency, Region 10 initiated the Puget Sound Scientific Studies and Technical Investigations Assistance Program, designed to support research in support of implementing the Puget Sound Action Agenda. The Action Agenda was created in response to Puget Sound having been designated as one of 28 estuaries of national significance under section 320 of the U.S. Clean Water Act, and its overall goal is to restore the Puget Sound Estuary's environment by 2020. The Southern Salish Sea Mapping Project was funded by the Assistance Program request for proposals process, which also supports a large number of coastal-zone- and ocean-management issues. The issues include the recommendations of the Marine Protected Areas Work Group to the Washington State Legislature (Van Cleve and others, 2009), which endorses a Puget Sound and coast-wide marine conservation needs assessment, gap analysis of existing Marine Protected Areas (MPA) and recommendations for action. This publication is the first of four U.S. Geological Survey Scientific Investigation Maps that make up the Southern Salish Sea Mapping Project. The remaining three map blocks to be published in the future, located south of Admiralty Inlet, are shown in figure 1. Puget Sound is a deep, fjord-type estuary covering an area of 2,330 km 2 in the Pacific Northwest region of the United States (fig. 1). It is connected to the ocean by the Strait of Juan de Fuca, a turbulent passage approximately 160 km in length and 22 km wide at its west end, expanding to over 40 km wide at its east end (Thomson, 1994). During the Pleistocene, the area was occupied several times by lobes of continental ice, resulting in a complex basin-fill of glacial and interglacial deposits that are locally as thick as 1100 m (Johnson and others, 2001). The last glaciation, called the Fraser glaciation, began after 28,800±740 14 C yr B.P. when ice started a slow expansion (Clague, 1981). At peak advance the westward Juan de Fuca lobe reached the edge of the continental shelf through the Juan de Fuca Strait shortly before 14,460±200 14 C yr B.P. (Herzer and Bornhold, 1982). The southward Puget lobe advanced to its terminal position in Puget Sound by around 14,150 14 C yr B.P. (Porter and Swanson, 1998). Ice retreated from its maximum to northern Whidbey Island by 13,650±350 14 C yr B.P. (Dethier and others, 1995). Retreating glaciers resulted in a thick sequence of ice-contact, glacial-marine sediment, and early post-glacial sediments (Linden and Schurrer, 1988). These deposits have experienced the effects of a marine transgression followed by regression, resulting in a sea-level several tens of meters lower than the present day (Linden and Schurrer, 1988). A second transgression brought sea level to about the present level by around 5,470±120 14 C yr B.P. (Clague and others, 1982) establishing the present oceanographic and geologic environment Puget Sound is separated into four interconnected basins; Whidbey, Central (Main), Hood Canal, and South (Thomson, 1994). The Whidbey, Central, and Hood Canal basins are the three main branches of the Puget Sound estuary and are separated from the Strait of Juan de Fuca by a double sill at Admiralty Inlet. The Admiralty Inlet map area includes the Inlet and a portion of the Whidbey Basin (fig. 1). The shallower South Basin is separated by a sill at Tacoma Narrows and is highly branched with numerous finger inlets. Flow within Puget Sound is dominated by tidal currents of as much as 1 m/s at Admiralty Inlet, reducing to approximately 0.5 m/s in the Central Basin (Lavelle and others, 1988). The lack of silt and clay-sized sediments in the Admiralty Inlet map area is likely a result of the strong currents (see Ground-Truth Studies for the Admiralty Inlet Map Area, sheet 3). The subtidal component of flow reaches approximately 0.1 m/s and is driven by density gradients arising from the contrast in salty ocean water at the entrance and freshwater inputs from stream flow (Lavelle and others, 1988). The total freshwater input to Puget Sound is approximately 3.4 x 10 6 m 3 /day, primarily from the Skagit River (Cannon, 1983). The subtidal circulation mostly consists of a two-layered flow in the basins with fresher water exiting at the surface and saltier water entering at depth (Ebbesmeyer and Cannon, 2001). In general, surface waters flow north and deeper waters flow south; variations arise from wind effects that can drive a surface current in the same direction as the wind, and a baroclinic response in the lower layer to about 100-m depth (Matsuura and Cannon, 1997). Oceanographic properties are influenced by temporal forcing parameters such as reduced stream flow during the 2000-01 drought that increased surface salinity and decreased differences between surface and bottom waters (Newton and others, 2003). On offshore seismic-reflection profiles, Pleistocene strata (excluding latest Pleistocene glacial and post-glacial deposits) form a distinct seismic unit, bounded below by pre-Tertiary or Tertiary basement and above by typically flat-lying latest Pleistocene to Holocene deposits that fill in erosional or depositional relief (Johnson and others, 2001). Cores from central Puget Sound have accumulation rates that range from 85 to 1200 mg/cm 2 /yr, or 0.12 to 2.4 cm/yr; the highest accumulation rates are near the southern end of central Puget Sound (Carpenter and others, 1985). Carpenter and others (1985) un-weighted arithmetic mean of accumulation rates for central Puget Sound deeper stations is 480±340 (± one standard deviation) mg/cm 2 /yr. Lavelle and others (1985) also found rates as high as 1200 mg/cm 2 /yr over the past approximately 70 years in cores in the Central Basin off of and north and south of Elliott Bay. Puget Sound basin rates are comparable to rates in midshelf silt deposits on the Washington coast north of the Columbia River (Nittrouer and others, 1979). The deep subtidal (in other words, below SCUBA depths) habitats of Puget Sound are relatively poorly known. A few subtidal surveys exist for several habitat types from the 1960s and 1970s (reviewed in Dethier, 1990), using grab and box core data. The Dethier (1990) review divides habitat up into Coast and Marine Ecological Classification Standard (CMECS) substrate, water column energy, and depth zones but does not attempt to map these habitats, rather it is an inventory of habitats found in the area and the flora and fauna associated with each habitat. The approach of the Southern Salish Sea Mapping project is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, and bottom-sediment sampling data. This approach is based in part on methods presented and data collection and product needs identified at the Washington State Seafloor Mapping Workshop (Washington State Seafloor Mapping Workshop Steering Committee, 2008), attended by coastal and marine managers and scientists. The map products display seafloor geomorphology and substrate, and identify potential marine benthic habitats. It is emphasized that the more interpretive habitat and geology maps rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. Oceanographic current and wave data is not included in this analysis, however, the accompanying geographic information system (GIS) data set is designed and intended to be combined with oceanographic and biologic data sets assembled by others in the future and some of the GIS data has already been incorporated in the unpublished Nature Conservancy Benthic Habitats of Puget Sound database. This publication includes four map sheets, explanatory text, and a descriptive pamphlet. Each map sheet is published as a portable document format (PDF) file. ESRI ArcGIS compatible geotiffs (for example, bathymetry) and shapefiles (for example video observation points) will be available for download in the data catalog associated with this publication (Cochrane, 2015). An ArcGIS Project File with the symbology used to generate the map sheets is also provided. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html .

Washington

Effects of the earthquake of March 27, 1964, at Seward, Alaska

Seward, in south-central Alaska, was one of the towns most devastated by the Alaska earthquake of March 27, 1964. The greater part of Seward is built on an alluvial fan-delta near the head of Resurrection Bay on the southeast coast of the Kenai Peninsula. It is one of the few ports in south-central Alaska that is ice free all year, and the town’s economy is almost entirely dependent upon its port facilities. The Alaska earthquake of March 27, 1964, magnitude approximately 8.3–8.4, began at 6:36 p.m. Its epicenter was in the northern part of the Prince William Sound area; focal depth was 20–50 km. Strong ground motion at Seward lasted 3–4 minutes. During the shaking, a strip of land 50–400 feet wide along the Seward waterfront, together with docks and other harbor facilities, slid into Resurrection Bay as a result of large-scale submarine landsliding. Fractures ruptured the ground for'severa1 hundred feet back from the landslide scarps. Additional ground was fractured in the Forest Acres subdivision and on the alluvial floor of the Resurrection River valley; fountaining and sand boils accompanied the ground fracturing. Slide-generated wares, possibly seiche waves, and seismic sea waves crashed onto shore; ware runup was as much as 30 feet above mean lower low water and caused tremendous damage; fire from burning oil tanks added to the destruction. Damage from strong ground motion itself was comparatively minor. Tectonic subsidence of about 3.6 feet resulted in low areas being inundated at high tide. Thirteen people were killed and five were injured as a result of the earthquake. Eighty-six houses were totally destroyed and 260 were heavily damaged. The harbor facilities were almost completely destroyed, and the entire economic base of the town was wiped out. The total cost to replace the destroyed public and private facilities was estimated at $22 million. Seward lies on the axis of the Chugach Mountains geosyncline. The main structural trend in the mapped area, where the rocks consist almost entirely of graywacke and phyllite, is from near north to N. 20° E. Beds and cleavage of the rocks commonly dip 70° W. or NW. to near vertical. Locally, the rocks are complexly folded or contorted. So major faults were found in the mapped area, but small faults, shear zones, and joints are common. Surficial deposits of the area hare been divided for mapping into the following units: drift deposits, alluvial fan deposits, valley alluvium, intertidal deposits, landslide deposits, and artificial fill. Most of these units intergrade and were deposited more or less contemporaneously. The drift deposits consist chiefly of till that forms moraines along the lower flanks of the Resurrection River valley and up tributary valleys. The till is predominantly silt and sand and lesser amounts of clay-size particles, gravel, cobbles, and boulders. Glacial outwash and stratified ice-contact deposits constitute the remainder of the drift deposits. Fans and fan-deltas have been deposited at the valley mouths of tributary streams. Some, including the one upon which Seward built, project into Resurrection Bay, and deltaic-type deposits form their distal edges. The larger fans—composed chiefly of loosely compacted and poorly sorted silt, sand, and gravel—form broad aprons having low gradients. The fan deposits range in thickness from about 100 feet to possibly several hundred feet and, at least in some places, lie on a platform of compact drift. Smaller fans at the mouths of several canyons have steep gradients and considerable local relief. Valley alluvium, deposited chiefly by the Resurrection River, consists mostly of coarse sand and fine to medium gravel. In the axial part of the valley it is probably more than 100 feet thick. Near the head of Resurrection Bay, the alluvium is underlain by at least 75 feet of marine deltaic sediments, which are in turn underlain by 600 or more feet of drift in the deepest part of the bedrock valley. Beach, deltaic, and estuarine sediments, deposited on intertidal flats at the head of the bay and along far1 margins that extend into the bay, arc mapped as intertidal deposits. They consist mostly of silt, sand, and fine gravel, and lesser amounts of clay-size particles. The earthquake reactivated old slides and trigged new ones in the mountains. Rock and snow avalanches, debris flows, and creep of talus deposits characterized slide activity on the steeper slops. The Seward waterfront had been extended before the earthquake by adding artificial fill consisting of loose sand and gravel; part of the lagoon area had been filled with refuse. After the earth- quake, fill, consisting of silt and sand dredged from the head of the bay, was pumped onto part of the lagoon area and also on land at the northwest corner of the bay. Response to the disaster was immediate and decisive. City, State, and Federal agencies, as well as other organizations and individuals, gave unstintingly of their time and facilities. Within a few days, there was temporary restoration of water, sewerage, and electrical facilities. The U.S. Army Corps of Engineers was authorized to select sites and construct a new dock for the Alaska Railroad, a new small-boat basin, and related facilities. The firm of Shannon and Wilson, Inc., under contract to the Corps of Engineers, investigated subsurface soils extensively to determine the factors responsible for the sliding along the Seward waterfront and to assist in site selection for reconstruction of the destroyed harbor facilities. Borings also made along the Seward waterfront and at the head of the bay, and laboratory tests were conducted on pertinent samples. These studies were augmented by geophysical studies both on land and in the bay. In addition, the Corps of Engineers made shallow borings on the intertidal flats at the head of the bay and performed pile-driving and load tests. Borings also were drilled and test pits were dug in the subdivision of Forest Acres. Sliding along the Seward waterfront markedly deepened the water along the former shoreline. Post-earthquake slopes of the bay floor immediately offshore also are steeper in places than before the earthquake. The strong ground motion of the earthquake triggered the landsliding, but several factors may have contributed to the magnitude and characteristics of the slides. These factors are: (1) the long duration of strong ground motion, (2) the grain size and texture of the material involved in the sliding, (3) the probability that the finer grained materials liquefied and flowed seaward, and (4) the added load of manmade facilities built on the edge of the shore, Secondary effects of the slides themselves—sudden drawdown of water, followed by the weight of returning waves—also may have contributed to the destruction. Submarine sliding at the northwest corner of the bay occurred in fine-grained deltaic deposits whose frontal slopes probably were in metastable equilibrium under static conditions. Uplift pressures from aquifers under hydrostatic head, combined with the probable liquefaction characteristics of the sediments when vibrated by strong ground motion, probably caused the material to slide and flow seaward as a heavy slurry. Under static conditions, no major shoreline or submarine landsliding is expected in the Seward area; in the event of another severe earthquake, however, additional sliding is likely along the Seward waterfront and also in the deltaic deposits at the northwest corner of the bay. Fractured ground in back of the present shoreline along the Seward waterfront is an area of incipient landslides that would be unstable under strong shaking. For this reason the Scientific and Engineering Task Force placed the area in a high-risk classification and recommended no repair, rehabilitation, or new construction in this area involving use of Federal funds; it was further recommended that the area should be reserved for park or other uses that do not involve large congregations of people. The deltaic deposits at the head of the bay probably also would be susceptible to sliding during another large earthquake. This sliding would result in further landward retreat of the present shoreline toward the new railroad dock. Specifications for the new dock, whose seaward end is now approximately 1,100 feet from the back scarp of the subaqueous landslide, require design pro- visions to withstand seismic shock up to certain limits. Earthquake-induced fracturing of the ground in the subdivision of Forest Acres was confined to the lower part of a broad alluvial fan. There, sewer and water lines were ruptured and the foundations of some homes were heavily damaged. Landsliding, such as occurred along the shoreline of the bay, was not a contributing cause of the fracturing. Two hypotheses are offered to explain the fracturing: 1. Seismic energy was transformed into visible surface waves of such amplitude that the strength of surface layer was exceeded and rupturing occurred; tensional and compressional stresses alternately opened and closed the fractures and forced out water and mud. 2. Compaction by vibration of the fine-grained deposits of the fan caused ground settlement and fracturing; ground water under temporary hydrostatic head was forced to the surface as fountains and carried the finer material with it. Water waves that crashed onto shore, while shaking was still continuing, were generated chiefly by onshore and offshore landsliding. Waves that overran the shores about 25 minutes after shaking stopped and that continued to arrive for the next several hours are believed to be seismic sea waves (tsunamis) that originated in an uplifted area in the Gulf of Alaska. During the time of seismic sea-wave activity and perhaps preceding it, seiche wares also may have been generated within Resurrection Bay and complicated the wave effects along the shoreline.

Alaska

Modeling flow and water quality in reservoir and river reaches of the Mahoning River Basin, Ohio

Executive Summary The U.S. Army Corps of Engineers (USACE) is considering changes to the management of water surface elevation in four lakes in the Mahoning River Basin. These changes would affect the timing and amounts of water released to the Mahoning River and could affect the water quality of those releases. To provide information on possible water-quality effects from these operational changes, flow and water-quality models were constructed for Berlin Lake, Lake Milton, Michael J Kirwan Reservoir, Mosquito Creek Lake, Mosquito Creek, and the Mahoning River from the dams downstream to Lowellville, Ohio. The models were calibrated for two calendar years each, with model years selected depending on the availability of water-quality data. Models were developed with CE-QUAL-W2 version 4.2 (Wells, S.A., 2020, CE-QUAL-W2—A two-dimensional, laterally averaged, hydrodynamic and water quality model [version 4.2]: Portland State University, variously paged), a two-dimensional, laterally averaged hydrodynamic and water-quality model. Modeled constituents included flow, velocity, ice cover, water temperature, total dissolved solids (TDS), sulfate, chloride, inorganic suspended sediment, nitrate, ammonia, total Kjeldahl nitrogen, orthophosphate, total phosphorus, dissolved and particulate organic matter, algae, and dissolved oxygen. Iron was included for the lake models, but not the river. A whole-basin model, with the four lake models and river model, was used to run model scenarios to examine the effects of altered lake water surface elevations on flow and water quality in the lakes, the lake outflows, and the Mahoning River. The initial whole-basin model, with calendar year 2013 hydrology and measured or typical water quality, was designated as scenario 0. Mahoning River flows for calendar year 2013 were close to a 20-year median flow. Four additional scenarios were constructed based on reservoir operations model (RES-SIM) model water surface elevations for the four lakes as provided by USACE. Scenario 1 was the RES-SIM base case, scenario 2 kept Berlin Lake water surface elevations higher in summer, scenario 3 allowed 25 percent of summer flood storage to extend the guide curve, and scenario 4 allowed more flexibility in lake management by removing any downstream Mahoning River minimum flow requirements. The Mahoning River model was not changed in any scenarios but received altered flows from the lakes. Significant findings from this study include the following: In two of the four lakes (Berlin and Mosquito Creek Lakes), development of lake model grids using recent bathymetric surveys suggests that sedimentation in these lakes has occurred since they were constructed, altering volume-elevation curves. Tests of model parameter sensitivity showed that modeled water temperature, TDS, and dissolved oxygen were relatively insensitive to model parameter values. Modeled chlorophyll a , a measure of algal concentration, was most sensitive to parameter values; nitrate and total phosphorus concentrations were affected by a few of the parameters tested. As a group, the lake model results were more sensitive to model parameter values compared to the Mahoning River model. Data gaps were identified for inflows, both for water quantity and water quality, that could be filled through future sampling programs. Ample data were available from within the waterbodies for model calibration. The model simulated the general spatial and temporal patterns of water temperature, TDS, chloride, sulfate, nutrients, suspended sediment, organic matter, chlorophyll a , and dissolved oxygen in the lakes and Mahoning River. From late spring to autumn in the years modeled (2006, 2013, 2017–19 depending on the lake), all lakes developed thermal stratification and periods of anoxia in bottom waters. Stratification was most stable in Michael J Kirwan Reservoir and least stable in Mosquito Creek Lake. The stratification and anoxia in Berlin Lake, Lake Milton, and Mosquito Creek Lake could be interrupted by high-flow inputs moving through those lakes. The model predicted the release of ammonia and iron during anoxic periods in the lake hypolimnions. Concentrations of TDS, nitrate, orthophosphate, and total phosphorus increased in the Mahoning River down to Lowellville, the end of the river model, in the years modeled. These concentrations were greater than those in upstream lake releases. Chloride and sulfate concentrations were underpredicted in the Mahoning River, suggesting the presence of unreported chloride and sulfate inputs to the river, at least in the years modeled. Model scenario 4 kept water surface elevations the highest in all lakes in the April to mid-December period, compared to scenarios 1–3. Model scenario 2 kept water surface elevations in Berlin Lake higher in summer and late autumn, compared to scenarios 1 and 3, but to satisfy downstream minimum flow requirements, water surface elevations in the other lakes had periods of lower water surface elevation. As a group, scenarios 1–3 had largely similar effects on flow and water surface elevation in the Mahoning River because the lake releases in those scenarios still met downstream Mahoning River flow targets. Modeling the removal of downstream flow targets, scenario 4 had periods of lower flow in the Mahoning River from April to mid-September as water was held in the lakes, and periods of higher Mahoning River flow from mid-September through November as the lakes were drawn down to prepare for winter flood-risk management. In the four scenarios, all the lakes and lake outflows had generally similar seasonal cycles of water quality, though some differences were predicted. For instance, higher concentrations of iron and ammonia in the Lake Milton hypolimnion were modeled during a period of both low inflows from Berlin Lake and low outflows at Lake Milton dam. It is possible that those changes could be minimized by maintaining more flow or pulses of higher flow through the lake. Compared to the scenario 1 base case, changes to Mahoning River water quality were relatively minor for scenarios 2 and 3, which maintained downstream flows but shifted the flow source among the upstream lakes. The largest changes in Mahoning River water quality were predicted between Leavittsburg and Lowellville for scenario 4. The periods of lower lake outflows between April and mid-September led to correspondingly higher concentrations of TDS, orthophosphate, total phosphorus, and nitrate in the river, compared to the base case scenario 1. Conversely, the overall greater lake outflows from mid-September through November in scenario 4 led to periods of lower concentrations of TDS and nutrients in that portion of the river, at that time of year.

Ohio