SEARCH · Geology Reports

Results for “Interfaces”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,315 records · Page 73Linked to original sources

USGS2009

Benthic flux of nutrients and trace metals in the northern component of San Francisco Bay, California

Two sets of sampling trips were coordinated in late summer 2008 (weeks of July 8 and August 6) to sample the interstitial and overlying bottom waters at 10 shallow locations (9 sites <3 meters in depth) within the northern component of the San Francisco Bay/Delta (herein referred to as North Bay). The work was performed to better understand sources of biologically reactive solutes (namely, dissolved macronutrients and trace metals) that may affect the base of the food web in this part of the estuary. A nonmetallic pore-water profiler was used to obtain the first centimeter-scale estimates of the vertical solute-concentration gradients for diffusive-flux determinations. This study, performed in collaboration with scientists from San Francisco State University?s Romberg Tiburon Center for Environmental Studies, provides information to assist in developing and refining management strategies for the Bay/Delta system and supports efforts to monitor changes in food-web structure associated with regional habitat modifications directed by the California Bay-Delta Authority. On July 7, 2008, and August 5, 2008, pore-water profilers were successfully deployed at six North Bay sites per trip to measure the concentration gradient of dissolved macronutrients and trace metals near the sediment-water interface. Only two of the sites (433 and SSB009 within Honker Bay) were sampled in both series of profiler deployments. At each sampling site, profilers were deployed in triplicate, while discrete samples and dataloggers were used to collect ancillary data from both the water column and benthos to help interpret diffusive-flux measurements. Benthic flux of dissolved (0.2-micron filtered) inorganic phosphate (that is, soluble reactive phosphorus (SRP)) ranged from negligible levels (-0.003?0.005 millimole per square meter per day (mmole m-2d-1) at Site 4.1 outside Honker Bay) to 0.060?0.006 mmole m-2d-1 near the northern coast of Brown?s Island. Except for the elevated flux at Browns Island, the benthic flux of soluble reactive phosphorus (SRP) was consistently: (1) lower than previously reported for South Bay sites, (2) an order of magnitude lower than oligotrophic Coeur d?Alene Lake, (3) two orders of magnitude lower than determined for eutrophic Upper Klamath Lake, and (4) an order of magnitude or more lower than the estimated summer riverine inputs for SRP (900 to 1,300 kilograms of phosphorous per day (kg-P d-1)). In contrast to fluxes reported for the South Bay, nitrate fluxes were consistently negative (that is, drawn from the water column into the sediment), except for one site with statistically insignificant nitrate fluxes (Site 409 within Suisun Bay). The most negative nitrate flux (-7.3?0.1 mmole m-2d-1) was observed within Grizzly Bay (Site 416). Observed nitrate fluxes bracketed the estimated summer fluvial flux of nitrate (3,500 to 5,000 kg-N d-1). With the exception of the two Grizzly Bay sites (416 and 417), the consistently positive benthic flux of ammonia generally counteracted the negative flux of nitrate to yield a net balance of dissolved inorganic nitrogen. Ammonia benthic fluxes extrapolated for Suisun Bay ranged from 320 kg-N d-1 (Site SSB009 near the entrance to Honker Bay) to 1,900 kg-N d-1 (Montezuma Island). These values represent a significant ammonia source to the water column relative to summer riverine inputs (approximately 400 to 600 kg-N d-1). Dissolved silica also displayed a consistently positive benthic flux, except for Site 409 within Suisun Bay, which showed insignificant fluxes (also insignificant for nitrate and SRP). As with the nitrate fluxes, Grizzly Bay and Browns Island sites yielded the highest dissolved silica fluxes (1.3?1.2 to 2.5?0.6 mmole m-2d-1, respectively). These initial diffusive-flux estimates are greater than those measured in the South Bay using core-incubation experiments, which include bioturbation and bioirrigation effects, but they are nevertheless probably one to t

Open-File Report↗

USGS2010

Simulation of Runoff and Reservoir Inflow for Use in a Flood-Analysis Model for the Delaware River, Pennsylvania, New Jersey, and New York, 2004-2006

A model was developed to simulate inflow to reservoirs and watershed runoff to streams during three high-flow events between September 2004 and June 2006 for the main-stem subbasin of the Delaware River draining to Trenton, N.J. The model software is a modified version of the U.S. Geological Survey (USGS) Precipitation-Runoff Modeling System (PRMS), a modular, physically based, distributed-parameter modeling system developed to evaluate the impacts of various combinations of precipitation, climate, and land use on surface-water runoff and general basin hydrology. The PRMS model simulates time periods associated with main-stem flooding that occurred in September 2004, April 2005, and June 2006 and uses both daily and hourly time steps. Output from the PRMS model was formatted for use as inflows to a separately documented reservoir and riverrouting model, the HEC-ResSim model, developed by the U.S. Army Corps of Engineers Hydrologic Engineering Center to evaluate flooding. The models were integrated through a graphical user interface. The study area is the 6,780 square-mile watershed of the Delaware River in the states of Pennsylvania, New Jersey, and New York that drains to Trenton, N.J. A geospatial database was created for use with a geographic information system to assist model discretization, determine land-surface characterization, and estimate model parameters. The USGS National Elevation Dataset at 100-meter resolution, a Digital Elevation Model (DEM), was used for model discretization into streams and hydrologic response units. In addition, geospatial processing was used to estimate initial model parameters from the DEM and other data layers, including land use. The model discretization represents the study area using 869 hydrologic response units and 452 stream segments. The model climate data for point stations were obtained from multiple sources. These sources included daily data for 22 National Weather Service (NWS) Cooperative Climate Station network stations, hourly data for 15 stations from the National Climatic Data Center, hourly data for 1 station from the NWS Middle Atlantic River Forecast Center records, and daily and hourly data for 7 stations operated by the New York City Department of Environmental Protection. The NWS Multisensor Precipitation Estimate data set for 2001-2007 was used for computing daily precipitation for the model and for computing hourly precipitation for storm simulation periods. Calibration of the PRMS model included regression and optimization algorithms, as well as manual adjustments of model parameters. The general goal of the calibration procedure was to minimize the difference between discharge measured at USGS streamgages and the corresponding discharge simulated by the model. Daily streamflow data from 35 USGS streamgages were used in model calibration. The streamflow data represent areas draining from 20.2 to 6,780 square miles. The PRMS model simulates reservoir inflow and watershed runoff for use as input into HECResSim for the purpose of evaluating and comparing the effects of different watershed conditions on main-stem flooding in the Delaware River watershed draining to Trenton, N.J. The PRMS model is useful as a planning tool to simulate the effects of land-use changes and different antecedent conditions on local runoff and reservoir inflow and, as input to the HEC-ResSim model, on flood flows in the main stem of the Delaware River.

Open-File Report↗

USGS2010

Documentation for initial seismic hazard maps for Haiti

In response to the urgent need for earthquake-hazard information after the tragic disaster caused by the moment magnitude (M) 7.0 January 12, 2010, earthquake, we have constructed initial probabilistic seismic hazard maps for Haiti. These maps are based on the current information we have on fault slip rates and historical and instrumental seismicity. These initial maps will be revised and improved as more data become available. In the short term, more extensive logic trees will be developed to better capture the uncertainty in key parameters. In the longer term, we will incorporate new information on fault parameters and previous large earthquakes obtained from geologic fieldwork. These seismic hazard maps are important for the management of the current crisis and the development of building codes and standards for the rebuilding effort. The boundary between the Caribbean and North American Plates in the Hispaniola region is a complex zone of deformation. The highly oblique ~20 mm/yr convergence between the two plates (DeMets and others, 2000) is partitioned between subduction zones off of the northern and southeastern coasts of Hispaniola and strike-slip faults that transect the northern and southern portions of the island. There are also thrust faults within the island that reflect the compressional component of motion caused by the geometry of the plate boundary. We follow the general methodology developed for the 1996 U.S. national seismic hazard maps and also as implemented in the 2002 and 2008 updates. This procedure consists of adding the seismic hazard calculated from crustal faults, subduction zones, and spatially smoothed seismicity for shallow earthquakes and Wadati-Benioff-zone earthquakes. Each one of these source classes will be described below. The lack of information on faults in Haiti requires many assumptions to be made. These assumptions will need to be revisited and reevaluated as more fieldwork and research are accomplished. We made two sets of maps using different assumptions about site conditions. One set of maps is for a firm-rock site condition (30-m averaged shear-wave velocity, Vs30, of 760 m/s). We also developed hazard maps that contain site amplification based on a grid of Vs30 values estimated from topographic slope. These maps take into account amplification from soils. We stress that these new maps are designed to quantify the hazard for Haiti; they do not consider all the sources of earthquake hazard that affect the Dominican Republic and therefore should not be considered as complete hazard maps for eastern Hispaniola. For example, we have not included hazard from earthquakes in the Mona Passage nor from large earthquakes on the subduction zone interface north of Puerto Rico. Furthermore, they do not capture all the earthquake hazards for eastern Cuba.

Open-File Report↗

USGS2010

Seismicity of the Earth 1900-2010, Aleutian arc and vicinity

This map shows details of the Aleutian arc not visible in an earlier publication. The Aleutian arc extends about 3,000 km from the Gulf of Alaska to the Kamchatka Peninsula. It marks the region where the Pacific plate subducts into the mantle beneath the North America plate. This subduction is responsible for the generation of the Aleutian Islands and the deep offshore Aleutian Trench. Relative to a fixed North America plate, the Pacific plate is moving northwest at a rate that increases from about 55 mm per year at the arc's eastern edge to 75 mm per year near its western terminus. In the east, the convergence of the plates is nearly perpendicular to the plate boundary. However, because of the boundary's curvature, as one travels westward along the arc, the subduction becomes more and more oblique to the boundary until the relative plate motion becomes parallel to the arc at the Near Islands near its western edge. Subduction zones such as the Aleutian arc are geologically complex and produce numerous earthquakes from multiple sources. Deformation of the overriding North America plate generates shallow crustal earthquakes, whereas slip at the interface of the plates generates interplate earthquakes that extend from near the base of the trench to depths of 40 to 60 km. At greater depths, Aleutian arc earthquakes occur within the subducting Pacific plate and can reach depths of 300 km. Since 1900, six great earthquakes have occurred along the Aleutian Trench, Alaska Peninsula, and Gulf of Alaska: M8.4 1906 Rat Islands; M8.6 1938 Shumagin Islands; M8.6 1946 Unimak Island; M8.6 1957 Andreanof Islands; M9.2 1964 Prince William Sound; and M8.7 1965 Rat Islands. Several relevant tectonic elements (plate boundaries and active volcanoes) provide a context for the seismicity presented on the main map panel. The plate boundaries are most accurate along the axis of the Aleutian Trench and more diffuse or speculative in extreme northeastern Russia. The active volcanoes parallel the Aleutian Trench from the Gulf of Alaska to the Rat Islands.

Alaska↗

USGS2013

Seismicity of the Earth 1900-2012 Philippine Sea plate and vicinity

The complex tectonics surrounding the Philippine Islands are dominated by the interactions of the Pacific, Sunda, and Eurasia plates with the Philippine Sea plate (PSP). The latter is unique because it is almost exclusively surrounded by zones of plate convergence. At its eastern and southeastern edges, the Pacific plate is subducted beneath the PSP at the Izu-Bonin, Mariana, and Yap trenches. Here, the subduction zone exhibits high rates of seismic activity to depths of over 600 km, though no great earthquakes (M>8.0) have been observed, likely because of weak coupling along the plate interface. In the northeast, the PSP subducts beneath Japan and the eastern margin of the Eurasia plate at the Nankai and Ryukyu trenches, extending westward to Taiwan. The Nankai portion of this subduction zone has hosted some of the largest earthquakes along the margins of the PSP, including a pair of Mw8.1 megathrust events in 1944 and 1946. Along its western margin, the convergence of the PSP and the Sunda plate is responsible for a broad and active plate boundary system extending along both sides of the Philippine Islands chain. The region is characterized by opposite-facing subduction systems on the east and west sides of the islands, and the archipelago is cut by a major transform structure: the Philippine Fault. Subduction of the Philippine Sea plate occurs at the eastern margin of the islands along the Philippine Trench and its northern extension, the East Luzon Trough. On the west side of Luzon, the Sunda Plate subducts eastward along a series of trenches, including the Manila Trench in the north, the smaller Negros Trench in the central Philippines, and the Sulu and Cotabato trenches in the south. Twentieth and early twentyfirst century seismic activity along the boundaries of the Philippine Sea plate has produced seven great (M>8.0) earthquakes and 250 large (M>7) events. Among the most destructive events were the 1923 Kanto, the 1948 Fukui, and the 1995 Kobe, Japan, earthquakes; the 1935 and the 1999 Chi-Chi, Taiwan, earthquakes; and the 1976 M7.6 Moro Gulf and 1990 M7.6 Luzon, Philippines, earthquakes.

Open-File Report↗

USGS2014

Seismicity of the Earth 1900-2012 Java and vicinity

The Sunda convergent margin extends for 5,600 km from the Bay of Bengal and the Andaman Sea, both located northwest of the map area, towards the island of Sumba in the southeast, and then continues eastward as the Banda arc system. This tectonically active margin is a result of the India and Australia plates converging with and subducting beneath the Sunda plate at a rate of approximately 50 to 70 mm/yr. The main physiographic feature associated with this convergent margin is the Sunda-Java Trench, which stretches for 3,000 km parallel to the Java and Sumatra land masses and terminates at 120° E. The convergence of the Indo-Australia and Sunda plates produces two active volcanic arcs: Sunda, which extends from 105 to 122° E and Banda, which extends from 122 to 128° E. The Sunda arc results solely from relatively simple oceanic plate subduction, while the Banda arc represents the transition from oceanic subduction to continental collision, where a complex, broad deforming zone is found. Based on modern activity, the Banda arc can be divided into three distinct zones: an inactive section, the Wetar Zone, bound by two active segments, the Flores Zone in the west and the Damar Zone in the east. The lack of volcanism in the Wetar Zone is attributed to the collision of Australia with the Sunda plate. The absence of gap in volcanic activity is underlain by a gap in intermediate depth seismicity, which is in contrast to nearly continuous, deep seismicity below all three sections of the arc. The Flores Zone is characterized by down-dip compression in the subducted slab at intermediate depths and late Quaternary uplift of the forearc. These unusual features, along with GPS data interpretations indicate that the Flores Zone marks the transition between subduction of oceanic crust in the west and the collision of continental crust in the east. The Java section of the Sunda arc is considered relatively aseismic historically when compared to the highly seismically active Sumatra section, despite both areas being located along the same active subduction margin. Shallow (0–20 km) events have occurred historically in the overlying Sunda plate, causing damage to local and regional communities. A recent example was the May 26, 2006 M6.3 left-lateral strike-slip event that occurred at a depth of 10 km in central Java, and caused over 5,700 fatalities. Intermediate depth (70–300 km) earthquakes frequently occur beneath Java as a result of intraplate faulting within the Australia slab. Deep (300–650 km) earthquakes occur beneath the Java Sea and the back-arc region to the north of Java. Similar to other intermediate depth events, these earthquakes are also associated with intraslab faulting. However, this subduction zone exhibits a gap in seismicity from 250 to 400 km, interpreted as the transition between extensional and compressional slab stresses. Historical examples of large intraplate events include: the 1903 M8.1 event, 1921 M7.5 event, 1977 M8.3 event, and August 2007 M7.5 event. Large thrust earthquakes close to the Java trench are typically interplate faulting events along the slab interface between the Australia and Sunda plates. These earthquakes also generally have high tsunamigenic potential due to their shallow hypocentral depths. In some cases, these events have demonstrated slow moment-release and have been defined as ‘tsunami’ earthquakes, where rupture is large in the weak crustal layers very close to the seafloor. These events are categorized by tsunamis that are significantly larger than predicted by the earthquake’s magnitude. The most notable tsunami earthquakes in the Java region occurred on June 2, 1994 (M7.8) and July 17, 2006 (M7.7). The 1994 event produced a tsunami with wave runup heights of 13 m, killing over 200 people. The 2006 event produced a tsunami of up to 15 m, and killed 730 people. Although both of these tsunami earthquakes were characterized by rupture along thrust faults, they were followed by an abundance of normal faulting aftershocks. These aftershocks are interpreted to result from extension within the subducting Australia plate, whereas the mainshocks represented interplate faulting between the Australia and Sunda plates.

Open-File Report↗

USGS2010

West Florida Shelf: A natural laboratory for the study of ocean acidificiation

Declining oceanic pH and carbonate-ion concentrations are well-known consequences of increased atmospheric and surface-ocean partial pressure of carbon dioxide (pCO2). The possible subject of shifts in seawater carbonate chemistry on biocalcification and survival rates of marine organisms provides questions amenable to both experimental and field study (Kleypas and Langdon, 2006). To date, limited quantitative data exist with which to formalize and test hypotheses regarding such impacts, particularly in continental-shelf settings. The continental shelves of Florida provide an ideal natural laboratory in which to test latitudinal (and temperature and depth) shifts in habitat ranges of calcifying organisms. Both the east and west Florida shelves extend from warm temperate to subtropical latitudes; additionally, the west Florida shelf has very little siliciclastic influx to mask the carbonate production. This study utilizes the natural laboratory of the west and southwest Florida shelf (fig 1.1) to examine the transition from foramol (predominately foraminifera and molluscan) carbonate sediments, characteristic of the west-central Florida shelf, to chlorozoan (algal and coral) sediments characteristic of the southwest Florida shelf. The west Florida shelf is a mixed siliciclastic carbonate ramp that to the south transitions to the carbonate-dominated southwest Florida shelf (Enos, 1977; Brooks and others, 2003). The west Florida shelf is a distally steepened carbonate ramp that is ~250 kilometers (km) wide (Read, 1985). It is covered by a veneer of unconsolidated sediment consisting of mainly biogenic carbonate and quartz in the near shore, with subordinate amounts of phosphate. The sediment-distribution pattern is largely a function of proximity to source, with physical processes playing a minor role in distribution. The carbonate sand-and-gravel fraction is produced by organisms within the depositional basin of the west Florida shelf (Brooks and others, 2003). The southwest Florida shelf is a rimmed carbonate margin where organisms produce virtually all of the substrate; it also exhibits a greater sediment thickness as compared to the west Florida shelf (Enos, 1977). Temperature, which is usually associated with latitude, plays a major role in locations of foramol versus chlorozoan assemblages, but other factors beyond latitude influence temperature on the west and southwest Florida shelves. The potential of cooler, deep-water upwelling and transport over the bottom waters of the shelf may have a significant role in the species assemblage at the sediment/water interface and ultimately on location of foramol versus chlorozoan production. Deep water transported onto and over the shelf may also have environmental ramifications beyond temperature by bringing in water of different chemistry.

Open-File Report↗

USGS2010

CO2calc: A User-Friendly Seawater Carbon Calculator for Windows, Mac OS X, and iOS (iPhone)

A user-friendly, stand-alone application for the calculation of carbonate system parameters was developed by the U.S. Geological Survey Florida Shelf Ecosystems Response to Climate Change Project in response to its Ocean Acidification Task. The application, by Mark Hansen and Lisa Robbins, USGS St. Petersburg, FL, Joanie Kleypas, NCAR, Boulder, CO, and Stephan Meylan, Jacobs Technology, St. Petersburg, FL, is intended as a follow-on to CO2SYS, originally developed by Lewis and Wallace (1998) and later modified for Microsoft Excel? by Denis Pierrot (Pierrot and others, 2006). Besides eliminating the need for using Microsoft Excel on the host system, CO2calc offers several improvements on CO2SYS, including: An improved graphical user interface for data entry and results Additional calculations of air-sea CO2 fluxes (for surface water calculations) The ability to tag data with sample name, comments, date, time, and latitude/longitude The ability to use the system time and date and latitude/ longitude (automatic retrieval of latitude and longitude available on iPhone? 3, 3GS, 4, and Windows? hosts with an attached National Marine Electronics Association (NMEA)-enabled GPS) The ability to process multiple files in a batch processing mode An option to save sample information, data input, and calculated results as a comma-separated value (CSV) file for use with Microsoft Excel, ArcGIS,? or other applications An option to export points with geographic coordinates as a KMZ file for viewing and editing in Google EarthTM

Open-File Report↗

USGS2011

PQLX: A seismic data quality control system description, applications, and users manual

We present a detailed description and users manual for a new tool to evaluate seismic station performance and characteristics by providing quick and easy transitions between visualizations of the frequency and time domains. The software is based on the probability density functions (PDF) of power spectral densities (PSD) (McNamara and Buland, 2004) and builds on the original development of the PDF stand-alone software system (McNamara and Boaz, 2005) and the seismological data viewer application PQL (IRIS-PASSCAL Quick Look) and PQLII (available through the IRIS PASSCAL program: http://www.passcal.nmt.edu/content/pql-ii-program-viewing-data). With PQLX (PQL eXtended), computed PSDs are stored in a MySQL database, allowing a user to access specific time periods of PSDs (PDF subsets) and time series segments through a GUI-driven interface. The power of the method and software lies in the fact that there is no need to screen the data for system transients, earthquakes, or general data artifacts, because they map into a background probability level. In fact, examination of artifacts related to station operation and episodic cultural noise allow us to estimate both the overall station quality and a baseline level of Earth noise at each site. The output of this analysis tool is useful for both operational and scientific applications. Operationally, it is useful for characterizing the current and past performance of existing broadband stations, for conducting tests on potential new seismic station locations, for evaluating station baseline noise levels (McNamara and others, 2009), for detecting problems with the recording system or sensors, and for evaluating the overall quality of data and metadata. Scientifically, the tool allows for mining of PSDs for investigations on the evolution of seismic noise (for example, Aster and others, 2008; and Aster and others, 2010) and other phenomena. Currently, PQLX is operational at several organizations including the USGS National Earthquake Information Center (NEIC), the USGS Albuquerque Seismological Laboratory (ASL), and the Incorporated Research Institutions in Seismology (IRIS) Data Management Center (DMC) for station monitoring and instrument response quality control. The PQLX system is available to the community at large through the U.S. Geological Survey (USGS) (http://ehpm-earthquake.wr.usgs.gov/research/software/pqlx.php) and IRIS (http://www.iris.edu/software/pqlx). Also provided is a fully searchable website for bug reporting and enhancement requests (http://wush.net/bugzilla/PQLX). The first part of this document aims to describe and illustrate some of the features and capabilities of the software. The second part of this document is a detailed users manual that covers installation procedures, system requirements, operations, bug reporting, and software components (Appendix).

Open-File Report↗

USGS2010

Geochemical data for core and bottom-sediment samples collected in 2007 from Grand Lake O' the Cherokees, northeast Oklahoma

Grand Lake O' the Cherokees is a large reservoir in northeast Oklahoma, below the confluence of the Neosho and Spring Rivers, both of which drain the Tri-State Mining District to the north. The Tri-State district covers an area of 1,200 mi 2 (3,100 km 2 ) and comprises Mississippi Valley-type lead-zinc deposits. A result of 120 years of mining activity is an estimated 75 million tons of processed mine tailings (chat) remaining in the district. Concerns of sediment quality and the possibility of human exposure to cadmium and lead through eating fish have led to several studies of the sediments in the Tri-State district. In order to record the transport and deposition of metals from the Tri-State district by the Spring and Neosho Rivers into Grand Lake O' the Cherokees, the U.S. Geological Survey collected 11 sediment cores and 15 bottom-sediment samples in September 2007. Subsamples from five selected cores and the bottom-sediment samples were analyzed for major and trace elements and forms of carbon. The sediment samples collected from the sediment-water interface had larger average concentrations of zinc, cadmium, and lead than local background. The core collected from the Spring River had the largest concentrations of mining-related elements. A core collected just south of Twin Bridges State Park, at the confluence of the Spring and Neosho Rivers, showed a mixing zone with more mining-related elements coming from the Spring River side. The element zinc showed the most definitive patterns in graphs depicting concentration-versus-depth profiles. A core collected from the main body of the reservoir showed affected sediment down to a depth of 85 cm (33 in). This core and two others appear to have penetrated to below mining-affected sediment.

Open-File Report↗

USGS2013

Continuous resistivity profiling data from Great South Bay, Long Island, New York

An investigation of submarine aquifers adjacent to the Fire Island National Seashore and Long Island, New York was conducted to assess the importance of submarine groundwater discharge as a potential nonpoint source of nitrogen delivery to Great South Bay. Over 200 kilometers of continuous resistivity profiling data were collected to image the fresh-saline groundwater interface in sediments beneath the bay. In addition, groundwater sampling was performed at sites (1) along the north shore of Great South Bay, particularly in Patchogue Bay, that were representative of the developed Long Island shoreline, and (2) at sites on and adjacent to Fire Island, a 50-kilometer-long barrier island on the south side of Great South Bay. Other field activities included sediment coring, stationary electrical resistivity profiling, and surveys of in situ pore water conductivity. Results of continuous resistivity profiling surveys are described in this report. The onshore and offshore shallow hydrostratigraphy of the Great South Bay shorelines, particularly the presence and nature of submarine confining units, appears to exert primary control on the dimensions and chemistry of the submarine groundwater flow and discharge zones. Sediment coring has shown that the confining units commonly consist of drowned and buried peat layers likely deposited in salt marshes. Low-salinity groundwater extends from 10 to 100 meters offshore along much of the north and south shores of Great South Bay based on continuous resistivity profiling data, especially off the mouths of tidal creeks and beneath shallow flats to the north of Fire Island adjacent to modern salt marshes. Human modifications of much of the shoreline and nearshore areas along the north shore of the bay, including filling of salt marshes, construction of bulkheads and piers, and dredging of navigation channels, has substantially altered the natural hydrogeology of the bay's shorelines by truncating confining units and increasing recharge near the shore in filled areas. Better understanding of the nature of submarine groundwater discharge along developed and undeveloped shorelines of embayments such as this could lead to improved models and mitigation strategies for nutrient overenrichment of estuaries.

New York↗

USGS2011

cloudPEST - A python module for cloud-computing deployment of PEST, a program for parameter estimation

This report documents cloudPEST-a Python module with functions to facilitate deployment of the model-independent parameter estimation code PEST on a cloud-computing environment. cloudPEST makes use of low-level, freely available command-line tools that interface with the Amazon Elastic Compute Cloud (EC2(TradeMark)) that are unlikely to change dramatically. This report describes the preliminary setup for both Python and EC2 tools and subsequently describes the functions themselves. The code and guidelines have been tested primarily on the Windows(Registered) operating system but are extensible to Linux(Registered).

Open-File Report↗

USGS2011

Description and testing of the Geo Data Portal: Data integration framework and Web processing services for environmental science collaboration

Interest in sharing interdisciplinary environmental modeling results and related data is increasing among scientists. The U.S. Geological Survey Geo Data Portal project enables data sharing by assembling open-standard Web services into an integrated data retrieval and analysis Web application design methodology that streamlines time-consuming and resource-intensive data management tasks. Data-serving Web services allow Web-based processing services to access Internet-available data sources. The Web processing services developed for the project create commonly needed derivatives of data in numerous formats. Coordinate reference system manipulation and spatial statistics calculation components implemented for the Web processing services were confirmed using ArcGIS 9.3.1, a geographic information science software package. Outcomes of the Geo Data Portal project support the rapid development of user interfaces for accessing and manipulating environmental data.

Open-File Report↗

USGS2011

Simulating daily water temperatures of the Klamath River under dam removal and climate change scenarios

A one-dimensional daily averaged water temperature model was used to simulate Klamath River temperatures for two management alternatives under historical climate conditions and six future climate scenarios. The analysis was conducted for the Secretarial Determination on removal of four hydroelectric dams on the Klamath River. In 2012, the Secretary of the Interior will determine if dam removal and implementation of the Klamath Basin Restoration Agreement (KBRA) (Klamath Basin Restoration Agreement, 2010) will advance restoration of salmonid fisheries and is in the public interest. If the Secretary decides dam removal is appropriate, then the four dams are scheduled for removal in 2020. Water temperature simulations were conducted to compare the effect of two management alternatives: the no-action alternative where dams remain in place, and the action alternative where dam removal occurs in 2020 along with habitat restoration. Each management alternative was simulated under historical climate conditions (1961-2010) and six 50-year (2012-2061) climate scenarios. The model selected for the study, River Basin Model-10 (RBM10), was used to simulate water temperatures over a 253-mile reach of the Klamath River located in south-central Oregon and northern California. RBM10 uses a simple equilibrium flow model, assuming discharge in each river segment on each day is transmitted downstream instantaneously. The model uses a heat budget formulation to quantify heat flux at the air-water interface. Inputs for the heat budget were calculated from daily-mean meteorological data, including net shortwave solar radiation, net longwave atmospheric radiation, air temperature, wind speed, vapor pressure, and a psychrometric constant needed to calculate the Bowen ratio. The modeling domain was divided into nine reaches ranging in length from 10.8 to 42.4 miles, which were calibrated and validated separately with measured water temperature data collected irregularly from 1961 to 2010. Calibration root mean square errors of observed versus simulated water temperatures for the nine reaches ranged from 0.8 to 1.5 degrees C. Mean absolute errors ranged from 0.6 to 1.2 degrees C. For model validation, a k-fold cross-validation technique was used. Validation root mean square error and mean absolute error for the nine reaches ranged from 0.8 to 1.4 degrees C and 0.8 to 1.2 degrees C, respectively. Input data for the six future climate scenarios (2012-2061) were derived from historical hydrological and meteorological data and simulated meteorological output from five Global Circulation Models. Total Maximum Daily Loads or other regulatory processes that might reduce future water temperatures were not included in the simulations. Under the current climate conditions scenario, impacts of dam removal on water temperatures were greatest near Iron Gate Dam (near Yreka, California) and were attenuated in the lower reaches of the Klamath River. May and October simulated mean water temperatures increased and decreased by approximately 1-2 degrees C and 2-4 degrees C, respectively, downstream of Iron Gate Dam after dam removal. Dam removal also resulted in an earlier annual temperature cycle shift of 18 days, 5 days, and 2 days, near Iron Gate Dam, Scott River, and Trinity River, respectively. Although the magnitude of precipitation and air temperature change predicted by the five Global Circulation Models varied, all five models resulted in progressive incremental increases in water temperatures with each decade from 2012 to 2061. However, dam removal under KBRA appeared to delay the effects of climate change to some extent near Iron Gate Dam. With dam removal under KBRA, annual-mean water temperatures exceeded the 49-year historical mean temperature beginning in 2045; whereas with dams, annual-mean temperatures exceeded the historical mean beginning in 2025. Potential changes in seasonal water temperatures resulting from dam removal, with or without future climat

Open-File Report↗

USGS2011

The Lake Tahoe Basin Land Use Simulation Model

This U.S. Geological Survey Open-File Report describes the final modeling product for the Tahoe Decision Support System project for the Lake Tahoe Basin funded by the Southern Nevada Public Land Management Act and the U.S. Geological Survey's Geographic Analysis and Monitoring Program. This research was conducted by the U.S. Geological Survey Western Geographic Science Center. The purpose of this report is to describe the basic elements of the novel Lake Tahoe Basin Land Use Simulation Model, publish samples of the data inputs, basic outputs of the model, and the details of the Python code. The results of this report include a basic description of the Land Use Simulation Model, descriptions and summary statistics of model inputs, two figures showing the graphical user interface from the web-based tool, samples of the two input files, seven tables of basic output results from the web-based tool and descriptions of their parameters, and the fully functional Python code.

Open-File Report↗

USGS2012

User's guide for mapIMG 3--Map image re-projection software package

Version 0.0 (1995), Dan Steinwand, U.S. Geological Survey (USGS)/Earth Resources Observation Systems (EROS) Data Center (EDC)--Version 0.0 was a command line version for UNIX that required four arguments: the input metadata, the output metadata, the input data file, and the output destination path. Version 1.0 (2003), Stephen Posch and Michael P. Finn, USGS/Mid-Continent Mapping Center (MCMC--Version 1.0 added a GUI interface that was built using the Qt library for cross platform development. Version 1.01 (2004), Jason Trent and Michael P. Finn, USGS/MCMC--Version 1.01 suggested bounds for the parameters of each projection. Support was added for larger input files, storage of the last used input and output folders, and for TIFF/ GeoTIFF input images. Version 2.0 (2005), Robert Buehler, Jason Trent, and Michael P. Finn, USGS/National Geospatial Technical Operations Center (NGTOC)--Version 2.0 added Resampling Methods (Mean, Mode, Min, Max, and Sum), updated the GUI design, and added the viewer/pre-viewer. The metadata style was changed to XML and was switched to a new naming convention. Version 3.0 (2009), David Mattli and Michael P. Finn, USGS/Center of Excellence for Geospatial Information Science (CEGIS)--Version 3.0 brings optimized resampling methods, an updated GUI, support for less than global datasets, UTM support and the whole codebase was ported to Qt4.

Open-File Report↗

USGS2011

Summary of November 2010 meeting to evaluate turbidite data for constraining the recurrence parameters of great Cascadia earthquakes for the update of national seismic hazard maps

This report summarizes a meeting of geologists, marine sedimentologists, geophysicists, and seismologists that was held on November 18–19, 2010 at Oregon State University in Corvallis, Oregon. The overall goal of the meeting was to evaluate observations of turbidite deposits to provide constraints on the recurrence time and rupture extent of great Cascadia subduction zone (CSZ) earthquakes for the next update of the U.S. national seismic hazard maps (NSHM). The meeting was convened at Oregon State University because this is the major center for collecting and evaluating turbidite evidence of great Cascadia earthquakes by Chris Goldfinger and his colleagues. We especially wanted the participants to see some of the numerous deep sea cores this group has collected that contain the turbidite deposits. Great earthquakes on the CSZ pose a major tsunami, ground-shaking, and ground-failure hazard to the Pacific Northwest. Figure 1 shows a map of the Pacific Northwest with a model for the rupture zone of a moment magnitude M w 9.0 earthquake on the CSZ and the ground shaking intensity (in ShakeMap format) expected from such an earthquake, based on empirical ground-motion prediction equations. The damaging effects of such an earthquake would occur over a wide swath of the Pacific Northwest and an accompanying tsunami would likely cause devastation along the Pacifc Northwest coast and possibly cause damage and loss of life in other areas of the Pacific. A magnitude 8 earthquake on the CSZ would cause damaging ground shaking and ground failure over a substantial area and could also generate a destructive tsunami. The recent tragic occurrence of the 2011 M w 9.0 Tohoku-Oki, Japan, earthquake highlights the importance of having accurate estimates of the recurrence times and magnitudes of great earthquakes on subduction zones. For the U.S. national seismic hazard maps, estimating the hazard from the Cascadia subduction zone has been based on coastal paleoseismic evidence of great earthquakes over the past 5,000 years. The instrumental catalog of earthquakes is of little use for constraining the hazard of the CSZ, because there are virtually no recorded earthquakes on most of the plate interface of the CSZ. There are no historical accounts in the past 150 years of large earthquakes on most of the CSZ. Until about 20 years ago, some interpreted this lack of recent and historical earthquakes as an indicator that the subduction zone was slipping aseismically and could not produce a great earthquake. The work of Brian Atwater and others, in the late 1980s and the 1990s (Atwater, 1987, 1992; Atwater and others, 1995; Nelson and others, 1996; Clague, 1997; Atwater and Hemphill-Haley, 1997; Atwater and others, 2004) demonstrated that submerged forests, buried soils, tsunami deposits, and liquefaction along and near the coast were compelling evidence of repeated great Cascadia earthquakes over at least the past 5,000 years. Atwater and Hemphill-Haley (1997) concluded from paleoseismic evidence at Willapa Bay, Washington, that great earthquakes ruptured the CSZ with an average recurrence time of about 500 years. The date of the last great CSZ earthquake, January 26, 1700, was established from historical records of the so-called orphan tsunami in Japan that is inferred to have been produced by this earthquake (Satake and others, 1996, 2003; Atwater and others, 2005) and is consistent with tree-ring data from drowned forests in Washington and Oregon. From modeling the observations of the tsunami, Satake and others (2003) estimated a moment magnitude of about 9.0 for this earthquake. Many other paleoseismic sites have been investigated along the Pacific Northwest coast from Vancouver Island to northern California and show evidence of great CSZ earthquakes. Nelson and others (2006) summarized the dates found from these studies and proposed correlations between sites indicating the extent of rupture for individual events. Dating of inferred tsunami deposits in Bradley Lake, Oregon by Kelsey and others (2005), as well as tsunami and subsidence evidence from Six Rivers, Oregon (Kelsey and others, 2002) and Coquille River (Witter and others, 2003), indicates that there were probably M w 8 ruptures in the southern portion of the CSZ in addition to the M w 9 events that rupture the whole length of the CSZ (Nelson and others, 2006). A parallel development over the past 20 years or more is the use of deep-sea turbidite deposits for identifying and dating great Cascadia earthquakes over the past 10,000 years (Adams, 1990; Goldfinger and others, 2003, 2008, in press; Goldfinger, 2011). Turbidites are sediment deposits in the deep ocean from turbidity currents, which are energetic flows of sediment and water along the continental shelf and slope. Adams (1990), using the counts of turbidites in deep-sea cores off the coast of Oregon and Washington collected and analyzed by Griggs (1969) and Griggs and others (1969), proposed that these turbidites were caused by the shaking of great Cascadia earthquakes. Part of his reasoning was that the number (13) of turbidite deposits that occurred since deposition of the Mazama Ash 7,000 years ago gave a recurrence time of about 500 years, consistent with that derived from the coastal submergence data. Adams (1990) also proposed the “confluence test” which evaluates the number of turbidites for submarine channels that form a confluence. He reported that the number of turbidites in the single downstream channel equaled the number in each of the tributary channels. He reasoned that this indicated that the turbidites in each tributary were simultaneously triggered and were, therefore, caused by a common forcing agent. He concluded that shaking from extended ruptures of great Cascadia earthquakes was the most likely cause of these turbidites. Based on the paleoseismic evidence of past great earthquakes, the hazard from the Cascadia subduction zone was included in the 1996 U.S. NSHM (Frankel and others, 1996), which were the basis for seismic provisions in the 2000 International Building Code. These hazard maps used the paleoseismic studies to constrain the recurrence rate of great CSZ earthquakes. Goldfinger and his colleagues have since collected many more deep ocean cores and done extensive analysis on the turbidite deposits that they identified in the cores (Goldfinger and others, 2003, 2008, in press; Goldfinger, 2011). Using their dating of the sediments and correlation of features in the logs of density and magnetic susceptibility between cores, they developed a detailed chronology of great earthquakes along the CSZ for the past 10,000 years (Goldfinger and others, in press). These correlations consist of attempting to match the peaks and valleys in logs of density and magnetic susceptibility between cores separated, in some cases, by hundreds of kilometers. Based on this work, Goldfinger and others (2003, 2008, in press) proposed that the turbidite evidence indicated the occurrence of great earthquakes (M w 8) that only ruptured the southern portion of the CSZ, as well as earthquakes with about M w 9 that ruptured the entire length of the CSZ. For the southernmost portion of the CSZ, Goldfinger and others (in press) proposed a recurrence time of M w 8 or larger earthquakes of about 230 years. This proposed recurrence time was shorter than the 500 year time that was incorporated in one scenario in the NSHM’s. It is important to note that the hazard maps of 1996 and later also included a scenario or set of scenarios with a shorter recurrence time for M w 8 earthquakes, using rupture zones that are distributed along the length of the CSZ (Frankel and others, 1996; Petersen and others, 2008). Originally, this scenario was meant to correspond to the idea that some of the 500-year averaged ruptures seen in the paleoseismic evidence could have been a series of M w 8 earthquakes that occurred over a short period of time (a few decades), rather than M w 9 earthquakes. Figure 2 shows the logic tree for the CSZ used in the 2008 NSHM’s (Petersen and others, 2008). This logic tree includes whole CSZ rupture earthquakes (M w 8.8–9.2) and partial CSZ rupture earthquakes (M w 8.0–8.7). In this latest version of the NSHM’s, the effective recurrence time of earthquakes on the CSZ with moment magnitudes greater than or equal to 8.0 over the various models is about 270 years (Petersen and others, 2008). This recurrence time applies to the entire CSZ, so that the hazard from great earthquakes was approximately equal along the whole zone, although the hazard estimates taper on the northern and southern ends of the CSZ, because of the way rupture zones of M w 8 earthquakes were distributed along the strike of the CSZ. The NSHM will be updated in 2013, as part of the standard update cycle that corresponds to the update cycle of the national model building codes that are based on the seismic hazard maps. A meeting was necessary to assemble a wide group of experts to hear Dr. Goldfinger explain his methodology for dating and correlating the turbidites and for developing the earthquake chronology. The overall goal of the workshop was to evaluate observations of turbidite deposits to provide constraints on the recurrence times and rupture extents of great Cascadia subduction zone earthquakes for the next update of the NSHM. Before the meeting, participants were supplied with the U.S. Geological Survey (USGS) Professional Paper of Goldfinger and others (in press), as well as material from Brian Atwater and Alan Nelson. The agenda of the meeting was developed by Art Frankel, with assistance from Chris Goldfinger, Brian Atwater, Alan Nelson, Mark Petersen, and Craig Weaver. The meeting was hosted by Chris Goldfinger of Oregon State University. We stress that it is difficult to evaluate in a two-day meeting the large amount of work that Goldfinger and his colleagues have done over the past 15 years or more. This meeting is the first step in a process that develops the inputs to the update of the national maps. The conclusions of this workshop will be discussed and possibly modified at the regional Pacific Northwest workshop for the hazard maps to be held in early 2012. Vetting new research results using informed expert opinion is an integral part of updating the national maps and does not reflect on the veracity of these results.

Cascadia↗

USGS2012

USGS Environmental health science strategy: providing environmental health science for a changing world: Public review release

America has an abundance of natural resources. We have bountiful clean water, fertile soil, and unrivaled national parks, wildlife refuges, and public lands. These resources enrich our lives and preserve our health and wellbeing. These resources have been maintained because of our history of respect for their value and an enduring commitment to their vigilant protection. Awareness of the social, economic, and personal value of the health of our environment is increasing. The emergence of environmentally driven diseases caused by environmental exposure to contaminants and pathogens is a growing concern worldwide. New health threats and patterns of established threats are affected by both natural and anthropogenic changes to the environment. Human activities are key drivers of emerging (new and re-emerging) health threats. Societal demands for land and natural resources, a better quality of life, improved economic prosperity, and the environmental impacts associated with these demands will continue to increase. Natural earth processes, climate trends, and related climatic events will add to the environmental impact of human activities. These environmental drivers will influence exposure to disease agents, including viral, bacterial, prion, and fungal pathogens, parasites, natural earth materials, toxins and other biogenic compounds, and synthetic chemicals and substances. The U.S. Geological Survey (USGS) defines environmental health science broadly as the interdisciplinary study of relations among the quality of the physical environment, the health of the living environment, and human health. The interactions among these three spheres are driven by human activities, ecological processes, and natural earth processes; the interactions affect exposure to contaminants and pathogens and the severity of environmentally driven diseases in animals and people. This definition provides USGS with a framework for synthesizing natural science information from across the Bureau and providing it to environmental, natural resource, agricultural, and public-health managers. The USGS is a Federal science agency with a broad range of natural science expertise relevant to environmental health. USGS provides scientific information and tools as a scientific basis for management and policy decision making. USGS specializes in science at the environment-health interface, by characterizing the processes that affect the interaction among the physical environment, the living environment, and people, and the resulting factors that affect ecological and human exposure to disease agents. This report describes a 10-year strategy that encompasses the portfolio of USGS environmental health science. It summarizes national environmental health priorities that USGS is best suited to address, and will serve as a strategic framework for USGS environmental health science goals, actions, and outcomes for the next decade. Implementation of this strategy is intended to aid coordination of USGS environmental health activities and to provide a focal point for disseminating information to stakeholders. The "One Health" paradigm advocated by the World Health Organization (WHO, 2011), and the American Veterinary Medicine Association (AVMA, 2008), among others, is based on a general recognition that the health of humans, animals, and the environment are inextricably linked. Thus, successful efforts to protect that health will require increased interdisciplinary research and increased communication and collaboration among the broader scientific and health community. This strategy is built upon that paradigm. The vision, mission, and five cornerstone goals of the USGS Environmental Health Science Strategy were developed with significant input from a wide range of stakeholders. Vision - The USGS is a premier source of the environmental health science needed to safeguard the health of the environment, fish, wildlife, and people. Mission - The mission of USGS in environmental health science is to contribute scientific information to environmental, natural resource, agricultural, and public-health managers, who use that science to support sound decision making. USGS provides the science to: - Goal 1: Identify, prioritize, and detect contaminants and pathogens of emerging environmental concern. - Goal 2: Reduce the impact of contaminants on the environment, fish, wildlife, and people. - Goal 3: Reduce the impact of pathogens on the environment, fish, wildlife, and people. - Goal 4: Discover the complex interactions and combined effects of exposure to contaminants and pathogens. - Goal 5: Prepare for and respond to environmental impacts and related health threats of natural and anthropogenic disasters. Goals 1 through 4 are intended to provide science to address environmental health threats in a logical order, from informing prevention and preparedness, to supporting systematic management response to environmental health issues. Goal 4 addresses the interaction among contaminants and pathogens, an issue of emerging concern in environmental health science. Goal 5 acknowledges the fact that natural and anthropogenic disasters can cause immediate and prolonged adverse environmental health threats. This strategy proposes that USGS take the following strategic science actions to achieve each of the five goals of this strategy: Goal 1: Identify, prioritize, and detect contaminants and pathogens of emerging environmental concern. - Strategic Science Action 1. - Prioritize contaminants and pathogens of emerging concern to guide research, detection, and management activities. - Strategic Science Action 2. - Conduct surveillance and monitoring to provide early warning of emerging health threats. - Strategic Science Action 3. - Develop approaches and tools that identify vulnerable environmental settings, ecosystems, and species. Goal 2: Reduce the impact of contaminants on the environment, fish, wildlife, and people. - Strategic Science Action 1. - Systematically characterize the sources, occurrence, transport and fate of environmental contaminants to guide efforts to manage and mitigate contamination. - Strategic Science Action 2. - Evaluate the threats of contamination on the health of the environment, fish, wildlife, and people, and inform the associated management and protection efforts. - Strategic Science Action 3. - Characterize potential human exposure to support establishment of health-based standards or guidelines and contamination-reduction efforts. Goal 3: Reduce the impact of pathogens on the environment, fish, wildlife, and people. - Strategic Science Action 1. - Determine the biotic and abiotic factors that control the ecology of infectious diseases affecting natural populations of aquatic and terrestrial species and potential transmission to other animals and humans. - Strategic Science Action 2. - Establish how natural and anthropogenic environmental changes affect the distribution and severity of infectious diseases in natural populations of aquatic and terrestrial species and potential transmission to other animals and humans. - Strategic Science Action 3. - Develop surveillance systems to identify changing patterns of disease activity in priority geographic areas. Goal 4: Discover the complex interactions and combined effects of exposure to contaminants and pathogens. - Strategic Science Action 1. - Identify how exposure to one class of disease agents (contaminants or pathogens) can make an organism more susceptible to effects from exposure to the other class of disease agents. - Strategic Science Action 2. - Implement interdisciplinary studies that characterize the effects of combined exposure to pathogens and contaminants. Goal 5: Prepare for and respond to the environmental impacts and related health threats of natural and anthropogenic disasters. - Strategic Science Action 1. - Establish a formal interdisciplinary science capability to rapidly assess the environmental health risks associated with disasters. - Strategic Science Action 2. - Enhance methods to anticipate, prepare for, and identify environmental, ecological, and related health impacts of future disasters. This strategy is one of seven USGS science strategies developed concurrently: - Climate and Land Use Change - Core Science Systems - Ecosystems - Energy and Mineral Resources - Environmental Health - Natural Hazards - Water. This strategy describes how USGS will address the highest priority environmental health issues facing the Nation. The ultimate intended outcome of this science strategy is prevention and reduction of adverse impacts to the quality of the environment, the health of our living resources, and human health. Communication with, and receiving input from, partners and stakeholders regarding their science needs is essential for successful implementation of this strategy. It is incumbent on USGS to reach out to all stakeholders to ensure that USGS efforts are focused on the highest priority environmental health issues and that products are provided in the most timely and usable form to all those who can use them. USGS must reach out to the scientific community, internally and externally, to ensure that our efforts are integrated with and take full advantage of the activities of others.

Open-File Report↗