SEARCH · Geology Reports

Results for “Algorithms”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,279 records · Page 71Linked to original sources

USGS1994

A comparison of change detection methods using multispectral scanner data

Change detection methods were investigated as a cooperative activity between the U.S. Geological Survey and the National Bureau of Surveying and Mapping, People's Republic of China. Subtraction of band 2, band 3, normalized difference vegetation index, and tasseled cap bands 1 and 2 data from two multispectral scanner images were tested using two sites in the United States and one in the People's Republic of China. A new statistical method also was tested. Band 2 subtraction gives the best results for detecting change from vegetative cover to urban development. The statistical method identifies areas that have changed and uses a fast classification algorithm to classify the original data of the changed areas by land cover type present for each image date.

Open-File Report↗

USGS1996

Catalog of earthquake hypocenters for Augustine, Redoubt, Iliamna, and Mount Spurr volcanoes, Alaska: January 1, 1991 - December 31, 1993

The Alaska Volcano Observatory (AVO), a cooperative program of the U.S. Geological Survey, the Geophysical Institute of the University of Alaska, Fairbanks, and the Alaska Division of Geological and Geophysical Surveys, has maintained a program of seismic monitoring at potentially active volcanoes in the Cook Inlet region since 1988. The principal objectives of this program include the seismic surveillance of the Cook Inlet volcanoes and the investigation of seismic processes associated with active volcanism. This catalog reflects the status and evolution of the seismic monitoring program, and presents the basic seismic data for the time interval January 1, 1991, to December 31, 1993. For an interpretation of these data the reader should refer to several recent articles on volcano related seismicity in the Cook Inlet region (e.g. Jolly and others, 1994; Power and others, 1995; and McNutt and others, 1995). A similar catalog covers the period from October 12, 1989 to December 31, 1991 (Power and others 1993). The AVO seismic monitoring program has undergone significant changes during the catalog period. The changes included 1) new seismic stations placed at Mount Spurr and Redoubt Volcano, resulting in increased earthquake detection capability and improved earthquake locations, 2) the addition of several regional stations to the seismic data acquisition system which improved location quality near the volcano and enhanced our ability to scale eruptions and measure magnitudes of the largest volcanic earthquakes, 3) installation of a new event detection algorithm XDETECT (Rogers, 1993), which offered increased data collection capabilities , 4) modifications to the earthquake location program HYPOELLIPSE (Lahr, 1989) which now allows distinct velocity models and station corrections at each volcanic center, and 5) the addition of seismic stations at Augustine and niamna volcanoes to the data acquisition/location system. The 1992 eruptions at Mount Spurr's Crater Peak vent provided the highlight of the catalog period. The crisis included three sub-plinian eruptions, which occurred on June 27, August 18, and September 16-17, 1992. The three eruptions punctuated a complex seismic sequence which included volcano-tectonic (VT) earthquakes, tremor, and both deep and shallow long period (LP) earthquakes. The seismic sequence began on August 18, 1991, with a small swarm of volcano-tectonic events beneath Crater Peak, and spread throughout the volcanic complex by November of the same year. Elevated levels of seismicity persisted at Mount Spurr beyond the catalog time period.

Alaska↗

USGS1997

A manual for a laboratory information management system (LIMS) for light stable isotopes

The reliability and accuracy of isotopic data can be improved by utilizing database software to (i) store information about samples, (ii) store the results of mass spectrometric isotope-ratio analyses of samples, (iii) calculate analytical results using standardized algorithms stored in a database, (iv) normalize stable isotopic data to international scales using isotopic reference materials, and (v) generate multi-sheet paper templates for convenient sample loading of automated mass-spectrometer sample preparation manifolds. Such a database program is presented herein. Major benefits of this system include (i) an increase in laboratory efficiency, (ii) reduction in the use of paper, (iii) reduction in workload due to the elimination or reduction of retyping of data by laboratory personnel, and (iv) decreased errors in data reported to sample submitters. Such a database provides a complete record of when and how often laboratory reference materials have been analyzed and provides a record of what correction factors have been used through time. It provides an audit trail for stable isotope laboratories. Since the original publication of the manual for LIMS for Light Stable Isotopes, the isotopes 3 H, 3 He, and 14 C, and the chlorofluorocarbons (CFCs), CFC-11, CFC-12, and CFC-113, have been added to this program.

Open-File Report↗

USGS1998

A manual for a Laboratory Information Management System (LIMS) for light stable isotopes

Open-File Report↗

USGS1999

SutraPlot, a graphical post-processor for SUTRA, a model for ground-water flow with solute or energy transport

This report documents a graphical display post-processor (SutraPlot) for the U.S. Geological Survey Saturated-Unsaturated flow and solute or energy TRAnsport simulation model SUTRA, Version 2D3D.1. This version of SutraPlot is an upgrade to SutraPlot for the 2D-only SUTRA model (Souza, 1987). It has been modified to add 3D functionality, a graphical user interface (GUI), and enhanced graphic output options. Graphical options for 2D SUTRA (2-dimension) simulations include: drawing the 2D finite-element mesh, mesh boundary, and velocity vectors; plots of contours for pressure, saturation, concentration, and temperature within the model region; 2D finite-element based gridding and interpolation; and 2D gridded data export files. Graphical options for 3D SUTRA (3-dimension) simulations include: drawing the 3D finite-element mesh; plots of contours for pressure, saturation, concentration, and temperature in 2D sections of the 3D model region; 3D finite-element based gridding and interpolation; drawing selected regions of velocity vectors (projected on principal coordinate planes); and 3D gridded data export files. Installation instructions and a description of all graphic options are presented. A sample SUTRA problem is described and three step-by-step SutraPlot applications are provided. In addition, the methodology and numerical algorithms for the 2D and 3D finite-element based gridding and interpolation, developed for SutraPlot, are described. 1

Open-File Report↗

USGS1982

Digital classification of Landsat data for vegetation and land-cover mapping in the Blackfoot River watershed, southeastern Idaho

This paper documents the procedures, results, and final products of a digital analysis of Landsat data used to produce a vegetation and landcover map of the Blackfoot River watershed in southeastern Idaho. Resource classes were identified at two levels of detail: generalized Level I classes (for example, forest land and wetland) and detailed Levels II and III classes (for example, conifer forest, aspen, wet meadow, and riparian hardwoods). Training set statistics were derived using a modified clustering approach. Environmental stratification that separated uplands from lowlands improved discrimination between resource classes having similar spectral signatures. Digital classification was performed using a maximum likelihood algorithm. Classification accuracy was determined on a single-pixel basis from a random sample of 25-pixel blocks. These blocks were transferred to small-scale color-infrared aerial photographs, and the image area corresponding to each pixel was interpreted. Classification accuracy, expressed as percent agreement of digital classification and photo-interpretation results, was 83.0:t 2.1 percent (0.95 probability level) for generalized (Level I) classes and 52.2:t 2.8 percent (0.95 probability level) for detailed (Levels II and III) classes. After the classified images were geometrically corrected, two types of maps were produced of Level I and Levels II and III resource classes: color-coded maps at a 1:250,000 scale, and flatbed-plotter overlays at a 1:24,000 scale. The overlays are more useful because of their larger scale, familiar format to users, and compatibility with other types of topographic and thematic maps of the same scale.

Professional Paper↗

USGS2006

Mapping ground water in three dimensions: An analysis of airborne geophysical surveys of the Upper San Pedro River Basin, Cochise County, southeastern Arizona

This report summarizes the results of two airborne geophysical surveys conducted in the upper San Pedro Valley of southeastern Arizona in 1997 and 1999. The combined surveys cover about 1,000 square kilometers and extend from the Huachuca Mountains on the west to the Mule Mountains and Tombstone Hills on the east and from north of the Babocomari River to near the Mexican border on the south. The surveys included the acquisition of high-resolution magnetic data, which were used to map depth to the crystalline basement rocks underlying the sediments filling the basin. The magnetic inversion results show a complex basement morphology, with sediment thickness in the center of the valley ranging from ~237 meters beneath the city of Sierra Vista to ~1,500 meters beneath Huachuca City and the Palominas area near the Mexican border. The surveys also included acquisition of 60-channel time-domain electromagnetic (EM) data. Extensive quality analyses of these data, including inversion to conductivity vs. depth (conductivity-depth-transform or CDT) profiles and comparisons with electrical well logs, show that the electrical conductor mapped represents the subsurface water-bearing sediments throughout most of the basin. In a few places (notably the mouth of Huachuca Canyon), the reported water table lies above where the electrical conductor places it. These exceptions appear to be due to a combination of outdated water-table information, significant horizontal displacement between the wells and the CDT profiles, and a subtle calibration issue with the CDT algorithm apparent only in areas of highly resistive (very dry) overburden. These occasional disparities appear in less than 5 percent of the surveyed area. Observations show, however, that wells drilled in the thick unsaturated zone along the Huachuca Mountain front eventually intersect water, at which point the water rapidly rises high into the unsaturated zone within the wellbore. This rising of water in a wellbore implies some sort of confinement below the thick unsaturated zone, a confinement that is not identified in the available literature. Occasional disparities notwithstanding, maps of the electrical conductor derived from the airborne EM system provide a synoptic view of the presence of water underlying the upper San Pedro Valley, including its three-dimensional distribution. The EM data even show faults previously only inferred from geologic mapping. The magnetic and electromagnetic data together appear to show the thickness of the sediments, the water in the saturated sediments down to a maximum of about 400 meters depth, and even places where the main ground-water body is not in direct contact with the San Pedro River. However, the geophysical data cannot reveal anything directly about hydraulic conductivity or ground-water flow. Estimating these characteristics requires new hydraulic modeling based in part on this report. One concern to reviewers of this report is the effect that clays may have on the electrical conductor mapped with the airborne geophysical system. Although the water in the basin is unusually conductive, averaging 338 microsiemens per centimeter, reasoning cited below suggests that the contribution of clays to the overall conductivity would be relatively small. Basic principles of sedimentary geology suggest that silts and clays should dominate the center of the basin, while sands and gravels would tend to dominate the margins. Although clay content may increase the amplitude of the observed electrical conductors somewhat, it will not affect the depths to the conductor derived from depth inversions. Further, fine-grained sediments generally have higher porosity and tend to lie toward a basin center, a fact in general agreement with the observed geophysical data.

Arizona↗

USGS2008

Evolving magma storage conditions beneath Mount St. Helens inferred from chemical variations in melt inclusions from the 1980-1986 and current (2004-2006) eruptions

Major element, trace element, and volatile concentrations in 187 glassy melt inclusions and 25 groundmass glasses from the 1980-86 eruption of Mount St. Helens are presented, together with 103 analyses of touching FE-Ti oxide pairs from the same samples. These data are used to evaluate the temporal evolution of the magmatic plumbing system beneath the volcano during 1980-86 and so provide a framework in which to interpret analyses of melt inclusions from the current (2004-2006) eruption. Major and trace element concentrations of all melt inclusions lie at the high SiO 2 end of the data array defined by eruptive products of the late Quaternary age from Mount St. Helens. For several major and trace elements, the glasses define a trend that is oblique to the whole-rock trend, indicating that different mineral assemblages were responsible for the two trends. The whole-rock trend can be ascribed to differentiation of hydrous basaltic parents in a deep-seated magma reservoir, probably at depths great enough to stabilize garnet. In contrast, the glass trends were generated by closed-system crystallization of the phenocryst and microlite mineral assemblages at low pressures. The dissolved H 2 O content of the melt inclusions from 1980-86, as measured by the ion microprobe, ranges from 0 to 6.7 wt. percent, with the highest values obtained from the plinian phase of May 18, 1980. Water contents decrease with increasing SiO 2 , consistent with decompression-driven crystallization. Preliminary data for dissolved CO 2 in melt inclusions from the May 18 plinian phase from August 7, 1980, indicate that X H 2 O in a vapor phase was approximately constant at 0.80, irrespective of H 2 O content, suggestive of closed-system degassing with a high bubble fraction or gas streaming through the subvolcanic system. Temperature and f O 2 estimates for touching Fe-Ti oxides show evidence for heating during crystallization owing to release of latent heat. Consequently, magmas with the highest microlite crystallinities record the highest temperatures. Magmas also become progressively reduced during ascent and degassing, probably as a result of redox equilibria between exsolving S-bearing gases and magmas. The lowest temperature oxides have f O 2 ≈ NNO, similar to high-temperature fumarole gases from the volcano. The temperature and f O 2 of the magma tapped by the plinian phase of May 18, 1980, are 870-875°C and NNO+0.8, respectively. The dissolved volatile contents of the melt inclusions have been used to calculate sealing pressures; that is, the pressure at which chemical exchange between inclusion and matrix melt ceased. These are greatest for the May 18 plinian magma (120 to 320 MPa); lower pressures are recorded by samples of the preplinian cryptodome and by all post-May 18 magmas. Magma crystallinity, calculated from melt-inclusion Rb contents, is negatively correlated with sealing pressure, consistent with decompression crystallization. Elevated contents of Li in melt inclusions from the cryptodome and post-May 18 samples are consistent with transfer of Li in a magmatic vapor phase from deeper parts of the magma system to magma stored at shallower levels. The Li enrichment attains its maximum extent at ~150 MPa, which is ascribed to separation of a single vapor phase into H 2 O-rich gas and dense Li-rich brine at the top of the magma column. There are striking correlations between melt-inclusion chemistry and monitoring data for the 1980-86 eruption. Dissolved SO 2 contents of melt inclusions from any given event, multiplied by the mass of magma erupted during that event, correlate with the measured flux of SO 2 at the surface, suggesting that magma degassing and melt-inclusion sealing are closely related in time and space. Textural and chemical evidence indicates that melt inclusions became effectively sealed (physically or kinetically) shortly before eruption. Thus by converting pressure to depth using a density model and edifice-loading algorithm for the volcano, changing depths of magma extraction with time can be tracked and compared to the seismic record. The plinian eruption of May 18, 1980, involved magma stored 5-11 km below sea level; this is inferred to be the subvolcanic magma chamber. The preceding eruptions, including the May 18, 1980, blast, involved magma withdrawal from the cryptodome and conduit down to 5 km below sea level. Subsequent 1980 eruptions tapped magma down to depths of ≤10 km below sea level. Tapping of magma stored deeper than 2 km below sea level stopped abruptly at the end of 1980, coincident with the onset of extensive shallow seismicity and a change from explosive to effusive eruption style from 1981 to 1986. Overall, the 1980-86 eruption is consistent with the evisceration of a thin, vertically extensive body of magma extending from 5 to at least 11 km below sea level and connected to the surface by a thin conduit. In the absence of sustained high magma-supply rates from depth, decompression crystallization of magma ascending through the system leads eventually to plugging of the conduit. The current eruption of Mount St. Helens shares some similarities with the 1981-86 dome-building phase of the previous eruption, in that there is extensive shallow seismicity and extrusion of highly crystalline material in the form of a sequence of flows and spines. Melt inclusions from the current eruption have low H 2 O contents, consistent with magma extraction from shallow depths. Highly enriched Li in melt inclusions suggests that vapor transport of Li is a characteristic feature of Mount St. Helens. Melt inclusions from the current eruption have subtly different trace-element chemistry from all but one of the 1980-86 melt inclusions, with steeper rareearth-element (REE) patterns and low U, Th, and high-fieldstrength elements (HFSE), indicating addition of a new melt component to the magma system. It is anticipated that increasing involvement of the new melt component will be evident as the current eruption proceeds.

Washington↗

USGS2008

Nutrient Concentrations and Their Relations to the Biotic Integrity of Nonwadeable Rivers in Wisconsin

Excessive nutrient [phosphorus (P) and nitrogen (N)] input from point and nonpoint sources is frequently associated with degraded water quality in streams and rivers. Point-source discharges of nutrients are fairly constant and are controlled by the U.S. Environmental Protection Agency's (USEPA) National Pollutant Discharge Elimination System. To reduce inputs from nonpoint sources, agricultural performance standards and regulations for croplands and livestock operations are being proposed by various States. In addition, the USEPA is establishing regionally based nutrient criteria that can be refined by each State to determine whether actions are needed to improve water quality. More confidence in the environmental benefits of the proposed performance standards and nutrient criteria would be possible with improved understanding of the biotic responses to a range of nutrient concentrations in different environmental settings. To achieve this general goal, the U.S. Geological Survey and the Wisconsin Department of Natural Resources collected data from 282 streams and rivers throughout Wisconsin during 2001 through 2003 to: (1) describe how nutrient concentrations and biotic-community structure differ throughout the State, (2) determine which environmental characteristics are most strongly related to the distribution of nutrient concentrations and biotic-community structure, (3) determine reference conditions for water quality and biotic indices for streams and rivers in the State, (4) determine how the biotic communities in streams and rivers in different areas of the State respond to differences in nutrient concentrations, (5) determine the best regionalization scheme to describe the patterns in reference conditions and the corresponding responses in water quality and the biotic communities (primarily for smaller streams), and (6) develop algorithms to estimate nutrient concentrations in streams and rivers from a combination of biotic indices. The ultimate goal of this study is to provide the information needed to guide the development of regionally based nutrient criteria for Wisconsin streams and rivers. In this report, data collected, primarily in 2003, from 42 nonwadeable rivers are used to describe nutrient concentrations and their relations to the biotic integrity of rivers in Wisconsin. In a separate report by Robertson and others (2006a), the data collected from 240 wadeable streams are used to describe these relations in streams in Wisconsin. Reference water-quality conditions for nonwadeable rivers were found to be similar throughout Wisconsin (approximately 0.035 milligrams per liter (mg/L) for total P (TP), 0.500 mg/L for total N (TN), 4 micrograms per liter for suspended chlorophyll a (SCHL), and greater than 110 centimeters for Secchi-tube depth (SD)). For each category of the biotic community (SCHL, macroinvertebrates, and fish), a few indices were more strongly related to differences in nutrient concentrations than were others. For the indices most strongly related to nutrient concentrations, reference conditions were obtained with a regression approach, from values corresponding to the worst 75th-percentile value from a subset of minimally impacted streams (streams having reference nutrient concentrations), and from the best 25th-percentile value of all the data. Concentrations of TP and TN in nonwadeable rivers increased as the percentage of agricultural land in the basin increased; these increases resulted in increased SCHL concentrations and decreased SDs. The responses in SDs and SCHL concentrations to changes in nutrient concentrations were similar throughout most of the State except in rivers in the southeastern part, where SCHL concentrations were lower than would be expected given their nutrient concentrations. Rivers in the southeastern part of the State had high concentrations of total suspended sediment compared to the SCHL concentrations. Many biotic indices responded to increases in nu

Professional Paper↗

USGS2010

Integrated satellite observations of the 2006 eruption of Augustine Volcano: Chapter 20 in The 2006 eruption of Augustine Volcano, Alaska

Satellite observations played an important role in monitoring the 2006 eruption of Augustine Volcano. It represented the first opportunity for observers to use, in an operational setting, new Web-based tools and techniques developed by the Alaska Volcano Observatory remote sensing group. The 'Okmok Algorithm' was used to analyze thermal infrared satellite data and highlight changes in the style and phases of activity. Temperature measurements were used to estimate ash cloud heights, which compared favorably to radar and ground-based observations, although larger discrepancies were seen when compared to pilot reports. Brightness temperature difference techniques were used to locate and track the 14 ash clouds produced during the explosive phase of the eruption. Stacking of these analyses allowed the creation of composite maps showing the distribution of airborne ash. The data from these maps were further combined with information from local reports and samples of ashfall to create a prototype of a concentration map that could be used to assess the potential hazard an eruption represents to aircraft, infrastructure and human health.

Alaska↗

USGS2023

Magnetotelluric investigation of northern Harrat Rahat, Kingdom of Saudi Arabia

Volcanism within the harrats (Arabic for “volcanic field”) of the Kingdom of Saudi Arabia includes at least one historical eruption occurring close to the holy city of Al Madīnah in 1256 C.E. As part of a volcanic- and seismic-hazard assessment of northern Harrat Rahat, magnetotelluric (MT) data were collected to investigate the structural setting of the area, the presence or absence of melt within the crust, and the mantle-derived magmatic source. Collected MT data were modeled in both two dimensions, where anisotropy can be estimated, and three dimensions. Interpretation of the preferred resistivity model includes a shallow sediment-filled graben beneath northern Harrat Rahat lavas, a melt-free upper crust, and a region of decompression melting in the asthenosphere below 60–70 kilometers depth. Models in two dimensions image the lower crust as anisotropic, demonstrating that a series of elongate conductivity anomalies with a strike of N. 10° E. within the lower crust of the three-dimensional model are artifacts of inverting anisotropic data with an isotropic modeling algorithm. Careful examination of the resistivity models, in combination with regional geological and geophysical data, suggests an anisotropic lower crust that is free of large zones of melt. Azimuthal anisotropy in the lower crust extends well beyond the limits of Harrat Rahat volcanic rocks, with a conductive direction oriented N. 10° E. and an anisotropy factor of 2–5 between the most and least conductive directions. Enhanced conductivity is likely caused by interconnected grain-boundary graphite, where the direction of anisotropy reflects either frozen-in fabric from the Neoproterozoic stabilization of the Arabian Shield or ductile deformation driven by channelized asthenospheric flow coupled with a thin rigid mantle lid. Asthenospheric melt is interpreted to transect the crust largely through diking, with limited melt storage and short residence times within the crustal column.

northern Harrat Rahat↗

USGS2023

Thickness of the Saudi Arabian crust

As part of a joint Saudi Geological Survey (SGS) and U.S. Geological Survey (USGS) project, we analyzed P-wave receiver functions from seismic stations covering most of the Kingdom of Saudi Arabia to map the thickness of the crust across the Arabia Plate. We present an update of crustal-thickness estimates and fill in gaps for the western Arabian Shield and the rifted margin at the Red Sea (the coastal plain), as well as the eastern Arabian Platform. We applied a conventional H-k stacking algorithm and included careful attention to stacking weights, two forms of sedimentary corrections for stations located on the Arabian Platform, and additional processing for noisy stations. We obtained useful results at 154 stations from 898 teleseismic events over a 2-year period from 1995–1997 (for non-SGS stations) and a 6-year period from 2008–2014 (for SGS stations). Average crustal thickness (that is, depth to the Mohorovičić discontinuity [Moho] below the surface) beneath the Red Sea coastal plain (the rift margin) is 29 kilometers (km), beneath the volcanic fields (known in Arabic as harra [plural] or harrat [singular]) is 35 km, beneath the Arabian Shield (excluding harrats) is 37 km, and beneath the Arabian Platform is 38 km. Crustal thinning appears not to extend east of the rift escarpment, suggesting uniform extension that is no broader at depth than at the surface. In contrast to some previous claims that the Arabian Platform crust is thicker than that of the Arabian Shield, we find no statistically significant difference between their whole crustal thicknesses. However, the average subsedimentary crustal thickness (that is, the crystalline crust) for stations on the Arabian Platform is 34 km, 3 km thinner than the crust of the Arabian Shield. Individual station P-wave (pressure) velocity and S-wave (shear) velocity ratios ( V P / V S ) are highly variable for the Arabia Plate, ranging from 1.60 to 1.97 and averaging 1.75, with a standard deviation of 0.07. There are no statistically significant differences between V P / V S ratios of the different geologic regions of Saudi Arabia. Similar V P / V S ratios, coupled with similar crustal thicknesses for harrats and the Arabian Shield, indicate that Cenozoic magmatism has contributed negligibly to crustal growth.

Professional Paper↗

USGS2023

Ambient seismic noise tomography of the Kingdom of Saudi Arabia

Harrat Rahat is a Cenozoic volcanic field in the west-central part of the Kingdom of Saudi Arabia, 150 kilometers east of the Red Sea, and is the site of the most recent eruption in the country (1256 C.E.; 654 in the year of the Hijra). The city of Al Madīnah lies at the north end of Harrat Rahat, and its volcanic and seismic risks are frequently reassessed. In 2009 C.E. an earthquake swarm at Harrat Lunayyir, 200 km west-northwest of Al Madīnah, produced significant shaking and some building damage in nearby communities, which prompted a revision of seismic hazard models for the region. A network of seismic stations on this volcanically active western side of the Arabia Plate was installed, and stations were also added in the tectonically active northern part of the country. Although regional earthquakes may be used to determine the crustal structure of the western Arabia Plate, such crustal models are hindered by insufficient numbers of earthquakes in the stable plate interior. Tomography studies can be used to infer material properties of the subsurface, such as presence of partial melt, and are beneficial for volcanic hazard assessments. We use ambient seismic noise to compute Rayleigh and Love surface-wave dispersion maps between 5 and 12 second periods for a subset of seismic stations within and near northern Harrat Rahat. The surface-wave maps are inverted to produce shear-wave velocities using a neighborhood algorithm and interpolated into a pseudo three-dimensional model. The distributions of surface-wave and shear-wave velocities are heterogenous, varying from ±3–8 percent. However, low velocities are not restricted to Harrat Rahat. We observed a difference between Rayleigh- and Love-wave velocities that extends north of the site of the 1256 C.E. eruption and coincides with a low gravity anomaly. We obtain a shear-wave velocity increase of 10–15 percent between 15- and 25-kilometers depth, which is consistent with the presence of a transition between the felsic upper crust and the mafic lower crust of the Arabian Shield. The average shear-wave velocities of the upper and lower crust are estimated to be 3.64 and 3.95 kilometers per second using Rayleigh waves and 3.53 and 4.16 kilometers per second using Love waves, which are in good agreement with the results of other geophysical surveys in this area. The modest magnitude of the low-velocity anomalies within the crust and their locations extending well beyond the limits of Harrat Rahat indicate that they are not caused by a crustal magma chamber. If magma chambers exist, they are smaller than can be imaged with our seismological method (resolution corresponds to a 15-kilometer wavelength), deeper than 30 kilometers, or shallower than 5 kilometers with a small velocity contrast. The last possibility is discounted by weak-to-absent surface geothermal activity and by a lack of shallow seismicity that characterizes other areas with known shallow magmas, such as Hawaiʻi. We then expanded our work to the entire Arabian Shield with principally the same methodology to look at large-scale regional patterns. We found modest shear-wave velocity deviations on the order of only ±3 percent, which are within expected ranges for lithological variation.

Professional Paper↗

USGS2012

Maps of estimated nitrate and arsenic concentrations in basin-fill aquifers of the southwestern United States

Human-health concerns and economic considerations associated with meeting drinking-water standards motivated a study of the vulnerability of basin-fill aquifers to nitrate contamination and arsenic enrichment in the southwestern United States. Statistical models were developed by using the random forest classifier algorithm to predict concentrations of nitrate and arsenic across a model grid representing about 190,600 square miles of basin-fill aquifers in parts of Arizona, California, Colorado, Nevada, New Mexico, and Utah. The statistical models, referred to as classifiers, reflect natural and human-related factors that affect aquifer vulnerability to contamination and relate nitrate and arsenic concentrations to explanatory variables representing local- and basin-scale measures of source and aquifer susceptibility conditions. Geochemical variables were not used in concentration predictions because they were not available for the entire study area. The models were calibrated to assess model accuracy on the basis of measured values. Only 2 percent of the area underlain by basin-fill aquifers in the study area was predicted to equal or exceed the U.S. Environmental Protection Agency drinking-water standard for nitrate as N (10 milligrams per liter), whereas 43 percent of the area was predicted to equal or exceed the standard for arsenic (10 micrograms per liter). Areas predicted to equal or exceed the drinking-water standard for nitrate include basins in central Arizona near Phoenix; the San Joaquin Valley, the Santa Ana Inland, and San Jacinto Basins of California; and the San Luis Valley of Colorado. Much of the area predicted to equal or exceed the drinking-water standard for arsenic is within a belt of basins along the western portion of the Basin and Range Physiographic Province that includes almost all of Nevada and parts of California and Arizona. Predicted nitrate and arsenic concentrations are substantially lower than the drinking-water standards in much of the study area-about 93 percent of the area underlain by basin-fill aquifers was less than one-half the standard for nitrate as N (5.0 milligrams per liter), and 50 percent was less than one-half the standard for arsenic (5.0 micrograms per liter). The predicted concentrations and the improved understanding of the susceptibility and vulnerability of southwestern basin-fill aquifers to nitrate contamination and arsenic enrichment can be used by water managers as a qualitative tool to assess and protect the quality of groundwater resources in the Southwest.

Arizona, California, Colorado, Idaho, Nevada, New ↗

USGS2018

Isostatic gravity map of Mountain Pass and vicinity, California and Nevada

Gravity investigations of Mountain Pass and vicinity were begun as part of an effort to study regional crustal structures as an aid to understanding the geologic framework and mineral resources of the eastern Mojave Desert. The study area, which straddles the state boundary between southeastern California and southern Nevada, encompasses Mountain Pass, which is host to one of the world’s largest rare earth element carbonatite deposits. The deposit is found along a north-northwest-trending, fault-bounded block that extends along the eastern parts of the Clark Mountain Range, Mescal Range, and Ivanpah Mountains. This Paleoproterozoic block is composed of a 1.7-Ga metamorphic complex of gneiss and schist that underwent widespread metamorphism and associated plutonism during the Ivanpah orogeny. The Paleoproterozoic rocks were intruded by a Mesoproterozoic (1.4 Ga) ultrapotassic alkaline intrusive suite and carbonatite body. The intrusive rocks include, from oldest to youngest, shonkinite, mesosyenite, syenite, quartz syenite, potassic granite, carbonatite, carbonatite dikes, and late shonkinite dikes. Generally speaking, gravity anomalies can be used to infer subsurface geologic structure, revealing variations in lithology and delineating features such as faults, plutons, volcanic centers, calderas, and deep sedimentary basins. As part of this study, gravity data from more than 2,400 stations were collected and processed to identify lateral changes in subsurface density. Gravity stations were distributed across parts of Shadow Valley, Clark Mountain Range, Mescal Range, Ivanpah Mountains, and Ivanpah Valley. The new gravity data were combined with preexisting gravity data from the surrounding areas in California and Nevada. All gravity data were gridded using a minimum curvature algorithm at an interval of 200 m, and the result is displayed as a color-contour isostatic gravity map.

California, Nevada↗

USGS2005

Analysis and mapping of post-fire hydrologic hazards for the 2002 Hayman, Coal Seam, and Missionary Ridge wildfires, Colorado

Wildfires caused extreme changes in the hydrologic, hydraulic, and geomorphologic characteristics of many Colorado drainage basins in the summer of 2002. Detailed assessments were made of the short-term effects of three wildfires on burned and adjacent unburned parts of drainage basins. These were the Hayman, Coal Seam, and Missionary Ridge wildfires. Longer term runoff characteristics that reflect post-fire drainage basin recovery expected to develop over a period of several years also were analyzed for two affected stream reaches: the South Platte River between Deckers and Trumbull, and Mitchell Creek in Glenwood Springs. The 10-, 50-, 100-, and 500-year flood-plain boundaries and water-surface profiles were computed in a detailed hydraulic study of the Deckers-to-Trumbull reach. The Hayman wildfire burned approximately 138,000 acres (216 square miles) in granitic terrain near Denver, and the predominant potential hazard in this area is flooding by sediment-laden water along the large tributaries to and the main stem of the South Platte River. The Coal Seam wildfire burned approximately 12,200 acres (19.1 square miles) near Glenwood Springs, and the Missionary Ridge wildfire burned approximately 70,500 acres (110 square miles) near Durango, both in areas underlain by marine shales where the predominant potential hazard is debris-flow inundation of low-lying areas. Hydrographs and peak discharges for pre-burn and post-burn scenarios were computed for each drainage basin and tributary subbasin by using rainfall-runoff models because streamflow data for most tributary subbasins were not available. An objective rainfall-runoff model calibration method based on nonlinear regression and referred to as the ?objective calibration method? was developed and applied to rainfall-runoff models for three burned areas. The HEC-1 rainfall-runoff model was used to simulate the pre-burn rainfall-runoff processes in response to the 100-year storm, and HEC-HMS was used for runoff hydrograph generation. Post-burn rainfall-runoff parameters were determined by adjusting the runoff-curve numbers on the basis of a weighting procedure derived from the U.S. Soil Conservation Service (now the National Resources Conservation Service) equation for precipitation excess and the effect of burn severity. This weighting procedure was determined to be more appropriate than simple area weighting because of the potentially marked effect of even small burned areas on the runoff hydrograph in individual drainage basins. Computed water-peak discharges from HEC-HMS models were increased volumetrically to account for increased sediment concentrations that are expected as a result of accelerated erosion after burning. Peak discharge estimates for potential floods in the South Platte River were increased by a factor that assumed a volumetric sediment concentration (Cv) of 20 percent. Flood hydrographs for the South Platte River and Mitchell Creek were routed down main-stem channels using watershed-routing algorithms included in the HEC-HMS rainfall-runoff model. In areas subject to debris flows in the Coal Seam and Missionary Ridge burned areas, debris-flow discharges were simulated by 100-year rainfall events, and the inflow hydrographs at tributary mouths were simulated by using the objective calibration method. Sediment concentrations (Cv) used in debris-flow simulations were varied through the event, and were initial Cv 20 percent, mean Cv approximately 31 percent, maximum Cv 48 percent, Cv 43 percent at the time of the water hydrograph peak, and Cv 20 percent for the duration of the event. The FLO-2D flood- and debris-flow routing model was used to delineate the area of unconfined debris-flow inundation on selected alluvial fan and valley floor areas. A method was developed to objectively determine the post-fire recovery period for the Hayman and Coal Seam burned areas using runoff-curve numbers (RCN) for all drainage basins for a 50-year period. A

Scientific Investigations Report↗

USGS2006

MODOPTIM: A general optimization program for ground-water flow model calibration and ground-water management with MODFLOW

MODOPTIM is a non-linear ground-water model calibration and management tool that simulates flow with MODFLOW-96 as a subroutine. A weighted sum-of-squares objective function defines optimal solutions for calibration and management problems. Water levels, discharges, water quality, subsidence, and pumping-lift costs are the five direct observation types that can be compared in MODOPTIM. Differences between direct observations of the same type can be compared to fit temporal changes and spatial gradients. Water levels in pumping wells, wellbore storage in the observation wells, and rotational translation of observation wells also can be compared. Negative and positive residuals can be weighted unequally so inequality constraints such as maximum chloride concentrations or minimum water levels can be incorporated in the objective function. Optimization parameters are defined with zones and parameter-weight matrices. Parameter change is estimated iteratively with a quasi-Newton algorithm and is constrained to a user-defined maximum parameter change per iteration. Parameters that are less sensitive than a user-defined threshold are not estimated. MODOPTIM facilitates testing more conceptual models by expediting calibration of each conceptual model. Examples of applying MODOPTIM to aquifer-test analysis, ground-water management, and parameter estimation problems are presented.

Scientific Investigations Report↗

USGS2008

Mobility statistics and automated hazard mapping for debris flows and rock avalanches

Power-law equations that are physically motivated and statistically tested and calibrated provide a basis for forecasting areas likely to be inundated by debris flows, rock avalanches, and lahars with diverse volumes. The equations A=α 1 V 2/3 and B=α 2 V 2/3 are based on the postulate that the maximum valley cross-sectional area (A) and total valley planimetric area (B) likely to be inundated by a flow depend only on its volume (V) and the topography of the flow path. Testing of these equations involves determining whether or not they fit data for documented flows satisfactorily, and calibration entails determining best-fit values of the coefficients α 1 and α 2 for debris flows, rock avalanches, and lahars. This report describes statistical testing and calibration of the equations by using field data compiled from many sources, and it describes application of the equations to delineation of debris-flow hazard zones. Statistical results show that for each type of flow (debris flows, rock avalanches, and lahars), the dependence of A and B on V is described well by power laws with exponents equal to 2/3. This value of the exponent produces fits that are effectively indistinguishable from the best fits obtained by using adjustable power-law exponents. Statistically calibrated values of the coefficients α 1 and α 2 provide scale-invariant indices of the relative mobilities of rock avalanches (α 1 = 0.2, α 2 = 20), nonvolcanic debris flows (α 1 = 0.1, α 2 = 20), and lahars (α 1 = 0.05, α 2 = 200). These values show, for example, that a lahar of specified volume can be expected to inundate a planimetric area ten times larger than that inundated by a rock avalanche or nonvolcanic debris flow of the same volume. The utility of the calibrated debris-flow inundation equations A=0.1V 2/3 and B=20V 2/3 is demonstrated by using them within the GIS program LAHARZ to delineate nested hazard zones for future debris flows in an area bordering the Umpqua River in the south-central Oregon Coast Range. This application requires use of high-resolution topographic data derived form LIDAR surveys, knowledge of local geology to specify a suitable range of prospective debris-flow volumes, and development and use of a new algorithm for identification of prospective debris-flow source areas in finely dissected terrain.

Oregon↗