Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geologic 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 1,207 records · Page 67Linked to original sources

Ground water in folded Cretaceous sandstone of the Bhachau area, Kutch, India, with reference to the Kandla Port water supply

This report is based on an investigation of the availability of ground-water supplies in the Bhachau area for the nearby Kandla Port and township development undertaken by the Government of India. This seaport lies on an estuary of the Gulf of Kutch in western India and in the eastern part of the State of Kutch. The fieldwork on the investigation was carried on from November 1952 through April 1953 with continuing hydrologic observations through 1954-55. The fieldwork included: geologic mapping and delimitation of the principal aquifers of the region; preparation of water-table maps; a detailed inventory of existing wells and springs; observations of significant water table fluctuations; preparation of isobicarb, isochlor and isosulf maps to show the areal distribution of ground-water salinity. The Bhachau area includes about 116 square miles in eastern Kutch and lies in a belt of semiarid low-latitude steppes. The mean annual rainfall is about 15 inches, most of which falls from late June to late September during the southwest monsoon. The area includes a central sandy upland ranging from about 100 to 250 feet above sea level ; a northern lowland of between about 50 to 125 feet altitude that slopes north to the Great Rann of Kutch; a belt of low buttes and discontinuous ridges ranging from about 200 to 275 feet above sea level; and southern lowland which slopes in a southerly to southeasterly direction from an altitude of about 125 feet to 25 feet or less near the Gulf of Kutch. The principal streams are Kageshwar Vokra and Kara Vokra which drain north to the Great Rann and Kotwala Vokra and Dalwala Vokra which drain south toward the Gulf of Kutch. The rocks of the Bhachau area include nonmarine and marine sediments of Mesozoic, Tertiary, and Quaternary age and volcanic rocks of late Mesozoic to early Tertiary age. The oldest rocks in the area are medium- to coarse-grained white to buff current-bedded friable sandstone with occasional partings of white silty shale of the Upper Bhuj series that has been assigned to the Early Cretaceous. The soft friable sandstone of the Upper Bhuj series constitutes the most productive ground-water reservoir in the Bhachau area. At present (1955) there are nine irrigated tracts for which water is obtained from dug wells less than 90 feet deep in the Upper Bhuj. These wells are worked by bullocks and 'motes' (leather bags) at withdrawal rates ranging from about 6,000 to 24,000imperial gallons per day; however, many existing individual wells if equipped with mechanical pumps are capable of yielding 100,000 gallons per day. The Deccan trap of Late Cretaceous to Eocene age occurs in a sequence of basaltic lava flows in the Bhachau area, but trap dikes, sills and plugs that are common in other parts of Kutch have not been observed in the area. Laterite of probable Eocene age is extensive at the top of the Deccan trap, and in places where the lava flows are thin the parent rock has been almost completely lateritized. The Deccan lava flows or the laterite, where the trap is absent, rest disconformably on the Upper Bhuj. No wells have been observed in the Deccan trap of the Bhachau area, but it is possible that locally small supplies of good water may be obtained from these rocks. The Tertiary sediments, which are assigned to the Manchhar series of Pliocene age, generally rest on the laterite or the Deccan trap ; but where both are absent, the Manchhar rests directly on the Upper Bhuj. The Manchhar series includes massive reddish-brown gypseous clay shales, laminated gray siltstones, some limestone, mottled sandstone, and laterite trap gravel. Only meager supplies of brackish water are obtained from wells in these sediments. Along the channels of Kotwala, Dalwala, Kageshwar, and Kara Vokras are narrow bands of unconsolidated coarse sand with fine gravel of Quaternary age. No wells were observed in these deposits, but it is possible that locally they may contain small supplies of bracki

Water Supply Paper↗

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↗

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↗

Archive of Side Scan Sonar and Swath Bathymetry Data collected during USGS Cruise 10CCT02 Offshore of Petit Bois Island Including Petit Bois Pass, Gulf Islands National Seashore, Mississippi, March 2010

In March of 2010, the U.S. Geological Survey (USGS) conducted geophysical surveys offshore of Petit Bois Island, Mississippi, and Dauphin Island, Alabama (fig. 1). These efforts were part of the USGS Gulf of Mexico Science Coordination partnership with the U.S. Army Corps of Engineers (USACE) to assist the Mississippi Coastal Improvements Program (MsCIP) and the Northern Gulf of Mexico (NGOM) Ecosystem Change and Hazards Susceptibility Project by mapping the shallow geologic stratigraphic framework of the Mississippi Barrier Island Complex. These geophysical surveys will provide the data necessary for scientists to define, interpret, and provide baseline bathymetry and seafloor habitat for this area and to aid scientists in predicting future geomorphological changes of the islands with respect to climate change, storm impact, and sea-level rise. Furthermore, these data will provide information for barrier island restoration, particularly in Camille Cut, and protection for the historical Fort Massachusetts on Ship Island, Mississippi. For more information please refer to http://ngom.usgs.gov/gomsc/mscip/index.html. This report serves as an archive of the processed swath bathymetry and side scan sonar data (SSS). Data products herein include gridded and interpolated surfaces, seabed backscatter images, and ASCII x,y,z data products for both swath bathymetry and side scan sonar imagery. Additional files include trackline maps, navigation files, GIS files, Field Activity Collection System (FACS) logs, and formal FGDC metadata. Scanned images of the handwritten and digital FACS logs are also provided as PDF files. Refer to the Acronyms page for expansion of acronyms and abbreviations used in this report.

Data Series↗

Cobalt-copper deposits of the Blackbird district, Lemhi County, Idaho

The Blackbird district is in east-central Idaho, about 20 miles west-southwest of Salmon. The area is one of deeply weathered, flat-topped upland surfaces cut by several steep-walled valleys, which are tributary to the canyon of Panther Creek. Most of the area has a heavy vegetative cover and outcrops are relatively scarce except in the walls of the steeper valleys. The rocks of the district consist mostly of metamorphosed sedimentary rocks of the pre-Cambrian Yellowjacket formation (Belt series); a part of the Cretaceous Idaho batholith cuts across the northern part of the district, and acid porphyry dikes and metamorphosed basic rocks cut the Yellowjacket rocks. Structurally, the sedimentary rocks are divided by faults into three, roughly north-south blocks. The center one (Blackbird structural block) appears to have been more tightly squeezed than the others into relatively tight folds, with the development of widespread schistosity (flow cleavage). The rocks of the two outside blocks are in more open folds. In general they are nonschistose, except for the north end of the western block, where there are schistose rocks cut by several north-dipping thrust faults. The northern parts of the central and western blocks contain considerable garnet, chloritoid, and cordierite. The Blackbird structural block is cut by a number of mineralized shear zones. Those dipping moderately northeast and striking northwest, and those dipping steeply and striking north and northeast; appear to be most important. The mineralized rock contains chalcopyrite, cobaltite, pyrite, and pyrrhotite into a gangue of quartz, biotite, tourmaline, ankerite, and muscovite; the deposits were formed mostly by replacement of the shear zones. In addition, the block is cut by north-dipping thrust faults of west to northwest strike, and a number of high-angle faults. The district was first prospected about 1893; considerable developmont was done at the Brown Bear mine in 1899-1902, at the Haynes-Stellite in l917-1920, and at the Uncle Sam mine in 1938-1941, During World War II, the U.S. Bureau of Mines explored in the district with bulldozer and diamond drill., and the U.S. Geological Survey mapped the district and logged the drill cores. The Howe Sound Co. also did some diamond drilling in the district, and in 1945 the Calera Mining Co., Blackbird Division (subsidiary of Howe Sound Co.), started underground development at the Calera adit. Although the district has had very little production to the present, it is believed that a large tonnage of copper-cobalt ore exists in the district which should permit mining to be continued over a long period. The report contains brief descriptions of all the accessible workings in the district, of which the most important are Calera, Brown Bear, Uncle Sam, and Hawkeye mines. In the Calera adit, about 1,700 feet of the mineralized zone, ranging in width from 3 feet to 40 feet and averaging about 15 feet; have been explored (August 1946); the zone lies on a wide northwest-striking shear zone dipping moderately ( 60° ±) northeast. The Brown Bear adit is in a wide, mineralized, north-south shear zone in which are higher-grade pods plunging 25° to 35° north. The Uncle Sam mine explores a relatively narrow north-south shear zone in which are two or three north-plunging ore shoots. The Hawkeye mine is in a broad zone of mineralized schist in which are several north-plunging lenses of ore.

Idaho↗

Terrestrial lidar datasets of New Orleans, Louisiana, levee failures from Hurricane Katrina, August 29, 2005

Hurricane Katrina made landfall with the northern Gulf Coast on August 29, 2005, as one of the strongest hurricanes on record. The storm damage incurred in Louisiana included a number of levee failures that led to the inundation of approximately 85 percent of the metropolitan New Orleans area. Whereas extreme levels of storm damage were expected from such an event, the catastrophic failure of the New Orleans levees prompted a quick mobilization of engineering experts to assess why and how particular levees failed. As part of this mobilization, civil engineering members of the United States Geological Survey (USGS) performed terrestrial lidar topographic surveys at major levee failures in the New Orleans area. The focus of the terrestrial lidar effort was to obtain precise measurements of the ground surface to map soil displacements at each levee site, the nonuniformity of levee height freeboard, depth of erosion where scour occurred, and distress in structures at incipient failure. In total, we investigated eight sites in the New Orleans region, including both earth and concrete floodwall levee breaks. The datasets extend from the 17th Street Canal in the Orleans East Bank area to the intersection of the Gulf Intracoastal Waterway (GIWW) with the Mississippi River Gulf Outlet (MRGO) in the New Orleans East area. The lidar scan data consists of electronic files containing millions of surveyed points. These points characterize the topography of each levee’s postfailure or incipient condition and are available for download through online hyperlinks. The data serve as a permanent archive of the catastrophic damage of Hurricane Katrina on the levee systems of New Orleans. Complete details of the data collection, processing, and georeferencing methodologies are provided in this report to assist in the visualization and analysis of the data by future users.

Louisiana↗

Hydrologic conditions and terrestrial laser scanning of post-firedebris flows in the San Gabriel Mountains, CA, U.S.A

To investigate rainfall-runoff conditions that generate post-wildfire debris flows, we instrumented and surveyed steep, small watersheds along the tectonically active front of the San Gabriel Mountains, California. Fortuitously, we recorded runoff-generated debris-flows triggered by one spatially restricted convective event with 28 mm of rainfall falling over 62 minutes. Our rain gages, nested hillslope overland-flow sensors and soil-moisture probes, as well as a time series of terrestrial laser scanning (TLS) revealed the effects of the storm. Hillslope overland-flow response, along two ~10-m long flow lines perpendicular to and originating from a drainage divide, displayed only a 10 to 20 minute delay from the onset of rainfall with accumulated totals of merely 5-10 mm. Depth-stratified soil-moisture probes displayed a greater time delay, roughly 20- 30 minutes, indicating that initial overland flow was Hortonian. Furthermore, a downstream channel-monitoring array recorded a pronounced discharge peak generated by the passage of a debris flow after 18 minutes of rainfall. At this time, only four of the eleven hillslope overland flow sensors confirmed the presence of surface-water flow. Repeat TLS and detailed field mapping using GPS document how patterns of rainsplash, overland-flow scour, and rilling contributed to the generation of metter-scale debris flows. In response to a single small storm, the debris flows deposited irregular levees and lobate terminal snouts on hillslopes and caused wide- spread erosion of the valley axis with ground surface lowering exceeding 1.5 m.

California↗

Candidate-penetrative-fracture mapping of the Grand Canyon area, Arizona, from spatial correlation of deep geophysical features and surficial lineaments

Some aquifers of the southwestern Colorado Plateaus Province are deeply buried and overlain by several impermeable shale layers, and so recharge to the aquifer probably is mainly by seepage down penetrative-fracture systems. The purpose of this 2-year study, sponsored by the U.S. National Park Service, was to map candidate deep penetrative fractures over a 120,000-km2 area, using gravity and aeromagnetic-anomaly data together with surficial-fracture data. The study area was on the Colorado Plateau south of the Grand Canyon and west of Black Mesa; mapping was carried out at a scale of 1:250,000. The resulting database constitutes a spatially registered estimate of deep-fracture locations. Candidate penetrative fractures were located by spatial correlation of horizontal- gradient and analytic-signal maximums of gravity and magnetic anomalies with major surficial lineaments obtained from geologic, topographic, side-looking-airborne-radar, and satellite imagery. The maps define a subset of candidate penetrative fractures because of limitations in the data coverage and the analytical technique. In particular, the data and analytical technique used cannot predict whether the fractures are open or closed. Correlations were carried out by using image-processing software, such that every pixel on the resulting images was coded to uniquely identify which datasets are correlated. The technique correctly identified known and many new deep fracture systems. The resulting penetrative-fracture-distribution maps constitute an objectively obtained, repeatable dataset and a benchmark from which additional studies can begin. The maps also define in detail the tectonic fabrics of the southwestern Colorado Plateaus Province. Overlaying the correlated lineaments on the normalized-density-of-vegetation-index image reveals that many of these lineaments correlate with the boundaries of vegetation zones in drainages and canyons and so may be controlling near-surface water availability in some places. Many derivative products can be produced from the database, such as fracture-density-estimate maps, and maps with the number of correlations color-coded to estimate the possible quality of correlation. The database contained in this report is designed to be used in a geographic information system and image-processing systems, and most data layers are in georeferenced tagged image format (Geotiff) or ARC grids. The report includes 163 map plates and various metadata, supporting, and statistical diagram files.

Data Series↗

Seismic history of the Middle America subduction zone along El Salvador, Guatemala, and Chiapas, Mexico: 1526–2000

We present a catalog of subduction zone earthquakes along the Pacific coast from central El Salvador to eastern Chiapas, Mexico, from 1526 to 2000. We estimate that the catalog is complete since 1690 for M S ≥7.4 thrust events and M ≥ 7.4 normal-faulting events within the upper 60 km of the down-going slab. New intensity maps were constructed for the 27 earthquakes since 1690, using mostly primary data sources. By calibrating with recent events we find that the long axis of the (MM) VII intensity contour for such large earthquakes well approximates the length and location of rupture along the subduction zone and can thus be used to estimate the locations and magnitudes of older events. The section from western El Salvador to Chiapas appears to have ruptured completely in a series of four to five earthquakes during each of the periods 1902–1915, 1743–1776, and possibly 1565–1577. Earthquakes of M W 7.75 ± 0.3 have caused major damage along the 200 km long section from San Salvador to Guatemala City every 71 ± 17 yr, apparently since at least 1575. Although the January 2001 El Salvador earthquake caused damage within part of this zone, no major thrust earthquake has occurred there since at least 1915. We find that much of this section has been relatively quiescent for moderate earthquakes shallower that 50 km since at least 1963. The conditional probability that an earthquake of M W 7.75 ± 0.3 will occur at this location in the next 20 yr is estimated at 50% (±30%).

Chiapas↗

Coastal Science Navigator companion guide—Discover the U.S. Geological Survey coastal science products you need

The Coastal Science Navigator is an online gateway to a wide variety of U.S. Geological Survey (USGS) coastal change hazards-related information, data, and tools relevant to stakeholders’ scientific and decision-making needs. The products within the Coastal Science Navigator provide data related to past, present, and future threats to our coastlines. The filter search allows users to see all available products and identify relevant options by time scale, geographic scope, coastal hazard theme, and other filters. The guided search suggests products based on users’ answers to a short series of questions. A comprehensive summary is available for each product. The idea for the Coastal Science Navigator arose in 2020 in response to stakeholder feedback identifying the need for a central source for USGS coastal science information. It was published in July 2023 and initially included 55 products. Regular updates are planned to integrate other existing and new products. This guide introduces some of the many coastal change hazards-related products available through the USGS. In it, we showcase the products included in the Coastal Science Navigator’s initial publication in July 2023. While it is not representative of all the information, tools, and data available, we hope it serves as a compelling snapshot of what the USGS has to offer and encourages you to explore the Coastal Science Navigator to discover more of the products you need. To navigate this guide, the products have been organized by the time scale they are best suited for—past, present, or future—although many products cover multiple time scales. An additional section features software, one of the many product types available as filters within the Coastal Science Navigator. Other products include downloadable data, websites, and geonarratives (web pages that combine text, images, and interactive maps into narratives you can scroll through). Featured geographic scopes are also highlighted within this guide, detailing some of the many regions in which the USGS conducts research and illustrating another way to filter products within the Coastal Science Navigator.

Circular↗

An Automated Cropland Classification Algorithm (ACCA) for Tajikistan by combining Landsat, MODIS, and secondary data

The overarching goal of this research was to develop and demonstrate an automated Cropland Classification Algorithm (ACCA) that will rapidly, routinely, and accurately classify agricultural cropland extent, areas, and characteristics (e.g., irrigated vs. rainfed) over large areas such as a country or a region through combination of multi-sensor remote sensing and secondary data. In this research, a rule-based ACCA was conceptualized, developed, and demonstrated for the country of Tajikistan using mega file data cubes (MFDCs) involving data from Landsat Global Land Survey (GLS), Landsat Enhanced Thematic Mapper Plus (ETM+) 30 m, Moderate Resolution Imaging Spectroradiometer (MODIS) 250 m time-series, a suite of secondary data (e.g., elevation, slope, precipitation, temperature), and in situ data. First, the process involved producing an accurate reference (or truth) cropland layer (TCL), consisting of cropland extent, areas, and irrigated vs. rainfed cropland areas, for the entire country of Tajikistan based on MFDC of year 2005 (MFDC2005). The methods involved in producing TCL included using ISOCLASS clustering, Tasseled Cap bi-spectral plots, spectro-temporal characteristics from MODIS 250 m monthly normalized difference vegetation index (NDVI) maximum value composites (MVC) time-series, and textural characteristics of higher resolution imagery. The TCL statistics accurately matched with the national statistics of Tajikistan for irrigated and rainfed croplands, where about 70% of croplands were irrigated and the rest rainfed. Second, a rule-based ACCA was developed to replicate the TCL accurately (∼80% producer’s and user’s accuracies or within 20% quantity disagreement involving about 10 million Landsat 30 m sized cropland pixels of Tajikistan). Development of ACCA was an iterative process involving series of rules that are coded, refined, tweaked, and re-coded till ACCA derived croplands (ACLs) match accurately with TCLs. Third, the ACCA derived cropland layers of Tajikistan were produced for year 2005 (ACL2005), same year as the year used for developing ACCA, using MFDC2005. Fourth, TCL for year 2010 (TCL2010), an independent year, was produced using MFDC2010 using the same methods and approaches as the one used to produce TCL2005. Fifth, the ACCA was applied on MFDC2010 to derive ACL2010. The ACLs were then compared with TCLs (ACL2005 vs. TCL2005 and ACL2010 vs. TCL2010). The resulting accuracies and errors from error matrices involving about 152 million Landsat (30 m) pixels of the country of Tajikistan (of which about 10 million Landsat size, 30 m, cropland pixels) showed an overall accuracy of 99.6% (k hat = 0.97) for ACL2005 vs. TCL2005. For the 3 classes (irrigated, rainfed, and others) mapped in ACL2005, the producer’s accuracy was >86.4% and users accuracy was >93.6%. For ACL2010 vs. TCL2010, the error matrix showed an overall accuracy on 96.2% (k hat = 0.96). For the 3 classes (irrigated, rainfed, and others) mapped in ACL2010, the producer’s and user’s accuracies for the irrigated areas were ≥82.9%. Any intermixing was overwhelmingly between irrigated and rainfed croplands, indicating that croplands (irrigated plus rainfed areas) as well as irrigated areas were mapped with high levels of accuracies (∼90% or higher) even for the independent year. The ACL2005 and ACL2010, each, were produced using ACCA algorithm in ∼30 min using a Dell Precision desktop T7400 computer for the entire country of Tajikistan once the MFDCs for the years were ready. The ACCA algorithm for Tajikistan is made available through US Geological Survey’s ScienceBase: http://www.sciencebase.gov/catalog/folder/4f79f1b7e4b0009bd827f548 or at: https://powellcenter.usgs.gov/globalcroplandwater/content/models-algorithms . The research contributes to the efforts of global food security through research on global croplands and their water use (e.g., https://powellcenter.usgs.gov/globalcroplandwater/ ). The above results clearly demonstrated the ability of a rule-based ACCA to rapidly and accurately produce cropland data layer year after year (hindcast, nowcast, forecast) for the country it was developed using MFDCs that consist of combining multiple sensor data and secondary data. It needs to be noted that the ACCA is applicable to the area (e.g., country, region) for which it is developed. In this case, ACCA is applicable for the Country of Tajikistan to hindcast, nowcast, and forecast agricultural cropland extent, areas, and irrigated vs. rainfed. The same fundamental concept of ACCA applies to other areas of the World where ACCA codes need to be modified to suite the area/region of interest. ACCA can also be expanded to compute other crop characteristics such as crop types, cropping intensities, and phenologies.

Remote Sensing↗

Digital map of the aquifer boundary for the High Plains aquifer in parts of Colorado, Kansas, Nebraska, New Mexico, Oklahoma, South Dakota, Texas, and Wyoming

This digital data set represents the extent of the High Plains aquifer in the central United States. The extent of the High Plains aquifer covers 174,000 square miles in eight states: Colorado, Kansas, Nebraska, New Mexico, Oklahoma, South Dakota, Texas, and Wyoming. This data set represents a compilation of information from digital and paper sources and personal communication. This boundary is an update to the boundary published in U.S. Geological Survey Professional Paper 1400-B, and this report supersedes Open-File Report 99-267. The purpose of this data set is to refine and update the extent of the High Plains aquifer based on currently available information. This data set represents a compilation of arcs from a variety of sources and scales that represent the 174,000 square-mile extent of the High Plains aquifer within the eight states. Where updated information was not available, the original boundary extent defined by OFR 99-267 was retained. The citations for the sources in each State are listed in the 00README.txt file. The boundary also contains internal polygons, or 'islands', that represent the areas within the aquifer boundary where the aquifer is not present due to erosion or non-deposition. The datasets that pertain to this report can be found on the U.S. Geological Survey's NSDI (National Spatial Data Infrastructure) Node, the links are provided on the sidebar.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗

Seafloor gas seeps triggered by a northern California earthquake

A series of marine surveys using high-resolution seismic-reflection and side-scan sonar profiling was conducted on the Klamath River delta from 1977 to 1985. In November 1980, near the middle of this 8-year time period, a major earthquake (magnitude 7) occurred in the area. The earthquake caused liquefaction of shelf deposits off the Klamath River, and while mapping the sediment failures, we observed an inordinate number of gas vents and acoustic anomalies in the sediment (gas charging) and water column (gas clouds). Our analyses indicate that a change in gas seepage occurred after the earthquake. There was a marked increase in gas abundance immediately after the earthquake, a gradual decrease over the next year, and stabilization at pre-earthquake gas levels within five years or less. The acoustic anomalies were observed in pre-Quaternary rocks as well as in surficial Holocene deposits, and gas seeps overlay truncated anticlines, faults, and syncline exposures. Thus, the gas released may have been in part thermogenic in origin. A large area (20 km 2 ) on the delta failed during the earthquake, and the abundance of gas seeps and craters in the area suggests that sediment degassing may have been an important process contributing to failure.

Marine Geology↗

Geology of the Ishmas gold district, Kingdom of Saudi Arabia

The Ishmas gold district was mapped at 1:25,000 scale to place auriferous mineralization into geologic perspective, to assist in creating an ore-deposit model, and to aid in devising a strategy for subsequent exploration elsewhere in the Jabal Ishmas-Wadi Tathlith gold belt. The precratonic evolution of the district began with the deposition of a tholeiitic mafic volcanic and volcaniclastic sequence that was intruded by diabase and gabbro. Basaltic to rhyodacitic volcanism following a calc-alkaline evolutionary trend ensued. Subsequent deposition of a thick wacke and sandstone unit represented the final phase in the volcano sedimentary accumulation. The emplacement of a large lopolithic layered-gabbro complex marked the end of the precratonic evolutionary cycle. The district coincides with the boundary of two allochthonous terranes. The collisional Nabitah orogeny represents the suturing of the two terranes. The effects of this event are manifested by numerous north-trending, steeply dipping faults, shear zones, and mylonite belts, as well as diapiric serpentinite. During the waning stages of the orogeny, auriferous quartz pods were precipitated in dilatant structures within the north-trending shear zones by deep-seated circulating fluids. The emplacement of a tonalite stock was closely followed by the formation of N. 20°-35° W. -trending faults. These faults influenced the emplacement of dacite porphyry stocks and associated auriferous quartz veins. The auriferous veins are massive, tabular open-fracture fillings that are spatially, temporally, and genetically related to the dacite porphyry. The emplacement of a quartz monzodiorite stock was responsible for additional auriferous quartz vein mineralization that is almost exclusively hosted by the intrusion. A nearly random orientation of the veins indicates that no regional structure influenced their formation. The formation of a series of N. 60°-80° W -trending faults represents the final episode in the district's geologic history and corresponds to the cratonic reactivation that affected a large part of the Arabian Shield (that is, the Najd faulting event).

Open-File Report↗

Geologic history of the Yosemite Valley

The problem of the origin of the Yosemite Valley inherently demands a solution in quantitative terms. Its essence is, To what extent is the valley a product of glacial action, to what extent a product of stream erosion? The principal result of the investigations upon which this report is based is the determination within narrow limits of the preglacial depth of the Yosemite Valley and of other facts concerning its preglacial development which permit fairly definite estimates of the proportionate shares of work performed by stream and by glacier. The investigations comprise a detailed survey of the glacial and geomorphologic features of the Yosemite region and an equally intensive study of its rock formations, supplemented by reconnaissance work of both kinds in adjoining parts of the Sierra Nevada. The petrologic studies were made by Frank C. Calkins; the glacial and geomorphologic survey by Francois E. Matthes. Detailed mapping of the morainal system of the ancient Yosemite Glacier has served not only to determine the farthest limits reached by that glacier but to throw new light on the significance of the hanging valleys of the Yosemite region. It is reasonably certain that the glacier never extended more than about a mile beyond the site of El Portal. The hanging side valleys of the Merced Canyon below El Portal, therefore, hang not because of any glacial deepening suffered by that canyon. The explanation is offered that they hang because their streamlets have been unable to trench as rapidly as the Merced River since the rejuvenation of the Merced by the last uptilting of the Sierra Nevada. The streamlets were handicapped not only by their comparatively small volume but also by the fact that their courses trend northwestward and southeastward, substantially at right angles to the direction of the tilting, and therefore have remained essentially unsteepened, whereas the Merced's course trends southwestward, directly down the slope of the Sierra block, and therefore has been appreciably steepened. Projection of the longitudinal profiles of these hanging valleys forward to the axis of the Merced Canyon shows that they are closely accordant in height. Their profiles indicate a series of points on a former profile of the Merced with respect to which the side streams had graded their courses prior to the last uplift. This old profile can be extended upward into the glaciated part of the Merced Canyon above El Portal and even into the profoundly glaciated Yosemite Valley, accordant points being furnished by a number of hanging side valleys (due allowance being made for glacial erosion suffered by those valleys). However, not all the hanging valleys of the Yosemite region are accordant with this set. Several of them, including the upland valley of Yosemite Creek, constitute a separate set indicating another old profile of the Merced at a level 600 to 1,000 feet higher than the first. Others, including the hanging gulch of lower Bridalveil Creek, point to an old profile of the Merced about 1,200 feet lower than the first. There are thus three distinct sets of hanging valleys produced in three cycles of stream erosion. The valleys of the upper set, like those of the middle set, were left hanging as a result of rapid trenching by the Merced induced by an uplift of the range, there having been two such uplifts. Only the valleys of the lower set hang because of glacial deepening and widening of the Yosemite Valley, the cycle in which they were cut having been interrupted by the advent of the Pleistocene glaciers. They consequently indicate the preglacial depth of the Yosemite Valley. That depth, measured from the brow of El Capitan, was about 2,400 feet; measured from the rim at Glacier Point it was about 2,000 feet.

California↗

Compilation of V S 30 Data for the United States

V S 30 , the time-averaged shear-wave velocity ( V S ) to a depth of 30 meters, is a key index adopted by the earthquake engineering community to account for seismic site conditions. V S 30 is typically based on geophysical measurements of V S derived from invasive and noninvasive techniques at sites of interest. Owing to cost considerations, as well as logistical and environmental concerns, V S 30 data are sparse or not readily available for most areas. Where data are available, V S 30 values are often assembled in assorted formats that are accessible from disparate and (or) impermanent Web sites. To help remedy this situation, we compiled V S 30 measurements obtained by studies funded by the U.S. Geological Survey (USGS) and other governmental agencies. Thus far, we have compiled V S 30 values for 2,997 sites in the United States, along with metadata for each measurement from government-sponsored reports, Web sites, and scientific and engineering journals. Most of the data in our V S 30 compilation originated from publications directly reporting the work of field investigators. A small subset (less than 20 percent) of V S 30 values was previously compiled by the USGS and other research institutions. Whenever possible, V S 30 originating from these earlier compilations were crosschecked against published reports. Both downhole and surface-based V S 30 estimates are represented in our V S 30 compilation. Most of the V S 30 data are for sites in the western contiguous United States (2,141 sites), whereas 786 V S 30 values are for sites in the Central and Eastern United States; 70 values are for sites in other parts of the United States, including Alaska (15 sites), Hawaii (30 sites), and Puerto Rico (25 sites). An interactive map is hosted on the primary USGS Web site for accessing V S 30 data ( https://earthquake.usgs.gov/data/vs30/us/ ).

Data Series↗

Spatial mapping and attribution of Wyoming wind turbines, 2012

These data represent locations of wind turbines found within Wyoming as of August 2012. We assigned each wind turbine to a wind farm and, in these data, provide information about each turbine’s potential megawatt output, rotor diameter, hub height, rotor height, the status of the land ownership where the turbine exists, the county each turbine is located in, wind farm power capacity, the number of units currently associated with each wind farm, the wind turbine manufacturer and model, the wind farm developer, the owner of the wind farm, the current purchaser of power from the wind farm, the year the wind farm went online, and the status of its operation. Some of the attributes are estimates based on the information we found via the American Wind Energy Association and other on-line reports. The locations are derived from National Agriculture Imagery Program (2009 and 2012) true color aerial photographs and have a positional accuracy of approximately +/-5 meters. These data will provide a planning tool for wildlife- and habitat-related projects underway at the U.S. Geological Survey’s Fort Collins Science Center and other government and non-government organizations. Specifically, we will use these data to support quantifying disturbances of the landscape as related to wind energy as well as to quantify indirect disturbances to flora and fauna. This data set represents an update to a previous version by O’Donnell and Fancher (2010).

Wyoming↗

Earthquake scenario and probabilistic ground-shaking hazard maps for the Albuquerque-Belen-Santa Fe, New Mexico, corridor

New Mexico's population is concentrated along the corridor that extends from Belen in the south to Española in the north and includes Albuquerque and Santa Fe. The Rio Grande rift, which encompasses the corridor, is a major tectonically, volcanically, and seismically active continental rift in the western U.S. Although only one large earthquake (moment magnitude (M) ≥ 6) has possibly occurred in the New Mexico portion of the rift since 1849, paleoseismic data indicate that prehistoric surface-faulting earthquakes of M 6.5 and greater have occurred on aver- age every 400 yrs on many faults throughout the Rio Grande rift. We have developed a series of nine scenario and probabilistic hazard maps that portray the ground shaking that could occur in the Albuquerque-Belen-Santa Fe corridor from future earthquakes in New Mexico. These maps, at a scale of 1:500,000, display color-contoured ground-motion values in terms of the parameters of peak horizontal acceleration and horizontal spectral accelerations at 0.2 and 1.0 second (sec) periods. The maps depict surficial ground shaking and incorporate the site-response effects at locations underlain by unconsolidated sediments. The scenario maps are for a M 7.0 earthquake rupturing the Sandia-Rincon faults, which are adjacent to and dip west beneath Albuquerque. The probabilistic maps are for the two annual exceedance probabilities of building code relevance, 10% and 2% exceedance probabilities in 50 yrs (corresponding to return periods of 500 and 2,500 yrs, respectively). We included 57 Quaternary faults, all located within the Rio Grande rift, in the probabilistic seismic hazard analysis. These faults were characterized in terms of their geometry, rupture behavior (including possible segmentation), maximum expected earthquake magnitude, recurrence model, probability of activity, and slip rate. Preferred maximum magnitude values for these faults ranged from M 6.1 to 7.4 and preferred slip rates from 0.01 to approximately 0.12 mm/yr. Regional source zones and Gaussian smoothing of the historical seismicity were also included in the probabilistic hazard analysis to account for the hazard from background earthquakes (M ≤ 6.5). A numerical ground-motion modeling approach and empirical attenuation relation- ships appropriate for extensional tectonic regimes were used to compute the scenario earthquake and probabilistic ground motions on rock. Amplification factors were then used to modify the rock motions and hence to incorporate site response into the hazard maps. These factors were based on three generalized geologic site-response categories (hard rock, soft rock, and firm/stiff soil) and were adopted from similar California-based categories because insufficient subsurface geologic and geotechnical data are available for the map area. The resulting hazard maps indicate that from both scenario and probabilistic perspectives, the ground-shaking hazard in the Albuquerque–Belen–Santa Fe corridor from future earthquakes could be severe, damaging, and potentially disastrous. In the event of a M 7.0 earthquake occurring on the Sandia–Rincon faults, ground shaking as characterized by peak ground acceleration could reach 0.7 g in much of the eastern half of the Albuquerque metropolitan area. (1 g = 980 cm/sec, the rate of gravitational acceleration.) These high ground motions will be attributable to the city’s location directly over the Sandia–Rincon faults and the amplifying effect of the unconsolidated sediments within the Albuquerque Basin. These levels of ground shaking will probably result in severe damage to traditional adobe construction and even to modern buildings. Long- period ground motions (> 1.0 sec), which are significant to long and tall structures (e.g., tall buildings, long bridges, and highway overpasses), will also be high (> 1.0 g). Injuries and loss of life will be likely. For the 500- and 2,500-yr return period maps, the highest peak accelerations are predicted to be at the damaging levels of 0.3 g and 0.6 g, respectively. All maps show dramatically the frequency-dependent amplification of unconsolidated sediments in the basins along the Rio Grande valley (e.g., Albuquerque Basin). The pattern of amplification and deamplification is clearly a function of the distribution of unconsolidated sediments. These maps are not intended to be a substitute for site-specific studies for engineering design nor to replace standard maps commonly referenced in building codes. Rather, we hope that these maps will be used as a guide by government agencies; the engineering, urban planning, emergency preparedness, and response communities; and the general public as part of an overall program to reduce earthquake risk and losses in New Mexico.

New Mexico↗