Search USGSSearch

SEARCH · Search USGS

Results for “Physics of Fluids”

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.

245 records · Page 14Linked to original sources

Can lava flow like water? Assessing applications of critical flow theory to channelized basaltic lava flows

Flowing lava and water have dramatically different physical properties but can form similar hydraulic structures, including undular hydraulic jumps, or standing wave trains. In water flows, undular hydraulic jumps are evidence of critical flow (Froude number ∼1) and open-channel hydraulic theory provides a powerful tool for estimating flow depth and velocity. Monitoring these parameters in an active lava channel is inherently challenging, but essential for calculating lava discharge (effusion rate), a primary control on the rate of flow front advance and ultimate flow runout distance. We analyze undular hydraulic jumps in both water and lava flows to assess the conditions under which they form and, by extension, the potential use of critical flow theory to estimate, in real time, lava flow velocity, depth, and discharge. Experimental data for water flows show that these structures mark the transition from supercritical to subcritical flow. Undular hydraulic jumps in the near-vent lava channel of the 2018 lower East Rift Zone eruption of Kīlauea, Hawaiʻi also reflect critical flow conditions; their wavelengths scale with flow depth and velocity, consistent with hydraulic theory. Calculated lava effusion rates are similar to estimates made using more traditional approaches (Jeffreys', 1925, https://doi.org/10.1080/14786442508634662 , equation based on lava viscosity, density, and channel slope) and with lava volumes derived from topographic-change mapping. From this we conclude that critical flow phenomena show great potential to track flow dynamics and inform hazard assessment for a wide range of geophysical fluids.

Journal of Geophysical Research - Earth Surface

Partially melted granodiorite and related rocks ejected from Crater Lake caldera, Oregon

Blocks of medium-grained granodiorite to 4 m, and minor diabase, quartz diorite, granite, aplite and granophyre, are common in ejecta of the ∼6,900 yr BP caldera-forming eruption of Mount Mazama. The blocks show degrees of melting from 0–50 vol%. Because very few have adhering juvenile magma, it is thought that the blocks are fragments of the Holocene magma chamber’s walls. Primary crystallisation of granodiorite produced phenocrystic pl + hyp + aug + mt + il + ap + zc, followed by qz + hb + bt + alkali feldspar (af). Presence of fluid inclusions in all samples implies complete crystallisation before melting. Subsolidus exchange with meteoric hydrothermal fluids before melting is evident in δ 18 O values of −3.4–+4.9‰ for quartz and plagioclase in partially melted granodiorites (fresh lavas from the region have δ 18 O values of +5.8–+7.0‰); δ 18 O values of unmelted granodiorites from preclimatic eruptive units suggest hydrothermal exchange began between ∼70 and 24 ka. Before eruption, the granitic rocks equilibrated at temperatures, estimated from Fe–Ti oxide compositions, of up to ∼1000°C for c . 10 2 –10 4 years at a minimum pressure of 100–180 MPa. Heating caused progressive breakdown or dissolution of hb, af, bt, and qz, so that samples with the highest melt fractions have residual pl + qz and new or re-equilibrated af + hyp + aug + mt + il in high-silica rhyolitic glass (75–77% SiO 2 ). Mineral compositions vary systematically with increasing temperature. Hornblende is absent in rocks with Fe–Ti oxide temperatures >870°C, and bt above 970°C. Oxygen isotope fractionation between qz, pl, and glass in partially fused granodiorite also is consistent with equilibration at T ≥900°C (Δ 18 O qz-pl = +0.7±0.5‰). Element partitioning between glass and crystals reflects the large fraction of refractory pl, re-equilibration of af and isolation or incomplete dissolution of accessory phases. Ba and REE contents of analysed glass separates can be successfully modelled by observed degrees of partial melting of granodiorite, but Rb, Sr and Sc concentrations cannot. Several samples have veins of microlite-free glass 1–5 mm thick that are compositionally and physically continuous with intergranular melt and which apparently formed after the climactic eruption began. Whole-rock H 2 O content, microprobe glass analysis sums near 100% and evidence for high temperature suggest liquids in the hotter samples were nearly anhydrous. The occurrence of similar granodiorite blocks at all azimuths around the 8 × 10 km caldera implies derivation from one pluton. Compositional similarity between granodiorite and pre-Mazama rhyodacites suggests that the pluton may have crystallised as recently as 0.4 Ma; compositional data preclude crystallisation from the Holocene chamber. The history of crystallisation, hydrothermal alteration, and remelting of the granitic rocks may be characteristic of shallow igneous systems in which the balance between hydrothermal cooling and magmatic input changes repeatedly over intervals of 10 4 –10 6 years.

Oregon

Partially melted granodiorite and related rocks ejected from Crater Lake caldera, Oregon

Blocks of medium-grained granodiorite to 4 m, and minor diabase, quartz diorite, granite, aplite and granophyre, are common in ejecta of the ∼6,900 yrBP calderaforming eruption of Mount Mazama. The blocks show degrees of melting from 0–50 vol%. Because very few have adhering juvenile magma, it is thought that the blocks are fragments of the Holocene magma chamber's walls. Primary crystallisation of granodiorite produced phenocrystic pl + hyp + aug + mt + il + ap + zc, followed by qz + hb + bt + alkali feldspar (af). Presence of fluid inclusions in all samples implies complete crystallisation before melting. Subsolidus exchange with meteoric hydrothermal fluids before melting is evident in δ 18 O values of −3·4+4·9‰ for quartz and plagioclase in partially melted granodiorites (fresh lavas from the region have δ 18 O values of +5·8−+7·0‰); δ 18 O values of unmelted granodiorites from preclimatic eruptive units suggest hydrothermal exchange began between ∼70 and 24 ka. Before eruption, the granitic rocks equilibrated at temperatures, estimated from Fe-Ti oxide compositions, of up to ∼1000°C for c. 10 2 –10 4 years at a minimum pressure of 100-180 MPa. Heating caused progressive breakdown or dissolution of hb, af, bt, and qz, so that samples with the highest melt fractions have residual pl + qz and new or re-equilibrated af + hyp + aug + mt + il in high-silica rhyolitic glass (75-77% SiO 2 ). Mineral compositions vary systematically with increasing temperature. Hornblende is absent in rocks with Fe-Ti oxide temperatures >870°C, and bt above 970°C. Oxygen isotope fractionation between qz, pl, and glass in partially fused granodiorite also is consistent with equilibration at T≥900°C (Δ 18 O qz.pl = +0·7 ± 0·5‰). Element partitioning between glass and crystals reflects the large fraction of refractory pl, re-equilibration of af and isolation or incomplete dissolution of accessory phases. Ba and REE contents of analysed glass separates can be successfully modelled by observed degrees of partial melting of granodiorite, but Rb, Sr and Sc concentrations cannot. Several samples have veins of microlite-free glass 1–5 mm thick that are compositionally and physically continuous with intergranular melt and which apparently formed after the climactic eruption began. Whole-rock H 2 O content, microprobe glass analysis sums near 100% and evidence for high temperature suggest liquids in the hotter samples were nearly anhydrous. The occurrence of similar granodiorite blocks at all azimuths around the 8 × 10 km caldera implies derivation from one pluton. Compositional similarity between granodiorite and pre-Mazama rhyodacites suggests that the pluton may have crystallised as recently as 0·4 Ma; compositional data preclude crystallisation from the Holocene chamber. The history of crystallisation, hydrothermal alteration, and remelting of the granitic rocks may be characteristic of shallow igneous systems in which the balance between hydrothermal cooling and magmatic input changes repeatedly over intervals of 10 4 -10 6 years.

Oregon

Teleseismic waves reveal anisotropic poroelastic response of wastewater disposal reservoir

Connecting earthquake nucleation in basement rock to fluid injection in basal, sedimentary reservoirs, depends heavily on choices related to the poroelastic properties of the fluid-rock system, thermo-chemical effects notwithstanding. Direct constraints on these parameters outside of laboratory settings are rare, and it is commonly assumed that the rock layers are isotropic. With the Arbuckle wastewater disposal reservoir in Osage County, Oklahoma, high-frequency formation pressure changes and collocated broadband ground velocities measured during the passing of large teleseismic waves show a poroelastic response of the reservoir that is both azimuthally variable and anisotropic; this includes evidence of static shifts in pressure that presumably relate to changes in local permeability. The azimuthal dependence in both the static response and shear coupling appears related to tectonic stress and strain indicators such as the orientations of the maximum horizontal stress and faults and fractures. Using dynamic strains from a nearby borehole strainmeter, we show that the ratio of shear to volumetric strain coupling is ∼ 0.41 "> ∼ 0.41 ∼0.41 which implies a mean Skempton's coefficient of A = 0.24 "> A = 0.24 A=0.24 over the plausible range of the undrained Poisson's ratio. Since these observations are made at relatively low confining pressure and differential stress, we suggest that the hydraulically conductive fracture network is a primary control on the coupling between pore pressure diffusion and elastic stresses in response to natural or anthropogenic sources.

Earth and Planetary Physics

Hydrogeologic setting and simulation of groundwater flow near the Canterbury and Leadville Mine Drainage Tunnels, Leadville, Colorado

The Leadville mining district is historically one of the most heavily mined regions in the world producing large quantities of gold, silver, lead, zinc, copper, and manganese since the 1860s. A multidisciplinary investigation was conducted by the U.S. Geological Survey, in cooperation with the Colorado Department of Public Health and Environment, to characterize large-scale groundwater flow in a 13 square-kilometer region encompassing the Canterbury Tunnel and the Leadville Mine Drainage Tunnel near Leadville, Colorado. The primary objective of the investigation was to evaluate whether a substantial hydraulic connection is present between the Canterbury Tunnel and Leadville Mine Drainage Tunnel for current (2008) hydrologic conditions. Altitude in the Leadville area ranges from about 3,018 m (9,900 ft) along the Arkansas River valley to about 4,270 m (14,000 ft) along the Continental Divide east of Leadville, and the high altitude of the area results in a moderate subpolar climate. Winter precipitation as snow was about three times greater than summer precipitation as rain, and in general, both winter and summer precipitation were greatest at higher altitudes. Winter and summer precipitation have increased since 2002 coinciding with the observed water-level rise near the Leadville Mine Drainage Tunnel that began in 2003. The weather patterns and hydrology exhibit strong seasonality with an annual cycle of cold winters with large snowfall, followed by spring snowmelt, runoff, and recharge (high-flow) conditions, and then base-flow (low-flow) conditions in the fall prior to the next winter. Groundwater occurs in the Paleozoic and Precambrian fractured-rock aquifers and in a Quaternary alluvial aquifer along the East Fork Arkansas River, and groundwater levels also exhibit seasonal, although delayed, patterns in response to the annual hydrologic cycle. A three-dimensional digital representation of the extensively faulted bedrock was developed and a geophysical direct-current resistivity field survey was performed to evaluate the geologic structure of the study area. The results show that the Canterbury Tunnel is located in a downthrown structural block that is not in direct physical connection with the Leadville Mine Drainage Tunnel. The presence of this structural discontinuity implies there is no direct groundwater pathway between the tunnels along a laterally continuous bedrock unit. Water-quality results for pH and major-ion concentrations near the Canterbury Tunnel showed that acid mine drainage has not affected groundwater quality. Stable-isotope ratios of hydrogen and oxygen in water indicate that snowmelt is the primary source of groundwater recharge. On the basis of chlorofluorocarbon and tritium concentrations and mixing ratios for groundwater samples, young groundwater (groundwater recharged after 1953) was indicated at well locations upgradient from and in a fault block separate from the Canterbury Tunnel. Samples from sites downgradient from the Canterbury Tunnel were mixtures of young and old (pre-1953) groundwater and likely represent snowmelt recharge mixed with older regional groundwater that discharges from the bedrock units to the Arkansas River valley. Discharge from the Canterbury Tunnel contained the greatest percentage of old (pre-1953) groundwater with a mixture of about 25 percent young water and about 75 percent old water. A calibrated three-dimensional groundwater model representing high-flow conditions was used to evaluate large-scale flow characteristics of the groundwater and to assess whether a substantial hydraulic connection was present between the Canterbury Tunnel and Leadville Mine Drainage Tunnel. As simulated, the faults restrict local flow in many areas, but the fracture-damage zones adjacent to the faults allow groundwater to move along faults. Water-budget results indicate that groundwater flow across the lateral edges of the model controlled the majority of flow in and out of the aquifer (79 percent and 63 percent of the total water budget, respectively). The largest contributions to the water budget were groundwater entering from the upper reaches of the watershed and the hydrologic interaction of the groundwater with the East Fork Arkansas River. Potentiometric surface maps of the simulated model results were generated for depths of 50, 100, and 250 m. The surfaces revealed a positive trend in hydraulic head with land-surface altitude and evidence of increased control on fluid movement by the fault network structure at progressively greater depths in the aquifer. Results of advective particle-tracking simulations indicate that the sets of simulated flow paths for the Canterbury Tunnel and the Leadville Mine Drainage Tunnel were mutually exclusive of one another, which also suggested that no major hydraulic connection was present between the tunnels. Particle-tracking simulations also revealed that although the fault network generally restricted groundwater movement locally, hydrologic conditions were such that groundwater did cross the fault network at many locations. This cross-fault movement indicates that the fault network controls regional groundwater flow to some degree but is not a complete barrier to flow. The cumulative distributions of adjusted age results for the watershed indicate that approximately 30 percent of the flow pathways transmit groundwater that was younger than 68 years old (post-1941) and that about 70 percent of the flow pathways transmit old groundwater. The particle-tracking results are consistent with the apparent ages and mixing ratios developed from the chlorofluorocarbon and tritium results. The model simulations also indicate that approximately 50 percent of the groundwater flowing through the study area was less than 200 years old and about 50 percent of the groundwater flowing through the study area is old water stored in low-permeability geologic units and fault blocks. As a final examination of model response, the conductance parameters of the Canterbury Tunnel and Leadville Mine Drainage Tunnel were manually adjusted from the calibrated values to determine if altering the flow discharge in one tunnel affects the hydraulic behavior in the other tunnel. The examination showed no substantial hydraulic connection. The multidisciplinary investigation yielded an improved understanding of groundwater characteristics near the Canterbury Tunnel and the Leadville Mine Drainage Tunnel. Movement of groundwater between the Canterbury Tunnel and Leadville Mine Drainage Tunnel that was central to this investigation could not be evaluated with strong certainty owing to the structural complexity of the region, study simplifications, and the absence of observation data within the upper sections of the Canterbury Tunnel and between the Canterbury Tunnel and the Leadville Mine Drainage Tunnel. There was, however, collaborative agreement between all of the analyses performed during this investigation that a substantial hydraulic connection did not exist between the Canterbury Tunnel and the Leadville Mine Drainage Tunnel under natural flow conditions near the time of this investigation.

Colorado

Effects of slip, slip rate, and shear heating on the friction of granite

The stability of fault slip is sensitive to the way in which frictional strength responds to changes in slip rate and in particular to the effective velocity dependence of steady state friction Δμ ss /Δ ln V . This quantity can vary substantially with displacement, temperature and slip rate. To investigate the physical basis for this behavior and the possible influence of shear heating, we slid initially bare granite surfaces in unconfined rotary shear to displacements of hundreds of millimeters at normal stresses, σ n of 10 and 25 MPa and at room temperature. We imposed step changes in slip rate within the range 10 −2 to 10 3.5 μm/s and also monitored frictional heating with thermistors embedded in the granite. The transient response of μ to slip rate steps was fit to a rate- and state-dependent friction law using two state variables to estimate the values of several parameters in the constitutive law. The first 20 mm of slip shows rising friction and falling Δμ ss /Δ ln V ; further slip shows roughly constant friction, Δμ ss /Δ ln V and parameter values, suggesting that a steady state condition is reached on the fault surface. At V ≤ 10 μm/s, Δμ ss /Δ ln V = −0.004 ± 0.001. At higher rates the response is sensitive to normal stress: At σ n = 25 MPa granite shows a transition to effective velocity strengthening (Δμ ss /Δ ln V = 0.008 ± 0.004) at the highest slip rates tested. At 10 MPa granite shows a less dramatic change to Δμ ss /Δ ln V ≈ 0 at the highest rates. The maximum temperature measured in the granite is ∼60°C at 25 MPa and 10 3.5 μm/s. Temperatures are in general agreement with a numerical model of heat conduction which assumes spatially homogeneous frictional heating over the sliding surface. The simplest interpretation of our measurements of Δμ ss /Δ ln V is that the granite is inherently velocity weakening (∂μ ss /∂ ln V < 0) and temperature strengthening (∂μ ss /∂ T −1 < 0) at all velocities. At high slip rates the response of μ to changes in temperature from shear heating may outweigh the response to changing velocity, such that the net effect Δμ ss /Δ ln V > 0 mimics velocity strengthening. These results have implications for the frictional behavior of faults during earthquakes. High slip rates may cause a switch to effective velocity strengthening which could limit peak coseismic slip rate and stress drop. For fluid-saturated faults, strengthening by this mechanism may be partly or fully offset by weakening due to thermal pressurization of a poorly drained pore fluid.

Journal of Geophysical Research B: Solid Earth

Results of an adaptive environmental assessment modeling workshop concerning potential impacts of drilling muds and cuttings on the marine environment

Drilling fluids or "muds" are essential components of modern drilling operations. They provide integrity for the well bore, a medium for removal of formation cuttings, and lubrication and cooling of the drill bit and pipe. The modeling workshop described in this report was conducted September 14-18, 1981 in Gulf Breeze, Florida to consider potential impacts of discharged drilling muds and cuttings on the marine environment. The broad goals of the workshop were synthesis of information on fate and effects, identification of general relationships between drilling fluids and the marine environment, and identification of site-specific variables likely to determine impacts of drilling muds and cuttings in various marine sites. The workshop was structured around construction of a model simulating fate and effects of discharges from a single rig into open water areas of the Gulf of Mexico, and discussion of factors that might produce different fate and effects in enclosed areas such as bays and estuaries. The simulation model was composed of four connected submodels. A Discharge/Fate submodel dealt with the discharge characteristics of the rig and the subsequent fate of discharged material. Three effects submodels then calculated biological responses at distances away from the rig for the water column, soft bottom benthos (assuming the rig was located over a soft bottom environment), and hard bottom benthos (assuming the rig was located over a hard bottom environment). The model focused on direct linkages between the discharge and various organisms rather than on how the marine ecosystem itself is interconnected. Behavior of the simulation model indicated relatively localized effects of drilling muds and cuttings discharged from a single platform into open water areas. Water column fate and effects were dominated by rapid dilution. Effects from deposition of spent mud and cuttings were spatially limited with relatively rapid recovery, especially in soft bottom benthic communities which were conceptualized as being adapted to frequent storms. This behavior was generated by the set of assumptions about linkages and functional relationships used to construct the model. Areas of uncertainty included methods for extrapolating 96-hr LC50 so results to exposures of varying lengths and concentrations; recovery rates of benthic communities; responses to various depths and rates of burial; fate and effects of the plume in relationship to stratification layers; and long-term and sub-lethal effects of slightly elevated concentrations of discharged materials. Evaluation of the assumptions of the Soft Bottom Submodel suggest that the assumptions used may have been relatively liberal estimates of resiliency of these communities. Discussion of "closed" water bodies such as bays and estuaries indicated several reasons to expect different and more complex fate and effects behavior in these areas. These factors included different species and communities (such as aquatic macrophytes and oyster beds), more complex circulation and stratification patterns, and potentially more active resuspension processes. Much of the possible difference in behavior in these areas centers around the extent to which they are “closed” or in the relative residence times of water and sediments in these areas as they determine the long-term dispersion of discharged material. Despite the complexity and variability of these areas, a large body of knowledge (such as that concerning fate and physical effects of dredge spoil) that could be effectively employed in analysis of potential fate and physical effects in enclosed areas was identified.

Report

Abrupt physical and chemical changes during 1992-1999, Anderson Springs, SE Geyser Geothermal Field, California

The Anderson Springs area is located about 90 miles (145 kilometers) north of San Francisco, California, in the southwestern part of Lake County. The area was first developed in the late 1800s as a health resort, which was active until the 1930s. Patrons drank a variety of cool to hot mineral waters from improved springs, swam in various baths and pools, and hiked in the rugged hills flanking Anderson Creek and its tributaries. In the bluffs to the south of the resort were four small mercury mines of the eastern Mayacmas quicksilver district. About 1,260 flasks of mercury were produced from these mines between 1909 and 1943. By the early 1970s, the higher ridges south and west of Anderson Springs became part of the southeast sector of the greater Geysers geothermal field. Today, several electric power plants are built on these ridges, producing energy from a vapor-dominated 240 °C reservoir. Only the main hot spring at Anderson Springs has maintained a recognizable identity since the 1930s. The hot spring is actually a cluster of seeps and springs that issue from a small fault in a ravine southwest of Anderson Creek. Published and unpublished records show that the maximum temperature (Tm) of this cluster fell gradually from 63°C in 1889 to 48°C in 1992. However, Tm of the cluster climbed to 77°C in 1995 and neared boiling (98°C) in 1998. A new cluster of boiling vents and small fumaroles (Tm = 99.3°C) formed in 1998 about 30 m north of the old spring cluster. Several evergreen trees on steep slopes immediately above these vents apparently were killed by the new activity. Thermal waters at Anderson Hot Springs are mostly composed of near-surface ground waters with some added gases and condensed steam from The Geysers geothermal system. Compared to gas samples from Southeast Geysers wells, the hot spring gases are higher in CO 2 and lower in H 2 S and NH 3 . As the springs increased in temperature, however, the gas composition became more like the mean composition of steam discharges from the Southeast Geysers. The hot spring waters are low in ions of Cl, B, and Li, but relatively high in HCO 3 , SO 4 and NH 4 . The stable-isotope compositions (deuterium and oxygen-18) of these waters plot near the global meteoric water line. Geochemical data through time reveal apparent maxima in the concentrations of SO 4 , Fe, and Mn in 1991 to 1992, before the cluster became hotter. The black-to-gray deposits from the new spring cluster are rich in pyrite and contain anomalous metals. About one-half mile to the east of the hot springs, mineralized water discharges intermittently from an old adit of the Schwartz (Anderson) mine, and enters a tributary of Anderson Creek. This drainage increased substantially in July 1998, and a slurry of mine water and precipitates were transported down the tributary and into Anderson Creek. In December 1998, the adit water was 22°C, and had a chemical composition that was similar to spring waters that once discharged in the ravines surrounding the old Anderson Springs resort. The cause for the abrupt changes that have occurred in thermal features at Anderson Springs is still not resolved. One possibility is that these changes are a response to withdrawal of steam from The Geysers geothermal field over more than 20 years of production. Pressure declines in the geothermal reservoir may have caused a "drying out" of the overlying condensation zone. Induced boiling in this zone and upflow of deep steam to shallower depths would cause heating and vaporization of shallow ground waters. In addition, earthquakes occurring in the vicinity of Anderson Springs have increased significantly after nearby geothermal power plants began operation. These earthquakes may have enhanced surface discharge of thermal fluids along fractures and faults.

California

Effects of the earthquake of March 27, 1964, on various communities

The 1964 earthquake caused wide-spread damage to inhabited places throughout more than 60,000 square miles of south-central Alaska. This report describes damage to all communities in the area except Anchorage, Whittier, Homer, Valdez, Seward, the communities of the Kodiak group of islands, and communities in the Copper River Basin; these were discussed in previous chapters of the Geological Survey's series of reports on the earthquake. At the communities discussed herein, damage resulted primarily from sea waves of diverse origins, displacements of the land relative to sea level, and seismic shaking. Waves took all of the 31 lives lost at those communities; physical damage was primarily from the waves and vertical displacements of the land relative to sea level. Destructive waves of local origin struck during or immediately after the earthquake throughout much of Prince William Sound, the southern Kenai Peninsula, and the shores of Kenai Lake. In Prince William Sound, waves demolished all but one home at the native village of Chenega, destroyed homesites at Point Nowell and Anderson Bay, and caused varying amounts of damage to waterfront facilities at Sawmill Bay, Latouche, Port Oceanic, Port Nellie Juan, Perry Island, and western Port Valdez. The local waves, which ran up as high as 70 feet above tide level at Chenega and more than 170 feet in several uninhabited parts of the Sound, took nearly all of the lives lost by drowning at these communities. Destructive local waves that devastated shores of Anderson Bay and adjacent parts of western Port Valdez probably were generated primarily by massive submarine slides of glacial and fluvioglacial deposits ; the origin of the waves that caused damage at most of the other communities and at extensive uninhabited segments of shoreline is not known. At these places the most probable generative mechanisms are: unidentified submarine slides of unconsolidated deposits, and (or) the horizontal tectonic displacements, of 20 to more than 60 feet, that occurred in the Prince William Sound region during the earthquake. A train of long-period seismic sea waves that began about 20 minutes after the start of the earthquake inundated shores along the Gulf of Alaska coast to a maximum height of 35 feet above tide level. At the communities described, they virtually destroyed two logging camps at Whidbey Bay and Puget Bay on the south coast of the Kenai Peninsula, caused moderate damage to boat harbors and docks at Seldovia and Cordova, floated away some beach cabins in the Cordova area, and drowned two people, one at Point Whitshed near Cordora and one at the Cape Saint Elias Light Station. The seismic sea waves were generated by regional tectonic uplift of the sea floor on the Continental Shelf. Vertical tectonic displacements of the land relative to sea level that accompanied the earthquake affected virtually all the coastal communities. Tectonic subsidence of 5 to 6 feet, augmented locally by surficial subsidence of unconsolidated deposits required either the relocation or raising of structures at Portage, Girdwood, and Hope on Turnagain Arm. Shoreline submergence resulting from about 3½ feet of tectonic subsidence at Seldovia necessitated raising all waterfront facilities and the airstrip above the level of high tides. On the other hand, tectonic uplift of the land in the Prince Williams Sound region required deepening of the small-boat harbors at Cordora and Tatitlek, dredging of the waterways in the Cordova area, and lengthening of some docks or piers at Cordova, the Cape Hinchinbrook Light Station, and in Sawmill Bay. Significant structural damage from direct seismic shaking was largely confined to fluid containers and a pier facility near Kenai. Indirect damage from fissuring and differential settling of foundation mterials in the vicinity of the Cordova airfield mused damage to a building, underground utilities, an airfield fill, and the highway. Minor amounts of direct and indirect damage from seismic vibrations were sustained by most of the communities situated on unconsolidated deposits as far east as Yakutat, north to Fairbanks, and west to King Salmon. Except for a few cracked or toppled chimney, all the damage from shaking was confined to areas of thick, unconsolidated deposits. Foundation damage was almost entirely restricted to water-saturated unconsolidated deposits which, when liquefied by seismic shaking, could spread laterally toward free faces and (or) settle differentially through compaction.

Alaska

Potential effects of energy development on environmental resources of the Williston Basin in Montana, North Dakota, and South Dakota—Water resources

The Williston Basin has been a leading oil and gas producing area for more than 50 years. While oil production initially peaked within the Williston Basin in the mid-1980s, production rapidly increased in the mid-2000s, largely because of improved horizontal (directional) drilling and hydraulic fracturing methods. In 2012, energy development associated with the Bakken Formation was identified as a priority requiring collaboration toward improved timeliness of issuing permits for new wells combined with reasonable measures to maintain environmental quality. Shortly thereafter, the Bakken Federal Executive Group was created to address common challenges associated with energy development. The Bakken Federal Executive Group partner agencies identified a gap in current understanding of the cumulative environmental challenges attributed to energy development throughout the area, resulting in an effort to aggregate scientific data and identify additional research and information needs related to natural resources within areas of energy development in the Williston Basin. As part of this effort, water resources in the area (including groundwater; streams and rivers; and lakes, reservoirs, and wetlands) were characterized and described in terms of physical occurrence, flow characteristics, recharge, water quality, and water use. Similarly, waters produced during energy-development activities also were characterized even though these waters are not considered usable resources within the area. Groundwater resources were characterized by the major hydrogeologic units, or aquifers, identifying the units that supply most groundwater used for domestic, stock, agricultural, and industrial purposes. The groundwater characterization included other deeper hydrogeologic units in the Williston Basin that may be a useable source of water with treatment, have utility as a reservoir for reinjection of produced waters, or be a source of minerals and energy resources. A generalized groundwater budget and flow system identifying the sources of recharge (stream infiltration, precipitation, and movement [leakage] from other aquifers) and the general groundwater flow direction is included for each of the major hydrogeologic units. Rivers and streams within the Williston Basin with 10 or more years of continuous streamflow data were identified. For a subset of these sites, streamflow characteristics, including the monthly and annual mean flow, were generated to identify seasonal and interannual changes in streamflow and thus provide information on the drivers and reliability of streamflow at the seasonal or multiyear scale. Daily streamflow and annual extreme flows (peak and low flow) also were estimated for the subset of sites. The daily streamflow and annual extreme flow values provide information on short-term or extreme events that are relevant to infrastructure design and evaluating spills, leaks, or accidental discharges of water or petroleum products. Surface-water features (lakes, ponds, and wetlands) were classified using the Cowardin system and identified on the National Wetlands Inventory maps generated by the U.S. Fish and Wildlife Service. The spatial distribution of the surface-water features was analyzed by State, county, and specifically in comparison to the Prairie Pothole Region. The proximity of the surface-water features to energy development infrastructure (specifically oil or gas well pads) was evaluated. It was determined that, although oil or gas wells are often near a surface-water feature, most surface-water features do not have wells nearby, with the exception of wells in the Prairie Pothole Region. Water-quality data were aggregated from two data sources: (1) the Water-Quality Portal, sponsored by the U.S. Geological Survey (USGS), U.S. Environmental Protection Agency (EPA), and National Water Quality Monitoring Council; and (2) a data compilation completed as part of the USGS National Water-Quality Assessment project. The Water-Quality Portal integrates publicly available water-quality data from databases maintained by the USGS, EPA, and U.S. Department of Agriculture, including water-quality data from Tribal, State, and local databases. Water-quality data for 15 commonly measured water-quality constituents were aggregated for groundwater, rivers and streams, and lakes and reservoirs. For each aggregated dataset (groundwater, rivers and streams, and lakes and reservoirs), analyses of the water-quality data included summary statistics, maps of spatial distribution of constituent values, boxplots of constituent values by timeframe or hydrogeologic unit, spatial comparisons of site locations and constituent values to petroleum well density, and comparisons of the constituent values measured to EPA drinking-water standards/guidelines. Produced water includes all fluids brought to the surface along with the targeted hydrocarbons as part of the oil and gas exploration and extraction processes. These fluids may include formation water (waters that co-exist with rock/oil/gas), hydraulic fracturing fluids, and other combinations of water and chemicals used during oil and gas well drilling, development, treatments, recompletions, and workovers. Produced water datasets were aggregated from two sources: the USGS National Produced Waters Geochemical database (ver. 2.1) and a series of projects focused specifically on sampling produced water in the Williston Basin from 2010 to 2014. The National Produced Waters Geochemical database was useful for a general understanding of produced-water chemistry. Produced waters are characterized by extreme salinity and contain elevated concentrations of other constituents (including arsenic, barium, cadmium, lead, zinc, radium-226/radium-228, and ammonium) that could negatively affect water and aquatic resources if released. Produced waters also have a generally unique chemical (isotopic) signature that may be useful in tracking water from different geologic units; for example, the oxygen/deuterium and strontium ratio values measured in brine waters from the Bakken Formation are distinct from brines collected from other geologic units in the Williston Basin. Water-use information related to energy production in the area also was aggregated and summarized. The summary of water use is not limited to oil and gas production but includes water used to produce all types of energy resources in the Williston Basin, including coal/lignite, thermoelectric power, oil and gas, hydropower, biomass and biofuels, wind, geothermal, and solar. Each State has its own methods for regulating and reporting water usage within its jurisdiction. These methods can introduce problems when examining water use from sources, such as the Missouri River or Fox Hills aquifer, that are shared across political boundaries. Without the one-to-one match for usage types and amounts used from a water source, it is difficult to develop a comprehensive water budget for the water source being evaluated. A large amount of freshwater is required to prepare a well for oil and gas well production; in some cases, 3 to 7 million gallons of water are needed per well. The EPA estimates that hydraulic fracturing in the Williston Basin uses between 70 to 140 billion gallons per year. Water also is used for myriad other purposes related to ancillary oil and gas extraction. In addition to water used for immediate energy development, the expanded human workforce migrating into the area and other support staff who have moved into the area during the development also use water. Research and information needs were identified that could be relevant in the evaluation of the effects of energy development on water resources. Information needs related to the evaluation of groundwater resources include the following: improved potentiometric-surface maps for glacial units; availability of a uniform stream network digital geographic coverage that spans the international boundary with Canada; enhanced surface-water use information with regards to the gain and loss of streamflow to shallow groundwater, which would increase understanding groundwater and surface-water interactions; and expanded geophysical assessments. Gaps in the availability of streamflow data include the lack of information on ice-jam flooding despite potential for effects to infrastructure (pipelines, roads, and facilities) and an understanding of the cumulative effects of largely undocumented stock and diversion dams. Although this study resulted in the aggregation of a large quantity of water-quality data, the availability of consistently collected, systematically processed and reported data over large parts of the Williston Basin is sparse. Few samples have been analyzed for constituents that may indicate the effect of energy development on water resources. Constituents that could be considered include boron, chloride, bromide, iodine, fluoride, manganese, lithium, radium, strontium isotopes, volatile organic compounds, and isotopes of inorganic ions (such as hydrogen and carbon). Collaboration between Tribal, Federal, State, and local entities to identify a common study design, common monitoring constituents, and consistent sampling locations would generate datasets with broad utility and would likely result in overall cost savings for monitoring over time. Similarly, there is a need for standardized sample collection, processing, laboratory analytical methods, and the collection of ancillary data for produced waters sampling. Additional characterization of the range of chemical, microbial, and isotopic compositions and quantities of “end-member” produced waters, and the collection of time-series datasets to document the changes in produced waters during and after well development also were needs identified during this study. Water-use estimates would be improved through the implementation of comprehensive studies of water use from groundwater and surface-water sources using consistent methodologies across the Williston Basin. The submission of chemical and water data related to hydraulic fracturing collected by the oil and gas industry would add to the quantity of available data. Consistent implementation of regulations and monitoring controls across political boundaries (State, county, and international) would further improve the consistency of data available for the estimates of water use.

Montana, North Dakota, South Dakota

Geology and regional metamorphism of some high-grade cordierite gneisses, Front Range, Colorado

Cordierite is common in regional metamorphic gneisses of Precambrian age in the central part of the Front Range. It occurs in discontinuous stratigraphic units that are structurally a minor component, except locally, of the thick succession of biotite gneisses that comprise the widespread Idaho Springs Formation. The rocks have mineral assemblages, that are characteristic of the sillimanite grade of metamorphism. The cordierite occurs in three principal rock types: (1) potassic feldspar-bearing cordierite-garnet-sillimanite-biotite gneiss, (2) cordierite-biotite gneiss, and (3) cordierite-gedrite-biotite gneiss; each type contains several characteristic mineral assemblages. The rock types are gradational and overlap in areal distribution, and mainly owe their diversity in mineralogy to differences in bulk chemical composition. The field relations are consistent with an interpretation that the diverse cordierite rocks were derived from original sedimentary rocks, largely pelitic sediments. The potassic feldspar-bearing cordierite-garnet gneisses were formed from shales that contained more MgO and FeO than the more abundant sedimentary facies that yielded sillimanitic biotite gneisses. Cordierite-gedrite-biotite gneisses contain much aluminum, iron, and magnesium and little sodium and potassium as compared to the other biotite gneisses; they have an extremely low content of minor elements. Although their chemical compositions are unlike those of known modern sediments, the cordierite-gedrite gneisses are considered also to have been derived from sedimentary rocks. The physical properties and chemical compositions of the mineral phases vary somewhat from one rock type to another. Biotite varies systematically in composition, and the changes are closely related to rock type and thus to bulk composition; the MgO/FeO ratios range from 1.7 in the more mafic cordierite-gedrite rocks to 0.49 in potassic feldspar-bearing cordierite-garnet gneisses. Cordierite is magnesium-rich and intermediate in the range of composition of all analyzed cordierites (Leake, 1960); its MgO/FeO ratio is higher in the gedrite-bearing gneisses than in the potassic feldspar-bearing gneisses. The garnets consist dominantly of the almandine and pyrope molecules, and range from 64 to 75 percent almandine and from 14 to 27 pyrope. These crystals are zoned; their rims are slightly more ferrous and less magnesian than their cores. Both monoclinic and triclinic alkali feldspars coexist in the potassic feldspar-bearing cordierite-garnet gneisses. The potassic feldspars contain from 18 to 27 weight percent NaAlSi 3 O 8 . Plagioclase (oligoclase-andesine) is uncommon in the rocks. Gedrite has an MgO/FeO ratio ranging from 1 to 1.2. Associated minor minerals include iron oxides, andalusite, spinel and its alteration product högbomite, and corundum. The mineral assemblages can be correlated imperfectly with episodes of deformation and metamorphism. Relict staurolite and associated garnet occur locally as remnants of an assemblage formed early in regional metamorphism, presumably early in the first period of deformation. The dominant assemblage biotite-cordierite-garnet-magnetite-plagioclase-potassic feldspar-quartz-sillimanite and associated assemblages having fewer phases, were formed during period one and period two deformations, the principal episodes of regional dynamothermal metamorphism in the central part of the Front Range. A minor assemblage andalusite-biotite-magnetite-plagioclase-quartz was formed later, possibly coincident with a third period of deformation, largely cataclastic in effects, which was more local than the earlier deformations and metamorphism. Phase equilibria studies of the assemblage biotite-cordierite-garnet-magnetite-plagioclase-potassic feldspar-quartz-sillimanite and associated assemblages are interpreted to indicate that the cordierite assemblages approach a state of chemical equilibrium. The scatter of points in a distribution diagram can be interpreted in terms of at least two sets of equilibrium conditions that prevailed during the major plastic deformations. Other discrepancies indicating departure from a homogeneous equilibrium can be explained as a result of mosaic equilibrium involving limited diffusion of iron and magnesium for short distances. The mineral assemblages and the compositions of the ferromagnesian minerals in the cordierite rocks of this region are dependent primarily on the bulk composition of the rocks and variations in the mineral species that comprise the rocks and, to a lesser degree, on the grade of metamorphism. Biotite and cordierite are markedly more magnesian in the more mafic cordierite-gedrite-biotite gneiss than in the potassic feldspar-bearing cordierite-garnet-sillimanite-biotite gneiss. Associated microcline gneiss and biotite-sillimanite gneiss that contains muscovite as a primary stable mineral provides a means to define the metamorphic grade in the area of study. It is concluded from analyses of the assemblages with respect to theoretical phase relations in the system SiO 2 -Al 2 O 3 -Na 2 O-K 2 O-H 2 O that at least some of the rocks in the Central City-Nederland area are above the sillimanite-potassic feldspar isograd as defined by Evans and Guidotti (1966). In rocks of appropriate composition, muscovite is a stable phase in assemblages containing potassic feldspar and sillimanite. The cordierite assemblages and associated rocks are inferred to have formed in an environment having a load pressure of 3–5 kilobars (fluid pressure equaled load pressure) and a temperature somewhat in excess of 620° C.

Colorado