USGS Science⌕ Search

SEARCH · USGS Science

Results for “Near Surface Geophysics”

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 937 records · Page 52Linked to original sources
USGS2012

Crustal strain near the Big Bend of the San Andreas Fault: Analysis of the Los Padres-Tehachapi Trilateration Networks, California

In the region of the Los Padres-Tehachapi geodetic network, the San Andreas fault (SAF) changes its orientation by over 30° from N40°W, close to that predicted by plate motion for a transform boundary, to N73°W. The strain orientation near the SAF is consistent with right-lateral shear along the fault, with maximum shear rate of 0.38±0.01 μrad/yr at N63°W. In contrast, away from the SAF the strain orientations on both sides of the fault are consistent with the plate motion direction, with maximum shear rate of 0.19±0.01 μrad/yr at N44°W. The strain rate does not drop off rapidly away from the fault, and thus the area is fit by either a broad shear zone below the SAF or a single fault with a relatively deep locking depth. The fit to the line length data is poor for locking depth d less than 25 km. For d of 25 km a buried slip rate of 30 ± 6 mm/yr is estimated. We also estimated buried slip for models that included the Garlock and Big Pine faults, in addition to the SAF. Slip rates on other faults are poorly constrained by the Los Padres-Tehachapi network. The best fitting Garlock fault model had computed left-lateral slip rate of 11±2 mm/yr below 10 km. Buried left-lateral slip of 15±6 mm/yr on the Big Pine fault, within the Western Transverse Ranges, provides significant reduction in line length residuals; however, deformation there may be more complicated than a single vertical fault. A subhorizontal detachment on the southern side of the SAF cannot be well constrained by these data. We investigated the location of the SAF and found that a vertical fault below the surface trace fits the data much better than either a dipping fault or a fault zone located south of the surface trace.

Journal of Geophysical Research Solid Earth↗
USGS2006

Surficial sediment character of the Louisiana offshore continental shelf region: A GIS compilation

The Louisiana coastal zone, comprising the Mississippi River delta plain stretching nearly 400 km from Sabine Pass at the Texas border east to the Chandeleur Islands at the Mississippi border, represents one of North America’s most important coastal ecosystems in terms of natural resources, human infrastructure, and cultural heritage. At the same time, this region has the highest rates of coastal erosion and wetland loss in the Nation due to a complex combination of natural processes and anthropogenic actions over the past century. Comparison of historical maps dating back to 1855 and recent aerial photography show the Louisiana coast undergoing net erosion at highly variable rates. Rates have increased significantly during the past several decades. Earlier published statewide average shoreline erosion rates were >6 m/yr; rates have increased recently to >10 m/yr. The increase is attributable to collective action of storms, rapid subsidence, and pervasive man-made alterations of the rivers and the coast. In response to the dramatic landloss, regional-scale restoration plans are being developed by a partnership of federal and state agencies for the delta plain that have the objectives of maintaining the barrier islands, reducing wetland loss, and enhancing the natural sediment delivery processes. There is growing awareness that the sustainability of coastal Louisiana's natural resources and human infrastructure depends on the successful restoration of natural geologic processes. Critical to the long term success of restoration is scientific understanding of the geologic history and processes of the coastal zone region, including interactions between the rivers, wetlands, coast, and inner shelf. A variety of geophysical studies and mapping of Late Quaternary sedimentary framework and coastal processes by U.S. Geological Survey and other scientists during the past 50 years document that the Louisiana delta plain is the product of a complex history of cyclic delta switching by the Mississippi River and its distributaries over the past ~10,000 years that resulted in laterally overlapping deltaic depocenters. The interactions among riverine, coastal, and inner shelf processes have been superimposed on the Holocene transgression resulting in distinctive landforms and sedimentary sequences. Four Holocene shelf-phase delta complexes have been identified using seismic reflection data and vibracores. Each delta complex is bounded by transgressive surfaces. Following each cycle of deposition and abandonment, the delta lobes undergo regional subsidence and marine reworking that forms transgressive coastal systems and barrier islands. Ultimately, the distal end of each of the abandoned delta lobes is marked by submerged marine sand bodies representing drowned barriers. These sand bodies (e.g. Ship Shoal, Outer Shoal, Trinity Shoal, Tiger Shoal, St. Bernard Shoal) offer the largest volumes and highest quality sand for beach nourishment and shoreline and wetlands restoration. These four large sand shoals on inner continental shelf, representing the reworked remnants of former prograded deltaic headlands that existed on the continental shelf at lower sea level, were generated in the retreat path of the Mississippi River delta plain during the Holocene transgression. Penland and others (1989) have shown these sand bodies represent former shoreline positions associated with lower still stands in sea level. Short periods of rapid relative sea-level rise led to the transgressive submergence of the shorelines which today can be recognized at the -10 m to -20 m isobaths on the Louisiana continental shelf. Trinity Shoal and Ship Shoal represent the -10 m middle-to-late Holocene shoreline trend, whereas Outer Shoal and the St. Bernard Shoals define the -20 m early Holocene shoreline trend (Penland and others, 1989). Collectively, these sand shoals constitute a large volume of high quality sandy sediment potentially suitable for barrier island nourishment and coastal restoration. The USGS has actively supported coastal and wetlands geologic research for the past two decades in partnership with universities (e.g., Louisiana State University, University of New Orleans), state agencies (e.g. Louisiana Geological Survey, Louisiana Department of Natural Resources), and private organizations (Williams and others, 1992a,b; Williams and Cichon, 1993; List and others, 1994). These studies have focused on regional-scale mapping of coastal and wetland change and developing a better understanding of the processes that cause coastal erosion and wetlands loss, particularly the rapid deterioration of Louisiana's barrier islands, estuaries, and wetlands environments. With a better understanding of these processes, the ability to model and predict erosion and wetlands loss will improve. More accurate predictions will, in turn, allow for proper management of coastal resources. Improved predictions will also allow for better assessments of the utility of different restoration alternatives.

Louisiana↗
USGS2015

Strain accumulation across the Prince William Sound asperity, Southcentral Alaska

The surface velocities predicted by the conventional subduction model are compared to velocities measured in a GPS array (surveyed in 1993, 1995, 1997, 2000, and 2004) spanning the Prince William Sound asperity. The observed velocities in the comparison have been corrected to remove the contributions from postseismic (1964 Alaska earthquake) mantle relaxation. Except at the most seaward monument (located on Middleton Island at the seaward edge of the continental shelf, just 50 km landward of the deformation front in the Aleutian Trench), the corrected velocities qualitatively agree with those predicted by an improved, two-dimensional, back slip, subduction model in which the locked megathrust coincides with the plate interface identified by seismic refraction surveys, and the back slip rate is equal to the plate convergence rate. A better fit to the corrected velocities is furnished by either a back slip rate 20% greater than the plate convergence rate or a 30% shallower megathrust. The shallow megathrust in the latter fit may be an artifact of the uniform half-space Earth model used in the inversion. Backslip at the plate convergence rate on the megathrust mapped by refraction surveys would fit the data as well if the rigidity of the underthrust plate was twice that of the overlying plate, a rigidity contrast higher than expected. The anomalous motion at Middleton Island is attributed to continuous slip at near the plate convergence rate on a postulated, listric fault that splays off the megathrust at depth of about 12 km and outcrops on the continental slope south-southeast of Middleton Island.

Journal of Geophysical Research B: Solid Earth↗
USGS1982

Use of microcomputer in mapping depth of stratigraphic horizons in National Petroleum Reserve in Alaska

REGIONAL MAPPER is a menu-driven system in the BASIC language for computing and plotting (1) time, depth, and average velocity to geologic horizons, (2) interval time, thickness, and interval velocity of stratigraphic intervals, and (3) subcropping and onlapping intervals at unconformities. The system consists of three programs: FILER, TRAVERSER, and PLOTTER. A control point is a shot point with velocity analysis or a shot point at or near a well with velocity check-shot survey. Reflection time to and code number of seismic horizons are filed by digitizing tablet from record sections. TRAVERSER starts at a point of geologic control and, in traversing to another, parallels seismic events, records loss of horizons by onlap and truncation, and stores reflection time for geologic horizons at traversed shot points. TRAVERSER is basically a phantoming procedure. Permafrost thickness and velocity variations, buried canyons with low-velocity fill, and error in seismically derived velocity cause velocity anomalies that complicate depth mapping. Two depths to the top of the pebble is based shale are computed for each control point. One depth, designated Zs on seismically derived velocity. The other (Zw) is based on interval velocity interpolated linearly between wells and multiplied by interval time (isochron) to give interval thickness. Z w is computed for all geologic horizons by downward summation of interval thickness. Unknown true depth (Z) to the pebble shale may be expressed as Z = Zs + es and Z = Zw + ew where the e terms represent error. Equating the two expressions gives the depth difference D = Zs + Zw = ew + es A plot of D for the top of the pebble shale is readily contourable but smoothing is required to produce a reasonably simple surface. Seismically derived velocity used in computing Zs includes the effect of velocity anomalies but is subject to some large randomly distributed errors resulting in depth errors (es). Well-derived velocity used in computing Zw does not include the effect of velocity anomalies, but the error (ew) should reflect these anomalies and should be contourable (non-random). The D surface as contoured with smoothing is assumed to represent ew, that is, the depth effect of variations in permafrost thickness and velocity and buried canyon depth. Estimated depth (Zest) to each geologic horizon is the sum of Z w for that horizon and a constant e w as contoured for the pebble shale, which is the first highly continuous seismic horizon below the zone of anomalous velocity. Results of this 'depthing' procedure are compared with those of Tetra Tech, Inc., the subcontractor responsible for geologic and geophysical interpretation and mapping.

Open-File Report↗
USGS2010

Rayleigh-wave phase-velocity maps and three-dimensional shear velocity structure of the western US from local non-plane surface wave tomography

We utilize two-and-three-quarter years of vertical-component recordings made by the Transportable Array (TA) component of Earthscope to constrain three-dimensional (3-D) seismic shear wave velocity structure in the upper 200 km of the western United States. Single-taper spectral estimation is used to compile measurements of complex spectral amplitudes from 44 317 seismograms generated by 123 teleseismic events. In the first step employed to determine the Rayleigh-wave phase-velocity structure, we implement a new tomographic method, which is simpler and more robust than scattering-based methods (e.g. multi-plane surface wave tomography). The TA is effectively implemented as a large number of local arrays by defining a horizontal Gaussian smoothing distance that weights observations near a given target point. The complex spectral-amplitude measurements are interpreted with the spherical Helmholtz equation using local observations about a succession of target points, resulting in Rayleigh-wave phase-velocity maps at periods over the range of 18–125 s. The derived maps depend on the form of local fits to the Helmholtz equation, which generally involve the nonplane-wave solutions of Friederich et al. In a second step, the phase-velocity maps are used to derive 3-D shear velocity structure. The 3-D velocity images confirm details witnessed in prior body-wave and surface-wave studies and reveal new structures, including a deep (>100 km deep) high-velocity lineament, of width ∼200 km, stretching from the southern Great Valley to northern Utah that may be a relic of plate subduction or, alternatively, either a remnant of the Mojave Precambrian Province or a mantle downwelling. Mantle seismic velocity is highly correlated with heat flow, Holocene volcanism, elastic plate thickness and seismicity. This suggests that shallow mantle structure provides the heat source for associated magmatism, as well as thinning of the thermal lithosphere, leading to relatively high stress concentration. Our images also confirm the presence of high-velocity mantle at 100 km depth beneath areas of suspected mantle delamination (southern Sierra Nevada; Grande Ronde uplift), low velocity mantle underlying active rift zones, and high velocity mantle associated with the subducting Juan de Fuca plate. Structure established during the Proterozoic appears to exert a lasting influence on subsequent volcanism and tectonism up to the Present.

Geophysical Journal International↗
USGS2012

Rapid intrusion of magma into wet rock: Groundwater flow due to pore pressure increases

Analytical and numerical solutions are developed to simulate the pressurization, expansion, and flow of groundwater contained within saturated, intact host rocks subject to sudden heating from the planar surface of an igneous intrusion. For most rocks, water diffuses more rapidly than heat, assuring that groundwater is not heated along a constant-volume pressure path and that thermal expansion and pressurization adjacent to the intrusion drives a flow that extends well beyond the heated region. The forcing parameter for pressurization and flow is α Δ T , where α is a thermal expansion coefficient reflecting the overall expansion of water heated through the temperature difference Δ T between the initial ambient and intrusive values. Pore pressure increases due to heating are greatest when the intrusion is emplaced rapidly and where the intrusive contact is impervious to groundwater contained in stiff, impermeable rocks with high thermal diffusivities and porosities. The maximum velocity of water flowing in pores decays with the inverse square root of time and is insensitive to hydraulic properties of the host rocks. Pressures are lessened and flow directions are reversed with the onset of hydrothermal convection. This occurs at times ranging from hours to weeks after onset of intrusion. As magma rises into near-surface rocks, steam can be generated. Solutions indicate that pressure increases and velocities are sensitive to the overall amount of expansion rather than the behavior of the water-steam transition. Both the overall thermal expansion coefficient α and the temperature difference Δ T are greater in shallow (<1 km) environments than in deep (∼5 km) ones. Thus, for rocks with similar transport properties, pressure increases due to heating are greatest in shallow environments. Although solutions can be applied to rocks with a wide variety of properties, pressure increases are calculated for compliant quartz-rich sedimentary rocks with a porosities between 1 and 20% and permeabilities between 1 darcy and 1 μdarcy, subject to temperature increases of 500 and 1000 K at depths ranging from 0.1 to 5 km in a region of hydrostatic pressures and normal geothermal gradient. Such rocks, with porosities greater than 5%, permeabilities less than a 0.1 mdarcy, and drained hydrostatic compressibilities of 10 −4 /MPa, undergo pressure increases greater than 10 MPa (100 bars)for conditions typical of water table depths of 2.5 km and heating to 500 K above ambient. Similar rocks, but with permeabilities less than 1 mdarcy, undergo pressure increases of 10 MPa for conditions typical of 1 km water table depth. Rocks commonly considered to be good aquifers undergo pressure increases of less than 1 MPa, primarily because of their high permeability. Although these estimates neglect the effects of fracturing and brecciation that may accompany such pressure increases, calculations indicate that pressure increases due to heating of cool groundwater can lead to failure of host rocks by a phreatic mechanism.

Journal of Geophysical Research Solid Earth↗
USGS2010

The Quaternary Silver Creek Fault Beneath the Santa Clara Valley, California

The northwest-trending Silver Creek Fault is a 40-km-long strike-slip fault in the eastern Santa Clara Valley, California, that has exhibited different behaviors within a changing San Andreas Fault system over the past 10-15 Ma. Quaternary alluvium several hundred meters thick that buries the northern half of the Silver Creek Fault, and that has been sampled by drilling and imaged in a detailed seismic reflection profile, provides a record of the Quaternary history of the fault. We assemble evidence from areal geology, stratigraphy, paleomagnetics, ground-water hydrology, potential-field geophysics, and reflection and earthquake seismology to determine the long history of the fault in order to evaluate its current behavior. The fault formed in the Miocene more than 100 km to the southeast, as the southwestern fault in a 5-km-wide right step to the Hayward Fault, within which the 40-km-long Evergreen pull-apart basin formed. Later, this basin was obliquely cut by the newly recognized Mt. Misery Fault to form a more direct connection to the Hayward Fault, although continued growth of the basin was sufficient to accommodate at least some late Pliocene alluvium. Large offset along the San Andreas-Calaveras-Mt Misery-Hayward Faults carried the basin northwestward almost to its present position when, about 2 Ma, the fault system was reorganized. This led to near abandonment of the faults bounding the pull-apart basin in favor of right slip extending the Calaveras Fault farther north before stepping west to the Hayward Fault, as it does today. Despite these changes, the Silver Creek Fault experienced a further 200 m of dip slip in the early Quaternary, from which we infer an associated 1.6 km or so of right slip, based on the ratio of the 40-km length of the strike-slip fault to a 5-km depth of the Evergreen Basin. This dip slip ends at a mid-Quaternary unconformity, above which the upper 300 m of alluvial cover exhibits a structural sag at the fault that we interpret as a negative flower structure. This structure implies some continuing strike slip on the Silver Creek Fault in the late Quaternary as well, with a transtensional component but no dip slip. Our only basis for estimating the rate of this later Quaternary strike slip on the Silver Creek Fault is to assume continuation of the inferred early Quaternary rate of less than 2 mm/yr. Faulting evident in a detailed seismic reflection profile across the Silver Creek Fault extends up to the limit of data at a depth of 50 m and age of about 140 ka, and the course of Coyote Creek suggests Holocene capture in a structural depression along the fault. No surface trace is evident on the alluvial plain, however, and convincing evidence of Holocene offset is lacking. Few instrumentally recorded earthquakes are located near the fault, and those that are near its southern end represent cross-fault shortening, not strike slip. The fault might have been responsible, however, for two poorly located moderate earthquakes that occurred in the area in 1903. Its southeastern end does mark an abrupt change in the pattern of abundant instrumentally recorded earthquakes along the Calaveras Fault-in both its strike and in the depth distribution of hypocenters-that could indicate continuing influence by the Silver Creek Fault. In the absence of convincing evidence to the contrary, and as a conservative estimate, we presume that the Silver Creek Fault has continued its strike-slip movement through the Holocene, but at a very slow rate. Such a slow rate would, at most, yield very infrequent damaging earthquakes. If the 1903 earthquakes did, in fact, occur on the Silver Creek Fault, they would have greatly reduced the short-term future potential for large earthquakes on the fault.

Open-File Report↗
USGS2008

Hydrologic models of modern and fossil geothermal systems in the Great Basin: Genetic implications for epithermal Au-Ag and Carlin-type gold deposits

The Great Basin region in the western United States contains active geothermal systems, large epithermal Au-Ag deposits, and world-class Carlin-type gold deposits. Temperature profiles, fluid inclusion studies, and isotopic evidence suggest that modern and fossil hydrothermal systems associated with gold mineralization share many common features, including the absence of a clear magmatic fluid source, discharge areas restricted to fault zones, and remarkably high temperatures (>200 ??C) at shallow depths (200-1500 m). While the plumbing of these systems varies, geochemical and isotopic data collected at the Dixie Valley and Beowawe geothermal systems suggest that fluid circulation along fault zones was relatively deep (>5 km) and comprised of relatively unexchanged Pleistocene meteoric water with small (<2.5%) shifts from the meteoric water line (MWL). Many fossil ore-forming systems were also dominated by meteoric water, but usually exhibit ??18O fluid-rock interactions with larger shifts of 5???-20??? from the MWL. Here we present a suite of two-dimensional regional (100 km) and local (40-50 km) scale hydrologic models that we have used to study the plumbing of modern and Tertiary hydrothermal systems of the Great Basin. Geologically and geophysically consistent cross sections were used to generate somewhat idealized hydrogeologic models for these systems that include the most important faults, aquifers, and confining units in their approximate configurations. Multiple constraints were used, including enthalpy, ??18O, silica compositions of fluids and/or rocks, groundwater residence times, fluid inclusion homogenization temperatures, and apatite fission track anomalies. Our results suggest that these hydrothermal systems were driven by natural thermal convection along anisotropic, subvertical faults connected in many cases at depth by permeable aquifers within favorable lithostratigraphic horizons. Those with minimal fluid ?? 18O shifts are restricted to high-permeability fault zones and relatively small-scale (???5 km), single-pass flow systems (e.g., Beowawe). Those with intermediate to large isotopic shifts (e.g., epithermal and Carlin-type Au) had larger-scale (???15 km) loop convection cells with a greater component of flow through marine sedimentary rocks at lower water/rock ratios and greater endowments of gold. Enthalpy calculations constrain the duration of Carlin-type gold systems to probably <200 k.y. Shallow heat flow gradients and fluid silica concentrations suggest that the duration of the modern Beowawe system is <5 k.y. However, fluid flow at Beowawe during the Quaternary must have been episodic with a net duration of ???200 k.y. to account for the amount of silica in the sinter deposits. In the Carlin trend, fluid circulation extended down into Paleozoic siliciclastic rocks, which afforded more mixing with isotopically enriched higher enthalpy fluids. Computed fission track ages along the Carlin trend included the convective effects, and ranged between 91.6 and 35.3 Ma. Older fission track ages occurred in zones of groundwater recharge, and the younger ages occurred in discharge areas. This is largely consistent with fission track ages reported in recent studies. We found that either an amagmatic system with more permeable faults (10-11 m2) or a magmatic system with less permeable faults (10-13 m2) could account for the published isotopic and thermal data along the Carlin trend systems. Localized high heat flow beneath the Muleshoe fault was needed to match fl uid inclusion temperatures at Mule Canyon. However, both magmatic and amagmatic scenarios require the existence of deep, permeable faults to bring hot fluids to the near surface. ?? 2008 Geological Society of America.

Geosphere↗
USGS1987

Summary mineral resource appraisal of the Richfield 1° x 2° quadrangle, west-central Utah

The mineral resource potential of the Richfield 1° x 2° quadrangle, Utah, has been appraised using geological, geophysical, geochemical, and remote-sensing techniques. These studies have led to many publications giving basic data and interpretations; of these, a series of 18 maps at 1:250,000 and 1:500,000 scales summarizing aspects of the geology, geophysics, geochemistry, and remote sensing is designated the CUSMAP (Conterminous United States Mineral Appraisal Program) folio. This circular uses the data shown on these maps to appraise the mineral resource potential of the quadrangle. The oldest rocks exposed in the Richfield quadrangle are small patches of Early Proterozoic (1.7 billion years old) gneiss and schist on the west side of the Mineral Mountains. These rocks presumably formed the basement on which many thousands of meters of Late Proterozoic, Paleozoic, and lower Mesozoic sedimentary strata were deposited. These rocks were deformed during the Late Cretaceous Sevier orogeny when Precambrian and Paleozoic strata in the western part of the quadrangle were thrust relatively eastward across Paleozoic and Mesozoic strata in the eastern part of the quadrangle. Late Cretaceous and early Tertiary highlands above the overthrust belt were eroded and much of the debris was deposited in broad basins east of the belt. Volcanism in Oligocene and earliest Miocene time formed an east-northeast-trending belt of calcalkalic volcanoes across the southern half of the quadrangle. In early Miocene time, the composition of the volcanic rocks changed to a bimodal assemblage of mafic rocks and high-silica alkali rhyolite that has been erupted episodically ever since. Syngenetic mineral resources developed during formation of both sedimentary and volcanic rocks. These include limestone and dolomite, silica-rich sandstone, metalliferous black shale, evaporite deposits, zeolite deposits, pumice, cinders and scoria, and evaporitic or diagenetic deposits in playa environments. Most of these deposits need to have markets established, or extraction and fabrication techniques developed, for them to be utilized. Most epigenetic deposits are of volcanogenic-hydrothermal origin. Deposits associated with calc-alkalic igneous activity largely contain Cu, Pb, Zn, Au, and Ag, and occur in a variety of types zoned around core intrusions. Younger deposits are mostly associated with silicic igneous centers belonging to the bimodal mafic-silicic igneous association. Resources associated with this latter group are likely to contain one or more of the elements Mo, W, U, Sn, Be, and F, as well as Pb, Zn, Au, and Ag. Alunite and kaolinite deposits are found at many mineralized centers. Most epigenetically mineralized areas expose only the upper, near-surface parts of the different hydrothermal systems; most of whatever mineral deposits formed in these systems probably still exist at depth, awaiting discovery. Our conclusion is that many mineralized areas have excellent possibilities for the occurrence of mineral resources. Each of the many identified centers of mineralization is discussed briefly in this report and an estimate made of its resource potential.

Utah↗
USGS1986

Application of the aerial profiling of terrain system

Test well Ch-Bf 146 was drilled to 1,650 ft below land surface to explore and evaluate the Patapsco Formation aquifers deeper than 1,000 ft in the Waldorf/La Plata area of Charles County, Maryland. The test hole penetrated two major aquifer systems above a depth of 1,000 ft and another major aquifer system and aquifer below 1,000 ft. The deepest aquifer system found in the test drilling, the Lower Patapsco aquifer system, is comprised of a group of sands which lie between 1,140 ft and the base of the Patapsco Formation at 1,417 ft. The Waldorf aquifer system, overlies the Upper Patapsco aquifer system and includes sands in the Waldorf/La Plata area that are referred to as the ' Magothy aquifer ' by other. At the test well site, the Waldorf aquifer system consists solely of sands in the Magothy and Monmouth Formations. The three major aquifer systems and the Middle Patapsco aquifer are correlative along the regional strike for at least 10 miles in the Waldorf/La Plata area. Correlations of geophysical logs indicate that individual sands 15 to 25 ft thick can be traced between wells. Geophysical log correlations also indicate that, at some well sites, the Waldorf aquifer system and the Upper Patapsco aquifer system are in sand-on-sand contact and function as one major hydrologic system at these sites. A production well (Ch-Bf 147) was drilled at the site of the test well (Ch-Bf 146). The initial static water level of these aquifers was 9.8 ft below sea level. Evaluation of the 24-hr pumping test of Ch-Bf 147 resulted in a transmissivity of 14,000 gpd/ft (gallons per day per foot) for the first 2 hours of the test and 7,300 gpd/ft for the remainder of the test. The decrease of 6,700 gpd/ft in the transmissivity indicates that the well 's expanding cone of depression intersected a transmissivity boundary after 2 hours of pumping. The Patapsco sands tested in Ch-Bf 147 have a storage coefficient of 0.0016. The hydraulic conductivity of the producing sands in Ch-Bf 147 is 127 gpd/sq ft. The quality of groundwater from the Lower Patapsco aquifer system and the Middle Patapsco aquifer is suitable for most purposes. The water is nearly neutral (pH = 7.4), low in dissolved solids (about 200 mg/L), low in specific conductance, slightly warm (22 C), very soft, naturally fluoridated, and moderately high in dissolved silica. (Lantz-PTT)

Open-File Report↗
USGS2007

Geochemical evidence for African dust inputs to soils of western Atlantic islands: Barbados, the Bahamas, and Florida

We studied soils on high-purity limestones of Quaternary age on the western Atlantic Ocean islands of Barbados, the Florida Keys, and the Bahamas. Potential soil parent materials in this region, external to the carbonate substrate, include volcanic ash from the island of St. Vincent (near Barbados), volcanic ash from the islands of Dominica and St. Lucia (somewhat farther from Barbados), the fine-grained component of distal loess from the lower Mississippi River Valley, and wind-transported dust from Africa. These four parent materials can be differentiated using trace elements (Sc, Cr, Th, and Zr) and rare earth elements that have minimal mobility in the soil-forming environment. Barbados soils have compositions that indicate a complex derivation. Volcanic ash from the island of St. Vincent appears to have been the most important influence, but African dust is a significant contributor, and even Mississippi River valley loess may be a very minor contributor to Barbados soils. Soils on the Florida Keys and islands in the Bahamas appear to have developed mostly from African dust, but Mississippi River valley loess may be a significant contributor. Our results indicate that inputs of African dust are more important to the genesis of soils on islands in the western Atlantic Ocean than previously supposed. We hypothesize that African dust may also be a major contributor to soils on other islands of the Caribbean and to soils in northern South America, central America, Mexico, and the southeastern United States. Dust inputs to subtropical and tropical soils in this region increase both nutrient-holding capacity and nutrient status and thus may be critical in sustaining vegetation. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2001

Rheology of ice I at low stress and elevated confining pressure

Triaxial compression testing of pure, polycrystalline water ice I at conditions relevant to planetary interiors and near-surface environments (differential stresses 0.45 to 10 MPa, temperatures 200 to 250 K, confining pressure 50 MPa) reveals that a complex variety of rheologies and grain structures may exist for ice and that rheology of ice appears to depend strongly on the grain structures. The creep of polycrystalline ice I with average grain size of 0.25 mm and larger is consistent with previously published dislocation creep laws, which are now extended to strain rates as low as 2×10 −8 s −1 . When ice I is reduced to very fine and uniform grain size by rapid pressure release from the ice II stability field, the rheology changes dramatically. At 200 and 220 K the rheology matches the grain-size-sensitive rheology measured by Goldsby and Kohlstedt [1997, this issue] at 1 atm. This finding dispels concerns that the Goldsby and Kohlstedt results were influenced by mechanisms such as microfracturing and cavitation, processes not expected to operate at elevated pressures in planetary interiors. At 233 K and above, grain growth causes the fine-grained ice to become more creep resistant. Scanning electron microscopy investigation of some of these deformed samples shows that grains have markedly coarsened and the strain hardening can be modeled by normal grain growth and the Goldsby and Kohlstedt rheology. Several samples also displayed very heterogeneous grain sizes and high aspect ratio grain shapes. Grain-size-sensitive creep and dislocation creep coincidentally contribute roughly equal amounts of strain rate at conditions of stress, temperature, and grain size that are typical of terrestrial and planetary settings, so modeling ice dynamics in these settings must include both mechanisms.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Seismicity induced by massive wastewater injection near Puerto Gaitán, Colombia

Seven years after the beginning of a massive wastewater injection project in eastern Colombia, local earthquake activity increased significantly. The field operator and the Colombian Geological Survey immediately reinforced the monitoring of the area. Our analysis of the temporal evolution of the seismic and injection data together with our knowledge of the geological parameters of the region indicate that the surge of seismicity is being induced by the re-injection of produced water into the same three producing reservoirs. Earthquake activity began on known faults once disposal rates had reached a threshold of ∼2 × 10 6 m 3 of water per month. The average reservoir pressure had remained constant at 7.6 MPa after several years of production, sustained by a large, active aquifer. Surface injection pressures in the seismically active areas remain below 8.3 MPa, a value large enough to activate some of the faults. Since faults are mapped throughout the region and many do not have seismicity on them, we conclude that the existence of known faults is not the only control on whether earthquakes are generated. Stress conditions of these faults are open to future studies. Earthquakes are primarily found in four clusters, located near faults mapped by the operator. The hypocentres reveal vertical planes with orientations consistent with focal mechanisms of these events. Stress inversion of the focal mechanisms gives a maximum compression in the direction ENE-WSW, which is in agreement with borehole breakout measurements. Since the focal mechanisms of the earthquakes are consistent with the tectonic stress regime, we can conclude that the seismicity is resulting from the activation of critically stressed faults. Slip was progressive and seismic activity reached a peak before declining to few events per month. The decline in seismicity suggests that most of the stress has been relieved on the main faults. The magnitude of a large majority of the recorded earthquakes was lower than 4, as the pore pressure disturbance did not reach the mapped large faults whose activation might have resulted in larger magnitude earthquakes. Our study shows that a good knowledge of the local fault network and conditions of stress is of paramount importance when planning a massive water disposal program. These earthquakes indicate that while faults provide an opportunity to dispose produced water at an economically attractive volume–pressure ratio, the possibility of induced seismicity must also be considered.

Puerto Gaitan↗
USGS1999

Nevada STORMS project: Measurement of mercury emissions from naturally enriched surfaces

Diffuse anthropogenic and naturally mercury-enriched areas represent long-lived sources of elemental mercury to the atmosphere. The Nevada Study and Tests of the Release of Mercury From Soils (STORMS) project focused on the measurement of mercury emissions from a naturally enriched area. During the project, concurrent measurements of mercury fluxes from naturally mercury-enriched substrate were made September 1–4, 1997, using four micrometeorological methods and seven field flux chambers. Ambient air mercury concentrations ranged from 2 to nearly 200 ng m −3 indicating that the field site is a source of atmospheric mercury. The mean daytime mercury fluxes, during conditions of no precipitation, measured with field chambers were 50 to 360 ng m −2 h −1 , and with the micrometeorological methods were 230 to 600 ng m −2 h −1 . This wide range in mercury emission rates reflects differences in method experimental designs and local source strengths. Mercury fluxes measured by many field chambers were significantly different ( p <0.05) but linearly correlated. This indicates that field chambers responded similarly to environmental conditions, but differences in experimental design and site heterogeneity had a significant influence on the magnitude of mercury fluxes. Data developed during the field study demonstrated that field flux chambers are ideal for assessment of the physicochemical processes driving mercury flux and development of an understanding of the magnitude of the influence of individual factors on flux. In general, mean mercury fluxes measured with micrometeorological methods during daytime periods were nearly 3 times higher than mean fluxes measured with field flux chambers. Micrometeorological methods allow for derivation of a representative mercury flux occurring from an unconstrained system and provide an assessment of the actual magnitude and variability of fluxes occurring from an area.

Journal of Geophysical Research D: Atmospheres↗
USGS1988

The Geology and Remarkable Thermal Activity of Norris Geyser Basin, Yellowstone National Park, Wyoming

Norris Geyser Basin, normally shortened to Norris Basin, is adjacent to the north rim of the Yellowstone caldera at the common intersection of the caldera rim and the Norris-Mammoth Corridor, a zone of faults, volcanic vents, and thermal activity that strikes north from the caldera rim to Mammoth Hot Springs. An east-west fault zone terminates the Gallatin Range at its southern end and extends from Hebgen Lake, west of the park, to Norris Basin. No local evidence exists at the surface in Norris Basin for the two oldest Yellowstone volcanic caldera cycles (~2.0 and 1.3 m.y.B.P.). The third and youngest cycle formed the Yellowstone caldera, which erupted the 600,000-year-old Lava Creek Tuff. No evidence is preserved of hydrothermal activity near Norris Basin during the first 300,000.years after the caldera collapse. Glaciation probably removed most of the early evidence, but erratics of hot-spring sinter that had been converted diagenetically to extremely hard, resistant chalcedonic sinter are present as cobbles in and on some moraines and till from the last two glacial stages, here correlated with the early and late stages of the Pinedale glaciation <150,000 years B.P.). Indirect evidence for the oldest hydrothermal system at Norris Basin indicates an age probably older than both stages of Pinedale glaciation. Stream deposits consisting mainly of rounded quartz phenocrysts of the Lava Creek Tuff were subaerial, perhaps in part windblown and redeposited by streams. A few small rounded pebbles are interpreted as chalcedonic sinter of a still older cycle. None of these are precisely dated but are unlikely to be more than 150,000 to 200,000 years old. ...Most studies of active hydrothermal areas have noted chemical differences in fluids and alteration products but have given little attention to differences and models to explain evolution in types. This report, in contrast, emphasizes the kinds of changes in vents and their changing chemical types of waters and then provides models for explaining these differences. Norris Basin is probably not an independent volcanic-hydrothermal system. The basin and nearby acid-leached areas (from oxidation of H2S-enriched vapor) are best considered as parts of the same system, extending from Norris Basin to Roaring Mountain and possibly to Mammoth. If so, are they parts of a single large system centered within the Yellowstone caldera, or are Norris Basin and the nearby altered areas both parts of one or more young independent corridor systems confined, at least in the shallow crust, to the Norris-Mammoth Corridor? Tentatively, we favor the latter relation, probably having evolved in the past ~300,000 years. A model for large, long-lived, volcanic-hydrothermal activity is also suggested, involving all of the crust and upper mantle and using much recent geophysical data bearing on crust-mantle interrelations. Our model for large systems is much superior to previous suggestions for explaining continuing hydrothermal activity over hundreds of thousands of years, but is less attractive for the smaller nonhomogenized volcanic system actually favored here for the Norris-Mammoth Corridor.

Professional Paper↗
USGS2002

35 - Strength and energetics of active fault zones

The strength of active fault zones, i.e., the shear stress level required to cause fault slip, is fundamental to understanding the physics of earthquakes and to assessing earthquake hazard. Although many researchers have concluded that fault zones are weak (shear stresses 10 MPa or less averaged between 0 and ˜20 km depth), others maintain that faults are strong (˜100 MPa average of an approximate linear increase with depth). Thus, despite 30 y of dedicated research, relevant data remain inconclusive and fault strength remains uncertain by an order of magnitude. In part, this is because the main source of energy release in earthquakes is at depths greater than 5 km, inaccessible to direct instrumental observation. Very large earthquakes rupture to the Earth's surface where direct observation of the shallow rupture process is possible. However, the rupture characteristics at shallow depth may differ from those at seismogenic depths. To date no great ( M ≥ 7.75) earthquakes have occurred within a network of modern strong-motion instruments, but the large, well-recorded 1999 earthquakes in Turkey and Taiwan, both M = 7.6, show that this data gap is rapidly being filled. Furthermore, shear stress changes at the earthquake source (and the resulting seismic waves) are nearly linear perturbations of the absolute stress field. Thus, an unknown absolute background stress does not greatly affect the basic characteristics of the observed low frequency seismic waves and the observed geodetic deformation. Other, less-direct data must then be used to infer the physical state and ambient stress levels on active faults. This chapter discusses available evidence and current ideas about fault zone strength and energetics. In our review we begin by outlining the general physical conditions prevailing in the Earth around active faults and summarize the generally agreed upon observational features of faulting and earthquake occurrence ( Section 2 ). In doing so we make as few assumptions as possible, because making uncertain assumptions may lead to logical inconsistencies and apparent paradoxes. For example, the absence of a measurable, frictionally generated heat flow anomaly near active faults, the so-called stress-heat flow paradox, rests upon assumptions that must somehow be incorrect. We wish to avoid such inconsistencies and to begin we introduce the observations and briefly state what they imply about fault zones. Subsequently, in Sections 3–8 , we place these observations in an interpretative context and show how they have been used to infer various measures of fault shear stress. Section 9 summarizes our assessment of the average shear stress state near faults, argues for the general importance of stress heterogeneity in faulting processes, and discusses its implications.

International Geophysics↗
USGS2000

The enigma of the Arthur's Pass, New Zealand, earthquake: 1. Reconciling a variety of data for an unusual earthquake sequence

The 1994 Arthur's Pass earthquake ( M W 6.7) is the largest in a recent sequence of earthquakes in the central South Island, New Zealand. No surface rupture was observed, the aftershock distribution was complex, and routine methods of obtaining the faulting orientation of this earthquake proved contradictory. We use a range of data and techniques to obtain our preferred solution, which has a centroid depth of 5 km, M 0 =1.3 × 10 19 N m, and a strike, dip, and rake of 221°, 47°, 112°, respectively. Discrepancies between this solution and the Harvard centroid moment tensor, together with the Global Positioning System (GPS) observations and unusual aftershock distribution, suggest that the rupture may not have occurred on a planar fault. A second, strike slip, subevent on a more northerly striking plane is suggested by these data but neither the body wave modeling nor regional broadband recordings show any complexity or late subevents. We relocate the aftershocks using both one-dimensional and three-dimensional velocity inversions. The depth range of the aftershocks (1–10 km) agrees well with the preferred mainshock centroid depth. The aftershocks near the hypocenter suggest a structure dipping toward the NW, which we interpret to be the mainshock fault plane. This structure and the Harper fault, ∼15 km to the south, appear to have acted as boundaries to the extensive aftershock zone trending NNW-SSE. Most of the M L ≥ 5 aftershocks, including the two largest ( M L 6.1 and M L 5.7), clustered near the Harper fault and have strike slip mechanisms consistent with motion on this fault and its conjugates. Forward modeling of the GPS data suggests that a reverse slip mainshock, combined with strike slip aftershock faulting in the south, is able to match the observed displacements. The occurrence of this earthquake sequence implies that the level of seismic hazard in the central South Island is greater than previous estimates.

Arthur's Pass↗
USGS2006

Distribution of stress drop, stiffness, and fracture energy over earthquake rupture zones

Using information provided by slip models and the methodology of McGarr and Fletcher (2002), we map static stress drop, stiffness (k = ????/u, where ???? is static stress drop and u is slip), and fracture energy over the slip surface to investigate the earthquake rupture process and energy budget. For the 1994 M6.7 Northridge, 1992 M7.3 Landers, and 1995 M6.9 Kobe earthquakes, the distributions of static stress drop show strong heterogeneity, emphasizing the importance of asperities in the rupture process. Average values of static stress drop are 17, 11, and 4 Mpa for Northridge, Landers, and Kobe, respectively. These values are substantially higher than estimates based on simple crack models, suggesting that the failure process involves the rupture of asperities within the larger fault zone. Stress drop as a function of depth for the Northridge and Landers earthquakes suggests that stress drops are limited by crustal strength. For these two earthquakes, regions of high slip are surrounded by high values of stiffness. Particularly for the Northridge earthquake, the prominent patch of high slip in the central part of the fault is bordered by a ring of high stiffness and is consistent with expectations based on the failure of an asperity loaded at its edge due to exterior slip. Stiffness within an asperity is inversely related to its dimensions. Estimates of fracture energy, based on static stress drop, slip, and rupture speed, were used to investigate the nature of slip weakening at four locations near the hypocenter of the Kobe earthquake for comparison with independent results based on a dynamic model of this earthquake. One subfault updip and to the NE of the hypocenter has a fracture energy of 1.1 MJ/m2 and a slip-weakening distance, Dc, of 0.66 m. Right triangles, whose base and height are Dc and the dynamic stress drop, respectively, approximately overlie the slip-dependent stress given by Ide and Takeo (1997) for the same locations near the hypocenter. The total fracture energy for the Kobe earthquake, 3.7 ?? 1014 J, is about the same as the seismic energy (Ea = 3.2 ?? 1014 J.

Journal of Geophysical Research B: Solid Earth↗