Search USGS⌕ Search

SEARCH · Search USGS

Results for “Water Use Map Series”

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 577 records · Page 32Linked to original sources

Seasonal and Spatial Distribution of Freshwater Flow and Salinity in the Ten Thousand Islands Estuary, Florida, 2007-2009

The watershed of the Ten Thousand Islands (TTI) estuary has been substantially altered through the construction of canals and roads for the Southern Golden Gate Estates (SGGE), Barron River Canal, and U.S. 41 (Tamiami Trail). Two restoration projects designed to improve freshwater delivery to the estuary are the Picayune Strand Restoration Project, which includes the Southern Golden Gate Estates, and the Tamiami Trail Culverts Project; both are part of the Comprehensive Everglades Restoration Plan. To address hydrologic information needs critical for monitoring the effects of these restoration projects, the U.S. Geological Survey initiated a study in October 2006 to characterize freshwater outflows from the rivers, internal circulation and mixing within the estuary, and surface-water exchange between the estuary and Gulf of Mexico. The effort is conducted in cooperation with the South Florida Water Management District and complemented by monitoring performed by the Rookery Bay National Estuarine Research Reserve. Surface salinity was measured during moving boat surveys using a flow-through system that operated at planing speeds averaging 20 miles per hour. The data were logged every 10 seconds by a data recorder that simultaneously logged location information from a Global Positioning System. The major rivers, bays, and nearshore Gulf of Mexico region of the TTI area were surveyed in approximately 5 hours by two boats traversing about 200 total miles. Salinity and coordinate data were processed using inverse distance weighted interpolation to create salinity contour maps of the entire TTI region. Ten maps were created from salinity surveys performed between May 2007 and May 2009 and illustrate the dry season, transitional, and wet season salinity patterns of the estuarine rivers, inner bays, mangrove islands, and Gulf of Mexico boundary. The effects of anthropogenic activities are indicated by exceptionally low salinities associated with point discharge into the estuary from the Faka Union Canal and Barron River during the wet season. Low salinities in Faka Union Bay may cause reduced diversity and density of submerged aquatic vegetation, fish, and benthic organisms compared with neighboring Fakahatchee Bay. The Faka Union Canal System reduced the size of the watershed for the western TTI estuary, resulting in increased wet season salinities compared to those for the eastern TTI estuary, the watershed of which is composed of the relatively pristine Fakahatchee Strand Preserve State Park. Minimal river discharge and high evaporation caused hypersaline conditions to develop throughout the entire TTI region during the dry season. The 2007-2008 drought and passage of Tropical Storm Fay on August 18-19, 2008, demonstrated the effects of seasonal rainfall on salinity patterns, with substantially higher salinities observed during the 2007 wet season compared to those for the 2008 wet season. The salinity maps, coupled with data from the monitoring stations, provide baseline information of seasonal and spatial distribution of freshwater flow and salinity in the TTI estuary, and a means of monitoring the effects of restoration in improving freshwater delivery to the estuary.

Data Series↗

Void-Filled SRTM Digital Elevation Model of Afghanistan

EXPLANATION The purpose of this data set is to provide a single consistent elevation model to be used for national scale mapping, GIS, remote sensing applications, and natural resource assessments for Afghanistan's reconstruction. For 11 days in February of 2000, the National Aeronautics and Space Administration (NASA), the National Geospatial-Intelligence Agency ian Space Agency (ASI) flew X-band and C-band radar interferometry onboard the Space Shuttle Endeavor. The mission covered the Earth between 60?N and 57?S and will provide interferometric digital elevation models (DEMs) of approximately 80% of the Earth's land mass when processing is complete. The radar-pointing angle was approximately 55? at scene center. Ascending and descending orbital passes generated multiple interferometric data scenes for nearly all areas. Up to eight passes of data were merged to form the final processed Shuttle Radar Topography Mission (SRTM) DEMs. The effect of merging scenes averages elevation values recorded in coincident scenes and reduces, but does not completely eliminate, the amount of area with layover and terrain shadow effects. The most significant form of data processing for the Afghanistan DEM was gap-filling areas where the SRTM data contained a data void. These void areas are as a result of radar shadow, layover, standing water, and other effects of terrain as well as technical radar interferometry phase unwrapping issues. To fill these gaps, topographic contours were digitized from 1:200,000 - scale Soviet General Staff Topographic Maps which date from the middle to late 1980's. Digital contours were gridded to form elevation models for void areas and subsequently were merged with the SRTM data through GIS and image processing techniques. The data contained in this publication includes SRTM DEM quadrangles projected and clipped in geographic coordinates for the entire country. An index of all available SRTM DEM quadrangles is displayed here: Index_Geo_DD.pdf. Also included are quadrangles projected into their appropriate Universal Transverse Mercator (UTM) projection. The country of Afghanistan spans three UTM Zones: Zone 41, Zone 42, and Zone 43. Maps are stored in their respective UTM Zone projection. Indexes of all available SRTM DEM quadrangles in their respective UTM zone are displayed here: Index_UTM_Z41.pdf, Index_UTM_Z42.pdf, Index_UTM_Z43.pdf.

Data Series↗

Statistical summaries of ground-water level data collected in the Suwannee River Water Management District, 1948 to 1994

Since 1948, ground-water level data have beensystematically collected from selected wells in theSuwannee River Water Management District (SRWMD) by the U.S. Geological Survey (USGS),the SRWMD, and other agencies. Records of waterlevels in the SRWMD (fig. 1), collected by the USGS and SRWMD through 1990, and by the SRWMD from 1990 to 1994, have been published for many years in the USGS annual report series "Water Resources Data for Florida." However, no systematic statistical summaries of water levels in the SRWMD have been previously published. The need for such statistical summary data forevaluations of drought severity, ground-water supplyavailability, and minimum water levels for regulatory purposes increases daily as demands for ground-water usage increase. Also, much of the base flow of the Suwannee River is dependent upon ground water. As the population and demand for ground water for drinking water and irrigation purposes increase, the ability to quickly and easily predict trends in ground-water availability will become paramount. In response to this need, the USGS, in cooperation with the SRWMD, compiled this report. Ground-water sta tistics for 136 sites are presented as well as figures showing water levels that were measured in wells from 1948 through September 1994. In 1994, the SRWMD and the USGS began a long- term program of cooperative studies designed tobetter understand minimum and maximum streamflows and ground-water levels in the SRWMD. Minimum and maximum flows and levels are needed by the district to manage the surface- and ground-water resources of the SRWMD and to maintain or improve the various ecosystems. Data evaluation was a necessary first step in the long- term SRWMD ground-water investigations program, because basic statistics for ground-water levels are not included in the USGS annual data reports such as "Water Resources Data for Florida, Water Year 1994" (Fran klin and others, 1995). Statistics included in this report were generated using the USGS computer pro gram ADAPS (Automatic Data Processing System) to characterize normal ground-water levels and depar tures from normal. The report has been organized so that the statisti cal analyses of water levels in the wells are presentedfollowing this introductory material, a description ofthe hydrogeology in the study area, and a description of the statistics used to present the water-level data. Specifically, the report presents statistical analyses for each well, as appropriate, in the following manner: Description of the well.Hydrographs of ground-water levels for the period of record, for the last 10 years of record, and for the last 5 years of record. Graphs of maximum, minimum, and mean of monthly mean ground-water levels for wells with 5 or more years of record.Frequency hydrographs (25, 50, and 75 percent) of monthly mean ground-water levels for wells with 5 or more years of record. Water-level data and statistical plots are grouped by county and sorted within the county by ascendingsite identification number. Well locations are plottedon county maps preceding the well descriptions andhydrographs.

Open-File Report↗

Ground-water-quality data in Pennsylvania: A compilation of computerized [electronic] databases, 1979-2004

This study, by the U.S. Geological Survey (USGS) in cooperation with the Pennsylvania Department of Environmental Protection (PADEP), provides a compilation of ground-water-quality data for a 25-year period (January 1, 1979, through August 11, 2004) based on water samples from wells. The data are from eight source agencies唯orough of Carroll Valley, Chester County Health Department, Pennsylvania Department of Environmental Protection-Ambient and Fixed Station Network, Montgomery County Health Department, Pennsylvania Drinking Water Information System, Pennsylvania Department of Agriculture, Susquehanna River Basin Commission, and the U.S. Geological Survey. The ground-water-quality data from the different source agencies varied in type and number of analyses; however, the analyses are represented by 12 major analyte groups:biological (bacteria and viruses), fungicides, herbicides, insecticides, major ions, minor ions (including trace elements), nutrients (dominantly nitrate and nitrite as nitrogen), pesticides, radiochemicals (dominantly radon or radium), volatile organic compounds, wastewater compounds, and water characteristics (dominantly field pH, field specific conductance, and hardness). A summary map shows the areal distribution of wells with ground-water-quality data statewide and by major watersheds and source agency. Maps of 35 watersheds within Pennsylvania are used to display the areal distribution of water-quality information. Additional maps emphasize the areal distribution with respect to 13 major geolithologic units in Pennsylvania and concentration ranges of nitrate (as nitrogen). Summary data tables by source agency provide information on the number of wells and samples collected for each of the 35 watersheds and analyte groups. The number of wells sampled for ground-water-quality data varies considerably across Pennsylvania. Of the 8,012 wells sampled, the greatest concentration of wells are in the southeast (Berks, Bucks, Chester, Delaware, Lancaster, Montgomery, and Philadelphia Counties), in the vicinity of Pittsburgh, and in the northwest (Erie County). The number of wells sampled is relatively sparse in south-central (Adams, Cambria, Cumberland, and Franklin Counties), central (Centre, Indiana, and Snyder Counties), and north-central (Bradford, Potter, and Tioga Counties) Pennsylvania. Little to no data are available for approximately one-third of the state. Water characteristics and nutrients were the most frequently sampled major analyte groups; approximately 21,000 samples were collected for each group. Major and minor ions were the next most-frequently sampled major analyte groups; approximately 17,000 and 12,000 samples were collected, respectively. For the remaining eight major analyte groups, the number of samples collected ranged from a low of 307 samples (wastewater compounds) to a high of approximately 3,000 samples (biological). The number of samples that exceeded a maximum contaminant level (MCL) or secondary maximum contaminant level (SMCL) by major analyte group also varied. Of the 2,988 samples in the biological analyte group, 53 percent had water that exceeded an MCL. Almost 2,500 samples were collected and analyzed for volatile organic compounds; 14 percent exceeded an MCL. Other major analyte groups that frequently exceeded MCLs or SMCLs included major ions (17,465 samples and a 33.9 percent exceedence), minor ions (11,905 samples and a 17.1 percent exceedence), and water characteristics (21,183 samples and a 20.3 percent exceedence). Samples collected and analyzed for fungicides, herbicides, insecticides, and pesticides (4,062 samples), radiochemicals (1,628 samples), wastewater compounds (307 samples), and nutrients (20,822 samples) had the lowest exceedences of 0.3, 8.4, 0.0, and 8.8 percent, respectively.

Pennsylvania↗

Riparian Land Use/Land Cover Data for Three Study Units in Group II of the Nutrient Enrichment Effects Topical Study of the National Water-Quality Assessment Program

This data set was developed as part of the National Water-Quality Assessment (NAWQA) Program, Nutrient Enrichment Effects Topical (NEET) study. This report is concerned with three of the eight NEET study units distributed across the United States: Ozark Plateaus, Upper Mississippi River Basin, and Upper Snake River Basin, collectively known as Group II of the NEET study. Ninety stream reaches were investigated during 2006-08 in these three study units. Stream segments, with lengths equal to the base-10 logarithm of the basin area, were delineated upstream from the stream reaches through the use of digital orthophoto quarter-quadrangle (DOQQ) imagery. The analysis area for each stream segment was defined by a streamside buffer extending laterally to 250 meters from the stream segment. Delineation of landuse and land-cover (LULC) map units within stream-segment buffers was completed using on-screen digitizing of riparian LULC classes interpreted from the DOQQ. LULC units were classified using a strategy consisting of nine classes. National Wetlands Inventory (NWI) data were used to aid in wetland classification. Longitudinal riparian transects (lines offset from the stream segments) were generated digitally, used to sample the LULC maps, and partitioned in accord with the intersected LULC map-unit types. These longitudinal samples yielded the relative linear extent and sequence of each LULC type within the riparian zone at the segment scale. The resulting areal and linear estimates of LULC extent filled in the spatial-scale gap between the 30-meter resolution of the 1990s National Land Cover Dataset and the reach-level habitat assessment data collected onsite routinely for NAWQA ecological sampling. The resulting data consisted of 12 geospatial data sets: LULC within 25 meters of the stream reach (polygon); LULC within 50 meters of the stream reach (polygon); LULC within 50 meters of the stream segment (polygon); LULC within 100 meters of the stream segment (polygon); LULC within 150 meters of the stream segment (polygon); LULC within 250 meters of the stream segment (polygon); frequency of gaps in woody vegetation at the reach scale (arc); stream reaches (arc); longitudinal LULC transect sample at the reach scale (arc); frequency of gaps in woody vegetation at the segment scale (arc); stream segments (arc); and longitudinal LULC transect sample at the segment scale (arc).

Data Series↗

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming↗

Riparian Land Use/Land Cover Data for Five Study Units in the Nutrient Enrichment Effects Topical Study of the National Water-Quality Assessment Program

This dataset was developed as part of the National Water-Quality Assessment (NAWQA) Program, Nutrient Enrichment Effects Topical (NEET) study for five study units distributed across the United States: Apalachicola-Chattahoochee-Flint River Basin, Central Columbia Plateau-Yakima River Basin, Central Nebraska Basins, Potomac River Basin and Delmarva Peninsula, and White, Great and Little Miami River Basins. One hundred forty-three stream reaches were examined as part of the NEET study conducted 2003-04. Stream segments, with lengths equal to the logarithm of the basin area, were delineated upstream from the downstream ends of the stream reaches with the use of digital orthophoto quarter quadrangles (DOQQ) or selected from the high-resolution National Hydrography Dataset (NHD). Use of the NHD was necessary when the stream was not distinguishable in the DOQQ because of dense tree canopy. The analysis area for each stream segment was defined by a buffer beginning at the segment extending to 250 meters lateral to the stream segment. Delineation of land use/land cover (LULC) map units within stream segment buffers was conducted using on-screen digitizing of riparian LULC classes interpreted from the DOQQ. LULC units were mapped using a classification strategy consisting of nine classes. National Wetlands Inventory (NWI) data were used to aid in wetland classification. Longitudinal transect sampling lines offset from the stream segments were generated and partitioned into the underlying LULC types. These longitudinal samples yielded the relative linear extent and sequence of each LULC type within the riparian zone at the segment scale. The resulting areal and linear LULC data filled in the spatial-scale gap between the 30-meter resolution of the National Land Cover Dataset and the reach-level habitat assessment data collected onsite routinely for NAWQA ecological sampling. The final data consisted of 12 geospatial datasets: LULC within 25 meters of the stream reach (polygon); LULC within 50 meters of the stream reach (polygon); LULC within 50 meters of the stream segment (polygon); LULC within 100 meters of the stream segment (polygon); LULC within 150 meters of the stream segment (polygon); LULC within 250 meters of the stream segment (polygon); frequency of gaps in woody vegetation LULC at the reach scale (arc); stream reaches (arc); longitudinal LULC at the reach scale (arc); frequency of gaps in woody vegetation LULC at the segment scale (arc); stream segments (arc); and longitudinal LULC at the segment scale (arc).

Data Series↗

National Assessment of Shoreline Change; historical shoreline change along the New England and Mid-Atlantic coasts

Beach erosion is a chronic problem along many open-ocean shores of the United States. As coastal populations continue to grow and community infrastructures are threatened by erosion, there is increased demand for accurate information regarding past and present trends and rates of shoreline movement. There is also a need for a comprehensive analysis of shoreline movement that is consistent from one coastal region to another. To meet these national needs, the U.S. Geological Survey (USGS) is conducting an analysis of historical shoreline changes along open-ocean sandy shores of the conterminous United States and parts of Hawaii, Alaska, and the Great Lakes. One purpose of this work is to develop standard, repeatable methods for mapping and analyzing shoreline movement so that periodic, systematic, internally consistent updates regarding coastal erosion and land loss can be made nationally. In the case of this study, the shoreline is the interpreted boundary between the ocean water surface and the sandy beach. This report on the New England and Mid-Atlantic coasts is the fifth in a series of reports on historical shoreline change. Previous investigations include analyses and descriptive reports of the Gulf of Mexico, the Southeast Atlantic, and, for California, the sandy shoreline and the coastal cliffs. The rates of change presented in this report represent conditions up to the date of the most recent shoreline data and therefore are not intended for predicting future shoreline positions or rates of change. Because of the geomorphology of the New England and Mid-Atlantic (rocky coastlines, large embayments and beaches) as well as data gaps in some areas, this report presents beach erosion rates for 78 percent of the 1,360 kilometers of the New England and Mid-Atlantic coasts. The New England and Mid-Atlantic shores were subdivided into a total of 10 analysis regions for the purpose of reporting regional trends in shoreline change rates. The average rate of long-term shoreline change for the New England and Mid-Atlantic coasts was -0.5 meters per year with an uncertainty in the long-term trend of plus or minus 0.09 meters per year. The rate is based on shoreline change rates averaged from 21,184 individual transects, of which 68 percent were eroding. In both the long and short term, the average rates of shoreline change for New England and the Mid-Atlantic were erosional. Long-term erosion rates were generally lower in New England than in the Mid-Atlantic. This is a function of the dominant coastal geomorphology; New England has a greater percentage of shore types that tend to erode more slowly (rocky coasts, pocket beaches, and mainland beaches), whereas the Mid-Atlantic is dominated by more vulnerable barrier islands and dynamic spit/inlet environments. However, the percentage of coastline eroding was higher in New England than in the Mid-Atlantic, highlighting that although rates of shoreline erosion may not be extreme, coastal erosion is still widespread along this region of the U.S. coastline. The average rate of short-term shoreline change for the New England and Mid-Atlantic coasts was also erosional but the rate of erosion decreased in comparison to long-term rates. The net short-term rate as averaged along 17,045 transects was -0.3 meters per year. Uncertainties for these rates range from 0.06 to 0.1 meters per year depending on the data sources used in the rate calculations. Of transects used to measure short-term change, 60 percent were erosional, as compared to 65 percent of coast eroding in the long term. The slight decrease (5 percent) in the amount of coastline eroding may be related to an increase in the frequency and extent of nourishment programs and (or) the effects of hardened structures during the more recent time period. The most stable (lower rates of erosion) beaches were more commonly found in New England. Despite an overall lowering of the average rates of erosion from long-term to short-term, the amount

Open-File Report↗

USGS Streamgages Linked to the Medium Resolution NHD

The locations of approximately 23,000 current and historical U.S. Geological Survey (USGS) streamgages in the United States and Puerto Rico (with the exception of Alaska) have been snapped to the medium resolution National Hydrography Dataset (NHD). The NHD contains geospatial information about mapped surface-water features, such as streams, lakes, and reservoirs, etc., creating a hydrologic network that can be used to determine what is upstream or downstream from a point of interest on the NHD network. An automated snapping process made the initial determination of the NHD location of each streamgage. These initial NHD locations were comprehensively reviewed by local USGS personnel to ensure that streamgages were snapped to the correct NHD reaches. About 75 percent of the streamgages snapped to the appropriate NHD reach location initially and 25 percent required adjustment and relocation. This process resulted in approximately 23,000 gages being successfully snapped to the NHD. This dataset contains the latitude and longitude coordinates of the point on the NHD to which the streamgage is snapped and the location of the gage house for each streamgage. A process known as indexing may be used to create reference points (event tables) to the NHD reaches, expressed as a reach code and measure (distance along the reach). Indexing is dependent on the version of NHD to which the indexing is referenced. These data are well suited for use in indexing because nearly all the streamgage NHD locations have been reviewed and adjusted if necessary, to ensure they will index to the appropriate NHD reach. Flow characteristics were computed from the daily streamflow data recorded at each streamgage for the period of record. The flow characteristics associated with each streamgage include: *First date (year, month, day) of streamflow data *Last date (year, month, day) of streamflow data *Number of days of streamflow data *Number of days of non-zero streamflow data *Minimum and maximum daily flow for the period of record (cubic feet per second) *Percentiles (1, 5, 10, 20, 25, 50, 75, 80, 90, 95, 99) of daily flow for the period of record (cubic feet per second) *Average and standard deviation of daily flow for the period of record (cubic feet per second) *Mean annual base-flow index (BFI) computed for the period of record (fraction, ranging from 0 to 1) *Year-to-year standard deviation of the annual base-flow index computed for the period of record (fraction) *Number of years of data used to compute the base-flow index (years) The streamflow data used to compute flow characteristics were copied from the NWIS-Web historical daily discharge archive (nadww01.er.usgs.gov:/www/htdocs/nwisweb/data/discharge) on June 15, 2005.

Data Series↗

Sediment distribution and hydrologic conditions of the Potomac aquifer in Virginia and parts of Maryland and North Carolina

Sediments of the heavily used Potomac aquifer broadly contrast across major structural features of the Atlantic Coastal Plain Physiographic Province in eastern Virginia and adjacent parts of Maryland and North Carolina. Thicknesses and relative dominance of the highly interbedded fluvial sediments vary regionally. Vertical intervals in boreholes of coarse-grained sediment commonly targeted for completion of water-supply wells are thickest and most widespread across the central and southern parts of the Virginia Coastal Plain. Designated as the Norfolk arch depositional subarea, the entire sediment thickness here functions hydraulically as a single interconnected aquifer. By contrast, coarse-grained sediment intervals are thinner and less widespread across the northern part of the Virginia Coastal Plain and into southern Maryland, designated as the Salisbury embayment depositional subarea. Fine-grained intervals that are generally avoided for completion of water-supply wells are increasingly thick and widespread northward. Fine-grained intervals collectively as thick as several hundred feet comprise two continuous confining units that hydraulically separate three vertically spaced subaquifers. The subaquifers are continuous northward but merge southward into the single undivided Potomac aquifer. Lastly, far southeastern Virginia and northeastern North Carolina are designated as the Albemarle embayment depositional subarea, where both coarse- and fine-grained intervals are of only moderate thickness. The entire sediment thickness functions hydraulically as a single interconnected aquifer. A substantial hydrologic separation from overlying aquifers is imposed by the upper Cenomanian confining unit. Potomac aquifer sediments were deposited by a fluvial depositional complex spanning the Virginia Coastal Plain approximately 100 to 145 million years ago. Westward, persistently uplifted granite and gneiss source rocks sustained a supply of coarse-grained sand and gravel. Immature, high-gradient braided streams deposited longitudinal bars and channel fills across the Norfolk arch subarea. By contrast, across the Salisbury and Albemarle embayment subareas, mature, medium- to low-gradient meandering streams deposited medium- to coarse-grained channel fills and point bars segregated from fine-grained overbank deposits. The Virginia depositional complex merged northward across the Salisbury embayment subarea with another complex in Maryland. Here, additional sediments were received from schist source rocks that underwent three cycles of initial uplift and rapid erosion followed by crustal stability and erosional leveling. Because of the predominance of coarse-grained sediments, transmissivity, hydraulic conductivity, and regional velocities of lateral flow through the Potomac aquifer are greatest across the Norfolk arch depositional subarea, but decrease progressively northward with increasingly fine-grained sediments. Confining units hydraulically separate the Potomac aquifer from overlying aquifers, as indicated by large vertical hydraulic gradients. By contrast, most of the Potomac aquifer internally functions hydraulically as a single interconnected aquifer, as indicated by uniformly small vertical gradients. Most fine-grained sediments within the aquifer do not hydraulically separate overlying and underlying coarse-grained sediments. Across the Salisbury embayment depositional subarea, however, hydraulic separation among the vertically spaced subaquifers is imposed by the intervening confining units. The Potomac aquifer is the largest and most heavily used source of groundwater in the Virginia Coastal Plain. Water-level declines as great as 200 feet create the potential for saltwater intrusion. Conventional stratigraphic correlation has been generally ineffective at accurately characterizing complexly distributed fluvial sediments that compose the Potomac aquifer. Consequently, the aquifer’s internal hydraulic connectivity and overall hydrologic function have not been well understood. Water-supply planning and development efforts have been hampered, and interpretations of regulatory criteria for allowable water-level declines have been ambiguous. An investigation undertaken during 2010–11 by the U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, provides a comprehensive regional description of the spatial distribution of Potomac aquifer sediments and their relation to hydrologic conditions. Altitudes and thicknesses of 2,725 vertical sediment intervals represent the spatial distribution of Potomac aquifer sediments in the Virginia Coastal Plain and adjacent parts of Maryland and North Carolina. Sediment intervals are designated as either dominantly coarse or fine grained and were determined by interpretation of geophysical logs and ancillary information from 456 boreholes. Sediment-interval and borehole summary statistical data indicate regional trends in sediment lithology and stratigraphic continuity, upon which three structurally based and hydrologically distinct sediment depositional subareas are designated. Broad patterns of sediment deposition over time are inferred from published sediment pollen-age data. Discrepancies in previously drawn hydrostratigraphic relations between southeastern Virginia and northeastern North Carolina are partly resolved based on borehole geophysical logs and a recently documented geologic map and corehole. A conceptual model theorizes the depositional history of the sediments and geologically accounts for their distribution. Documented pumping tests of the Potomac aquifer at 197 locations produced 336 values of transmissivity and 127 values of storativity. Based on effective aquifer thicknesses, 296 values of sediment hydraulic conductivity and 113 values of sediment specific storage are calculated. Vertical hydraulic gradients are calculated from 9,479 pairs of water levels measured between November 17, 1953, and October 4, 2011, in 129 closely spaced pairs of wells. Borehole sediment-interval and related data provide a means to achieve high yielding production wells in the Potomac aquifer by site-specific targeting of drilling operations toward water-bearing coarse-grained sand and gravel. Advance knowledge of the potential of different parts of the aquifer also aids in planning optimal groundwater-development areas. Depositional subareas further provide a possible context for resource management. Current (2013) regulatory limits on water-level declines are relative to top surfaces of subdivided upper, middle, and lower Potomac aquifers across the entire Virginia Coastal Plain, but have the potential to exceed the same limit relative to a single undivided Potomac aquifer. By contrast, designation of the sediments as a single aquifer in the Norfolk arch and Albemarle embayment subareas—and as a series of vertically spaced subaquifers and intervening confining units in the Salisbury embayment subarea—best reflects understanding of the Potomac aquifer and can avoid the potential for excessive water-level declines. Simulation modeling to evaluate effects of groundwater withdrawals could be designed similarly, including vertical discretization and (or) zonation of the Potomac aquifer based on depositional subareas and a geostatistical distribution of aquifer properties derived from borehole sediment-interval data. Further resource-management information needs extend beyond the developed part of the Potomac aquifer, particularly across the Northern Neck and Middle Peninsula where only the shallowest part of the aquifer is known, and include structural aspects such as faults, basement bedrock, and the Chesapeake Bay impact crater.

Maryland, North Carolina, Virginia↗

Estimating Water Storage Capacity of Existing and Potentially Restorable Wetland Depressions in a Subbasin of the Red River of the North

Executive Summary Concern over flooding along rivers in the Prairie Pothole Region has stimulated interest in developing spatially distributed hydrologic models to simulate the effects of wetland water storage on peak river flows. Such models require spatial data on the storage volume and interception area of existing and restorable wetlands in the watershed of interest. In most cases, information on these model inputs is lacking because resolution of existing topographic maps is inadequate to estimate volume and areas of existing and restorable wetlands. Consequently, most studies have relied on wetland area to volume or interception area relationships to estimate wetland basin storage characteristics by using available surface area data obtained as a product from remotely sensed data (e.g., National Wetlands Inventory). Though application of areal input data to estimate volume and interception areas is widely used, a drawback is that there is little information available to provide guidance regarding the application, limitations, and biases associated with such approaches. Another limitation of previous modeling efforts is that water stored by wetlands within a watershed is treated as a simple lump storage component that is filled prior to routing overflow to a pour point or gaging station. This approach does not account for dynamic wetland processes that influence water stored in prairie wetlands. Further, most models have not considered the influence of human-induced hydrologic changes, such as land use, that greatly influence quantity of surface water inputs and, ultimately, the rate that a wetland basin fills and spills. The goals of this study were to (1) develop and improve methodologies for estimating and spatially depicting wetland storage volumes and interceptions areas and (2) develop models and approaches for estimating/simulating the water storage capacity of potentially restorable and existing wetlands under various restoration, land use, and climatic scenarios. To address these goals, we developed models and approaches to spatially represent storage volumes and interception areas of existing and potentially restorable wetlands in the upper Mustinka subbasin within Grant County, Minn. We then developed and applied a model to simulate wetland water storage increases that would result from restoring 25 and 50 percent of the farmed and drained wetlands in the upper Mustinka subbasin. The model simulations were performed during the growing season (May-October) for relatively wet (1993; 0.79 m of precipitation) and dry (1987; 0.40 m of precipitation) years. Results from the simulations indicated that the 25 percent restoration scenario would increase water storage by 21-24 percent and that a 50 percent scenario would increase storage by 34-38 percent. Additionally, we estimated that wetlands in the subbasin have potential to store 11.57-20.98 percent of the total precipitation that fell over the entire subbasin area (52,758 ha). Our simulation results indicated that there is considerable potential to enhance water storage in the subbasin; however, evaluation and calibration of the model is necessary before simulation results can be applied to management and planning decisions. In this report we present guidance for the development and application of models (e.g., surface area-volume predictive models, hydrology simulation model) to simulate wetland water storage to provide a basis from which to understand and predict the effects of natural or human-induced hydrologic alterations. In developing these approaches, we tried to use simple and widely available input data to simulate wetland hydrology and predict wetland water storage for a specific precipitation event or a series of events. Further, the hydrology simulation model accounted for land use and soil type, which influence surface water inputs to wetlands. Although information presented in this report is specific to the Mustinka subbasin, the approaches

Open-File Report↗

Magnitude, frequency, and trends of floods at gaged and ungaged sites in Washington, based on data through water year 2014

An investigation into the magnitude and frequency of floods in Washington State computed the annual exceedance probability (AEP) statistics for 648 U.S. Geological Survey unregulated streamgages in and near the borders of Washington using the recorded annual peak flows through water year 2014. This is an updated report from a previous report published in 1998 that used annual peak flows through the water year 1996. New in this report, a regional skew coefficient was developed for the Pacific Northwest region that includes areas in Oregon, Washington, Idaho and western Montana within the Columbia River drainage basin south of the United States-Canada border, the coastal areas of Oregon and western Washington, and watersheds draining into Puget Sound, Washington. The skew coefficient is an important term in the Log Pearson Type III equation used to define the distribution of the log-transformed annual peaks. The Expected Moments Algorithm was used to fit historical and censored peak-flow data to the log Pearson Type III distribution. A Multiple Grubb-Beck test was employed to censor low outliers of annual peak flows to improve on the frequency distribution. This investigation also includes a section on observed trends in annual peak flows that showed significant trends ( p -value < 0.05) in 21 of 83 long-term sites, but with small magnitude Kendall tau values suggesting a limited monotonic trend in the time series of annual peaks. Most of the sites with a significant trend in western Washington were positive and all the sites with significant trends (three sites) in eastern Washington were negative. Multivariate regression analysis with measured basin characteristics and the AEP statistics at long-term, unregulated, and un-urbanized (defined as drainage basins with less than 5 percent impervious land cover for this investigation) streamgages within Washington and some in Idaho and Oregon that are near the Washington border was used to develop equations to estimate AEP statistics at ungaged basins. Washington was divided into four regions to improve the accuracy of the regression equations; a set of equations for eight selected AEPs and for each region were constructed. Selected AEP statistics included the annual peak flows that equaled or exceeded 50, 20, 10, 4, 2, 1, 0.5 and 0.2 percent of the time equivalent to peak flows for peaks with a 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals, respectively. Annual precipitation and drainage area were the significant basin characteristics in the regression equations for all four regression regions in Washington and forest cover was significant for the two regression regions in eastern Washington. Average standard error of prediction for the regional regression equations ranged from 70.19 to 125.72 percent for Regression Regions 1 and 2 on the eastern side of the Cascade Mountains and from 43.22 to 58.04 percent for Regression Regions 3 and 4 on the western side of the Cascade Mountains. The pseudo coefficient of determination (where a value of 100 signifies a perfect regression model) ranged from 68.39 to 90.68 for Regression Regions 1 and 2, and 92.35 to 95.44 for Regions 3 and 4. The calculated AEP statistics for the streamgages and the regional regression equations are expected to be incorporated into StreamStats after the publication of this report. StreamStats is the interactive Web-based map tool created by the U.S. Geological Survey to allow the user to choose a streamgage and obtain published statistics or choose ungaged locations where the program automatically applies the regional regression equations and computes the estimates of the AEP statistics.

Washington↗

Natural-color and color-infrared image mosaics of the Colorado River corridor in Arizona derived from the May 2009 airborne image collection

The Grand Canyon Monitoring and Research Center (GCMRC) of the U.S. Geological Survey (USGS) periodically collects airborne image data for the Colorado River corridor within Arizona (fig. 1) to allow scientists to study the impacts of Glen Canyon Dam water release on the corridor’s natural and cultural resources. These data are collected from just above Glen Canyon Dam (in Lake Powell) down to the entrance of Lake Mead, for a total distance of 450 kilometers (km) and within a 500-meter (m) swath centered on the river’s mainstem and its seven main tributaries (fig. 1). The most recent airborne data collection in 2009 acquired image data in four wavelength bands (blue, green, red, and near infrared) at a spatial resolution of 20 centimeters (cm). The image collection used the latest model of the Leica ADS40 airborne digital sensor (the SH52), which uses a single optic for all four bands and collects and stores band radiance in 12-bits. Davis (2012) reported on the performance of the SH52 sensor and on the processing steps required to produce the nearly flawless four-band image mosaic (sectioned into map tiles) for the river corridor. The final image mosaic has a total of only 3 km of surface defects in addition to some areas of cloud shadow because of persistent inclement weather during data collection. The 2009 four-band image mosaic is perhaps the best image dataset that exists for the entire Arizona part of the Colorado River. Some analyses of these image mosaics do not require the full 12-bit dynamic range or all four bands of the calibrated image database, in which atmospheric scattering (or haze) had not been removed from the four bands. To provide scientists and the general public with image products that are more useful for visual interpretation, the 12-bit image data were converted to 8-bit natural-color and color-infrared images, which also removed atmospheric scattering within each wavelength-band image. The conversion required an evaluation of the histograms of each band’s digital-number population within each map tile throughout the corridor and the determination of the digital numbers corresponding to the lower and upper one percent of the picture-element population within each map tile. Visual examination of the image tiles that were given a 1-percent stretch (whereby the lower 1- percent 12-bit digital number is assigned an 8-bit value of zero and the upper 1-percent 12-bit digital number is assigned an 8-bit value of 255) indicated that this stretch sufficiently removed atmospheric scattering, which provided improved image clarity and true natural colors for all surface materials. The lower and upper 1-percent, 12-bit digital numbers for each wavelength-band image in the image tiles exhibit erratic variations along the river corridor; the variations exhibited similar trends in both the lower and upper 1-percent digital numbers for all four wavelength-band images (figs. 2–5). The erratic variations are attributed to (1) daily variations in atmospheric water-vapor content due to monsoonal storms, (2) variations in channel water color due to variable sediment input from tributaries, and (3) variations in the amount of topographic shadows within each image tile, in which reflectance is dominated by atmospheric scattering. To make the surface colors of the stretched, 8-bit images consistent among adjacent image tiles, it was necessary to average both the lower and upper 1-percent digital values for each wavelength-band image over 20 river miles to subdue the erratic variations. The average lower and upper 1-percent digital numbers for each image tile (figs. 2–5) were used to convert the 12-bit image values to 8-bit values and the resulting 8-bit four-band images were stored as natural-color (red, green, and blue wavelength bands) and color-infrared (near-infrared, red, and green wavelength bands) images in embedded geotiff format, which can be read and used by most geographic information system (GIS) and image-processing software. The tiff world files (tfw) are provided, even though they are generally not needed for most software to read an embedded geotiff image. All image data are projected in the State Plane (SP) map projection using the central Arizona zone (202) and the North American Datum of 1983 (NAD83). The map-tile scheme used to segment the corridor image mosaic followed the standard USGS quarter-quadrangle (QQ) map borders, but the high resolution (20 cm) of the images required further quarter segmentation (QQQ) of the standard QQ tiles, where the image mosaic covered a large fraction of a QQ map tile (segmentation shown in (figure 6), where QQ_1 to QQ_4 shows the number convention used to designate a quarter of a QQ tile). To minimize the size of each image tile, each image or map tile was subset to only include that part of the tile that had image data. In addition, some QQQ image tiles within a QQ tile were combined when adjacent QQQ map tiles were small. Thus, some image tiles consist of combinations of QQQ map tiles, some consist of an entire QQ map tile, and some consist of two adjoining QQ map tiles. The final image tiles number 143, which is a large number of files to list on the Internet for both the natural-color and color-infrared images. Thus, the image tiles were placed in seven file folders based on the one-half-degree geographic boundaries within the study area (fig. 7). The map tiles in each file folder were compressed to minimize folder size for more efficient downloading. The file folders are sequentially referred to as zone 1 through zone 7, proceeding down river (fig. 7). The QQ designations of the image tiles contained within each folder or zone are shown on the index map for each respective zone (figs. 8–14).

Arizona↗

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&ndash;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 &ldquo;confluence test&rdquo; 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&rsquo;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&rsquo;s (Petersen and others, 2008). This logic tree includes whole CSZ rupture earthquakes (M w 8.8&ndash;9.2) and partial CSZ rupture earthquakes (M w 8.0&ndash;8.7). In this latest version of the NSHM&rsquo;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↗

EAARL coastal topography-western Florida, post-Hurricane Charley, 2004: seamless (bare earth and submerged.

Project Description These remotely sensed, geographically referenced elevation measurements of lidar-derived seamless (bare-earth and submerged) topography were produced as a collaborative effort between the U.S. Geological Survey (USGS), Coastal and Marine Geology Program (CMGP), St. Petersburg, FL, and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, VA. This project provides highly detailed and accurate datasets of a portion of the western Florida coastline beachface, acquired post-Hurricane Charley on August 17 and 18, 2004. The datasets are made available for use as a management tool to research scientists and natural-resource managers. An innovative airborne lidar instrument originally developed at the NASA Wallops Flight Facility, and known as the Experimental Advanced Airborne Research Lidar (EAARL), was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532-nanometer) lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive lidar, a down-looking red-green-blue (RGB) digital camera, a high-resolution multispectral color infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers, and an integrated miniature digital inertial measurement unit, which provide for sub-meter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine Cessna 310 aircraft, but the instrument may be deployed on a range of light aircraft. A single pilot, a lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system, and the resulting data were then processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of lidar data in an interactive or batch mode. Modules for presurvey flight-line definition, flight-path plotting, lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is used routinely to create maps that represent submerged or sub-aerial topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation. For more information about similar projects, please visit the Decision Support for Coastal Science and Management website. Selected References Brock, J.C., Wright, C.W., Sallenger, A.H., Krabill, W.B., and Swift, R.N., 2002, Basis and methods of NASA airborne topographic mapper Lidar surveys for coastal studies: Journal of Coastal Research, v. 18, no. 1, p. 1-13. Crane, Michael, Clayton, Tonya, Raabe, Ellen, Stoker, Jason, Handley, Larry, Bawden, Gerald, Morgan, Karen, and Queija, Vivian, 2004, Report of the U.S. Geological Survey Lidar workshop sponsored by the Land Remote Sensing Program and held in St. Petersburg, FL, November 2002: U.S. Geological Survey Open-File Report 2004-1456, 72 p. Nayegandhi, Amar, Brock, J.C., and Wright, C.W., 2009, Small-footprint, waveform-resolving Lidar estimation of submerged and sub-canopy topography in coastal environments: International Journal of Remote Sensing, v. 30, no. 4, p. 861-878. Sallenger, A.H., Wright, C.W., and Lillycrop, Jeff, 2005, Coastal impacts of the 2004 hurricanes measured with airborne Lidar; initial results: Shore and Beach, v. 73, nos. 2-3, p. 10-14. Resources Included Readme.txt File

Data Series↗

Topobathymetric model of Mobile Bay, Alabama

Topobathymetric Digital Elevation Models (DEMs) are a merged rendering of both topography (land elevation) and bathymetry (water depth) that provides a seamless elevation product useful for inundation mapping, as well as for other earth science applications, such as the development of sediment-transport, sea-level rise, and storm-surge models. This 1/9-arc-second (approximately 3 meters) resolution model of Mobile Bay, Alabama was developed using multiple topographic and bathymetric datasets, collected on different dates. The topographic data were obtained primarily from the U.S. Geological Survey (USGS) National Elevation Dataset (NED) (http://ned.usgs.gov/) at 1/9-arc-second resolution; USGS Experimental Advanced Airborne Research Lidar (EAARL) data (2 meters) (http://pubs.usgs.gov/ds/400/); and topographic lidar data (2 meters) and Compact Hydrographic Airborne Rapid Total Survey (CHARTS) lidar data (2 meters) from the U.S. Army Corps of Engineers (USACE) (http://www.csc.noaa.gov/digitalcoast/data/coastallidar/). Bathymetry was derived from digital soundings obtained from the National Oceanic and Atmospheric Administration’s (NOAA) National Geophysical Data Center (NGDC) (http://www.ngdc.noaa.gov/mgg/geodas/geodas.html) and from water-penetrating lidar sources, such as EAARL and CHARTS. Mobile Bay is ecologically important as it is the fourth largest estuary in the United States. The Mobile and Tensaw Rivers drain into the bay at the northern end with the bay emptying into the Gulf of Mexico at the southern end. Dauphin Island (a barrier island) and the Fort Morgan Peninsula form the mouth of Mobile Bay. Mobile Bay is 31 miles (50 kilometers) long by a maximum width of 24 miles (39 kilometers) with a total area of 413 square miles (1,070 square kilometers). The vertical datum of the Mobile Bay topobathymetric model is the North American Vertical Datum of 1988 (NAVD 88). All the topographic datasets were originally referenced to NAVD 88 and no transformations were made to these input data. The NGDC hydrographic, multibeam, and trackline surveys were transformed from mean low water (MLW) or mean lower low water (MLLW) to NAVD 88 using VDatum (http://vdatum.noaa.gov). VDatum is a tool developed by the National Geodetic Survey (NGS) that performs transformations among tidal, ellipsoid-based, geoid-based, and orthometric datums using calibrated hydrodynamic models. The vertical accuracy of the input topographic data varied depending on the input source. Because the input elevation data were derived primarily from lidar, the vertical accuracy ranges from 6 to 20 centimeters in root mean square error (RMSE). he horizontal datum of the Mobile Bay topobathymetric model is the North American Datum of 1983 (NAD 83), geographic coordinates. All the topographic and bathymetric datasets were originally referenced to NAD 83, and no transformations were made to the input data. The bathymetric surveys were downloaded referenced to NAD 83 geographic, and therefore no horizontal transformations were required. The topbathymetric model of Mobile Bay and detailed metadata can be obtained from the USGS Web sites: http://nationalmap.gov/.

Alabama↗

Archive of digital boomer seismic reflection data collected during USGS field activity 02LCA02 in Lakes Ada, Crystal, Jennie, Mary, Rice, and Sylvan, Central Florida, July 2002

In July of 2002, the U.S. Geological Survey and St. Johns River Water Management District (SJRWMD) conducted geophysical surveys in Lakes Ada, Crystal, Jennie, Mary, Rice, and Sylvan, central Florida, as part of the USGS Lakes and Coastal Aquifers (LCA) study. This report serves as an archive of unprocessed digital boomer seismic reflection data, trackline maps, navigation files, Geographic Information System (GIS) files, and formal Federal Geographic Data Committee (FGDC) metadata. Filtered and gained (a relative increase in signal amplitude) digital images of the seismic profiles are also provided. Refer to the Acronyms page for expansions of acronyms and abbreviations used in this report. The archived trace data are in standard Society of Exploration Geophysicists (SEG) SEG-Y format (Barry and others, 1975) and may be downloaded and processed with commercial or public domain software such as Seismic Unix (SU). Example SU processing scripts and USGS software for viewing the SEG-Y files (Zihlman, 1992) are also provided. The USGS Florida Integrated Science Center (FISC) - St. Petersburg assigns a unique identifier to each cruise or field activity. For example, 02LCA02 tells us the data were collected in 2002 for the Lakes and Coastal Aquifers (LCA) study and the data were collected during the second field activity for that study in that calendar year. Refer to http://walrus.wr.usgs.gov/infobank/programs/html/definition/activity.html for a detailed description of the method used to assign the field activity ID. The boomer plate is an acoustic energy source that consists of capacitors charged to a high voltage and discharged through a transducer in the water. The transducer is towed on a sled floating on the water surface and when discharged emits a short acoustic pulse, or shot, which propagates through the water, sediment column, or rock beneath. The acoustic energy is reflected at density boundaries (such as the seafloor, sediment, or rock layers beneath the seafloor), detected by the receiver, and recorded by a PC-based seismic acquisition system. This process is repeated at timed intervals (for example, 0.5 s) and recorded for specific intervals of time (for example, 100 ms). In this way, a two-dimensional (2-D) vertical profile of the shallow geologic structure beneath the ship track is produced. Figure 1 displays the acquisition geometry. Refer to table 1 for a summary of acquisition parameters. Table 2 lists trackline statistics. The unprocessed seismic data are stored in SEG-Y format (Barry and others, 1975). For a detailed description of the data format, refer to the SEG-Y Format page. See the How To Download SEG-Y Data page for download instructions. The printable profiles provided here are GIF images that were filtered and gained using Seismic Unix software. Refer to the Software page for details about the processing and examples of the processing scripts. The processed SEG-Y data were exported to Chesapeake Technology, Inc. (CTI) SonarWeb software to produce an interactive Web page of the profile, which allows the user to obtain a geographic location and depth from the profile for a curser position. This information is displayed in the status bar of the browser.

Florida↗

California State Waters map series — Offshore of Aptos, California

Introduction In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California&rsquo;s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar bathymetric data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow subsurface geology. The Offshore of Aptos map area is located on the Pacific Coast, on the north side of Monterey Bay, about 105 km southeast of San Francisco. The largest incorporated city in the map area, Capitola, and numerous unincorporated towns including Aptos, lie on uplifted marine terraces between the shoreline, and the northwest-trending Santa Cruz Mountains, part of California&rsquo;s Coast Ranges. The map area includes the northernmost part of Santa Cruz Harbor, and Moss Landing Harbor is located about 10 km south of the map area. The offshore part of the map area is entirely within California&rsquo;s State Waters and is also part of the Monterey Bay National Marine Sanctuary. In the southern part of the map area, the Soquel Canyon State Marine Conservation Area extends eastward from the limit of California&rsquo;s State Waters across Soquel Canyon. The Offshore of Aptos map area is on the western margin of North American Plate&mdash;the only continental margin in the world delineated largely by transform faults. The San Andreas Fault Zone cuts through the Santa Cruz Mountains just 3.5 km northeast of the map area. The San Gregorio Fault Zone, another major plate-boundary structure, cuts through Monterey Canyon about 13 km southwest of the map area. Ongoing deformation associated with and between these major fault zones has uplifted the Santa Cruz Mountains and formed well-developed sets of marine terraces that characterize almost the entire coastal zone of the map area. The offshore part of the map area consists of relatively flat and shallow continental shelf which is underlain by variable amounts (0 to 30 m) of upper Quaternary shelf, estuarine, and fluvial sediments deposited as sea level fluctuated in the late Pleistocene. Along the southern edge of the map area, the shelf is incised by the head of Soquel Canyon, a northeast-trending tributary to the west-trending Monterey Canyon. During the last sea-level lowstand (the Last Glacial Maximum [LGM], about 21,000 years ago) Soquel Creek flowed through a paleochannel across the emergent shelf into the head of Soquel Canyon. This canyon was disconnected from its onshore watershed during the post-LGM sea-level rise of about 125 m; the abandoned paleochannel was subsequently filled with marine sediment. Coastal sediment in the Offshore of Aptos map area is supplied by coastal watersheds and bluff erosion. Sediment transport in this part of the Santa Cruz littoral cell is primarily from the northwest to the southeast and terminates in the submarine Monterey Canyon. Longshore drift is impeded by jetties at Santa Cruz Harbor, resulting in high beach erosion rates east of the harbor. Sediment dredged from the harbor mouth (estimated 300,000 yds 3 /yr) is currently being used to nourish beaches directly to the east. Farther downcoast, the rapidly eroding beach at Capitola was stabilized by construction of an about 75-m-long groin. This part of central California is exposed to large North Pacific swells from the northwest throughout the year. North Pacific swell heights range from 2 to 10 meters, with larger swells occurring from October to May. During El Ni&ntilde;o-Southern Oscillation (ENSO) events, winter storms track farther south than they do in normal (non-ENSO) years, thereby impacting the map area more frequently and with waves of larger heights. Bedrock exposed in coastal cliffs is relatively erosion-resistant, and significant erosional events primarily are restricted to storm-wave activity that also erodes the overlying unconsolidated marine-terrace sediments. The Offshore of Aptos map area lies within the cold-temperate biogeographic zone that is called either the &ldquo;Oregonian&rdquo; province or the &ldquo;northern California ecoregion.&rdquo; This biogeographic province is maintained by the long-term stability of the southward-flowing California Current, the eastern limb of the North Pacific subtropical gyre that flows from southern British Columbia to Baja California. At its midpoint off central California, the California Current transports subarctic surface (0&ndash;500 m deep) waters southward, about 150 to 1,300 km from shore. Seasonal northwesterly winds that are, in part, responsible for the California Current, generate coastal upwelling. The south end of the Oregonian province is at Point Conception (about 310 km southeast of the map area), although its associated phylogeographic group of marine fauna may extend beyond to the area offshore of Los Angeles in southern California. The ocean off of central California has experienced a warming over the last 50 years that is driving an ecosystem shift away from the productive subarctic regime towards a depopulated subtropical environment. Seafloor habitats in the Offshore of Aptos map area lie within the &ldquo;Shelf&rdquo; Megahabitat class and include patchy rocky habitats, gravel-rich scour depressions, minor bedrock habitat, and predominantly sandy inner shelf. Sandy shelf habitats grade offshore to mud-dominated habitats on the midshelf in deeper water. Biological productivity resulting from coastal upwelling supports populations of Sooty Shearwater, Western Gull, Common Murre, Cassin&rsquo;s Auklet, and many other less populous bird species. In addition, an observable recovery of Humpback and Blue Whales has occurred in the area; both species are dependent on coastal upwelling to provide nutrients. California sea lions and Pacific harbor seals are abundant in the map area. Common bottlenose dolphins are often observed very close to shore and in the surf zone. The large extent of exposed inner shelf bedrock supports large forests of &ldquo;bull kelp,&rdquo; which is well adapted for high-wave-energy environments. The kelp beds are the northernmost known habitat for the population of southern sea otters. Common fish species found in the kelp beds and rocky reefs include blue rockfish, black rockfish, olive rockfish, kelp rockfish, gopher rockfish, black-and-yellow rockfish, painted greenling, kelp greenling, and lingcod.

California↗