USGS Science⌕ Search

SEARCH · USGS Science

Results for “Interpretation”

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 973 records · Page 54Linked to original sources

Saudi Arabian seismic deep-refraction profiles; final project report

In February 1978 a seismic deep-refraction profile was recorded by the U.S. Geological Survey along a 1000-km line across the Arabian Shield in western Saudi Arabia. The line begins in Mesozoic cover rocks near Riyadh on the Arabian Platform, leads southwesterly across three major Precambrian tectonic provinces, traverses Cenozoic rocks of the coastal plain near Jizan (Tihamat-Asir), and terminates at the outer edge of the Farasan Bank in the southern Red Sea. More than 500 surveyed recording sites were occupied, including 19 in the Farasan Islands. Six shot points were used: five on land, with most charges placed below the water table in drill holes, and one at sea, with charges placed on the sea floor and detonated from a ship. Slightly more than 61 metric tons of explosives were used in 19 discrete firings. Seismic energy was recorded by 100 newly-developed portable seismic stations deployed in approximately 200 km-long arrays for each firing. Each station consisted of a standard 2-Hz vertical component geophone coupled to a self-contained analog recording instrument equipped with a magnetic-tape cassette. In this final report, we fully document the field and data-processing procedures and present the final seismogram data set as both a digital magnetic tape and as record sections for each shot point. Record sections include a normalized set of seismograms, reduced at 6 km/s, and a true-amplitude set, reduced at 8 km/s, which have been adjusted for amplifier gain, individual shot size, and distance from the shot point. Appendices give recorder station and shot information, digital data set descriptions, computer program listings, arrival times used in the interpretation, and a bibliography of reports published as a result of this project. We used two-dimensional ray-tracing techniques in the data analysis, and our interpretation is based primarily on horizontally layered models. The Arabian Shield is composed, to first-order, of two layers, each about 20 km thick, with average velocities of 6.3 km/s and 7.0 km/s, respectively. At the western shield margin the crust thins to less than 20 km total thickness, beyond which the Red Sea shelf and coastal plain are interpreted to be underlain by oceanic crust. A major crustal lateral velocity inhomogeneity northeast of Sabhah in the Shammar Tectonic Province is interpreted as the suture zone of two crustal blocks of different composition. Several high-velocity anomalies in the upper crust correlate with mapped gneissic dome structures. Two intra-crustal reflectors at13 km depth are interpreted as the tops of mafic intrusives. The Mohorovicic discontinuity beneath the shield varies from 43 km depth in the northeast with 8.2 km/s mantle velocity to 38 km depth in the southwest with 8.0 km/s mantle velocity. Two velocity discontinuities are identified in the upper mantle, at 59 and 70 km depth. We suggest further work, including refined analyses of the data employing filtering and synthetic seismogram techniques, as well as consideration of attenuation properties. Extension of the seismic refraction profile to the Arabian Gulf and some short profiles perpendicular to the existing profile would be fruitful areas for future field work.

Open-File Report↗

Two hundred years of magma transport and storage at Kīlauea Volcano, Hawai'i, 1790-2008

This publication summarizes the evolution of the internal plumbing of Kīlauea Volcano on the Island of Hawaiʻi from the first documented eruption in 1790 to the explosive eruption of March 2008 in Halemaʻumaʻu Crater. For the period before the founding of the Hawaiian Volcano Observatory in 1912, we rely on written observations of eruptive activity, earthquake swarms, and periodic draining of magma from the lava lake present in Kīlauea Caldera. After 1912 the written observations are supplemented by continuous measurement of tilting of the ground at Kīlauea’s summit and by a continuous instrumental record of earthquakes, both measurements made during 1912–56 by a single pendulum seismometer housed on the northeast edge of Kīlauea’s summit. Interpretations become more robust following the installation of seismic and deformation networks in the 1960s. A major advance in the 1990s was the ability to continuously record and telemeter ground deformation to allow its precise correlation with seismic activity before and after eruptions, intrusions, and large earthquakes. We interpret specific events in Kīlauea’s 200- year written history as steps in a broad transition from summit lava-lake activity in Kīlauea Caldera to shield building on the east rift zone. The ability of the magmatic plumbing to deliver magma to eruption is critical to the history of eruption and intrusion. When the rate of magma supply equals the rate of eruption, there is little ground deformation or intrusion. When the magma supply rate is greater than the rate of eruption, then the edifice responds through any or all of summit inflation, intrusion, increased spreading rate, and large flank earthquakes. In Kīlauea’s 200-year history we identify three regions of the volcano in which magma is stored and supplied from below. Source 1 is at 1-km depth or less beneath Kīlauea’s summit and fed Kīlauea’s summit lava lakes throughout most of the 19th century and again from 1907 to 1924. Source 1 was used up in the series of small Halemaʻumaʻu eruptions following the end of lava-lake activity in the summit collapse of 1924. Source 2 is the magma reservoir at a depth of 2–6 km beneath Kīlauea’s summit that has been imaged by seismic and deformation measurements beginning in the 1960s. This source was first identified in the summit collapses of 1922 and 1924. Source 3 is a diffuse volume of magma-permeated rock between 5 and 11 km depth beneath the east rift zone and above the near-horizontal decollement at the base of the Kīlauea edifice. Magma distribution within source 2 has been derived by combining petrologic study of the three chemically uniform summit eruptions of 1952, 1961, and 1967–68 and the east rift eruptions within this interval with both observation of migrating centers of inflation determined from leveling surveys conducted before the 1967–68 eruption and with published models of expected deformation from different source geometries. We adopt a model of concatenated magmatic plugs with nodes beneath the inflation centers. Addition of erupted and intruded volumes of the three summit magma batches yields a liquid magma volume of about 0.2 km3, with dimensions of ~1 km by 1 km by 200 m centered at about 3-km depth within source 2. Following the Halemaʻumaʻu eruption of 1967–68, the chemistry of magma coming into Kīlauea’s summit reservoir has changed frequently, and during the eruption that began in 1983, chemical changes have been subtle and continuous. In this period we interpret changes in chemistry as related to an increase in magma supply resulting from increased partial melting in an expanding mantle source volume. We know from instrumental recording of eruptions since the long Halemaʻumaʻu eruption in 1952 that stress in the edifice accumulates as magma is added underground and is relieved by eruption and by dilation of the rift zones associated with seaward movement (spreading) of Kīlauea’s south flank. During and after the last half of the 20th century, magma transfer to the rift zone has dominantly occurred from source 2. High rates of flank motion have been correlated with high rates of endogenous growth; alternatively, lower rates of motion have characterized periods when the underground magmatic plumbing was being refilled following lateral removal of magma, as well as periods when a more open magmatic plumbing favored continuous eruption. Since at least 1952, source 3 has not drained during deflations, which was apparently not the case before 1924. Triangulation and leveling conducted in 1912, 1921, and 1926, combined with post-1912 tilt measurements, identified a broad regional uplift in 1918–19 and an equally broad collapse in 1924, neither of which has been seen since. We associate these elevation changes with addition or subtraction of magma from all three magma sources, dominantly source 3. We interpret the intrusion beneath the east rift zone during the 1924 collapse to have stabilized the rift zone-south flank relationship, preventing loss of magma from source 3 in subsequent collapses. Rates of seaward spreading were low until 1952, when earthquakes in 1950 and 1951 associated with surges of magma from the hotspot triggered a large offshore south flank earthquake swarm that unlocked the south flank and enabled a greatly increased rate of seaward spreading. Magma supply rates have been derived for the entire period of study. Between 1823 and 1840, magma was supplied from source 1 at a very high rate of more than 0.2 km 3 /yr, which we interpret as recovery from a substantial draining of magma from beneath Kīlauea in 1790. Inferred magma supply rates diminished to one-tenth of that value after 1840, in part because of increase in the activity of Mauna Loa beginning in 1843. Magma supply rates between 1918 and 1924 were about 0.024 km 3 /yr, matching that of the period from 1840 to 1894. During 1950–52 the magma supply rate increased to about 0.06 km 3 /yr, in part because of the great reduction in Mauna Loa activity following its large eruption in June 1950. Following the summit eruption of 1967–68, magma supply increased further to ~0.1 km 3 /yr, and further increases to more than 0.2 km 3 /yr occurred during the east rift eruption that began in 1983. Eruption at Kīlauea’s summit took place in 1952, and eruptive activity steadily increased as increased magma supply also drove increased spreading rates. The inability of magma supply to be accommodated by a combination of eruption and spreading during the 1969–74 Mauna Ulu period stressed Kīlauea’s south flank. The stress was relieved in part by the M7.2 earthquake of 29 November 1975. That earthquake, in turn, dilated Kīlauea’s east rift zone as the south flank moved seaward, producing a favorable condition for continuous east rift eruption, which began in 1983. The 1975 earthquake also resulted in the ability of the south flank to move independently under the influence of gravity, effectively decoupling the spreading rate from changes in the magma supply rate. The continuing increase in magma supply after 1983 was instead manifested in rift dilation, increased intrusion, and ultimately in the launching of a second eruption in Halemaʻumaʻu in March 2008, the first instance in Kīlauea’s recorded history of simultaneous eruption at the summit and on the east rift zone. Kīlauea’s history can be considered in cycles of equilibrium, crisis, and recovery. The approach of a crisis is driven by a magma supply rate that greatly exceeds the capacity of the plumbing to deliver magma to the surface. Crises can be anticipated by inflation measured at Kīlauea’s summit coupled with an increase in overall seismicity, particularly manifest by intrusion and eruption in the southwest sector of the volcano. Unfortunately the nature of the crisis—for example, large earthquake, new eruption, or edifice-changing intrusion—cannot be specified ahead of time. We conclude that Kīlauea’s cycles are controlled by nonlinear dynamics, which underscores the difficulty in predicting eruptions and earthquakes.

Hawai'i↗

Geophysical Characterization of the Quaternary-Cretaceous Contact Using Surface Resistivity Methods in Franklin and Webster Counties, South-Central Nebraska

To help manage and understand the Platte River system in Nebraska, the Platte River Cooperative Hydrology Study (COHYST), a group of state and local governmental agencies, developed a regional ground-water model. The southern boundary of this model lies along the Republican River, where an area with insufficient geologic data immediately north of the Republican River led to problems in the conceptualization of the simulated flow system and to potential problems with calibration of the simulation. Geologic descriptions from a group of test holes drilled in south-central Nebraska during 2001 and 2002 indicated a possible hydrologic disconnection between the Quaternary-age alluvial deposits in the uplands and those in the Republican River lowland. This disconnection was observed near a topographic high in the Cretaceous-age Niobrara Formation, which is the local bedrock. In 2003, the U.S. Geological Survey, in cooperation with the COHYST, collected surface geophysical data near these test holes to better define this discontinuity. Two-dimensional imaging methods for direct-current resistivity and capacitively coupled resistivity were used to define the subsurface distribution of resistivity along several county roads near Riverton and Inavale, Nebraska. The relation between the subsurface distribution of resistivity and geology was defined by comparing existing geologic descriptions of test holes to surface-geophysical resistivity data along two profiles and using the information gained from these comparisons to interpret the remaining four profiles. In all of the resistivity profile sections, there was generally a three-layer subsurface interpretation, with a resistor located between two conductors. Further comparison of geologic data with the geophysical data and with surficial features was used to identify a topographic high in the Niobrara Formation near the Franklin Canal which was coincident with a resistivity high. Electrical properties of the Niobrara Formation made accurate interpretation of the resistivity profile sections difficult and less confident because of similar resistivity of this formation and that of the coarser-grained sediment of the Quaternary-age deposits. However, distinct conductive features were identified within the resistivity profile sections that aided in delineating the contact between the resistive Quaternary-age deposits and the resistive Niobrara Formation. Using this information, an interpretive boundary was drawn on the resistivity profile sections to represent the contact between the Quaternary-age alluvial deposits and the Cretaceous-age Niobrara Formation. A digital elevation model (DEM) of the top of the Niobrara Formation was constructed using the altitudes from the interpreted contact lines. This DEM showed that the general trend of top of the Niobrara Formation dips to the southeast. At the north edge of the study site, the Niobrara Formation topographic high trends east-west with an altitude range of 559 meters in the west to 543 meters in the east. Based on the land-surface elevation and the Niobrara Formation DEM, the estimated thickness of the Quaternary-age alluvial deposits throughout the study area was mapped and showed a thinning of the Quaternary-age alluvial deposits to the north, approximately where the topographic high of the Niobrara Formation is located. This topographic high in the Niobrara Formation has the potential to act as a barrier to ground-water flow from the uplands alluvial aquifer to the Republican River alluvial aquifer as shown in the resistivity profile sections. The Quaternary-age alluvial deposits in the uplands and those in the Republican River Valley are not fully represented as disconnected because it is possible that there are ground-water flow paths that were not mapped during this study.

Scientific Investigations Report↗

Three-dimensional geologic mapping of the Cenozoic basin fill, Amargosa Desert basin, Nevada and California

Understanding the subsurface geologic framework of the Cenozoic basin fill that underlies the Amargosa Desert in southern Nevada and southeastern California has been improved by using borehole data to construct three-dimensional lithologic and interpreted facies models. Lithologic data from 210 boreholes from a 20-kilometer (km) by 90-km area were reduced to a limited suite of descriptors based on geologic knowledge of the basin and distributed in three-dimensional space using interpolation methods. The resulting lithologic model of the Amargosa Desert basin portrays a complex system of interfingered coarse- to fine-grained alluvium, playa and palustrine deposits, eolian sands, and interbedded volcanic units. Lithologic units could not be represented in the model as a stacked stratigraphic sequence due to the complex interfingering of lithologic units and the absence of available time-stratigraphic markers. Instead, lithologic units were grouped into interpreted genetic classes, such as playa or alluvial fan, to create a three-dimensional model of the interpreted facies data. Three-dimensional facies models computed from these data portray the alluvial infilling of a tectonically formed basin with intermittent internal drainage and localized regional groundwater discharge. The lithologic and interpreted facies models compare favorably to resistivity, aeromagnetic, and geologic map data, lending confidence to the interpretation.

California;Nevada↗

Mineralogical applications of electron diffraction. 1. Theory and techniques

The small wavelengths used in electron-diffraction experiments and the thinness of the crystals necessary for the transmission of the electron beam combine to require a somewhat different diffraction geometry for the interpretation of electron-diffraction patterns than is used in the interpretation of X-ray diffraction patterns. This geometry, based on the reciprocal lattice concept and geometrical construction of Ewald, needed for the interpretation. of transmission electron-diffraction single-crystal patterns is here reviewed. Transmission electron-diffraction single-crystal patterns of two monoclinic substances, colemanite [CaB 3 O 3 (OH) 3 •H 2 O] and potassium chlorate (KC10 3 ), are examined and the .theory necessary for their interpretation is given in detail. The study of these patterns furnishes a basis for the interpretation of single-crystal patterns of materials belonging to any crystal system. It is shown that useful unit-cell data, accurate to a few tenths of a percent, can be obtained from the patterns of colemanite and KClO 3 . A method of evaluating unit-cell data from measurements of such single-crystal patterns is given. The transmission electron-diffraction powder pattern obtained from an oriented aggregate of thin crystals gives the same unit-cell data as are given by the electron-diffraction single-crystal pattern obtained from one crystal of the aggregate., A graphical method is given for precisely evaluating unit-cell constants from measurements of such a powder pattern.

Trace Elements Investigations↗

Borehole-geophysical investigation of the University of Connecticut landfill, Storrs, Connecticut

A borehole-geophysical investigation was conducted to help characterize the hydrogeology of the fractured-rock aquifer and the distribution of unconsolidated glacial deposits near the former landfill and chemical waste-disposal pits at the University of Connecticut in Storrs, Connecticut. Eight bedrock boreholes near the landfill and three abandoned domestic wells located nearby were logged using conventional and advanced borehole-geophysical methods from June to October 1999. The conventional geophysical-logging methods included caliper, gamma, fluid temperature, fluid resistivity, and electromagnetic induction. The advanced methods included deviation, optical and acoustic imaging of the borehole wall, heat-pulse flowmeter, and directional radar reflection. Twenty-one shallow piezometers (less than 50-feet deep) were logged with gamma and electromagnetic induction tools to delineate unconsolidated glacial deposits. Five additional shallow bedrock wells were logged with conventional video camera, caliper, electromagnetic induction, and fluid resistivity and temperature tools. The rock type, foliation, and fracturing of the site were characterized from high-resolution optical-televiewer (OTV) images of rocks penetrated by the boreholes. The rocks are interpreted as fine- to medium-grained quartz-feldsparbiotite-garnet gneiss and schist with local intrusions of quartz diorite and pegmatite and minor concentrations of sulfide mineralization similar to rocks described as the Bigelow Brook Formation on regional geologic maps. Layers containing high concentrations of sulfide minerals appear as high electrical conductivity zones on electromagneticinduction and borehole-radar logs. Foliation in the rocks generally strikes to the northeast-southwest and dips to the west, consistent with local outcrop observations. The orientation of foliation and small-scale gneissic layering in the rocks, however, varies locally and with depth in some of the boreholes. In two of the boreholes, the foliation strikes predominantly to the northwest and dips to the northeast. Although small-scale faults and lithologic discontinuities were observed in the OTV data, no large-scale faults were observed that appear on regional geologic maps. Fractures were located and characterized through the use of conventional geophysical, OTV, acoustic-televiewer (ATV), and borehole-radar logs. The orientation of fractures varies considerably across the site; some fractures are parallel to the foliation, whereas others cross-cut the foliation. Many of the transmissive fractures in the bedrock boreholes strike about N170°E and N320°E with dips of less than 45°. Other transmissive fractures strike about N60°E with dips of more than 60°. Most of the transmissive fractures in the domestic wells strike about N60°E and N22°E with dips of more than 45°. The strike of N60°E is parallel to the trend of a thrust fault that appears on regional geologic maps. Vertical flow in the boreholes was measured with the heat-pulse flowmeter under ambient and (or) pumping conditions. Results of ATV, OTV, and conventional logs were used to locate specific zones for flowmeter testing. Ambient downflow was measured in three boreholes, ambient upflow was measured in two other boreholes, and both ambient downflow and upflow were measured in a sixth borehole. The other five bedrock boreholes and domestic wells did not have measurable vertical flow. The highest rate of ambient flow was measured in the background borehole in which upflow and downflow converged and exited the borehole at a fracture zone near a depth of 62 feet. Ambient flow of about 340 gallons per day was measured. In the other five wells, ambient flow of about 20 to 35 gallons per day was measured. Under low-rate pumping (0.25 to 1 gallon per minute), one to six inflow zones were identified in each well. Usually the fractures that are active under ambient conditions contribute to the well under pumping conditions. To prevent ambient vertical flow and the potential for cross-contamination, temporary borehole liners were installed in five of the boreholes. Specific-capacity and open-hole transmissivity values were determined in eight boreholes completed in bedrock. The specific capacity estimated for these boreholes ranges from 0.14 to 1.6 gallons per minute per foot. The values for open-hole transmissivity range over two orders of magnitude and when proportioned to individual fracture transmissivity, range from 23 to 340 feet squared per day. Two boreholes had been drilled to intersect electrically conductive zones identified by previous surface-geophysical investigations. The borehole-geophysical results indicate that the boreholes penetrate electrically conductive structures consistent with the anomalies interpreted from the surface-geophysical data. Borehole MW121R was located to intersect a dipping electrically conductive anomaly at about 60 feet, interpreted from the two-dimensional direct currentresistivity survey conducted on the western side of the landfill. The electromagnetic-conductivity log in the borehole contains a high electrical conductivity anomaly at a depth of 69 feet. The magnitude of this anomaly is nearly 10,000 millisiemens per meter and is coincident with a layer containing sulfide mineralization, rather than fractures. The other borehole, MW105R, was located to intersect another anomaly south of the landfill. This anomaly was interpreted as a north-south striking, westward dipping feature. In the borehole, two south-striking, westward dipping fractures were identified in the ATV, OTV, and radar logs. The specific conductance of the fluid measured near these fractures was as high as 1,250 microsiemens per centimeter. Water-quality samples collected in October 1999 from an isolated zone from 71.5 to 76.5 feet indicated high specific conductance (810 microsiemens per centimeter), high concentrations of iron and cadmium, negative oxidation-reduction potential, and chlorobenzene. Collectively, these parameters indicate that the high specific conductance in the borehole logs for MW105R was caused by landfill leachate. Therefore, the anomaly identified by boreholeand surface-geophysical surveys is interpreted as a conductive lithologic feature and a permeable fracture zone that contains landfill leachate.

Connecticut↗

Integrated use of surface geophysical methods for site characterization — A case study in North Kingstown, Rhode Island

A suite of complementary, non‐invasive surface geophysical methods was used to assess their utility for site characterization in a pilot investigation at a former defense site in North Kingstown, Rhode Island. The methods included frequency‐domain electromagnetics (FDEM), ground‐penetrating radar (GPR), electrical resistivity tomography (ERT), and multi‐channel analysis of surface‐wave (MASW) seismic. The results of each method were compared to each other and to drive‐point data from the site. FDEM was used as a reconnaissance method to assess buried utilities and anthropogenic structures; to identify near‐surface changes in water chemistry related to conductive leachate from road‐salt storage; and to investigate a resistive signature possibly caused by groundwater discharge. Shallow anomalies observed in the GPR and ERT data were caused by near‐surface infrastructure and were consistent with anomalies observed in the FDEM data. Several parabolic reflectors were observed in the upper part of the GPR profiles, and a fairly continuous reflector that was interpreted as bedrock could be traced across the lower part of the profiles. MASW seismic data showed a sharp break in shear wave velocity at depth, which was interpreted as the overburden/bedrock interface. The MASW profile indicates the presence of a trough in the bedrock surface in the same location where the ERT data indicate lateral variations in resistivity. Depths to bedrock interpreted from the ERT, MASW, and GPR profiles were similar and consistent with the depths of refusal identified in the direct‐push wells. The interpretations of data collected using the individual methods yielded non‐unique solutions with considerable uncertainty. Integrated interpretation of the electrical, electromagnetic, and seismic geophysical profiles produced a more consistent and unique estimation of depth to bedrock that is consistent with ground‐truth data at the site. This test case shows that using complementary techniques that measure different properties can be more effective for site characterization than a single‐method investigation.

Rhode Island↗

The South Fork detachment fault, Park County, Wyoming: discussion and reply ( USA).

Blackstone (1985) published an interpretation of South form detachment fault and related features. His interpretation of the area between Castle and Hardpan transverse faults is identical to mine of 1941. Subsequent detailed mapping has shown that the structure between the transverse faults is more complicated than originally envisioned and resurrected by Blackstone. The present paper describes and discusses geologic features that are the basis for my interpretations; also discussed are differences between my interpretations and those of Blackstone. Most data are shown on the geologic map of the Wapiti Quadrangle (Pierce and Nelson, 1969). Blackstone's 'allochthonous' masses are part of the South Form fault. Occurrences of Sundance Formation, which he interpreted as the upper plate of his 'North Fork fault', are related to Heart Mountain fault. Volcaniclastic rocks south of Jim Mountain mapped as Aycross Formation by Torres and Gingerich may be Cathedral Cliffs Formation, emplaced by movement of the Heart Mountain fault. - Author

Contributions to Geology - University of Wyoming, ↗

Oceans on Mars: An assessment of the observational evidence and possible fate

If the large Late Hesperian outflow channels were eroded by extensive floods, as appears likely, then large bodies of water must have once occupied the northern plains during that period. Previous estimates of the sizes of bodies of water in the northern lowlands range up to 3 ?? 108 km3. Several contacts have been previously mapped around the edges of the northern plains and interpreted to be shorelines remaining from these former standing bodies of water. We examine the elevations and geologic relations along these contacts in detail and find little support for their interpretation as shorelines. Some contacts are clearly of volcanic origin, and all have significant variations in elevation. Better support for the former presence of water over large parts of the northern plains is provided by the Vastitas Borealis Formation (VBF). Most of the post-Noachian fill within the northern basin is ridged plains of Lower Hesperian age, interpreted to be volcanic in origin. Overlying the ridged plains is the VBF, a thin veneer of material of Upper Hesperian age. The VBF may have been deposited from large floods. Support for this interpretation is the similarity in age between the outflow channels and the VBF, the presence of the VBF at the lower ends of the outflow channels, and identification of numerous features in the outcrop areas of the VBF that are suggestive of basal melting of an ice sheet. To cover all the area over which the VBF is exposed would require ???2.3 ?? 107 km3 of water. Spread over the entire surface of Mars, this volume is equal to a global layer (global equivalent layer, or GEL) ???156 m deep. We find no support for the larger estimates of ocean volumes that range up to 3 ?? 108 km3 and which imply comparable amounts of water per unit area as are currently present on the surface of the Earth. Under present climatic conditions on Mars an ocean would freeze in a geologically short time period (???104 years), then would sublimate away at rates strongly dependent on the presence or absence of debris on the ice surface. The present VBF is interpreted as a sublimation residue from the ponded outflow channel effluents. The fate of a volume of water thought to have been emplaced by the outflow channels (???2.3 ?? 107 km3) is largely accounted for by the presence of other existing reservoirs on the planet. An approximately 20-30 m GEL of water is estimated to be in the present polar caps, and a 50 m GEL may have escaped to space since the Hesperian, leaving ???80 m GEL unaccounted for. This amount may be partly trapped in other volatile-rich deposits on the surface, and a significant amount could have reentered the groundwater system by south polar basal melting and been progressively cold-trapped at the base of a growing cryosphere. On the basis of our assessment of the Hesperian-aged deposits, we predict that testing of the Clifford and Parker [2001] hypothesis that a Noachian-aged ocean covered up to one third of the surface of Mars will be made very difficult by the enhanced degradation rates in the Noachian and subsequent geological events in the northern lowlands.

Journal of Geophysical Research E: Planets↗

Modeling the mesozoic-cenozoic structural evolution of east texas

The U.S. Geological Survey (USGS) recently assessed the undiscovered technically recoverable oil and gas resources within Jurassic and Cretaceous strata of the onshore coastal plain and State waters of the U.S. Gulf Coast. Regional 2D seismic lines for key parts of the Gulf Coast basin were interpreted in order to examine the evolution of structural traps and the burial history of petroleum source rocks. Interpretation and structural modeling of seismic lines from eastern Texas provide insights into the structural evolution of this part of the Gulf of Mexico basin. Since completing the assessment, the USGS has acquired additional regional seismic lines in east Texas; interpretation of these new lines, which extend from the Texas-Oklahoma state line to the Gulf Coast shoreline, show how some of the region's prominent structural elements (e.g., the Talco and Mount Enterprise fault zones, the East Texas salt basin, and the Houston diapir province) vary along strike. The interpretations also indicate that unexplored structures may lie beneath the current drilling floor. Structural restorations based upon interpretation of these lines illustrate the evolution of key structures and show the genetic relation between structural growth and movement of the Jurassic Louann Salt. 1D thermal models that integrate kinetics and burial histories were also created for the region's two primary petroleum source rocks, the Oxfordian Smackover Formation and the Cenomanian-Turonian Eagle Ford Shale. Integrating results from the thermal models with the structural restorations provides insights into the distribution and timing of petroleum expulsion from the Smackover Formation and Eagle Ford Shale in eastern Texas.

Texas↗

Seismic source characterization in central and eastern United States

This report documents earthquake geology studies conducted in the central and eastern United States (CEUS) between 2018 and 2023. The overarching goal of the earthquake geology studies was to identify and characterize fault sources in the New Madrid seismic zone, the eastern Tennessee seismic zone, and the Charleston seismic zone. These studies primarily relied on the interpretation of high-resolution topography, morphotectonic analysis of 10-meter (m) digital elevation models, field reconnaissance, and integration with new and existing subsurface datasets. Results from each study are placed into a seismic hazard framework. The first two chapters focus on the New Madrid seismic zone, which is the most seismically active region in the CEUS, to better characterize the southern extent of the Reelfoot fault in the Obion River valley, and identify and characterize proposed faults along the margins of Crowleys Ridge in the center of the Mississippi River embayment. Crowleys Ridge is a ~320-kilometer (km) long landform that has been proposed to be formed by intermittent faulting or erosion from the Mississippi River. Legacy seismic reflection data documented Eocene, Mesozoic, and older strata offset by faults on the margins of the ridge (Van Arsdale et al., 1999). A morphotectonic analysis of the topography along Crowleys Ridge indicated southward increases in catchment hypsometric integral and slope. Neotectonic mapping on high-resolution lidar data revealed a series of subparallel linear fault and fold scarps along the margins of Crowleys Ridge. These scarps offset previously mapped fluvial and alluvial surfaces that were <56 thousand years ago (ka) (Rittenour et al., 2007). Integration of new seismic reflection and airborne electromagnetic data, legacy seismic data, with the landscape morphotectonic analysis and neotectonic mapping supports the interpretation that Crowleys Ridge is bound by faults that have been active in the late Quaternary but do not experience modern seismicity. By placing these newly characterized faults into the regional framework, the Crowleys Ridge faults can be interpreted as thrust faults associated with stepovers in the overall dextral fault system, similar to the Reelfoot fault. The results of these studies indicate that seismicity has migrated eastward during the late Quaternary, with the Reelfoot fault accommodating much of the modern-day strain. The Reelfoot fault ruptured in the 1811-1812 New Madrid earthquake sequence (Fuller, 1912) and has a record of past large earthquakes (Tuttle et al., 2002; 2019). Whereas the northern extent of the fault intersects the Mississippi River and has well-expressed surface deformation, questions remained regarding the southern limit of seismicity, deformation, and overall fault length. Analysis of <25 ka river terraces along the Obion River valley in western Tennessee revealed subtle fold scarps, with increased magnitudes of deformation on progressively older terraces, indicative of a longer record of fault movement than has been documented by paleoseismic and paleoliquefaction studies alone (e.g., Kelson et al., 1996; Gold et al., 2019; Tuttle et al., 2019). The third chapter integrates surface and subsurface data to assess the landscape record of surface deformation in the eastern Tennessee seismic zone (ETSZ), which is the second most seismically active region in the CEUS. Seismicity in the ETSZ is deep (5-26 km) and focused in the Proterozoic rock below the Paleozoic detachment. Although the region has experienced moderate (≤ Mw 4.8) seismicity in the instrumental record, the lack of a large historical earthquake and limited paleoseismic and paleoliquefaction evidence of prehistoric large ground shaking events have made seismic hazard characterization of the region challenging. Two models have emerged to predict the style and orientation of active surface deformation in the ETSZ. One suggests that surface rupturing earthquakes would exploit pre-existing faults along the regional structural grain (Cox et al., 2022). The second model advocates for primarily strike-slip motion on east-west or north-south faults, following analyses of recent seismicity (Chapman et al., 1997; Dunn and Chapman, 2006; Daniels and Peng, 2022) and the modern stress field (Levandowski et al., 2018). In our analysis, neotectonic mapping of high-resolution lidar data reveals a concentration of lineaments with east-west orientations, with some corresponding to previously mapped east-west faults that crosscut the northeast-southwest regional structural grain from the Paleozoic orogeny. A morphotectonic analysis of catchments and river segments within the same lithology indicates subtle differences that may be indicative of a longer-term tectonic uplift signal. These changes in morphotectonic metrics spatially correspond to newly mapped lineaments and previously mapped east-west trending faults. Within a regional framework, we suggest that diffuse surface deformation associated with deep seismicity is accommodated on a network on east-west faults. However, further work is needed to better understand the potential late Pleistocene fault activity of these lineaments, and we propose that the ETSZ is still best characterized as an area source in seismic hazard models. The final chapter presents preliminary work to create a geographic information system (GIS) database of recent studies analyzing surface and subsurface datasets and interpretations in the Charleston seismic zone, South Carolina. This GIS database will serve as a foundation for future work to analyze new quality level 1 (QL1; <0.5 m resolution) lidar data over the 1886 magnitude (M) 7 Charleston epicentral region to identify possible fault source sources responsible for the 1886 or other surface-rupturing events. This chapter briefly summarizes three new subsurface datasets, including seismic reflection (Pratt et al., 2022; Liberty, 2022) and seismicity data (Chapman et al., 2016), and one new surface dataset (Marple and Hurd, 2020), and their interpretations of potentially active lineaments and faults in the region.

central and eastern United States↗

Ida and Dactyl: Spectral reflectance and color variations

Galileo SSI color data between 0.4 and 1.0 ??m demonstrate that both Ida and Dactyl are S-type asteroids with similar, but distinct spectra. Small but definite color variations are also observed on Ida itself and involve both the blue part of the spectrum and the depth of the 1-??m pyroxene-olivine band. Ida's surface can be classified into two color terrains: Terrain A has a shallower 1-??m absorption and a steeper visible red slope than does Terrain B. Qualitatively, the color-albedo systematics of these two terrains follow those noted for color units on Gaspra and the variations in 1-??m band depth with weathering described by Gaffey et al. (Gaffey, M. J., J. F. Bell, R. H. Brown, T. H. Burbine, J. Piatek, K. L. Reed, and D. A. Chaky 1993. Icarus 106, 573-602). Terrain A, with its slightly lower albedo, its shallower 1-??m band, and its slightly steeper visible red slope relative to Terrain B could be interpreted as the "more processed," "more mature," or the "more weathered" of the two terrains. Consistent with this interpretation is that Terrain A appears to be the ubiquitous background on most of Ida, while Terrain B is correlated with some small craters as well as with possible ejecta from the 10-km Azzurra impact structure. Because of these trends, it is less likely that differences between Terrains A and B are caused by an original compositional inhomogeneity within the body of Ida, although they do fall within the range known to occur within the Koronis family. The spectrum of Dactyl is similar to, but definitely different from, that of Terrain B on Ida. It does not conform to the pattern that obtains between the colors and albedos of Terrains A and B: the satellite's 1-??m band is deeper than that of Terrain B, but its albedo is lower, rather than higher. By itself, the deeper band depth could be interpreted, following Gaffey et al., to mean that Dactyl is a less weathered version of Terrain B on Ida, but such an interpretation is at odds with Dactyl's redder spectral slope. Thus, the explanation for the color difference between Dactyl and Ida is likely to be different from that which accounts for the differences between the two terrains on Ida. Given that Dactyl and Ida have very similar photometric properties (Helfenstein, P., J. Veverka, P. C. Thomas, D. P. Simonelli, K. Klassen, T. V. Johnson, F. Fanale, J. Granahan, A. S. McEwen, M. J. S. Belton, and C. R. Chapman 1996 Icarus 120, 48-65), thus ruling out any dramatic texture differences between the two surfaces, the most likely explanation is that the satellite has a slightly different composition (more pyroxene?) than Ida. The spectral difference is within the range reported by Binzel et al. (Binzel, R. P., S. Xu, and S. J. Bus 1993. Icarus 106, 608-611.) for members of the Koronis family, and could be caused by compositional inhomogeneities of the Koronis parent body rather than by post-breakup regolith processes. ?? 1996 Academic Press, Inc.

Icarus↗

Interactive visualization to advance earthquake simulation

The geological sciences are challenged to manage and interpret increasing volumes of data as observations and simulations increase in size and complexity. For example, simulations of earthquake-related processes typically generate complex, time-varying data sets in two or more dimensions. To facilitate interpretation and analysis of these data sets, evaluate the underlying models, and to drive future calculations, we have developed methods of interactive visualization with a special focus on using immersive virtual reality (VR) environments to interact with models of Earth's surface and interior. Virtual mapping tools allow virtual "field studies" in inaccessible regions. Interactive tools allow us to manipulate shapes in order to construct models of geological features for geodynamic models, while feature extraction tools support quantitative measurement of structures that emerge from numerical simulation or field observations, thereby enabling us to improve our interpretation of the dynamical processes that drive earthquakes. VR has traditionally been used primarily as a presentation tool, albeit with active navigation through data. Reaping the full intellectual benefits of immersive VR as a tool for scientific analysis requires building on the method's strengths, that is, using both 3D perception and interaction with observed or simulated data. This approach also takes advantage of the specialized skills of geological scientists who are trained to interpret, the often limited, geological and geophysical data available from field observations. ?? Birkhaueser 2008.

Pure and Applied Geophysics↗

Patch-reef morphology as a proxy for Holocene sea-level variability, Northern Florida Keys, USA

A portion of the northern Florida Keys reef tract was mapped with the NASA Experimental Advanced Airborne Research Lidar (EAARL) and the morphology of patch reefs was related to variations in Holocene sea level. Following creation of a lidar digital elevation model (DEM), geospatial analyses delineated morphologic attributes of 1,034 patch reefs (reef depth, basal area, height, volume, and topographic complexity). Morphometric analysis revealed two morphologically different populations of patch reefs associated with two distinct depth intervals above and below a water depth of 7.7 m. Compared to shallow reefs, the deep reefs were smaller in area and volume and showed no trend in topographic complexity relative to water depth. Shallow reefs were more variable in area and volume and became flatter and less topographically complex with decreasing water depth. The knoll-like morphology of deep reefs was interpreted as consistent with steady and relatively rapidly rising early Holocene sea level that restricted the lateral growth of reefs. The morphology of shallow 'pancake-shaped' reefs at the highest platform elevations was interpreted as consistent with fluctuating sea level during the late Holocene. Although the ultimate cause for the morphometric depth trends remains open to interpretation, these interpretations are compatible with a recent eustatic sea-level curve that hindcasts fluctuating late Holocene sea level. Thus it is suggested that the morphologic differences represent two stages of reef accretion that occurred during different sea-level conditions. ?? 2008 Springer-Verlag.

Coral Reefs↗

Geoelectrical structure of the central zone of Piton de la Fournaise volcano (Reunion)

A study of the geoelectrical structure of the central part of Piton de la Fournaise volcano (Reunion, Indian Ocean) was made using direct current electrical (DC) and transient electromagnetic soundings (TEM). Piton de la Fournaise is a highly active oceanic basaltic shield and has been active for more than half a million years. Joint interpretation of the DC and TEM data allows us to obtain reliable 1D models of the resistivity distribution. The depth of investigation is of the order of 1.5 km but varies with the resistivity pattern encountered at each sounding. Two-dimensional resistivity cross sections were constructed by interpolation between the soundings of the 1D interpreted models. Conductors with resistivities less than 100 ohm-m are present at depth beneath all of the soundings and are located high in the volcanic edifice at elevations between 2000 and 1200 m. The deepest conductor has a resistivity less than 20 ohm-m for soundings located inside the Enclos and less than 60-100 ohm-m for soundings outside the Enclos. From the resistivity distributions, two zones are distinguished: (a) the central zone of the Enclos; and (b) the outer zone beyond the Enclos. Beneath the highly active summit area, the conductor rises to within a few hundred meters of the surface. This bulge coincides with a 2000-mV self-potential anomaly. Low-resistivity zones are inferred to show the presence of a hydrothermal system where alteration by steam and hot water has lowered the resistivity of the rocks. Farther from the summit, but inside the Enclos the depth to the conductive layers increases to approximately 1 km and is inferred to be a deepening of the hydrothermally altered zone. Outside of the Enclos, the nature of the deep, conductive layers is not established. The observed resistivities suggest the presence of hydrated minerals, which could be found in landslide breccias, in hydrothermally altered zones, or in thick pyroclastic layers. Such formations often create perched water tables. The known occurrence of large eastward-moving landslides in the evolution of Piton de la Fournaise strongly suggests that large volumes of breccias should exist in the interior of the volcano; however, extensive breccia deposits are not observed at the bottom of the deep valleys that incise the volcano to elevations lower than those determined for the top of the conductors. The presence of the center of Piton de la Fournaise beneath the Plaine des Sables area during earlier volcanic stages (ca. 0.5 to 0.150 Ma) may have resulted in broad hydrothermal alteration of this zone. However, this interpretation cannot account for the low resistivities in peripheral zones. It is not presently possible to discriminate between these general interpretations. In addition, the nature of the deep conductors may be different in each zone. Whatever the geologic nature of these conductive layers, their presence indicates a major change of lithology at depth, unexpected for a shield volcano such as Piton de la Fournaise.

Bulletin of Volcanology↗

Geochemistry and the understanding of ground-water systems

Geochemistry has contributed significantly to the understanding of ground-water systems over the last 50 years. Historic advances include development of the hydrochemical facies concept, application of equilibrium theory, investigation of redox processes, and radiocarbon dating. Other hydrochemical concepts, tools, and techniques have helped elucidate mechanisms of flow and transport in ground-water systems, and have helped unlock an archive of paleoenvironmental information. Hydrochemical and isotopic information can be used to interpret the origin and mode of ground-water recharge, refine estimates of time scales of recharge and ground-water flow, decipher reactive processes, provide paleohydrological information, and calibrate ground-water flow models. Progress needs to be made in obtaining representative samples. Improvements are needed in the interpretation of the information obtained, and in the construction and interpretation of numerical models utilizing hydrochemical data. The best approach will ensure an optimized iterative process between field data collection and analysis, interpretation, and the application of forward, inverse, and statistical modeling tools. Advances are anticipated from microbiological investigations, the characterization of natural organics, isotopic fingerprinting, applications of dissolved gas measurements, and the fields of reaction kinetics and coupled processes. A thermodynamic perspective is offered that could facilitate the comparison and understanding of the multiple physical, chemical, and biological processes affecting ground-water systems.

Hydrogeology Journal↗

Geochemistry of a Tertiary sedimentary phosphate deposit: Baja California Sur, Mexico

The San Gregorio Formation in Baja California Sur, a phosphate-enriched sedimentary unit of late Oligocene to early Miocene age, has been analyzed in two areas (La Purisima and San Hilario) for its chemical composition (major oxides, Cu, Cd, Cr, Co, V, and rare-earth elements - REE) and isotopic composition (??18O and ??13C). A detrital and a marine component were determined from major oxides. The detrital component consists of an unaltered volcanic-ash fraction and a terrigenous clay-silt fraction. The marine component, which accumulated initially as biogenic and hydrogenous material, is now present as opal-A, opal-CT, CaCO3, organic matter, and an authigenic phosphate fraction, mostly pelletal and composed of the carbonate-fluorapatite mineral francolite. The minor elements have been partitioned into these components by assuming a constant composition for the two detrital fractions. The composition of the marine component of minor elements can then be interpreted by assuming that the stoichiometry of the original accumulating organic matter was equal to that of modern plankton. The Cu and Cd contents in the marine component of all rocks require that the seawater-derived fractions of these two metals were supplied to the seafloor solely by organic matter. Enrichments of Cr and V at both sites required an additional marine input. On the basis of their geochemistry in the modern ocean, Cr and V could have precipitated, or been adsorbed, onto settling particles from an O2 minimum zone in which the O2 content was low enough to promote denitrification rather than oxygen respiration. An enrichment of the REE, now within the apatite fraction, resulted from their adsorption onto particulates also in the O2 minimum zone and to the dissolution and alteration of biogenic phases (predominantly silica) within the sediment. Co and Fe2O3 show no enrichment above a detrital contribution. The ??18O-values of apatites from the La Purisima site are heavier than those of apatites from the San Hilario site, whereas the ??13C-values show the opposite trend. One possible interpretation of these variations is that ??18O reflects seawater values and ??13C sediment pore water values. This interpretation suggests that upwelling rates and primary productivity within the water column were greater at La Purisima, an interpretation that is corroborated by a greater abundance of apatite measured in outcrop at La Purisima. The Ce anomalies of the phosphate-enriched samples also differ between the two sites, indicating that they also recorded water masses, similar to the ??18O-values. ?? 1991.

Chemical Geology↗

Further evidence for a low U/Pb source in the moon: U-Th-Pb, Sm-Nd, and Ar-Ar, isotopic systematics of lunar meteorite Yamato-793169

The coarse-grained lunar meteorites, Yamato-793169 and Asuka-881757, represent a new type of low-Ti mare basalt. This paper reports the results of a U-Th-Pb, Sm-Nd, and Ar-Ar isotopic study of Yamato-793169 performed as part of a consortium studies of lunar basaltic meteorites. The isotopic study was carried out on a small sample (100 mg) so that only three density fractions could be separated. These fractions were leached with dilute acid in order to eliminate terrestrial Pb contamination. However, the leaching procedure did not completely remove this contamination in some fractions and also apparently caused a fractionation of elements (U, Th, Pb) due to preferential leaching effects, producing a secondary disturbance of the systematics. Furthermore, the Ar-Ar analyses indicate that the isotopic systematics in this meteorite might have been disturbed sometime later than 750 Ma. For these reasons, the ages obtained using different isotopic systems disagree with each other and a precise formation age could not be obtained for this meteorite. However, using what results there are, two reasonable interpretations can be made: (1) the Sm-Nd system yielded an age of 3.4 Ga that could be interpreted as the formation age, assuming that this system is possibly the least disturbed during the metamorphic event(s) that this meteorite experienced at least once, and that all other isotopic systems and their corresponding ages were disturbed, and (2) the U-Pb system yielded a nearly concordant age of 3.8 Ga that could be interpreted as the formation age and the Sm-Nd isotopic systematics were somehow disturbed. We prefer the second interpretation for the reasons discussed below. On a concordia diagram, the CDT (Cañon Diablo troilite)-corrected U-Pb isotopic data yield a discordia line similar to the lunar catastrophic array, indicating that the source of the meteorite formed during early lunar differentiation (∼4.4 Ga) and that the basalt was generated near 3.9 Ga. Total Ar-Ar age on plagioclase is 3.26 Ga, which seems to be too old if the formation age is 3.4 Ga, because low temperature fractions lost large amounts of radiogenic 40 Ar during the late thermal event. If we assume a formation age of 3.8 Ga, the estimated source 238 U/ 204 Pb (μ) is 21.6 ± 3.5 and ϵ Nd is 3.9 ± 0.3. These results indicate that the source of Yamato-793169 is more depleted than Apollo 12 and 15 LT basalts, but less depleted than Asuka-881757. Therefore, Yamato-793169 may represent a new type of LT- or VLT-like mare basalt that is different from Asuka-881757. The wide variety of lead and neodymium isotopic characteristics among LT and VLT mare basalts indicate that the lunar mantle was very heterogeneous with respect to trace element abundances.

Geochimica et Cosmochimica Acta↗