Search USGS⌕ Search

SEARCH · Search USGS

Results for “Geophysical Investigations Map”

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

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

At least 1,351 records · Page 75Linked to original sources

Inversion of multi-frequency electromagnetic induction data for 3D characterization of hydraulic conductivity

Electromagnetic induction (EMI) instruments provide rapid, noninvasive, and spatially dense data for characterization of soil and groundwater properties. Data from multi-frequency EMI tools can be inverted to provide quantitative electrical conductivity estimates as a function of depth. In this study, multi-frequency EMI data collected across an abandoned uranium mill site near Naturita, Colorado, USA, are inverted to produce vertical distribution of electrical conductivity ( EC ) across the site. The relation between measured apparent electrical conductivity ( EC a ) and hydraulic conductivity ( K ) is weak (correlation coefficient of 0.20), whereas the correlation between the depth dependent EC obtained from the inversions, and K is sufficiently strong to be used for hydrologic estimation (correlation coefficient of − 0.62). Depth-specific EC values were correlated with co-located K measurements to develop a site-specific ln( EC )–ln( K ) relation. This petrophysical relation was applied to produce a spatially detailed map of K across the study area. A synthetic example based on EC a values at the site was used to assess model resolution and correlation loss given variations in depth and/or measurement error. Results from synthetic modeling indicate that optimum correlation with K occurs at ~ 0.5 m followed by a gradual correlation loss of 90% at 2.3 m. These results are consistent with an analysis of depth of investigation (DOI) given the range of frequencies, transmitter–receiver separation, and measurement errors for the field data. DOIs were estimated at 2.0 ± 0.5 m depending on the soil conductivities. A 4-layer model, with varying thicknesses, was used to invert the EC a to maximize available information within the aquifer region for improved correlations with K . Results show improved correlation between K and the corresponding inverted EC at similar depths, underscoring the importance of inversion in using multi-frequency EMI data for hydrologic estimation.

Journal of Applied Geophysics↗

Transient electromagnetic sounding for groundwater

The feasibility of using the transient electromagnetic sounding (TS or TDEM) method for groundwater exploration can be studied by means of numerical models. As examples of its applicability to groundwater exploration, we study four groundwater exploration problems: (1) mapping of alluvial fill and gravel zones over bedrock; (2) mapping of sand and gravel lenses in till; (3) detection of salt or brackish water interfaces in freshwater aquifers; and (4) determination of hydrostratigraphy. These groundwater problems require determination of the depth to bedrock; location of resistive, high‐porosity zones associated with fresh water; determination of formation resistivity to assess water quality; and determination of lithology and geometry, respectively. The TS method is best suited for locating conductive targets, and has very good vertical resolution. Unlike other sounding techniques where the receiver‐transmitter array must be expanded to sound more deeply, the depth of investigation for the TS method is a function of the length of time the transient is recorded. Present equipment limitations require that exploration targets with resistivities of 50 Ω ⋅ m or more be at least 50 m deep to determine their resistivity. The maximum depth of exploration is controlled by the geoelectrical section and background electromagnetic (EM) noise. For a particular exploration problem, numerical studies are recommended to determine if the target is detectable.

Geophysics↗

Three-Dimensional Geologic Model of Complex Fault Structures in the Upper Seco Creek Area, Medina and Uvalde Counties, South-Central Texas

This multimedia report shows and describes digital three-dimensional faulted geologic surfaces and volumes of the lithologic units of the Edwards aquifer in the upper Seco Creek area of Medina and Uvalde Counties in south-central Texas. This geologic framework model was produced using (1) geologic maps and interpretations of depositional environments and paleogeography; (2) lithologic descriptions, interpretations, and geophysical logs from 31 drill holes; (3) rock core and detailed lithologic descriptions from one drill hole; (4) helicopter electromagnetic geophysical data; and (5) known major and minor faults in the study area. These faults were used because of their individual and collective effects on the continuity of the aquifer-forming units in the Edwards Group. Data and information were compared and validated with each other and reflect the complex relationships of structures in the Seco Creek area of the Balcones fault zone. This geologic framework model can be used as a tool to visually explore and study geologic structures within the Seco Creek area of the Balcones fault zone and to show the connectivity of hydrologic units of high and low permeability between and across faults. The software can be used to display other data and information, such as drill-hole data, on this geologic framework model in three-dimensional space.

Scientific Investigations Report↗

Investigating landslides caused by earthquakes - A historical review

Post-earthquake field investigations of landslide occurrence have provided a basis for understanding, evaluating, and mapping the hazard and risk associated with earthquake-induced landslides. This paper traces the historical development of knowledge derived from these investigations. Before 1783, historical accounts of the occurrence of landslides in earthquake are typically so incomplete and vague that conclusions based on these accounts are of limited usefulness. For example, the number of landslides triggered by a given event is almost always greatly underestimated. The first formal, scientific post-earthquake investigation that included systematic documentation of the landslides was undertaken in the Calabria region of Italy after the 1783 earthquake swarm. From then until the mid-twentieth century, the best information on earthquake-induced landslides came from a succession of post-earthquake investigations largely carried out by formal commissions that undertook extensive ground-based field studies. Beginning in the mid-twentieth century, when the use of aerial photography became widespread, comprehensive inventories of landslide occurrence have been made for several earthquakes in the United States, Peru, Guatemala, Italy, El Salvador, Japan, and Taiwan. Techniques have also been developed for performing "retrospective" analyses years or decades after an earthquake that attempt to reconstruct the distribution of landslides triggered by the event. The additional use of Geographic Information System (GIS) processing and digital mapping since about 1989 has greatly facilitated the level of analysis that can applied to mapped distributions of landslides. Beginning in 1984, synthesis of worldwide and national data on earthquake-induced landslides have defined their general characteristics and relations between their occurrence and various geologic and seismic parameters. However, the number of comprehensive post-earthquake studies of landslides is still relatively small, and one of the most pressing needs in this area of research is for the complete documentation of landslides triggered by many more earthquakes in a wider variety of environments.

Surveys in Geophysics↗

Fluctuations in groundwater levels related to regional and local withdrawals in the fractured-bedrock groundwater system in northern Wake County, North Carolina, March 2008-February 2009

A study of dewatering of the fractured-bedrock aquifer in a localized area of east-central North Carolina was conducted from March 2008 through February 2009 to gain an understanding of why some privately owned wells and monitoring wells were intermittently dry. Although the study itself was localized in nature, the resulting water-resources data and information produced from the study will help enable resource managers to make sound water-supply and water-use decisions in similar crystalline-rock aquifer setting in parts of the Piedmont and Blue Ridge Physiographic Provinces. In June 2005, homeowners in a subdivision of approximately 11 homes on lots approximately 1 to 2 acres in size in an unincorporated area of Wake County, North Carolina, reported extremely low water pressure and temporarily dry wells during a brief period. This area of the State, which is in the Piedmont Physiographic Province, is undergoing rapid growth and development. Similar well conditions were reported again in July 2007. In an effort to evaluate aquifer conditions in the area of intermittent water loss, a study was begun in March 2008 to measure and monitor water levels and groundwater use. During the study period from March 2008 through February 2009, regular dewatering of the fractured-bedrock aquifer was documented with water levels in many wells ranging between 100 and 200 feet below land surface. Prior to this period, water levels from the 1980s through the late 1990s were reported to range from 15 to 50 feet below land surface. The study area includes three community wells and more than 30 private wells within a 2,000-foot radius of the dewatered private wells. Although groundwater levels were low, recovery was observed during periods of heavy rainfall, most likely a result of decreased withdrawals owing to less demand for irrigation purposes. Similar areal patterns of low groundwater levels were delineated during nine water-level measurement periods from March 2008 through February 2009. Correlation of groundwater-level distribution patterns with orientations of geologic structures obtained from surficial mapping, borehole geophysical measurements, and interpretation of fracture traces suggests two dominant trends striking north-south and N. 65 degrees W. A variation in overall response to groundwater withdrawals was noted in the continuous groundwater-level records for the monitored observation wells and dewatered private wells. The largest overall declines during the study period were observed in an observation well in which the water-level declined as much as 247 feet from mid-July through early August 2008, during a period of heavy usage. A private well had a water-level decline of about 94 feet during the same monitoring period. The large declines recorded in the observation well and the private well indicated a substantial temporary loss of storage in the fractured-bedrock aquifer near the wells, thus reducing the amount of water available to shallow wells in the area (those wells with total depths of about 300 feet), and resulting in temporary well failures until such time as the aquifer recovered.

North Carolina↗

Delineating the Pierre Shale from geophysical surveys east and southeast of Ellsworth Air Force Base, South Dakota, 2021

The U.S. Geological Survey, in cooperation with the U.S. Air Force Civil Engineer Center, used surface-geophysical methods to delineate the top of Cretaceous Pierre Shale along survey transects in selected areas east and southeast of Ellsworth Air Force Base, South Dakota, from April to September 2021. Two complementary geophysical methods—electrical resistivity and passive seismic—were used along 21 colocated transect surveys east and southeast of Ellsworth Air Force Base for a total of 24.7 line-kilometers. Electrical resistivity results were analyzed using EarthImager2D electrical resistivity tomography processing and inversion software. Two-dimensional earth models showing the electrical properties of the subsurface were evaluated by directly comparing the high and low subsurface resistivity values to a surficial-geologic map and nearby wells with drillers logs. Passive seismic data were analyzed using the horizontal-to-vertical spectral ratio method to determine the depth to the Cretaceous Pierre Shale at each survey point. The depth to the Pierre Shale along the transects ranged from 0.0 to about 19.8 meters, and the mean and median depths were about 6.1 and 5.6 meters, respectively. The elevation of the Pierre Shale and thickness of unconsolidated deposits generally increased with land-surface elevation from south to north; however, some transects displayed topographically high and low areas that did not correlate with land-surface topography.

South Dakota↗

Investigation of the high-frequency attenuation parameter, κ (kappa), from aftershocks of the 2010 M w 8.8 Maule, Chile earthquake

The Bío Bío region of Chile experienced a vigorous aftershock sequence following the 2010 February 27 M w 8.8 Maule earthquake. The immediate aftershock sequence was captured by two temporary seismic deployments: the Quake Catcher Network Rapid Aftershock Mobilization Program (QCN RAMP) and the Incorporated Research Institutions for Seismology CHile Aftershock Mobilization Program (IRIS CHAMP). Here, we use moderate to large aftershocks ( M L ≥ 4.0) occurring between 2010 March 1 and June 30 recorded by QCN RAMP and IRIS CHAMP stations to determine the spectral decay parameter, kappa ( κ ). First, we compare waveforms and κ estimates from the lower-resolution QCN stations to the IRIS CHAMP stations to ensure the QCN data are of sufficient quality. We find that QCN stations provide reasonable estimates of κ in comparison to traditional seismic sensors and provide valuable additional observations of local ground motion variation. Using data from both deployments, we investigate the variation in κ for the region to determine if κ is influenced primarily by local geological structure, path attenuation, or source properties (e.g. magnitude, mechanism and depth). Estimates of κ for the Bío Bío region range from 0.0022 to 0.0704 s with a mean of 0.0295 s and are in good agreement with κ values previously reported for similar tectonic environments. κ correlates with epicentral distance and, to a lesser degree, with source magnitude. We find little to no correlation between the site kappa, κ 0 , and mapped geology, although we were only able to compare the data to a low-resolution map of surficial geology. These results support an increasing number of studies that suggest κ observations can be attributed to a combination of source, path and site properties; additionally, measured κ are often highly scattered making it difficult to separate the contribution from each of these factors. Thus, our results suggest that contributions from the site, path and source should be carefully considered when interpreting κ values.

Bío Bío Region↗

Archive of bathymetry and backscatter data collected in 2014 nearshore Breton and Gosier Islands, Breton National Wildlife Refuge, Louisiana

As part of the Barrier Island Monitoring Project, scientists from the U.S. Geological Survey (USGS) St. Petersburg Coastal and Marine Science Center conducted nearshore geophysical surveys off Breton and Gosier Islands, Louisiana, in July and August of 2014. To assist the United States Fish and Wildlife Service (USFWS) with restoration planning efforts, the USGS was tasked with answering fundamental questions about the physical environment of the southern Chandeleur Islands, including the geology, morphology, and oceanography. Baseline data needed to answer these questions were either insufficient or missing. The USGS conducted a comprehensive geologic investigation in the summer of 2014, collecting geophysical and sedimentological data. Breton Island, located at the southern end of the Chandeleur Island chain in southeastern Louisiana, was recognized as a natural, globally significant nesting sanctuary for several bird species and was established as the Breton National Wildlife Refuge (NWR) in 1904. The areal extent of Breton Island has diminished 90 percent since 1920. Land loss is attributed to ongoing relative sea-level rise, diminished sediment supply, and storm impacts. The bird population on Breton Island has also declined over the years, most notably after Hurricane George in 1998 and after Hurricane Katrina in 2015; the latter completely submerged the island. Despite decreasing habitable acreage, migratory seabirds continue to return and nest on Breton Island. To prevent the island from being submerged in the future, and to protect, stabilize, and provide more nesting and foraging areas for the bird population, the USFWS proposed a restoration effort to rebuild Breton Island to its pre-Katrina footprint. This data series serves as an archive of processed interferometric swath and single-beam bathymetry data, and side-scan sonar data, collected in the nearshore of Breton and Gosier Islands, NWR, Louisiana. The data were collected during two USGS cruises (USGS Field Activity Numbers 2014-314-FA and 2014-317-FA) in July and August 2014. Geographic information system data products include a 100-meter-cell-size interpolated bathymetry grid, trackline maps, and point data files. Additional files include error analysis maps, Field Activity Collection System logs, and formal Federal Geographic Data Committee metadata. NOTE: These data are scientific in nature and are not to be used for navigation. Any use of trade names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Louisiana↗

Contributions to a shallow aquifer study by reprocessed seismic sections from petroleum exploration surveys, eastern Abu Dhabi, United Arab Emirates

The US Geological Survey, in cooperation with the National Drilling Company of Abu Dhabi, is conducting a 4-year study of the fresh and slightly saline groundwater resources of the eastern Abu Dhabi Emirate. Most of this water occurs in a shallow aquifer, generally less than 150 m deep, in the Al Ain area. A critical part of the Al Ain area coincides with a former petroleum concession area where about 2780 km of vibroseis data were collected along 94 seismic lines during 1981–1983. Field methods, acquistion parameters, and section processing were originally designed to enhance reflections expected at depths ranging from 5000 to 6000 m, and subsurface features directly associated with the shallow aquifer system were deleted from the original seismic sections. The original field tapes from the vibroseis survey were reprocessed in an attempt to extract shallow subsurface information (depths less than 550 m) for investigating the shallow aquifer. A unique sequence of reproccessing parameters was established after reviewing the results from many experimental tests. Many enhancements to the resolution of shallow seismic reflections resulted from: (1) application of a 20-Hz, low-cut filter; (2) recomputation of static corrections to a datum nearer the land surface; (3) intensive velocity analyses; and (4) near-trace muting analyses. The number, resolution, and lateral continuity of shallow reflections were greatly enhanced on the reprocessed sections, as was the delineation of shallow, major faults. Reflections on a synthetic seismogram, created from a borehole drilled to a depth of 786 m on seismic line IQS-11, matched precisely with shallow reflections on the reprocessed section. The 33 reprocessed sections were instrumental in preparing a map showing the major structural features that affect the shallow aquifer system. Analysis of the map provides a better understanding of the effect of these shallow features on the regional occurrence, movement, and quality of groundwater in the concession area. Results from this study demonstrate that original seismic field tapes collected for deep petroleum exploration can be reprocessed to explore for groundwater.

Journal of Applied Geophysics↗

Postaudit of head and transmissivity estimates and ground-water flow models of Avra Valley, Arizona

Ground water from regional alluvial-aquifer systems is the main source of water in the alluvial basins of Arizona, such as Avra Valley. Ground-water flow models are used to assess ground-water availability and the effects of development on the regional ground-water resources. A postaudit of regional-head and transmissivity estimates and the ground-water flow models of Avra Valley was used to evaluate potential errors in the distribution of aquifer properties and recharge that can cause predictive errors in ground-water models. Simulations of predevelopment conditions in 1940 and historical development conditions for 1960-79 provided the basis of comparison for assessing predictive errors of historical conditions for two regional ground-water flow models. Potential errors in the estimation of the regional-head and transmissivity and alternate conceptual models were compared with an existing calibrated two-layer flow model for predevelopment (1940) and developed conditions (1940-85). Measured heads can be subdivided into a north-central region and a region south of the basin constriction. A more variable regional-head surface typical of developed aquifer systems was indicated by kriged developed heads (1985) with about 50 percent more uncertainty than predevelopment heads (1940). Incorporating heads from adjacent basins at the ground-water inflow and outflow regions reduced uncertainty in kriged heads for these boundary areas. Universal cokriging of heads with the strongly correlated land-surface altitudes may improve regional-head estimates and model comparisons where head data are sparse. Local transmissivity estimates can be subdivided into northern and south-central regions that are distributed along the valley axis and the Santa Cruz River. Regional geostatistical estimates of transmissivity, which are based solely on local estimates, are low in the northern part of Avra Valley and are high in the south-central part when compared with the head-conditioned model-derived estimates. These differences may be related to a systematic bias between aquifer-test conditions and methods of aquifer-test analysis. Cokriging transmissivities with specific capacity and silt-and-clay content provided the least uncertainty of all kriged estimates. Predictive errors for the Avra Valley model are the result of a different combination of factors that become significant in the simulation of ground-water flow for the periods representing predevelopment, historical development, and future development conditions. Predictive errors for simulation of predevelopment conditions are caused by potential systematic errors in estimates of local transmissivity, uncertainty in long-term mountain-front recharge, and uncertainty in predevelopment heads along the margins of the basin where recharge and transmissivity estimates are constrained by heads during model calibration. Analog-model historical predictions of future development indicate changes to 1985 were as much as 50 to 100 feet different from actual declines that were caused by errors in the spatial distribution and not the total amount of estimated future pumpage. Predictive errors for simulation of historical development (1960-79) appear to be caused to a greater extent by combined errors in estimates of transmissivity and storage properties and to a lesser extent by estimates of net withdrawal and subsidence. Comparison of the two digital models resulted in differences in transmissivity of as much as 30,000 feet squared per day and differences in specific yield of as much as 0.1. In combination with some differences in net withdrawal, these model-parameter differences resulted in local differences in change in storage of as much as 4,000 acre-feet per square mile and are equivalent to historical predictive errors in water levels of as much as 40 feet. Areas with no differences in model parameters yield comparable simulated water-level declines that are similar to measured declines. The pattern of differences in transmissivity and storage parameters are similar to differences between model-derived estimates conditioned on heads and related geostatistical estimates derived from aquifer-test estimates. A postaudit analysis of alternate conceptual models was explored on the basis of well-by-well comparisons of reductions in mean error and variance, and through the use of standardized calibration-error maps for predevelopment heads (1940) and developed heads (1985). Calibration-error maps provide a useful tool for exploring the spatial structure of model errors and the relative adequacy of model fit that is not available from traditional methods of model comparison. Calibration-error maps indicate estimated heads were too high in the southern part of Avra Valley, and estimated heads were too low in the northern part. Increased transmissivity in the southern part of the lower model layer; decreased hydraulic conductivity in the southwestern part of the upper layer; reduced ground-water inflow from Altar Valley; and increased recharge along the Tortolita Mountains, Tucson Mountains, and Brawley Wash yielded a significantly better model for predevelopment but not for developed conditions (1940-85). This may indicate that alternate conceptual models are different for different time periods or require analysis of time-varying model parameters for developed conditions, such as climatically variable recharge. Predictive errors for future simulations (1986-2025) also could potentially include errors of more than 40 ft from omission of subsidence from the simulation of regional ground-water flow in Avra Valley. Further refinement of the changing conceptual model of an aquifer system under continuing development and variable climate, such as Avra Valley, will require a variety of additional geophysical, geochemical, and hydraulic field data.

Arizona↗

Magnetotelluric and audiomagnetotelluric groundwater survey along the Humu'ula portion of Saddle Road near and around the Pohakuloa Training Area, Hawaii

The Pohakuloa Training Area (PTA), operated by the U.S. Army on the Big Island of Hawaii, is in need of a reliable potable water supply to sustain ongoing operations by staff and trainees. In an effort to acquire baseline hydrologic data with which to develop a plan for providing that water, a series of magnetotelluric (MT) geophysical surveys was performed that spanned the Mauna Loa/Mauna Kea Saddle region of Hawaii Island. These surveys provided electrical resistivity profiles and resistivity maps at several elevations along the axis of the field measurements that can be interpreted to yield information on the depth to the water table. In 2004 a preliminary sequence of 23 audiomagnetotelluric (AMT) soundings was collected along Saddle Road extending from the Waikii Ranch area, west of the PTA, to Department of Hawaiian Home Lands Humu'ula properties east of the Mauna Kea access road. The results of those soundings showed that highly resistive rocks, consistent with dry basalts, were present to depths of at least one kilometer, the maximum depth to which the AMT technique can reliably reach in Hawaii's rocks. A second survey was conducted in 2008 using MT instruments capable of recovering resistivity data to depths of several kilometers below sea level where saturated formations are known to exist. A total of 30 MT soundings was performed along a roughly east to west transect that extended from the (recently acquired) Keamuku PTA lands on the west to as far as the County of Hawaii's upper Kaumana water supply well to the east. Inversion and processing of the field data yielded an electrical cross-section following the Saddle that roughly parallels the geologic contact between the Mauna Kea and Mauna Loa lavas. Several additional electrical sections were constructed normal to the main transect to investigate the three-dimensional nature of the contact. These resistivity data and models suggest that the elevation of saturated rock in places are 400 to 600 meters above mean sea level beneath the surveyed region. Highest elevations for water-saturated zones based upon preferred electrical models are located between training area 3 and training area 6 southwest of training area 4.

Hawai'i↗

Spatial heterogeneity of within-stream methane concentrations

Streams, rivers, and other freshwater features may be significant sources of CH 4 to the atmosphere. However, high spatial and temporal variabilities hinder our ability to understand the underlying processes of CH 4 production and delivery to streams and also challenge the use of scaling approaches across large areas. We studied a stream having high geomorphic variability to assess the underlying scale of CH 4 spatial variability and to examine whether the physical structure of a stream can explain the variation in surface CH 4 . A combination of high-resolution CH 4 mapping, a survey of groundwater CH 4 concentrations, quantitative analysis of methanogen DNA, and sediment CH 4 production potentials illustrates the spatial and geomorphic controls on CH 4 emissions to the atmosphere. We observed significant spatial clustering with high CH 4 concentrations in organic-rich stream reaches and lake transitions. These sites were also enriched in the methane-producing mcrA gene and had highest CH 4 production rates in the laboratory. In contrast, mineral-rich reaches had significantly lower concentrations and had lesser abundances of mcrA . Strong relationships between CH 4 and the physical structure of this aquatic system, along with high spatial variability, suggest that future investigations will benefit from viewing streams as landscapes, as opposed to ecosystems simply embedded in larger terrestrial mosaics. In light of such high spatial variability, we recommend that future workers evaluate stream networks first by using similar spatial tools in order to build effective sampling programs.

Journal of Geophysical Research G: Biogeosciences↗

Large surface-rupture gaps and low surface fault slip of the 2021 Mw 7.4 Maduo earthquake along a low-activity strike-slip fault, Tibetan Plateau

Based on field investigations, interpretations of high-resolution UAV images, and analyses of available InSAR data, we mapped the fault geometry and surface ruptures of the 2021 M w 7.4 Maduo earthquake that occurred on a low-activity strike-slip fault within the Tibetan Plateau. The results indicate that (a) the earthquake activated a fault that is ∼161 km long and has complicated structural geometry; (b) the surface rupture occurs over a distance of 148 km, but is separated into three distinct segments by two large gaps (38 and 20 km, respectively); (c) within the surface-rupture segments, the horizontal and vertical displacements are typically 0.2–2.6 m (much lower than the InSAR-based slip maximum of 2–6 m at depth) and ≤0.4 m, respectively. The two large gaps of the Maduo surface rupture represent the two largest surface-rupture discontinuities of strike-slip earthquakes ever documented, and coincide with structurally complicated fault portions and near-surface soft sediments.

Tibetan Plateau↗

Hydrogeologic setting and ground water flow beneath a section of Indian River Bay, Delaware

The small bays along the Atlantic coast of the Delmarva Peninsula (Delaware, Maryland, and Virginia) are a valuable natural resource, and an asset for commerce and recreation. These coastal bays also are vulnerable to eutrophication from the input of excess nutrients derived from agriculture and other human activities in the watersheds. Ground water discharge may be an appreciable source of fresh water and a transport pathway for nutrients entering the bays. This paper presents results from an investigation of the physical properties of the surficial aquifer and the processes associated with ground water flow beneath Indian River Bay, Delaware. A key aspect of the project was the deployment of a new technology, streaming horizontal resistivity, to map the subsurface distribution of fresh and saline ground water beneath the bay. The resistivity profiles showed complex patterns of ground water flow, modes of mixing, and submarine ground water discharge. Cores, gamma and electromagnetic-induction logs, and in situ ground water samples collected during a coring operation in Indian River Bay verified the interpretation of the resistivity profiles. The shore-parallel resistivity lines show subsurface zones of fresh ground water alternating with zones dominated by the flow of salt water from the estuary down into the aquifer. Advective flow produces plumes of fresh ground water 400 to 600 m wide and 20 m thick that may extend more than 1 km beneath the estuary. Zones of dispersive mixing between fresh and saline ground water develop on the upper, lower, and lateral boundaries of the the plume. the plumes generally underlie small incised valleys that can be traced landward to stream draining the upland. The incised valleys are filled with 1 to 2 m of silt and peat that act as a semiconfining layer to restrict the downward flow of salt water from the estuary. Active circulation of both the fresh and saline ground water masses beneath the bay is inferred from the geophysical results and supported by geochemical data.

Delaware↗

Composition and Distribution of Streambed Sediments in the Penobscot River, Maine, May 1999

Sediment samples were collected and geophysical surveys were run along 50 miles of the Penobscot River, Maine, in the spring of 1999 to produce maps that describe the composition and distribution of streambed sediments for selected areas in the river channel. The objective of the sediment survey was to locate areas along the river where fine-grained, easily transportable sediment types were deposited between Old Town and Medway, Maine. These data can be used to design future sediment-sampling programs to assess the quality of streambed sediments and evaluate the health of the Penobscot River. This report describes the results of the sediment survey and the methods used to collect, analyze, and interpret the data used to create maps of streambed-sediment types in the study area. Deposits of fine-grained sediments (mud and sand) are scattered along the shorelines of the mainland and the islands and at the downstream ends of islands and at the mouths of brooks and streams. The most extensive depositional areas were found in the Mattaseunk Dam impoundment near Medway. The main areas of the river channel consist primarily of gravel, sand, and rock.

Water-Resources Investigations Report↗

A field guide to the central, creeping section of the San Andreas fault and the San Andreas Fault Observatory at Depth

This field trip is along the central section of the San Andreas fault and consists of eight stops that illustrate surface evidence of faulting, in general, and features associated with active fault creep, in particular. Fault creep is slippage along a fault that occurs either in association with small-magnitude earthquakes or without any associated large-magnitude earthquakes. Another aspect of the trip is to highlight where there are multiple fault traces along this section of the San Andreas fault zone in order to gain a better understanding of plate-boundary processes. The first stop is along the Calaveras fault, part of the San Andreas fault system, at a location where evidence of active fault creep is abundant and readily accessible. The stops that follow are along the San Andreas fault and at convenient locations to present and discuss rock types juxtaposed across the fault that have been transported tens to hundreds of kilometers by right-lateral motion along the San Andreas fault. Stops 6 and 7 are examples of recent studies of different aspects of the fault: drilling into the fault at the depth of repeating magnitude (M) 2 earthquakes with the San Andreas Fault Observatory at Depth (SAFOD) and the geological, geophysical, and seismological study of M 6 earthquakes near the town of Parkfield. Along with the eight official stops on this field trip are 12 “rolling stops”—sites of geologic interest that add to the understanding of features and processes in the creeping section of the fault. Many of the rolling stops are located where stopping is difficult to dangerous; some of these sites are not appropriate for large vehicles (buses) or groups; some sites are not appropriate for people at all. We include photographs of or from many of these sites to add to the reader's experience without adding too many stops or hazards to the trip. An extensive set of literature is available for those interested in the San Andreas fault or in the creeping section, in particular. For more scientifically oriented overviews of the fault, see Wallace ( 1990 ) and Irwin ( 1990 ); for a more generalized overview with abundant, colorful illustrations, see Collier ( 1999 ). Although the presence of small sections of the San Andreas fault was known before the great 1906 San Francisco earthquake, it was only after that event and subsequent geologic investigations reported in Lawson ( 1908 ) that showed the fault as a long structure, extending all the way from east of Los Angeles into northern California. Prentice ( 1999 ) described the importance of the 1908 “Lawson report” and how it pivotally influenced the understanding of the San Andreas. Hill ( 1981 ) presented a wonderful introduction to the evolution of thought on the San Andreas. Geologic maps and maps of the most recently active fault trace in the creeping section, or large parts of it, include those by Brown ( 1970 ), Dibblee ( 1971 , 1980 ), and Wagner et al. ( 2002 ); detailed geologic maps are discussed at various stops in this guide. Various aspects of the creeping section of the San Andreas fault have been the focus of many geologic field trips in the past few decades. Guidebooks for some of those trips include those by Gribi ( 1963a , 1963b ), Brabb et al. ( 1966 ), Rogers ( 1969 ), Bucknam and Haller ( 1989 ), Harden et al. ( 2001 ), and Stoffer ( 2005 ). The creeping section of the San Andreas fault zone lies between areas that experienced large-displacement surface breakage during great earthquakes in 1857 and 1906 (Fig. 1 inset). Burford and Harsh ( 1980 ) divided the creeping section into three segments: (1) a northwest section where the creep rate increases to the southeast in step-like increments, (2) a central section where the creep rate is relatively constant at a maximum value of ∼30 mm/yr (∼1.2 in/yr), and (3) a southeast section where the creep rate decreases to the southeast (Fig. 2 ). The rate of slip along the creeping section of the fault zone has been measured using creepmeters, alignment arrays, and laser distance-measuring devices. The aperture of measurements over which these measurements are made ranges from 10 m (∼33 ft) (creepmeters) to 100 m (∼330 ft) (alignment arrays) to kilometers and tens of kilometers (laser measuring devices). Creepmeter and alignment-array measurements are here termed “near-fault” measurements; laser measurements over distances of 1–2 km (∼0.6–1.2 mi) are termed “intermediate-scale” measurements; laser measurements over tens of kilometers (miles) are termed “broadscale” measurements. Comparisons among near-fault, intermediate-scale, and broadscale measurements and geologic maps show that the northwest part of the creeping section of the fault is composed of two narrow zones of active deformation, one along the San Andreas fault and one along the Calaveras-Paicines fault, whereas the central and southeast sections are both composed of a single relatively narrow zone of deformation. The southeast section is transitional to a locked zone southeast of Cholame; a locked fault is one that slips only in association with a moderate to large earthquake. Throughout the creeping section of the San Andreas fault zone, broadscale measurements generally indicate more deformation than near-fault and intermediate-scale measurements, which are in reasonably close agreement except at Monarch Peak (Mustang Ridge), near the center of the creeping section and our Stop 5 ((Figs. 1 ) and 2 ). Features that we see on this trip include offset street curbs, closed depressions (sag ponds), fault scarps (steep slopes formed by movement along a fault), a split and displaced tree, offset fence lines, fresh fractures, and offset road lines (Fig. 3 is a sketch showing some of the landforms that represent deformation by an active fault). We also see evidence of long-term maturity of the San Andreas fault, as indicated by fault features and displaced rock types (Fig. 4 ). Finally, we will visit sites of ongoing research into the processes associated with earthquakes and their effects. Discussions include drilling into the San Andreas fault at the SAFOD drill site and the 2004 Parkfield earthquake and its effects and implications.

California↗

Map of the approximate inland extent of saltwater at the base of the Biscayne aquifer in Miami-Dade County, Florida, 2018

The inland extent of saltwater at the base of the Biscayne aquifer in eastern Miami-Dade County, Florida, was mapped in 2011, and it was mapped in the Model Land Area in 2016. The saltwater interface has continued to move inland in some areas and is now near several active well fields. An updated approximation of the inland extent of saltwater has been created by using data collected during March 8–December 13, 2018, from 111 monitoring wells open to the Biscayne aquifer near its base. Chloride concentrations in water samples from the monitoring wells and bulk conductivity from geophysical logs and measurements of the specific conductance of groundwater were used to approximate the position of the isochlor representing a chloride concentration of 1,000 milligrams per liter (mg/L) at the base of the Biscayne aquifer. An average rate of saltwater interface movement of about 102 meters per year in the Model Land Area along SW 360 Street was estimated from the approximated dates of arrival of the 250-, 500-, and 1,000-mg/L isochlors at wells TPGW-7L (2013–2014) and ACI-MW-05-FS (2017–2018). This estimate assumes that the interface is traveling in a path parallel to an imaginary line connecting the two monitoring wells. Of the 111 wells from which data were used, 80 wells have open intervals of ≤ 4 meters, 20 of the wells have open intervals that range from 4.3 to 39.6 meters, and the lengths of the open intervals could not be determined in 11 wells. Studies have shown that long open intervals might allow water from various depths to mix under ambient or pumped conditions, which in turn could alter the maximum chloride concentration sampled in the well, or it might change the depth at which the maximum specific conductance is measured within a well, relative to its depth in the aquifer. The approximation of the inland extent of the saltwater interface and the estimated rate of movement of the interface are dependent on the quality of existing data. Improved estimates could be obtained by installing uniformly designed monitoring wells in systematic transects extending landward of the advancing saltwater interface. To achieve this goal, Miami-Dade County and some other organizations are routinely adding new monitoring wells with short open intervals and replacing poorly designed or positioned monitoring wells to improve spatial coverage of the network.

Florida↗

The North American upper mantle: Density, composition, and evolution

The upper mantle of North America has been well studied using various seismic methods. Here we investigate the density structure of the North American (NA) upper mantle based on the integrative use of the gravity field and seismic data. The basis of our study is the removal of the gravitational effect of the crust to determine the mantle gravity anomalies. The effect of the crust is removed in three steps by subtracting the gravitational contributions of (1) topography and bathymetry, (2) low-density sedimentary accumulations, and (3) the three-dimensional density structure of the crystalline crust as determined by seismic observations. Information regarding sedimentary accumulations, including thickness and density, are taken from published maps and summaries of borehole measurements of densities; the seismic structure of the crust is based on a recent compilation, with layer densities estimated from P-wave velocities. The resultant mantle gravity anomaly map shows a pronounced negative anomaly (−50 to −400 mGal) beneath western North America and the adjacent oceanic region and positive anomalies (+50 to +350 mGal) east of the NA Cordillera. This pattern reflects the well-known division of North America into the stable eastern region and the tectonically active western region. The close correlation of large-scale features of the mantle anomaly map with those of the topographic map indicates that a significant amount of the topographic uplift in western NA is due to buoyancy in the hot upper mantle, a conclusion supported by previous investigations. To separate the contributions of mantle temperature anomalies from mantle compositional anomalies, we apply an additional correction to the mantle anomaly map for the thermal structure of the uppermost mantle. The thermal model is based on the conversion of seismic shear-wave velocities to temperature and is consistent with mantle temperatures that are independently estimated from heat flow and heat production data. The thermally corrected mantle density map reveals density anomalies that are chiefly due to compositional variations. These compositional density anomalies cause gravitational anomalies that reach ~250 mGal. A pronounced negative anomaly (−50 to −200 mGal) is found over the Canadian shield, which is consistent with chemical depletion and a corresponding low density of the lithospheric mantle, also referred to as the mantle tectosphere. The strongest positive anomaly is coincident with the Gulf of Mexico and indicates a positive density anomaly in the upper mantle, possibly an eclogite layer that has caused subsidence in the Gulf. Two linear positive anomalies are also seen south of 40°N: one with a NE-SW trend in the eastern United States, roughly coincident with the Grenville-Appalachians, and a second with a NW-SE trend beneath the states of Texas, New Mexico, and Colorado. These anomalies are interpreted as being due to (1) the presence of remnants of an oceanic slab in the upper mantle beneath the Grenville-Appalachian suture and (2) mantle thickening caused by a period of shallow, flat subduction during the Laramie orogeny, respectively. Based on these geophysical results, the evolution of the NA upper mantle is depicted in a series of maps and cartoons that display the primary processes that have formed and modified the NA crust and lithospheric upper mantle.

Journal of Geophysical Research B: Solid Earth↗