USGS Science⌕ Search

SEARCH · USGS Science

Results for “Final Report”

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

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

At least 1,711 records · Page 95Linked to original sources

Lake-level frequency analysis for Devils Lake, North Dakota

Two approaches were used to estimate future lake-level probabilities for Devils Lake. The first approach is based on an annual lake-volume model, and the second approach is based on a statistical water mass-balance model that generates seasonal lake volumes on the basis of seasonal precipitation, evaporation, and inflow. Autoregressive moving average models were used to model the annual mean lake volume and the difference between the annual maximum lake volume and the annual mean lake volume. Residuals from both models were determined to be uncorrelated with zero mean and constant variance. However, a nonlinear relation between the residuals of the two models was included in the final annual lake-volume model. Because of high autocorrelation in the annual lake levels of Devils Lake, the annual lakevolume model was verified using annual lake-level changes. The annual lake-volume model closely reproduced the statistics of the recorded lake-level changes for 1901-93 except for the skewness coefficient However, the model output is less skewed than the data indicate because of some unrealistically large lake-level declines. The statistical water mass-balance model requires as inputs seasonal precipitation, evaporation, and inflow data for Devils Lake. Analysis of annual precipitation, evaporation, and inflow data for 1950-93 revealed no significant trends or long-range dependence so the input time series were assumed to be stationary and short-range dependent. Normality transformations were used to approximately maintain the marginal probability distributions; and a multivariate, periodic autoregressive model was used to reproduce the correlation structure. Each of the coefficients in the model is significantly different from zero at the 5-percent significance level. Coefficients relating spring inflow from one year to spring and fall inflows from the previous year had the largest effect on the lake-level frequency analysis. Inclusion of parameter uncertainty in the model for generating precipitation, evaporation, and inflow indicates that the upper lake-level exceedance levels from the water mass-balance model are particularly sensitive to parameter uncertainty. The sensitivity in the upper exceedance levels was caused almost entirely by uncertainty in the fitted probability distributions of the quarterly inflows. A method was developed for using long-term streamflow data for the Red River of the North at Grand Forks to reduce the variance in the estimated mean. Comparison of the annual lake-volume model and the water mass-balance model indicates the upper exceedance levels of the water mass-balance model increase much more rapidly than those of the annual lake-volume model. As an example, for simulation year 5, the 99-percent exceedance for the lake level is 1,417.6 feet above sea level for the annual lake-volume model and 1,423.2 feet above sea level for the water mass-balance model. The rapid increase is caused largely by the record precipitation and inflow in the summer and fall of 1993. Because the water mass-balance model produces lake-level traces that closely match the hydrology of Devils Lake, the water mass-balance model is superior to the annual lake-volume model for computing exceedance levels for the 50-year planning horizon.

Open-File Report↗

Lake-level frequency analysis for Devils Lake, North Dakota

Two approaches were used to estimate future lake-level probabilities for Devils Lake. The first approach is based on an annual lake-volume model, and the second approach is based on a statistical water mass-balance model that generates seasonal lake volumes on the basis of seasonal precipitation, evaporation, and inflow. Autoregressive moving average models were used to model the annual mean lake volume and the difference between the annual maximum lake volume and the annual mean lake volume. Residuals from both models were determined to be uncorrelated with zero mean and constant variance. However, a nonlinear relation between the residuals of the two models was included in the final annual lakevolume model. Because of high autocorrelation in the annual lake levels of Devils Lake, the annual lake-volume model was verified using annual lake-level changes. The annual lake-volume model closely reproduced the statistics of the recorded lake-level changes for 1901-93 except for the skewness coefficient. However, the model output is less skewed than the data indicate because of some unrealistically large lake-level declines. The statistical water mass-balance model requires as inputs seasonal precipitation, evaporation, and inflow data for Devils Lake. Analysis of annual precipitation, evaporation, and inflow data for 1950-93 revealed no significant trends or long-range dependence so the input time series were assumed to be stationary and short-range dependent. Normality transformations were used to approximately maintain the marginal probability distributions; and a multivariate, periodic autoregressive model was used to reproduce the correlation structure. Each of the coefficients in the model is significantly different from zero at the 5-percent significance level. Coefficients relating spring inflow from one year to spring and fall inflows from the previous year had the largest effect on the lake-level frequency analysis. Inclusion of parameter uncertainty in the model for generating precipitation, evaporation, and inflow indicates that the upper lake-level exceedance levels from the water mass-balance model are particularly sensitive to parameter uncertainty. The sensitivity in the upper exceedance levels was caused almost entirely by uncertainty in the fitted probability distributions of the quarterly inflows. A method was developed for using long-term streamflow data for the Red River of the North at Grand Forks to reduce the variance in the estimated mean. Comparison of the annual lake-volume model and the water mass-balance model indicates the upper exceedance levels of the water mass-balance model increase much more rapidly than those of the annual lake-volume model. As an example, for simulation year 5, the 99-percent exceedance for the lake level is 1,417.6 feet above sea level for the annual lake-volume model and 1,423.2 feet above sea level for the water mass-balance model. The rapid increase is caused largely by the record precipitation and inflow in the summer and fall of 1993. Because the water mass-balance model produces lake-level traces that closely match the hydrology of Devils Lake, the water mass-balance model is superior to the annual lake-volume model for computing exceedance levels for the 50-year planning horizon.

Water Supply Paper↗

Intrusions and intrusive complexes in an ophiolite near San Luis Obispo; a chemical and petrologic study

The San Luis Obispo ophiolite is a sequence of mafic and ultramafic rocks that lie as a thrust-like slice on top of Franciscan melange in the southern California Coast Range. The ophiolite comprises the classic Steinmann Trinity of basal serpentinized ultramafic rocks overlain sequentially by spilitized pillow lavas and cherts. The ophiolite is intruded by numerous diabasic dikes and several large intrusive complexes that contain peridotite, gabbro, diabase and felsic rocks. Field relations indicate that this ophiolite originated as a sequence of submarine volcanic flows and breccia units extruded onto an ultramafic basement. Later the volcanic and ultramafic rocks were intruded along their contact by mafic magma that crystallized as a stratiform complex of cumulate peridotite and gabbro. Next, this stratiform complex was intruded by diorite and albite-granite (trondhjemite) sills, which may have been produced by partial melting of gabbros in the upper levels of the stratiform complex. Finally all earlier units were intruded by diabasic sills and dikes that probably were late feeders for the pre-existing volcanic unit. All the ophiolite rocks underwent subsequent highly localized brecciation, under greenschist-facies conditions, in which felsic rocks were incipiently melted by frictional heating. The host rocks and intrusions were metamorphosed at low temperature and the mafic rocks were altered by Na- and Si-metasomatism. Despite this pervasive alteration, evidences of the formational stages of the ophiolite are preserved in disequilibrium mineral assemblages, primarily within mafic and ultramafic rocks of the intrusive complexes. Whole-rock analyses of the ophiolite suite indicate that alteration has obscured but not obliterated the original chemical character of the mafic rocks. The ophiolite diabases and basalts contain up to 6 weight percent Na 2 0 and 58 percent Si0 2 , in contrast to normal Si0 2 contents of 50 percent or less and Na 2 0 contents of less than 3 percent for unaltered oceanic and continental tholeiite basalts. However, the low K 2 0 contents (less than 0.8 percent) of all units in the ophiolite suite indicate that these rocks are not part of an alkalic kindred, but rather are altered tholeiites. Compositions of relict amphibole and plagioclase in the mafic rocks indicate that the ophiolite formed in a high-temperature environment and the deep seated rocks were metamorphosed under amphibolite-facies conditions. A later change in the environmental conditions allowed subsequent low-temperature (greenschist facies) metamorphism and metasomatism. The structure and chemistry of the San Luis Obispo ophiolite support the hypothesis that this body may represent a fragment of oceanic crust and mantle. The volcanic rocks and the intrusive complexes probably formed at a mid-ocean ridge, where conditions of high heat-flow provided the requisite high-temperature environment and promoted the production of mafic magmas. The change to low-temperature conditions probably was due to inception or acceleration of crustal spreading at that portion of the ridge where this ophiolite formed. The nature of brecciation in the ophiolite indicates that the low-temperature environment also was one of greater tectonic disturbance, relative to conditions at the ridge. This tectonically active environment probably corresponds to an island arc where the ophiolite-bearing plate was subducted in an adjacent trench. The absence of a high-pressure mineral assemblage in the San Luis Obispo ophiolite suggests that this segment of ocean-crust and upper mantle was not buried in a subduction zone, but rather was rafted on top of subduction melange (Franciscan Formation) into its present position.

California↗

Geology and paleontology of Canal Zone and adjoining parts of Panama: Description of Tertiary mollusks (additions to gastropods, scaphopods, pelecypods Nuculidae to Malleidae)

Chapter E adds 112 described species and subspecies (a few briefly described) to the some 440 covered in preceding chapters : 27 additional gastropods, 18 scaphopods, and 67 pelecypods in 10 families. It is estimated that about 125 species are to be added in chapter F, the final chapter. The Arcidae is by far the largest family in the pelecypods so far studied : 30 species in six genera. The genus Anadara is the largest in that family : 25 species, representing the subgenera Hawaiarca? , Rasia, Tosarca, Grandiarca, Potiarca , and Cunearca . Potiarca , based on a living western Pacific species, is adopted for many American species heretofore assigned to Cunearca . Eighteen of the 27 additional gastropods were found at a new locality in the late Eocene part of the Gatuncillo formation. They 'are assigned to 16 genera, 14 of which are not represented at other Gatuncillo localities. The genera include Faunus, Turnpanotonos , and Bezaneonia , which are rare in America, and the rare endemic American genus Harrisianella . The marine member of the Bohio (?) formation is now assigned to the late Eocene, instead of late Eocene or early Oligocene—the only change in age assignment. Eocene affinities are especially strengthened by the presence of Samanoetia samanensis and Volvariella? , if indeed the unnamed species represents that genus, as seems likely. Both the marine member of the Bohio( ?) and the late Oligocene part of the Bohio contain early small species of Anadara , subgenus Rasia . The fossils from the moderately deep-water facies of the late Oligocene Caimito formation include the earliest hexagonal Dentalium of the subgenus Dentalium s.s., from the American mainland, and a widespread and abundant species of the essentially deep-water subgenus Fissidentalium . Anadara chiriquiensis chiriquiensis , a relatively small form of the subgenus Grandiarca that is widely distributed in the Miocene Caribbean province, occurs in the Culebra formation. For the most part it indicates brackish water. Among the fossils from the La Boca formation proper, Cyphoma aff. C. intermedia is the earliest representative of an endemic American genus, and Isognomon mimeticus , a remarkable, narrowly mytiliform species is the type of the new subgenus Mimonion . Both the Culebra and La Boca, both of early Miocene age, contain predecessors of Gatun species. As in other chapters, the largest number of species is from the middle Miocene Gatun formation: 50 species. Bailya crossata is the first Tertiary species of the genus from the Caribbean region, and the earliest now known. A bizarre, incomplete gastropod is described as xancid?, genus?. The Gatun contains 13 of the 24 species of Anadara and half of the 12 species of the subgenus Rasia of that genus. In consideration of two battered and worn specimens of a large form of the subgenus Grandiarca , identified as Anadara grandis patricia? , the similar Tertiary forms in the Caribbean region are reviewed and divided into a brackish-water group and a marine group. The Gatun fauna now totals some 330 species.

Professional Paper↗

Assessment of water availability in the Osage Nation using an integrated hydrologic-flow model

The Osage Nation of northeastern Oklahoma, conterminous with Osage County, covers about 2,900 square miles. The area is primarily rural with 62 percent of the land being native prairie grass, and much of the area is used for cattle ranching and extraction of petroleum and natural gas. Protection of water rights are important to the Osage Nation because of its reliance on cattle ranching and the potential for impairment of water quality by petroleum extraction. Additionally, the potential for future population increases, demands for water from neighboring areas such as the Tulsa metropolitan area, and expansion of petroleum and natural-gas extraction on water resources of this area further the need for the Osage Nation to better understand its water availability. Therefore, the U.S. Geological Survey, in cooperation with the Osage Nation, completed a hydrologic investigation to assess the status and availability of surface-water and groundwater resources in the Osage Nation. A transient integrated hydrologic-flow model was constructed using the U.S. Geological Survey fully integrated hydrologic-flow model called the MODFLOW One-Water Hydrologic Model. The integrated hydrologic-flow model, called the Osage Nation Integrated Hydrologic Model (ONIHM), was constructed and uses an orthogonal grid of 276 rows and 289 columns, and each grid cell measures 1,312.34 feet (ft; 400 meters) per side, with eight variably thick vertical layers that represented the alluvial and bedrock aquifers within the study area, including the alluvial aquifer, the Vamoosa-Ada aquifer, and the minor Pennsylvanian bedrock aquifers, and the confining units. Landscape and groundwater-flow processes were simulated for two periods: (1) the 1950–2014 period from January 1950 through September 2014 and (2) the forecast period from October 2014 through December 2099. The 1950–2014 period ONIHM simulated past conditions using measured or estimated inputs, and the forecast-period ONIHM simulated three separate potential forecast conditions under constant dry, average, or wet climate conditions using calibrated input values from the 1950–2014 period ONIHM. The 1950–2014 period ONIHM was calibrated by linking the Parameter Estimation software (PEST) with the MODFLOW One-Water Hydrologic Model. PEST uses statistical parameter estimation techniques to identify the best set of parameter values to minimize the difference between measured or estimated calibration targets and their simulated equivalent values (residuals). Tikhonov regularization and singular-value decomposition-assist features of PEST were used during the calibration process. The 1950–2014 period ONIHM was calibrated to 713 measured groundwater levels at 195 wells; 95,636 estimated monthly mean groundwater levels at 124 wells; 5,307 measured streamflows at 13 streamgages; and 8,679 simulated mean monthly streamflows at 10 streamgages extracted from a surface-water model by adjusting 231 parameters. The estimated groundwater-level observations and streamflows were included as observations to improve the spatial and temporal density of observation targets during calibration. The best set of parameter values obtained during the calibration process of the 1950–2014 model was then used as the input parameter values for the forecast model simulations. A comparison of the calibration targets to their corresponding simulated values indicated that the model adequately reproduced streamflows and groundwater levels for some streamgages and wells and underestimated streamflows and groundwater levels at other locations. Measured and simulated streamflows correlated adequately with a coefficient of determination of 0.938, as did water levels with a coefficient of determination of 0.795. The 1950–2014 period ONIHM underestimated certain groundwater levels and streamflows, but generally measured or estimated calibration targets correlated well with simulated equivalents, which indicated that the model can adequately simulate the response of the hydrologic system to stresses in the 1950–2014 and forecast periods. In the 1950–2014 period ONIHM, the calibrated mean horizontal hydraulic conductivity for layer 1 alluvial aquifer was 30.7 feet per day, and the seven lower layers had a calibrated mean horizontal hydraulic conductivity of less than 3.3 feet per day. The mean calibrated groundwater-level residual was 16.6 ft, and the mean calibrated streamflow residual of the Arkansas River at Ralston, Oklahoma, streamgage (U.S. Geological Survey station 07152500) was within 6 percent (373 cubic feet per second) of mean measured streamflow for the 1950–2014 period ONIHM. The ONIHM simulated landscape fluxes of precipitation; groundwater applied by irrigation wells; evapotranspiration from precipitation, groundwater, and irrigation; runoff from precipitation; and deep percolation from precipitation. The largest loss of water from the landscape was evapotranspiration from precipitation with a calibrated mean annual outflow of 32 inches (in.): mean annual precipitation was about 36 in. Calibrated mean annual runoff and deep percolation (recharge to the water table) rates were 4.7 inches per year (in/yr) and 0.70 in/yr, respectively, for the 1950–2014 period ONIHM. The calibrated 1950–2014 period ONIHM groundwater fluxes included net farm net recharge (calculated as the difference between the inflow of recharge to the water table and the outflow of evapotranspiration from the water table such that negative values indicate that evapotranspiration from the water table was greater than deep percolation [recharge to the water table] and vice versa). Net farm net recharge was the largest flux from the groundwater system with a mean annual net outflow of 153.4 cubic feet per second. Stream leakage was the largest flux to the groundwater system with a mean annual net inflow of 152.5 cubic feet per second, indicating that, on average, the groundwater/surface-water interaction was a “losing” system where stream water leaked into the subsurface and recharged the water table. Simulated monthly trends demonstrated that net stream leakage was the largest inflow to the groundwater-flow system for 10 of the 12 months; for the other 2 months (January and March), farm net recharge (January) and net storage (March) were the largest inflow to the groundwater-flow system. A saline groundwater interface map was created for the study and compared to the water levels from the final stress period of the 1950–2014 model to identify the presence of fresh/marginal groundwater throughout the study area. Fresh/marginal groundwater was characterized as groundwater with less than 1,500 milligrams per liter of total dissolved solids. Fresh/marginal groundwater thickness ranged from 0 to 438.2 ft within the study area. The thickest regions of fresh/marginal groundwater were in the eastern part of the study area near Sand Creek, Bird Creek, and Hominy Creek and in the Arkansas River alluvial aquifer in the region downstream from the Arkansas River at Ralston, Okla. Like the 1950–2014 model, forecast model results for the landscape indicated that transpiration from precipitation was the largest flux out of the landscape for all three forecasts, constituting 77, 73, and 58 percent of precipitation for the dry, average, and wet forecasts, respectively. The dry and average forecast landscape fluxes demonstrated similar trends and magnitudes, whereas the wet forecast landscape fluxes indicated the largest changes compared to the average forecast fluxes. Most notably, runoff increased from a mean of 1.1 and 1.6 in/yr for the dry and average forecasts, respectively, to 10 in/yr for the wet forecast. Similar changes occurred for the other wet forecast landscape fluxes. The calibrated 1950–2014 period ONIHM simulated three forecasts to assess the effects of potential climatic changes on the hydrologic system from October 2014 to December 2099. The three forecasts simulated theoretical dry, average, and wet conditions using precipitation and potential evapotranspiration datasets from selected years in the calibrated 1950–2014 period ONIHM. Annual precipitation amounts were 26.89, 35.47, and 50.73 in. for the dry, average, and wet forecasts, respectively. Groundwater-flow component forecast results indicated that stream leakage is always a net inflow to the groundwater-flow system for dry, average, and wet conditions, meaning the study area stream network is always predominantly a “losing” regime where stream water infiltrates into the underlying aquifer. Storage was only a net outflow from the groundwater-flow system and indicated a replenishment to groundwater storage that resulted in an increase in groundwater levels only during the wet forecast. Further, these gains in groundwater storage for the wet forecast occurred only during February through June. Mean fresh/marginal groundwater saturated thicknesses were 125 and 126 ft for the dry and average forecast conditions, respectively, and wet forecast average thickness was 145 ft and ranged from 0 to 443 ft. The spatial extents of fresh/marginal groundwater at the end of the dry, average, and wet forecast model periods (December 2099) did not change substantially from the end of the 1950–2014 model period (September 2014).

Kansas, Oklahoma↗

Enhancements to the Mississippi Embayment Regional Aquifer Study (MERAS) groundwater-flow model and simulations of sustainable water-level scenarios

Arkansas continues to be one of the largest users of groundwater in the Nation. As such, long-term planning and management are essential to ensure continued availability of groundwater and surface water for years to come. The Mississippi Embayment Regional Aquifer Study (MERAS) model was developed previously as a tool to evaluate groundwater availability within the Mississippi embayment, which encompasses much of eastern Arkansas where the majority of groundwater is used. The Arkansas Water Plan is being updated for the first time since 1990 and serves as the State’s primary, comprehensive water-resources planning and guidance document. The MERAS model was selected as the best available tool for evaluation of specific water-use pumping scenarios that are currently being considered by the State of Arkansas. The model, developed as part of the U.S. Geological Survey Groundwater Resources Program’s assessment of the Nation’s groundwater availability, is proving to be invaluable to the State as it works toward development of a sustained yield pumping strategy. One aspect of this investigation was to evaluate multiple methods to improve the match of observed to simulated groundwater levels within the Mississippi River Valley alluvial and middle Claiborne (Sparta) aquifers in the MERAS model. Five primary methods were evaluated: (1) explicit simulation of evapotranspiration (ET), (2) upgrade of the Multi-Node Well (MNW2) Package, (3) geometry improvement within the Streamflow Routing (SFR) Package, (4) parameter estimation of select aquifer properties with pilot points, and (5) modification of water-use estimates. For the planning purposes of the Arkansas Water Plan, three scenarios were developed to evaluate potential future conditions: (1) simulation of previously optimized pumping values within the Mississippi River Valley alluvial and the middle Claiborne aquifers, (2) simulated prolonged effects of pumping at average recent (2000–5) rates, and (3) simulation of drawdown constraints on most pumping wells. The explicit simulation of ET indicated little, if any, improvement of model fit at the expense of much longer simulation time and was not included in further simulations. Numerous attempts to fully utilize the MNW2 Package were unsuccessful in achieving model stability, though modifications made to the water-use dataset remained intact. Final improvements in the residual statistics may be attributed to a single method, or a cumulative effect of all other methods (geometry improvement with the SFR Package, parameter estimation with pilot points, and modification of water-use estimates) attempted. The root mean squared error (RMSE) for all observations in the model is 22.65 feet (ft) over a range in observed hydraulic head of 741.66 ft. The RMSE for water-level observations in the Mississippi River Valley alluvial aquifer is 14.14 ft (an improvement of almost 3 ft) over a range in observed hydraulic head of 297.25 ft. The RMSE for the Sparta aquifer is 32.02 ft (an improvement of approximately 3 ft) over a range in observed hydraulic head of 634.94 ft. Three scenarios were developed to utilize a steady-state version of the MERAS model. Scenario 1 was developed to use pumping values resulting from the optimization of baseline rates (typically 1997 pumping rates) from previous optimization modeling of the alluvial aquifer and the Sparta aquifer. Scenario 2 was developed to evaluate the prolonged effects of pumping from the alluvial aquifer at recent pumping rates. Scenario 3A was designed to evaluate withdrawal limits from the alluvial aquifer by utilizing drawdown constraints equal to an altitude of approximately 50 percent of the predevelopment saturated thickness of the alluvial aquifer or 30 ft above the bottom of the alluvial aquifer, whichever was greater. The results of scenario 1 indicate large water-level declines throughout the area of the alluvial aquifer, regardless of the substitution of the optimized pumping values from earlier model simulations. The results of scenario 2 also indicate large areas of water-level decline, as compared to half of the saturated thickness, throughout the alluvial aquifer. The results of scenario 3A reveal some effects from the inclusion of multiple aquifers in a single simulation. The initial configuration of scenario 3A resulted in water levels well below the defined drawdown constraint, and some areas of depleted aquifer (water levels that are near or below the bottom of the aquifer) in east-central Arkansas. A fourth simulation (scenario 3B) was configured to apply the same drawdown constraints from the alluvial aquifer wells to the Sparta aquifer wells in the depleted area. These drawdown constraints reduce leakage from the alluvial aquifer to the underlying Sparta aquifer. This configuration did not produce depleted areas within the alluvial aquifer. Scenarios 3A and 3B indicate that even when pumping is limited in the alluvial aquifer, water levels in the alluvial aquifer may continue to decline in some areas because of pumping in the underlying Sparta aquifer.

Arkansas↗

Effects of urbanization on stream ecosystems along an agriculture-to-urban land-use gradient, Milwaukee to Green Bay, Wisconsin, 2003-2004

In 2003 and 2004, 30 streams near Milwaukee and Green Bay, Wisconsin, were part of a national study by the U.S. Geological Survey to assess urbanization effects on physical, chemical, and biological characteristics along an agriculture-to-urban land-use gradient. A geographic information system was used to characterize natural landscape features that define the environmental setting and the degree of urbanization within each stream watershed. A combination of land cover, socioeconomic, and infrastructure variables were integrated into a multi-metric urban intensity index, scaled from 0 to 100, and assigned to each stream site to identify a gradient of urbanization within relatively homogeneous environmental settings. The 35 variables used to develop the final urban intensity index characterized the degree of urbanization and included road infrastructure (road area and road traffic index), 100-meter riparian land cover (percentage of impervious surface, shrubland, and agriculture), watershed land cover (percentage of impervious surface, developed/urban land, shrubland, and agriculture), and 26 socioeconomic variables (U.S. Census Bureau, 2001). Characteristics examined as part of this study included: habitat, hydrology, stream temperature, water chemistry (chloride, sulfate, nutrients, dissolved and particulate organic and inorganic carbon, pesticides, and suspended sediment), benthic algae, benthic invertebrates, and fish. Semipermeable membrane devices (SPMDs) were used to assess the potential for bioconcentration of hydrophobic organic contaminants (specifically polycyclic aromatic hydrocarbons, polychlorinated biphenyls, and organochlorine and pyrethroid insecticides) in biological membranes, such as the gills of fish. Physical habitat measurements reflective of channel enlargement, including bankfull channel size and bank erosion, increased with increasing urbanization within the watershed. In this study, percentage of riffles and streambed substrate size were more strongly related to local geologic setting, slope, watershed topography, and river-engineering practices than to urbanization. Historical local river-engineering features such as channelization, bank stabilization, and grade controls may have confounded relations among habitat characteristics and urbanization. A number of hydrologic-condition metrics (including flashiness and duration of high flow during pre- or post-ice periods) showed strong relations to the urban intensity index. Hydrologic-condition metrics cannot be used alone to predict habitat or geomorphic change. Chloride and SPMD measures of potential toxicity and polycyclic aromatic hydrocarbon concentrations showed the strongest positive correlations to urbanization including increases in road infrastructure, percentage of impervious surface in the watershed, urban land cover, and land-distribution related to urban land cover. This suggests that automobiles and the infrastructure required to support automobiles are a significant source of these compounds in this study area. Chloride in spring and summer showed a significant positive correlation with the urban intensity index; concentrations increased with increasing road infrastructure, urban land cover, and a number of landscape variables related to urbanization. Spring concentrations of sulfate, prometon, and diazinon correlated to fewer urban characteristics than chloride, including increases in road infrastructure, percentage of impervious surface, and urban land cover. Changes in biological communities correlated to the urban intensity index or individual urban-associated variables. Decreased percentages of pollution-sensitive diatoms and diatoms requiring high dissolved-oxygen saturation correlated to increases in the percentage of developed urban land, total impervious surface, stream flashiness, population density, road-area density, and decreases in the percentage of wetland in the watershed. Invertebrate taxa richness and Coleop

Scientific Investigations Report↗

Performance evaluation testing of wells in the gradient control system at a federally operated Confined Disposal Facility using single well aquifer tests, East Chicago, Indiana

The U.S. Geological Survey (USGS) performed tests to evaluate the hydrologic connection between the open interval of the well and the surrounding Calumet aquifer in response to fouling of extraction well pumps onsite. Two rounds of air slug testing were performed on seven monitoring wells and step drawdown and subsequent recovery tests on three extraction wells on a U.S. Army Corps of Engineers Confined Disposal Facility (CDF) in East Chicago, Indiana. The wells were tested in 2014 and again in 2015. The extraction and monitoring wells are part of the gradient control system that establishes an inward gradient around the perimeter of the facility. The testing established a set of protocols that site personnel can use to evaluate onsite well integrity and develop a maintenance procedure to evaluate future well performance. The results of the slug test analysis data indicate that the hydraulic connection of the well screen to the surrounding aquifer material in monitoring wells on the CDF and the reliability of hydraulic conductivity estimates of the surrounding geologic media could be increased by implementing well development maintenance. Repeated air slug tests showed increasing hydraulic conductivity until, in the case of the monitoring wells located outside of the groundwater cutoff wall (MW–4B, MW–11B, MW–14B), the difference in hydraulic conductivity from test to test decreased, indicating the results were approaching the optimal hydraulic connection between the aquifer and the well screen. Hydraulic conductivity values derived from successive tests in monitoring well D40, approximately 0.25 mile south of the CDF, were substantially higher than those derived from wells on the CDF property. Also, values did not vary from test to test like those measured in monitoring wells located on the CDF property, which indicated that a process may be affecting the connectivity of the wells on the CDF property to the Calumet aquifer. Derived hydraulic conductivity values from the initial air slug test during the 2015 testing period for MW–11A and MW–14A are an order of magnitude less than those derived from the final test during the 2014 testing period indicating the development of a low conductivity skin between the final test of the 2014 testing period and the beginning of the 2015 testing period that created a decrease in the connection of the monitoring well screen to the surrounding aquifer material. Repeated step drawdown and recovery testing of the extraction wells tested during this study provided results that indicate a slight increase in the development of a skin and a decrease in the connectivity of the extraction wells with the Calumet aquifer. Hydraulic conductivity values obtained from the test results were relatively similar in EW–4B and EW–14A but were substantially lower for EW–11C. This difference may be due to the presence of finer grained silt deposits in the area surrounding well nest 11. Skin factors calculated during the step drawdown and recovery analysis were lowest in EW–11C and relatively similar in EW–4B and EW–14A. Calculated skin factors increased slightly in the analysis of data collected in 2015 from that collected in 2014. Comparisons of the specific-capacity values calculated from well development data collected following extraction well installation to those calculated during the single well aquifer tests at EW–4B, EW–14A and EW–11C indicate that the productivity of extraction wells on the CDF property has diminished since 2008. Values calculated for monitoring wells MW–4A, MW–11A, and MW–14A were used to evaluate the decrease in air slug derived hydraulic conductivity for monitoring wells within the groundwater cutoff wall between testing in 2014 and 2015. Results from testing by this study indicate that implementation of an air slug testing regimen of the monitoring wells that control the gradient control system at the CDF throughout the course of a year may help sustain the connectivity between the monitoring wells and the surrounding aquifer and provide data to evaluate the need for different types of well development approaches to address chemical or biological fouling issues. Repeated step drawdown and recovery testing of the extraction wells tested during this study provided results that indicate a slight increase in the development of a skin and a decrease in the connectivity of the extraction wells with the Calumet aquifer. Implementation of a specific capacity testing regimen can provide data to record and track well condition through time for individual extraction wells. Results from aquifer testing by this study indicate that specific capacity test results, when paired with recovery testing, provide useful data to measure the development of any low conductivity wellbore skin through the skin factors derived for the individual extraction wells. An initial annual schedule of specific capacity and recovery tests would provide sufficient data to identify substantial short-term changes in the operating condition of the extraction wells.

Indiana↗

Updating data inputs, assessing trends, and evaluating a method to estimate probable high groundwater levels in selected areas of Massachusetts

A method to estimate the probable high groundwater level in Massachusetts, excluding Cape Cod and the islands, was developed in 1981. The method uses a groundwater measurement from a test site, groundwater measurements from an index well, and a distribution of high groundwater levels from wells in similar geologic and topographic settings. The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, conducted an update to the Frimpter method for estimating the probable high groundwater levels in Massachusetts. The study evaluated the potential changes to the method resulting from four decades of additional groundwater-level data and the expansion of the network of wells for monitoring groundwater levels. The differences and potential benefits of daily, as opposed to monthly, measurements in the application of the method were examined because of the increased availability of high-frequency (subdaily) groundwater-level data. The study also considered long-term trends in groundwater levels that may alter the accuracy of the method. Finally, the accuracy of the estimated high groundwater levels was evaluated, and improved implementation guidance was prepared. For this study, groundwater levels in 153 wells in Massachusetts and surrounding States with records with lengths of 16 to 78 years were analyzed. The highest recorded groundwater levels ranged from 1.2 feet (ft) above land surface (flooded conditions) to 45.8 ft below land surface, with a median of 4.6 ft below land surface. The maximum annual groundwater-level range was 1.4 to 17.9 ft, with a median of 5.5 ft. The within-month variation, maximum annual groundwater-level range, and highest recorded groundwater level were computed using daily mean groundwater-level values from 28 wells with continuous records. The use of daily data resulted in larger maximum annual groundwater-level ranges (0.02 to 2.94 ft larger, with a median of 0.58 ft larger) and shallower highest-recorded groundwater levels (0.0 to 1.60 ft shallower, with a median of 0.18 ft shallower) than computations based on monthly measurements in the same wells. Statistical tests showed moderate to strong evidence of trends in measurements of both high and low groundwater levels within most of the periods during which water levels were analyzed. High groundwater levels rose beneath the land surface at most sites during four of the six periods used for analysis (1966–2015, 1986–2015, 1991–2010, and 1981–2010). Low groundwater levels also increased at many sites during most of the periods evaluated, but this trend was less widespread than the similar trends in high groundwater levels, and the trend was to deeper low groundwater levels at more sites than the trend to deeper high groundwater levels. There was no clear trend in annual groundwater-level ranges at most sites during the six periods analyzed. In general, the Frimpter method predicted shallower (higher) high groundwater levels than were observed but correctly classified sites according to their suitabilities for unmounded septic systems. The mean error of the predictions (difference between the estimated and observed groundwater levels) ranged from −3.23 ft to −1.40 ft for various approaches to estimating the groundwater-level range and selecting an index well. The method correctly classified 83 to 86 percent of monitoring-well sites according to their suitability for an unmounded septic system for many approaches to estimating the annual groundwater-level range and selecting an index well. The approach selected for estimating the annual groundwater-level range and selecting an index well will depend upon the importance of an accurate estimate of the high groundwater level as compared to the importance of an estimated high groundwater level that is less likely to be exceeded.

Connecticut, Massachusetts, New Hampshire, Rhode I↗

Porphyry copper assessment of western Central Asia

The U.S. Geological Survey conducted an assessment of resources associated with porphyry copper deposits in the western Central Asia countries of Kyrgyzstan, Uzbekistan, Kazakhstan, and Tajikistan and the southern Urals of Kazakhstan and Russia as part of a global mineral resource assessment. The purpose of the study was to (1) delineate permissive areas (tracts) for undiscovered porphyry copper deposits; (2) compile a database of known porphyry copper deposits and significant prospects; (3) where data permit, estimate numbers of undiscovered deposits within those permissive tracts; and (4) provide probabilistic estimates the amounts of copper (Cu), molybdenum (Mo), gold (Au), and silver (Ag) that could be contained in those undiscovered deposits. Western Central Asia, a region diverse in its geologic complexity, is situated north of the Tarim and North China tectonic blocks and sandwiched between the East European and Siberian cratons. The Ural Mountains form the western margin of the region; the southern margin is formed by the high-standing ranges that make up the Tian Shan mountain range in the border regions of Kazakhstan, Kyrgyzstan, and western China, where the effects of collisional tectonics are well displayed. The tectonic collage that makes up the core of western Central Asia is perhaps the geologically least understood part of the region. There is broad agreement that the early Paleozoic is made up of tectonically juxtaposed blocks that vary from Precambrian-cored microcontinents to magmatic arc and related complexes to subduction-related accretionary complexes. The rudiments of an incipient, contiguous single Kazakhstan block were formed by the end of the Silurian. In the middle to late Paleozoic, the block was unconformably superposed by two large, nested magmatic-arc belts, one Devonian, the other Carboniferous. Both magmatic-arc complexes were folded into a horseshoe-shaped, southeast-opening orocline in response to the final collisions of the various surrounding cratonic blocks with the Kazakhstan block. Additional deformation in the upper Cenozoic derived from the collision of India and China significantly redistributed fragments of the various mosaicked blocks, particularly in the central and southern parts of the western Central Asian region. Porphyry copper deposits are associated with many of the magmatic-arc fragments and belts throughout the geologically complex region, and economically important deposits are found in arc sequences of all Paleozoic Periods. The economically most productive arcs are Carboniferous. The assessment includes a discussion of the tectonic and geologic setting of porphyry copper deposits in western Central Asia (chapter 1), an application of remote sensing data for hydrothermal alteration mapping as a tool for porphyry copper assessment in the region (chapter 2), and a probabilistic assessment of undiscovered porphyry copper resources in four areas that represent Ordovician and Late Paleozoic (Carboniferous-Permian) magmatic arcs (chapter 3). The principal litho-tectonic terrane concept used to delineate permissive tracts was that of a magmatic arc that formed in the subduction boundary zone above a subducting plate. Eight permissive tracts are delineated on the basis of mapped and inferred subsurface distributions of igneous rocks assigned to magmatic arcs of specified age ranges that define areas where the occurrence of porphyry copper deposits within 1 kilometer of the Earth’s surface is possible. These tracts range in area from about 8,000 to 200,000 square kilometers and host 18 known porphyry copper deposits that contain about 54 million metric tons of copper. Available data included geologic maps, the distribution of significant porphyry copper occurrences and potentially related deposit types, the distribution of hydrothermal alteration patterns that are consistent with porphyry copper mineralization, and information on possible subsurface extensions of permissive rocks. On the basis of analyses of these data, the assessment team estimated a mean of 25 undiscovered porphyry copper deposits for the study area. Estimates of numbers of undiscovered deposits were combined with grade and tonnage models in a Monte Carlo simulation to yield a mean estimate of 95 million metric tons of copper in undiscovered porphyry copper deposits; this represents about twice the amount of identified porphyry copper resources (54 million metric tons). Detailed descriptions of each permissive tract, including the rationales for delineation and assessment, are given in appendixes, along with a geographic information system (GIS) that includes permissive tract boundaries, point locations of known porphyry copper deposits and significant occurrences, and hydrothermal alteration data based on analysis of remote sensing data.

Ural Mountains↗

Distribution, facies, ages, and proposed tectonic associations of regionally metamorphosed rocks in northern Alaska

Approximately half of the exposed bedrock in northern Alaska has been regionally metamorphosed. The most widespread metamorphic episode that affected northern Alaska occurred under low-grade, initially high-pressure (blueschist-facies) conditions during Mesozoic time. This episode is thought to have been related to the obduction of one or more oceanic terranes onto the continental margin of North America. Rocks whose metamorphism is considered to have been part of this major episode have an aerial distribution of approximately 10,000 km 2 in the southern Brooks Range, 5,000 km 2 across much of the Seward Peninsula, and 800 km 2 in the Ruby geanticline within the southeastern borderlands of the Yukon-Koyukuk basin. In the southern Brooks Range and on the Seward Peninsula, continental rocks experienced a clockwise pressure-temperature path that evolved during Middle Jurassic to late Early Cretaceous time from the low- to high-temperature subfacies of the blueschist facies and, finally, due to decreasing pressure, evolved to the greenschist facies. Metamorphism in the southern Brooks Range was associated with north-vergent compression along a south-dipping subduction zone that emplaced the oceanic rocks of the Angayucham terrane (represented by klippen of ultramafic rocks and prehnite-pumpellyite-facies metabasite, metatuff, metachert, and metasedimentary rocks) onto the continental margin. The present structural and metamorphic relation between the continental blueschistand greenschist-facies rocks and the structurally overlying lower temperature and pressure oceanic prehnite-pumpellyite-facies rocks to the south indicates that postmetamorphic or late metamorphic down-to-the-south, low-angle extensional faulting has dismembered the upper plate and removed much of the section that originally buried the blueschists. High-pressure metamorphism on the Seward Peninsula probably had a similar origin to that in the southern Brooks Range, but remnants of the overriding plate have not been identified, and the mechanism (such as left-lateral strike-slip faulting or oroclinal bending) by which the high-pressure rocks in the two areas were separated is not known. A significant difference between the thermal histories of these two metamorphic belts also is enigmatic: low-grade metamorphism on the Seward Peninsula was followed by intermediate-pressure amphibolite-facies and, locally, granulite-facies metamorphism as well as plutonism in mid-Cretaceous time, but no such high-temperature events have been documented within the high-pressure belt of the southern Brooks Range. In the Ruby geanticline, glaucophane, attesting to high-pressure metamorphism, is sporadically developed both within the continental rocks of the lower plate and, less commonly, near the base of the overlying oceanic thrust sheets. The direction from which the oceanic rocks were thrust and determination of which oceanic sheets were involved is unclear. Although the majority of the metamorphic episodes that affected northern Alaska occurred during the Mesozoic, older episodes have been documented or are suspected in a few areas. Late Proterozoic medium-grade metamorphism has been documented in a small area in the southwestern Brooks Range, Middle to Late Proterozoic or early Paleozoic metamorphism just south of the Ruby geanticline, and middle Paleozoic metamorphism in the eastern Brooks Range. High-temperature, high-pressure metamorphism of klippen in the Yukon-Tanana upland is tentatively considered to have occurred in early Paleozoic time. Metamorphic ages of units in several areas, particularly in the Ruby geanticline and the Yukon-Tanana upland, can be bracketed only between the probable Paleozoic age of their protoliths and the late Early Cretaceous age of postmetamorphic granitoids that intrude them.

Alaska↗

Behavior patterns and fates of adult steelhead, Chinook salmon, and coho salmon released into the upper Cowlitz River Basin, 2005–09 and 2012, Washington

A multiyear radiotelemetry evaluation was conducted to monitor adult steelhead ( Oncorhynchus mykiss ), Chinook salmon ( O. tshawytscha ), and coho salmon ( O. kisutch ) behavior and movement patterns in the upper Cowlitz River Basin. Volitional passage to this area was eliminated by dam construction in the mid-1960s, and a reintroduction program began in the mid-1990s. Fish are transported around the dams using a trap-and-haul program, and adult release sites are located in Lake Scanewa, the uppermost reservoir in the system, and in the Cowlitz and Cispus Rivers. Our goal was to estimate the proportion of tagged fish that fell back downstream of Cowlitz Falls Dam before the spawning period and to determine the proportion that were present in the Cowlitz and Cispus Rivers during the spawning period. Fallback is important because Cowlitz Falls Dam does not have upstream fish passage, so fish that pass the dam are unable to move back upstream and spawn. A total of 2,051 steelhead and salmon were tagged for the study, which was conducted during 2005–09 and 2012, and 173 (8.4 percent) of these regurgitated their transmitter prior to, or shortly after release. Once these fish were removed from the dataset, the final number of fish that was monitored totaled 1,878 fish, including 647 steelhead, 770 Chinook salmon, and 461 coho salmon. Hatchery-origin (HOR) and natural-origin (NOR) steelhead, Chinook salmon, and coho salmon behaved differently following release into Lake Scanewa. Detection records showed that the percentage of HOR fish that moved upstream and entered the Cowlitz River or Cispus River after release was relatively low (steelhead = 38 percent; Chinook salmon = 67 percent; coho salmon = 41 percent) compared to NOR fish (steelhead = 84 percent; Chinook salmon = 82 percent; coho salmon = 76 percent). The elapsed time from release to river entry was significantly lower for NOR fish than for HOR fish for all three species. Tagged fish entered the Cowlitz River in greater proportions than the Cispus River, regardless of origin. We found that 23–47 percent of the HOR fish entered the Cowlitz River and 12–38 percent entered the Cispus River. Similarly, 67–70 percent of the NOR fish entered the Cowlitz River and 38–66 percent entered the Cispus River. These behavioral differences translated into similar differences in fates during the spawning periods as higher percentages of tagged fish were assigned Cowlitz River fates than Cispus River fates. Fallback rates were affected by fish origin and release site. Overall, 12 percent of steelhead, 19 percent of Chinook salmon, and 8 percent of coho salmon fell back downstream of Cowlitz Falls Dam prior to spawning. Fallback rates were lower for fish that were released in the Cowlitz River or the Cispus River than for reservoir-released fish, but statistical comparisons were not robust because of small sample sizes at the river release sites. Fallback rates for fish released at the river release sites were 10 percent lower for steelhead, 4 percent lower for Chinook salmon, and 9 percent lower for coho salmon than for reservoir-released fish. However, fallback rates also were different between HOR and NOR fish. Fallback rates were significantly higher for HOR reservoir-released fish than for NOR reservoir-released fish. This study provided data that were insightful for understanding behavior and movement patterns in the upper Cowlitz River Basin and yielded estimates of fallback rates and fish fates that may be useful for fishery managers in the years to come. Studies from other systems have shown that factors such as prespawn mortality and fallback have resulted in substantial losses to spawning populations where trap-and-haul programs are being used as a restoration tool. Future research in the upper Cowlitz River Basin may use additional telemetry studies, genetic analyses, and spawning ground surveys to provide answers for new questions and to continue to monitor the progress of the reintroduction effort.

Washington↗

Hydrology-driven chemical loads transported by the Green River to the Lower Duwamish Waterway near Seattle, Washington, 2013–17

The sediments in the Lower Duwamish Waterway Superfund site in Seattle, Washington, are contaminated with chemicals including metals such as arsenic, polychlorinated biphenyls (PCBs), carcinogenic polycyclic aromatic hydrocarbons (cPAHs), and dioxins/furans from decades of intense anthropogenic activities. The U.S. Geological Survey, in cooperation with the Washington State Department of Ecology, collected new data from 2013 to 2017 to estimate sediment and chemical loads transported by the Green/Duwamish River to the Lower Duwamish Waterway Superfund site (the final 8-kilometer reach of the river) in support of sediment remediation within the site. Chemical loads were calculated as the product of river suspended-sediment loads and suspended sediment-bound chemical concentrations measured at river kilometer 16.7. Using four different approaches, annual suspended sediment-bound chemical load estimates transported by the river to the Lower Duwamish Waterway were in the range of 1,120–1,470 kilograms arsenic, 2,810–8,200 grams (g) toxic equivalent cPAHs, 205–407 milligrams toxic equivalent dioxins/furans, and 340–1,180 g PCBs. Storm events contributed a disproportionately large amount of the load of anthropogenic organic compounds such as cPAHs (54 percent), dioxins/furans (44 percent), and PCBs (52 percent) as compared to overall time (17 percent). Chemical concentrations and load estimates often were underestimated using results from unfiltered water analysis only, especially in samples with high suspended-sediment concentrations and for hydrophobic organic chemicals such as cPAHs that prefer to sorb to particulates and are at low concentrations near or below the analytical limits of water methods. For metals and PCBs, the dissolved concentration was relatively low and consistent between sampling events, whereas the suspended sediment-bound chemical concentrations contributed most of the chemical concentration in the water column during periods of high river suspended-sediment concentrations. However, the dissolved fraction, on average, contributed more than one-third of the estimated total chemical load in the river system for arsenic and PCBs, even given the hydrophobic nature of the chemicals. These results suggest that the sum of the chemical concentrations measured on two separate fractions—the particulate fraction and the dissolved fraction—more fully represents the total chemical concentration as compared to analysis of an unfiltered water sample, especially in samples with high suspended-sediment concentrations. Most of the suspended-sediment load (97 percent) and sediment-bound chemical load (92–94 percent) occurred during the wet winter half of the year from October 15 to April 14. However, the highest sediment-bound chemical concentrations often occurred during short intense storms or “first flush” autumn runoff events during the dry summer half of the year from April 15 to October 14. Because of the highly variable and dynamic river system characteristics (including precipitation, discharge, sediment concentration, and tidal fluctuations), it is critical to characterize the occurrence, frequency, concentrations, and loads during extreme conditions (for example, when the river is affected by storm-derived runoff) rather than time-averaged conditions. These short extreme events have a high potential for acute effects on ecological and human health, and may have a great influence on the effectiveness of the sediment remediation activities that are underway in the Lower Duwamish Waterway.

Washington↗

Tectonic map of the Circum-Pacific region, Pacific basin sheet

Circum-Pacific Map Project : The Circum-Pacific Map Project was a cooperative international effort designed to show the relationship of known energy and mineral resources to the major geologic features of the Pacific basin and surrounding continental areas. Available geologic, mineral, and energy-resource data are being complemented by new, project-developed data sets such as magnetic lineations, seafloor mineral deposits, and seafloor sediment. Earth scientists representing some 180 organizations from more than 40 Pacific-region countries are involved in this work. Six overlapping equal-area regional maps at a scale of 1:10,000,000 form the cartographic base for the project: the four Circum-Pacific Quadrants (Northwest, Southwest, Southeast, and Northeast), and the Antarctic and Arctic Sheets. There is also a Pacific Basin Sheet at a scale of 1:17,000,000. The Base Map Series and the Geographic Series (published from 1977 to 1990), the Plate-Tectonic Series (published in 1981 and 1982), the Geodynamic Series (published in 1984 and 1985), and the Geologic Series (published from 1984 to 1989) all include six map sheets. Other thematic map series in preparation include Mineral-Resources, Energy-Resources and Tectonic Maps. Altogether, more than 50 map sheets are planned. The maps were prepared cooperatively by the Circum-Pacific Council for Energy and Mineral Resources and the U.S. Geological Survey and are available from the Branch of Distribution, U. S. Geological Survey, Box 25286, Federal Center, Denver, Colorado 80225, U.S.A. The Circum-Pacific Map Project is organized under six panels of geoscientists representing national earth-science organizations, universities, and natural-resource companies. The six panels correspond to the basic map areas. Current panel chairmen are Tomoyuki Moritani (Northwest Quadrant), R. Wally Johnson (Southwest Quadrant), Ian W.D. Dalziel (Antarctic Region), vacant. (Southeast Quadrant), Kenneth J. Drummond (Northeast Quadrant), and George W. Moore (Arctic Region). Project coordination and final cartography was being carried out through the cooperation of the Office of the Chief Geologist of the U.S. Geological Survey, under the direction of General Chairman, George Gryc of Menlo Park, California. Project headquarters were located at 345 Middlefield Road, MS 952, Menlo Park, California 94025, U.S.A. The framework for the Circum-Pacific Map Project was developed in 1973 by a specially convened group of 12 North American geoscientists meeting in California. The project was officially launched at the First Circum-Pacific Conference on Energy and Mineral Resources, which met in Honolulu, Hawaii, in August 1974. Sponsors of the conference were the AAPG, Pacific Science Association (PSA), and the Coordinating Committee for Offshore Prospecting for Mineral Resources in Offshore Asian Areas (CCOP). The Circum-Pacific Map Project operates as an activity of the Circum-Pacific Council for Energy and Mineral Resources, a nonprofit organization that promotes cooperation among Circum-Pacific countries in the study of energy and mineral resources of the Pacific basin. Founded by Michel T. Halbouty in 1972, the Council also sponsors conferences, topical symposia, workshops and the Earth Science Series books. Tectonic Map Series : The tectonic maps distinguish areas of oceanic and continental crust. Symbols in red mark active plate boundaries, and colored patterns show tectonic units (volcanic or magmatic arcs, arc-trench gaps, and interarc basins) associated with active plate margins. Well-documented inactive plate boundaries are shown by symbols in black. The tectonic development of oceanic crust is shown by episodes of seafloor spreading. These correlate with the rift and drift sequences at passive continental margins and episodes of tectonic activity at active plate margins. The recognized episodes of seafloor spreading seem to reflect major changes in plate kinematics. Oceanic plateaus and other prominences of greater than normal oceanic crustal thickness such as hotspot traces are also shown. Colored areas on the continents show the ages of deformation and metamorphism of basement rocks and the emplacement of igneous rocks. Transitional tectonic (molassic) and reactivation basins are shown by a colored boundary, and if they are deformed, a colored horizontal line pattern indicates the age of deformation. Colored bands along basin boundaries indicate age of inception, and isopachs indicate thickness of platform strata on continental crust and cover on oceanic crust. Colored patterns at separated continental margins show the age of inception of rift and drift (breakup) sequences. Symbols mark folds and faults, and special symbols show volcanoes and other structural features. Affiliations are as of compilation of the data. This map was created in quadrants and then compiled together. They are the Northwest land, Northwest Marine (different compilers), Northeast, Southwest and Southeast, and parts in plate-boundary sections.

Circum-Pacific Map↗

Estimating peak-flow frequency statistics for selected gaged and ungaged sites in naturally flowing streams and rivers in Idaho

The U.S. Geological Survey, in cooperation with the Idaho Transportation Department, updated regional regression equations to estimate peak-flow statistics at ungaged sites on Idaho streams using recent streamflow (flow) data and new statistical techniques. Peak-flow statistics with 80-, 67-, 50-, 43-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities (1.25-, 1.50-, 2.00-, 2.33-, 5.00-, 10.0-, 25.0-, 50.0-, 100-, 200-, and 500-year recurrence intervals, respectively) were estimated for 192 streamgages in Idaho and bordering States with at least 10 years of annual peak-flow record through water year 2013. The streamgages were selected from drainage basins with little or no flow diversion or regulation. The peak-flow statistics were estimated by fitting a log-Pearson type III distribution to records of annual peak flows and applying two additional statistical methods: (1) the Expected Moments Algorithm to help describe uncertainty in annual peak flows and to better represent missing and historical record; and (2) the generalized Multiple Grubbs Beck Test to screen out potentially influential low outliers and to better fit the upper end of the peak-flow distribution. Additionally, a new regional skew was estimated for the Pacific Northwest and used to weight at-station skew at most streamgages. The streamgages were grouped into six regions (numbered 1_2, 3, 4, 5, 6_8, and 7, to maintain consistency in region numbering with a previous study), and the estimated peak-flow statistics were related to basin and climatic characteristics to develop regional regression equations using a generalized least squares procedure. Four out of 24 evaluated basin and climatic characteristics were selected for use in the final regional peak-flow regression equations. Overall, the standard error of prediction for the regional peak-flow regression equations ranged from 22 to 132 percent. Among all regions, regression model fit was best for region 4 in west-central Idaho (average standard error of prediction=46.4 percent; pseudo-R 2 >92 percent) and region 5 in central Idaho (average standard error of prediction=30.3 percent; pseudo-R 2 >95 percent). Regression model fit was poor for region 7 in southern Idaho (average standard error of prediction=103 percent; pseudo-R 2 <78 percent) compared to other regions because few streamgages in region 7 met the criteria for inclusion in the study, and the region’s semi-arid climate and associated variability in precipitation patterns causes substantial variability in peak flows. A drainage area ratio-adjustment method, using ratio exponents estimated using generalized least-squares regression, was presented as an alternative to the regional regression equations if peak-flow estimates are desired at an ungaged site that is close to a streamgage selected for inclusion in this study. The alternative drainage area ratio-adjustment method is appropriate for use when the drainage area ratio between the ungaged and gaged sites is between 0.5 and 1.5. The updated regional peak-flow regression equations had lower total error (standard error of prediction) than all regression equations presented in a 1982 study and in four of six regions presented in 2002 and 2003 studies in Idaho. A more extensive streamgage screening process used in the current study resulted in fewer streamgages used in the current study than in the 1982, 2002, and 2003 studies. Fewer streamgages used and the selection of different explanatory variables were likely causes of increased error in some regions compared to previous studies, but overall, regional peak‑flow regression model fit was generally improved for Idaho. The revised statistical procedures and increased streamgage screening applied in the current study most likely resulted in a more accurate representation of natural peak-flow conditions. The updated, regional peak-flow regression equations will be integrated in the U.S. Geological Survey StreamStats program to allow users to estimate basin and climatic characteristics and peak-flow statistics at ungaged locations of interest. StreamStats estimates peak-flow statistics with quantifiable certainty only when used at sites with basin and climatic characteristics within the range of input variables used to develop the regional regression equations. Both the regional regression equations and StreamStats should be used to estimate peak-flow statistics only in naturally flowing, relatively unregulated streams without substantial local influences to flow, such as large seeps, springs, or other groundwater-surface water interactions that are not widespread or characteristic of the respective region.

Idaho↗

High-temperature, large-volume, lavalike ash-flow tuffs without calderas in southwestern Idaho

Rhyolitic rocks were erupted from vents in and adjacent to the Owyhee Mountains and Owyhee Plateau of southwestern Idaho from 16 m.y. ago to about 10 m.y. ago. They were deposited on a highly irregular surface developed on a variety of basement rocks that include granitic rocks of Cretaceous age, quartz latite and rhyodacite tuffs and lava flows of Eocene age, andesitic and basaltic lava flows of Oligocene age, and latitic and basaltic lava flows of early Miocene age. The rhyolitic rocks are principally welded tuffs that, regardless of their source, have one feature in common-namely internal characteristics indicating en-masse, viscous lavalike flowage. The flowage features commonly include considerable thicknesses of flow breccia at the bases of various cooling units. On the basis of the tabular nature of the rhyolitic deposits, their broad areal extents, and the local preservation of pyroclastic textures at the bases, tops, and distal ends of some of the deposits, we have concluded that the rocks were emplaced as ash flows at extremely high temperatures and that they coalesced to liquids before final emplacement and cooling. Temperatures of l090?C and higher are indicated by iron-titanium oxide compositions. Rhyolites that are about 16 m.y. old are preserved mostly in the downdropped eastern and western flanks of the Silver City Range and they are inferred to have been erupted from the Silver City Range. They rarely contain more than about 2 percent phenocrysts that consist of quartz and subequal amounts of plagioclase and alkali feldspar; commonly, they contain biotite, and they are the only rhyolitic rocks in the area to do so. The several rhyolitic units that are 14 m.y. to about 10 m.y. old contain only pyroxene-principally ferriferous and intermediate pigeonites-as mafic constituents. The rhyolites of the Silver City Range comprise many cooling units, none of which can be traced for great distances. Rocks erupted from the Owyhee Plateau include two sequences that were traced over areas having diameters of about 100 km. These two sheets are the herein-named Swisher Mountain Tuff, which is about 13.8 m.y. old, and the Little Jacks Tuff, which is about 10 m.y. old. The Swisher Mountain Tuff was erupted from the Juniper Mountain volcanic center, a gentle dome that is not bounded by arcuate faults indicative of cauldron subsidence. The tuff is 200 m thick over a considerable area in and adjacent to its source. It apparently thins gradually toward its distal edges, and it is inferred to be uniformly distributed around its source at Juniper Mountain. The unit contains vitrophyres at various intervals from base to top, and, although the vitrophyres are, in general, flow layered and commonly flow brecciated, they occasionally contain well-defined pumice clasts. The vitrophyres indicate compound cooling, and, near the distal edges of the sheet, some of them probably represent complete cooling breaks. The Little Jacks Tuff onlaps the Swisher Mountain Tuff in expo sures east of Juniper Mountain, and it is inferred to have been erupted from a source on the part of the Owyhee Plateau that lies just east of the area studied. This inferred source area, like that at Juniper Mountain, is also expressed today as a gentle dome without structural features indicative of cauldron subsidence. The Little Jacks Tuff, in most exposures in the deep canyons of the Plateau, consists of at least four cooling units, and, in places in the eastern part of the studied area near the source area, it possibly comprises as many as six. Although there is no obvious evidence of erosion between the various cooling units, magnetic polarity measurements indicate that there were at least two magnetic reversals during the eruption interval of the Little Jacks Tuff. Like the Swisher Mountain Tuff, the Little Jacks has flattened pumice clasts in a few outcrops-principally at the bases of the various cooling units. The two tuff sequences are calc-a

Professional Paper↗

Water Quality of Combined Sewer Overflows, Stormwater, and Streams, Omaha, Nebraska, 2006-07

The U.S. Geological Survey, in cooperation with the City of Omaha, investigated the water quality of combined sewer overflows, stormwater, and streams in the Omaha, Nebraska, area by collecting and analyzing 1,175 water samples from August 2006 through October 2007. The study area included the drainage area of Papillion Creek at Capeheart Road near Bellevue, Nebraska, which encompasses the tributary drainages of the Big and Little Papillion Creeks and Cole Creek, along with the Missouri River reach that is adjacent to Omaha. Of the 101 constituents analyzed during the study, 100 were detected in at least 1 sample during the study. Spatial and seasonal comparisons were completed for environmental samples. Measured concentrations in stream samples were compared to water-quality criteria for pollutants of concern. Finally, the mass loads of water-quality constituents in the combined sewer overflow discharges, stormwater outfalls, and streams were computed and compared. The results of the study indicate that combined sewer overflow and stormwater discharges are affecting the water quality of the streams in the Omaha area. At the Papillion Creek Basin sites, Escherichia coli densities were greater than 126 units per 100 milliliters in 99 percent of the samples (212 of 213 samples analyzed for Escherichia coli) collected during the recreational-use season from May through September (in 2006 and 2007). Escherichia coli densities in 76 percent of Missouri River samples (39 of 51 samples) were greater than 126 units per 100 milliliters in samples collected from May through September (in 2006 and 2007). None of the constituents with human health criteria for consumption of water, fish, and other aquatic organisms were detected at levels greater than the criteria in any of the samples collected during this study. Total phosphorus concentrations in water samples collected in the Papillion Creek Basin were in excess of the U.S. Environmental Protection Agency's proposed criterion in all but four stream samples (266 of 270). Similarly, only 2 of 84 Missouri River samples had total phosphorus concentrations less than the proposed criterion. The proposed total nitrogen criterion for the Corn Belt and Northern Great Plains ecoregion was surpassed in 80 percent of the water samples collected from the stream sites. Samples with total nitrogen concentrations greater than the proposed criterion were most common at Papillion Creek and Big Papillion Creek sites, where the proposed criterion was surpassed in 90 and 96 percent of the samples collected, respectively. Elevated concentrations of total nitrogen were less common at the Missouri River sites, with 33 percent of the samples analyzed having concentrations that surpassed the proposed nutrient criterion for total nitrogen. The three constituents with measured concentrations greater than their respective health-based screening levels were nickel, zinc, and dichlorvos. Differences in water quality during the beginning, middle, and end of the combined sewer overflow discharge and the stream hydrograph rise, peak, and recession were investigated. Concentrations from the ending part of the combined sewer overflow hydrograph were significantly different than those from the beginning and middle parts for 3 and 11 constituents, respectively. No constituents were significantly different between the beginning and middle parts of the combined sewer overflow discharge hydrograph. For the stream site upstream from combined sewer overflow outfalls on Cole Creek, the constituents with geometric mean values for the hydrograph rise that were at least twice those for the values of the peak and recession were specific conductance, magnesium, nitrite, N,N-diethyl-meta-toluamide (DEET), methyl salicylate, p-cresol, and Escherichia coli. Similarly, the constituents where the hydrograph peak was at least twice that for the rise and recession at the upstream Cole Creek site were total suspended solids, silver, an

Scientific Investigations Report↗

Coastal sensitivity to sea level rise— A focus on the Mid-Atlantic Region

This Synthesis and Assessment Product (SAP), developed as part of the U.S. Climate Change Science Program, examines potential effects of sea-level rise from climate change during the twenty-first century, with a focus on the mid-Atlantic coast of the United States. Using scientific literature and policy-related documents, the SAP describes the physical environments; potential changes to coastal environments, wetlands, and vulnerable species; societal impacts and implications of sea-level rise; decisions that may be sensitive to sea-level rise; opportunities for adaptation; and institutional barriers to adaptation. The SAP also outlines the policy context in the mid-Atlantic region and describes the implications of sea-level rise impacts for other regions of the United States. Finally, this SAP discusses ways natural and social science research can improve understanding and prediction of potential impacts to aid planning and decision making. Projections of sea-level rise for the twenty-first century vary widely, ranging from several centimeters to more than a meter. Rising sea level can inundate low areas and increase flooding, coastal erosion, wetland loss, and saltwater intrusion into estuaries and freshwater aquifers. Existing elevation data for the mid-Atlantic United States do not provide the degree of confidence needed for local decision making. Systematic nationwide collection of high-resolution elevation data would improve the ability to conduct detailed assessments in support of planning. The coastal zone is dynamic and the response of coastal areas to sea-level rise is more complex than simple inundation. Much of the United States consists of coastal environments and landforms such as barrier islands and wetlands that will respond to sea-level rise by changing shape, size, or position. The combined effects of sea-level rise and other climate change factors such as storms may cause rapid and irreversible coastal change. All these changes will affect coastal habitats and species. Increasing population and development in coastal areas also affects the ability of natural ecosystems to adjust to sea-level rise. Coastal communities and property owners have responded to coastal hazards by erecting shore protection structures, elevating land and buildings, or relocating inland. Accelerated sea-level rise would increase the costs and environmental impacts of these responses. Shoreline armoring can eliminate the land along the shore to which the public has access; beach nourishment projects often increase access to the shore. Preparing for sea-level rise can be justified in many cases, because the cost of preparing now is small compared to the cost of reacting later. Examples include wetland protection, flood insurance, longlived infrastructure, and coastal land-use planning. Nevertheless, preparing for sea-level rise has been the exception rather than the rule. Most coastal institutions were based on the implicit assumption that sea level and shorelines are stable. Efforts to plan for sea-level rise can be thwarted by several institutional biases, including government policies that encourage coastal development, flood insurance maps that do not consider sea-level rise, federal policies that prefer shoreline armoring over soft shore protection, and lack of plans delineating which areas would be protected or not as sea level rises. The prospect of accelerated sea-level rise and increased vulnerability in coastal regions underscores the immediate need for improving our scientific understanding of and ability to predict the effects of sea-level rise on natural systems and society. These actions, combined with development of decision support tools for taking adaptive actions and an effective public education program, can lessen the economic and environmental impacts of sea-level rise.

Report↗