SEARCH · Geology Reports

Results for “Journal of Geophysical Research – Earth Surface”

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

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23Linked to original sources

USGS2017

Refining the formation and early evolution of the Eastern North American Margin: New insights from multiscale magnetic anomaly analyses

To investigate the oceanic lithosphere formation and early seafloor spreading history of the North Atlantic Ocean, we examine multiscale magnetic anomaly data from the Jurassic/Early Cretaceous age Eastern North American Margin (ENAM) between 31 and 40°N. We integrate newly acquired sea surface magnetic anomaly and seismic reflection data with publicly available aeromagnetic and composite magnetic anomaly grids, satellite-derived gravity anomaly, and satellite-derived and shipboard bathymetry data. We evaluate these data sets to (1) refine magnetic anomaly correlations throughout the ENAM and assign updated ages and chron numbers to M0–M25 and eight pre-M25 anomalies; (2) identify five correlatable magnetic anomalies between the East Coast Magnetic Anomaly (ECMA) and Blake Spur Magnetic Anomaly (BSMA), which may document the earliest Atlantic seafloor spreading or synrift magmatism; (3) suggest preexisting margin structure and rifting segmentation may have influenced the seafloor spreading regimes in the Atlantic Jurassic Quiet Zone (JQZ); (4) suggest that, if the BSMA source is oceanic crust, the BSMA may be M series magnetic anomaly M42 (~168.5 Ma); (5) examine the along and across margin variation in seafloor spreading rates and spreading center orientations from the BSMA to M25, suggesting asymmetric crustal accretion accommodated the straightening of the ridge from the bend in the ECMA to the more linear M25; and (6) observe anomalously high-amplitude magnetic anomalies near the Hudson Fan, which may be related to a short-lived propagating rift segment that could have helped accommodate the crustal alignment during the early Atlantic opening.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

A dynamic balance between magma supply and eruption rate at Kilauea volcano, Hawaii

The dynamic balance between magma supply and vent output at Kilauea volcano is used to estimate both the volume of magma stored within Kilauea volcano and its magma supply rate. Throughout most of 1991 a linear decline in volume flux from the Kupaianaha vent on Kilauea's east rift zone was associated with a parabolic variation in the elevation of Kilauea's summit as vent output initially exceeded then lagged behind the magma supply to the volcano. The correspondence between summit elevation and tilt established with over 30 years of data provided daily estimates of summit elevation in terms of summit tilt. The minimum in the parabolic variation in summit tilt and elevation (or zero elevation change) occurs when the magma supply to the reservoir from below the volcano equals the magma output from the reservoir to the surface, so that the magma supply rate is given by vent flux on that day. The measurements of vent flux and tilt establish that the magma supply rate to Kilauea volcano on June 19, 1991, was 217,000±10,000 m 3 /d (or 0.079±0.004 km 3 /yr). This is close to the average eruptive rate of 0.08 km 3 /yr between 1958 and 1984. In addition, the predictable response of summit elevation and tilt to each east rift zone eruption near Puu Oo since 1983 shows that summit deformation is also a measure of magma reservoir pressure. Given this, the correlation between the elevation of the Puu Oo lava lake (4 km uprift of Kupaianaha and 18 km from the summit) and summit tilt provides an estimate for magma pressure changes corresponding to summit tilt changes. The ratio of the change in volume to the change in reservoir pressure ( dV / dP ) during vent activity may be determined by dividing the ratio of volume erupted to change in summit tilt ( dV / d tilt ) by the ratio of pressure change to change in summit tilt ( dP / d tilt ). This measure of dV / dP , when combined with laboratory measurements of the bulk modulus of tholeitic melt, provides an estimate of 240±50 km 3 for the volume of Kilauea's magma reservoir. This estimate is much larger than traditional estimates but consistent with seismic tomographic imaging and geophysical modeling of Kilauea's magma system.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

The Moon Mineralogy Mapper (M3) imaging spectrometer for lunar science: Instrument description, calibration, on‐orbit measurements, science data calibration and on‐orbit validation

The NASA Discovery Moon Mineralogy Mapper imaging spectrometer was selected to pursue a wide range of science objectives requiring measurement of composition at fine spatial scales over the full lunar surface. To pursue these objectives, a broad spectral range imaging spectrometer with high uniformity and high signal-to-noise ratio capable of measuring compositionally diagnostic spectral absorption features from a wide variety of known and possible lunar materials was required. For this purpose the Moon Mineralogy Mapper imaging spectrometer was designed and developed that measures the spectral range from 430 to 3000 nm with 10 nm spectral sampling through a 24 degree field of view with 0.7 milliradian spatial sampling. The instrument has a signal-to-noise ratio of greater than 400 for the specified equatorial reference radiance and greater than 100 for the polar reference radiance. The spectral cross-track uniformity is >90% and spectral instantaneous field-of-view uniformity is >90%. The Moon Mineralogy Mapper was launched on Chandrayaan-1 on the 22nd of October. On the 18th of November 2008 the Moon Mineralogy Mapper was turned on and collected a first light data set within 24 h. During this early checkout period and throughout the mission the spacecraft thermal environment and orbital parameters varied more than expected and placed operational and data quality constraints on the measurements. On the 29th of August 2009, spacecraft communication was lost. Over the course of the flight mission 1542 downlinked data sets were acquired that provide coverage of more than 95% of the lunar surface. An end-to-end science data calibration system was developed and all measurements have been passed through this system and delivered to the Planetary Data System (PDS.NASA.GOV). An extensive effort has been undertaken by the science team to validate the Moon Mineralogy Mapper science measurements in the context of the mission objectives. A focused spectral, radiometric, spatial, and uniformity validation effort has been pursued with selected data sets including an Earth-view data set. With this effort an initial validation of the on-orbit performance of the imaging spectrometer has been achieved, including validation of the cross-track spectral uniformity and spectral instantaneous field of view uniformity. The Moon Mineralogy Mapper is the first imaging spectrometer to measure a data set of this kind at the Moon. These calibrated science measurements are being used to address the full set of science goals and objectives for this mission.

Journal of Geophysical Research E: Planets↗

USGS2022

Thermophysical and compositional properties of paleobedforms on Mars

Bedforms on Earth and Mars are often preserved in the rock record in the form of sedimentary rock with distinct cross-bedding. On rare occasions, the full-surface geometry of a bedform can be preserved through burial and lithification. These features, known as paleobedforms, are found in a variety of geographic locations on Mars. Evidence in the morphology of paleobedforms, such as the retention of impact craters and steep erosional scarps, suggests that these features are well-lithified and capable of withstanding prolonged weathering and erosion. Here, we present results from thermophysical and compositional analyses on a subset of the best preserved paleobedform candidate fields on Mars. Thermophysical modeling elucidates the changes these bedforms underwent from their unconsolidated, particulate nature to their currently observed properties. Certain paleobedforms have elevated thermal inertias (e.g., ∼300–500 J·m −2 ·s −1/2 ·K −1 ) when compared with modern bedforms (∼250 J·m −2 ·s −1/2 ·K −1 ), and modeling indicates that they have cement volumes of 0.8%–1.5% even as high as 30%. However, most paleobedform candidates have unexpectedly low thermal inertia when compared with modern dunes. Additionally, compositional analyses reveal a range of spectral characteristics within paleobedforms (e.g., primary and secondary alteration products). These features add to the already existing class of Martian surfaces in which thermal inertia does not seem to correspond to erodibility, cohesion, or mechanical strength. Studying paleobedforms with both raised and nonraised thermal inertia has provided new insights into lithification on Mars and constrained the environmental conditions leading to the formation of these enigmatic features.

Journal of Geophysical Research - Planets↗

USGS2012

Rapid extension recorded by cooling‐age patterns and brittle deformation, Naxos, Greece

The metamorphic core complex exposed as the island of Naxos in the Aegean Sea records an unusually complete sequence of structures developed as a result of continental extension. The structures formed during Miocene rise and cooling from ductile, upper amphibolite facies and anatectic conditions to brittle near‐surface conditions beneath the Naxos detachment fault. Top‐to‐the‐north ductile fabrics in the footwall, which initially developed during amphibolite facies prograde metamorphism, were overprinted by a succession of north directed, normal sense lower‐temperature brittle structures as the footwall was tectonically unloaded and unroofed. Pseudotachylite and cataclasite formation, brittle faulting, alteration, and erosion of the footwall occurred during continued slip and tectonic denudation. Neogene conglomerate and megabreccia, in part derived from exhumation of this footwall, lie structurally above the peripheral Naxos fault. Published K‐Ar and 40 Ar/ 39 Ar ages for hornblende, white mica, and biotite in the footwall decrease northwestward; apparent ages 17–50 Ma in the southeast correspond to areas of low metamorphic grade where preextension argon was partially retained. Published ages 16–10 Ma in higher‐grade rocks of the domal core in the north are cooling ages that for each of the three minerals show a component of younging in the NNE direction of extension. Assuming this is the direction of unroofing, we interpret the rate of this younging as the fault slip rate as the footwall rocks moved >20 km SSW relative to their hanging wall along the base of the Naxos detachment fault. The calculated rates of slip average 5–8 mm/yr, comparable to maximum rates reported in the Basin and Range province.

Naxos↗

USGS2012

Lead and strontium isotopes and related trace elements as genetic tracers in the Upper Cenozoic rhyolite-basalt association of the Yellowstone Plateau volcanic field

Supported by various field geologic and petrologic data, the contents of Pb, U, Th, Rb, and Sr and the isotopic compositions of Pb and Sr for upper Cenozoic volcanic rocks of the Yellowstone Plateau volcanic field are consistent with the hypothesis of derivation of the basaltic and rhyolitic magmas by partial melting of distinct source regions in the upper mantle and lower crust, respectively. All the basalt samples analyzed but one have systematically lower values of 207 Pb/ 204 Pb and 87 Sr/ 86 Sr than the rhyolites. The values of 206 Pb/ 204 Pb are smaller, and 87 Sr/ 86 Sr are mostly larger than known values in oceanic basalts. In all but one case, the values of 207 Pb/ 204 Pb are higher than expected from an extrapolation of known values in oceanic basalts to less radiogenic values of 206 Pb/ 204 Pb. Because there are no xenoliths, phenocrysts are only moderate to sparse in abundance, REE patterns are low and flat at the radiogenic end of lead isotopic compositions, several values of Rb/Sr are low, and 80% of the basalt samples form a well-developed secondary isochron separate from the rhyolites, we favor an interpretation for basalt genesis wherein isotopic signatures of most mafic magmas were attained in a continental ‘keel’ of mantlelike character about 2.6 b.y. old or somewhat older attached to the crust, and these signatures were unaltered by magma passage through the crust. At the very least, the current data continue to cast serious doubt as to the inevitability of crustal contamination for basaltic magma intruding the continental environment and postulate that much can be learned about the mantle under continents through the study of continental basalts. One basalt unit with an unusually low value of 207 Pb/ 204 Pb and an 87 Ar/ 86 Ar less than 0.704 may represent subcontinental ‘keel’-derived magma that rose unaltered to the surface. Our data also are not consistent with formation of this rhyolite-basalt association primarily by such processes as crystal fractionation, separation of immiscible silicate liquids from a common parental magma, or fractional melting of a homogeneous source. Rather as a conceptual model, we envision large mafic intrusions to have been injected into the lower crust resulting in rhyolite generation through partial anatexis of the adjacent wall rocks which probably had a 206 Pb/ 204 Pb < 17 and 87 Sr/ 86 Sr > 0.709; a model that has much in common with that proposed by Holmes (1931). All the other hypotheses listed have the necessary added complication that either the basalt or the rhyolite or both become contaminated after the two magma types separated, have problems accounting for the lack of igneous rocks of intermediate compositions or production of such large volumes of rhyolitic material (∼5000 km 3 ), and fail to explain why rhyolitic magma is not a more common occurrence in the ocean basin. We appeal to bouyancy of rhyolites to generate a barrier for basalt magma migration and account for the great preponderance of rhyolite relative to basalt at the surface. Furthermore, the complex isotopic picture in the rhyolites indicates that many of these magmas interacted with the upper crustal geologic units that they traversed. The interactions involved diverse processes, probably including reacton with hydrothermal fluids or hydrothermally altered rocks at high levels as well as by contamination with Phanerozoic sedimentary and Precambrian crystalline rocks at deeper levels. At the very least, we feel our study adds a cautionary note to the currently increasingly popular hypothesis that differentiation of basalt or gabbro magmas to rhyolite or granite (as distinct from tonalite or dacite) is a common occurrence and is therefore an important continential building process. Models for formation of rhyolite and granite predominantly by reworking of crust (anatexis) must still be considered. The primitive Archean mantle of the region was characterized by higher Rb/Sr, U/Pb, and Th/U values than are typical of modern suboceanic mantle. The mantle residuum within the continental subcrustal lithosperic ‘keel’ that resulted from the Archean crustal differentiation event probably was depleted in Rb/Sr and U/Pb, and the crust was correspondingly enriched in these ratios. The crust probably was further differentiated by an Archean high-grade metamorphism, during or after the primary event, into a granulitic lower crust depleted in U/Pb and Rb/Sr and a lower-grade upper crust enriched in these ratios.

Journal of Geophysical Research Solid Earth↗

USGS2007

Hydrogeologic controls imposed by mechanical stratigraphy in layered rocks of the Chateauguay River Basin, a U.S.-Canada transborder aquifer

[1] The Châteauguay River Basin delineates a transborder watershed with roughly half of its surface area located in northern New York State and half in southern Québec Province, Canada. As part of a multidisciplinary study designed to characterize the hydrogeologic properties of this basin, geophysical logs were obtained in 12 wells strategically located to penetrate the four major sedimentary rock formations that constitute the regional aquifers. The layered rocks were classified according to their elastic properties into three primary units: soft sandstone, hard sandstone, and dolostone. Downhole measurements were analyzed to identify fracture patterns associated with each unit and to evaluate their role in controlling groundwater flow. Fracture networks are composed of orthogonal sets of laterally extensive, subhorizontal bedding plane partings and bed-delimited, subvertical joints with spacings that are consistent with rock mechanics principles and stress models. The vertical distribution of transmissive zones is confined to a few select bedding plane fractures, with soft sandstone having the fewest (one per 70-m depth) and hard sandstone the most (five per 70-m depth). Bed-normal permeability is examined using a probabilistic model that considers the lengths of flow paths winding along joints and bedding plane fractures. Soft sandstone has the smallest bed-normal permeability primarily because of its wide, geomechanically undersaturated joint spacing. Results indicate that the three formations have similar values of bulk transmissivity, within roughly an order of magnitude, but that each rock unit has its own unique system of groundwater flow paths that constitute that transmissivity.

Chateauguay River Basin↗

USGS2005

Slicing up the San Francisco Bay Area: Block kinematics and fault slip rates from GPS-derived surface velocities

Observations of surface deformation allow us to determine the kinematics of faults in the San Francisco Bay Area. We present the Bay Area velocity unification (BA??VU??, "bay view"), a compilation of over 200 horizontal surface velocities computed from campaign-style and continuous Global Positioning System (GPS) observations from 1993 to 2003. We interpret this interseismic velocity field using a three-dimensional block model to determine the relative contributions of block motion, elastic strain accumulation, and shallow aseismic creep. The total relative motion between the Pacific plate and the rigid Sierra Nevada/Great Valley (SNGV) microplate is 37.9 ?? 0.6 mm yr-1 directed toward N30.4??W ?? 0.8?? at San Francisco (??2??). Fault slip rates from our preferred model are typically within the error bounds of geologic estimates but provide a better fit to geodetic data (notable right-lateral slip rates in mm yr-1: San Gregorio fault, 2.4 ?? 1.0; West Napa fault, 4.0 ?? 3.0; zone of faulting along the eastern margin of the Coast Range, 5.4 ?? 1.0; and Mount Diablo thrust, 3.9 ?? 1.0 of reverse slip and 4.0 ?? 0.2 of right-lateral strike slip). Slip on the northern Calaveras is partitioned between both the West Napa and Concord/ Green Valley fault systems. The total convergence across the Bay Area is negligible. Poles of rotation for Bay Area blocks progress systematically from the North America-Pacific to North America-SNGV poles. The resulting present-day relative motion cannot explain the strike of most Bay Area faults, but fault strike does loosely correlate with inferred plate motions at the time each fault initiated. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Strain accumulation across the Prince William Sound asperity, Southcentral Alaska

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

Journal of Geophysical Research B: Solid Earth↗

USGS2018

Findings and lessons learned from the assessment of the Mexico-United States transboundary San Pedro and Santa Cruz aquifers: The utility of social science in applied hydrologic research

Study Region This study region encompasses the Transboundary San Pedro and Santa Cruz aquifers which are shared between the states of Sonora (Mexico) and Arizona (US). Special regional considerations include a semi-arid climate, basin-fill aquifers with predominantly montane recharge areas, economic drivers in the mining, trade, and military sectors, groundwater-dependent cities with expanding cones of depression, interbasin groundwater transfers, ground- and surface-water contamination, and protected aquatic and riparian habitats that act as significant migration corridors for hundreds of species, including some that are threatened and endangered. Study Focus We focus on lessons learned from the hydrologic assessment of the Transboundary San Pedro and Santa Cruz aquifers. We conducted the work, in two phases: (1) laying the groundwork and (2) implementation. The “laying the groundwork” phase consisted of binational meetings with stakeholders and key actors (agencies and individuals), and the development of an understanding of the physical, institutional, historical, and socio-political context. This led to signing of the binational Transboundary Aquifer Assessment Program (TAAP) agreement in 2009 and detailed the process for cooperation and coordination in the assessment of shared aquifers. The implementation phase began with an agreement to proceed with the study of four “focus” aquifers (Santa Cruz, San Pedro, Mesilla (Conejos-Médanos in Mexico), and Hueco Bolson (Bolsón del Hueco in Mexico)) and development of associated technical teams. Though we do include a brief discussion of the lessons learned from the physical science portion of the study, the results have been described and published elsewhere. The bulk of the paper instead focuses on the findings and lessons learned from the integration of social-science perspectives into a largely physical-science based program, since there is a growing recognition of the need for this type of approach especially in the management and assessment of transboundary aquifers. New Hydrological Insights for the Region The Sonora-Arizona effort succeeded because both countries were adequately represented, and because of flexibility of skills and ability of teams comprising both university and government scientists. Teams included social and earth scientists. Including the social sciences was critical to research design and implementation, and to addressing the cultural, institutional, and socio-political contexts of transboundary aquifer assessment. Significant components of the continuing implementation phase include strategic planning, data compilation and analysis, cross-border integration of datasets, geophysical and geochemical surveys, and internal, peer, and stakeholder engagement.

Arizona, Sonora↗

USGS2022

Quantifying large-scale surface change using SAR amplitude images: Crater morphology changes during the 2019-2020 Shishaldin Volcano eruption

Morphological processes often induce meter-scale elevation changes. When a volcano erupts, tracking such processes provides insights into the style and evolution of eruptive activity and related hazards. Compared to optical remote-sensing products, synthetic aperture radar (SAR) observes surface change during inclement weather and at night. Differential SAR interferometry estimates phase change between SAR acquisitions and is commonly applied to quantify deformation. However, large deformation or other coherence loss can limit its use. We develop a new approach applicable when repeated digital elevation models (DEMs) cannot be otherwise retrieved. Assuming an isotropic radar cross-section, we estimate meter-scale vertical morphological change directly from SAR amplitude images via an optimization method that utilizes a high-quality DEM. We verify our implementation through simulation of a collapse feature that we modulate onto topography. We simulate radar effects and recover the simulated collapse. To validate our method, we estimate elevation changes from TerraSAR-X stripmap images for the 2011–2012 eruption of Mount Cleveland. Our results reproduce those from two previous studies; one that used the same dataset, and another based on thermal satellite data. By applying this method to the 2019–2020 eruption of Shishaldin Volcano, Alaska, we generate elevation change time series from dozens of co-registered TerraSAR-X high-resolution spotlight images. Our results quantify previously unresolved cone growth in November 2019, collapses associated with explosions in December–January, and further changes in crater elevations into spring 2020. This method can be used to track meter-scale morphology changes for ongoing eruptions with low latency as SAR imagery becomes available.

Alaska↗

USGS2012

Fluid expulsion sites on the Cascadia accretionary prism: mapping diagenetic deposits with processed GLORIA imagery

Point-discharge fluid expulsion on accretionary prisms is commonly indicated by diagenetic deposition of calcium carbonate cements and gas hydrates in near-surface (<10 m below seafloor; mbsf) hemipelagic sediment. The contrasting clastic and diagenetic lithologies should be apparent in side scan images. However, sonar also responds to variations in bottom slope, so unprocessed images mix topographic and lithologic information. We have processed GLORIA imagery from the Oregon continental margin to remove topographic effects. A synthetic side scan image was created initially from Sea Beam bathymetric data and then was subtracted iteratively from the original GLORIA data until topographic features disappeared. The residual image contains high-amplitude backscattering that we attribute to diagenetic deposits associated with fluid discharge, based on submersible mapping, Ocean Drilling Program drilling, and collected samples. Diagenetic deposits are concentrated (1) near an out-of-sequence thrust fault on the second ridge landward of the base of the continental slope, (2) along zones characterized by deep-seated strikeslip faults that cut transversely across the margin, and (3) in undeformed Cascadia Basin deposits which overlie incipient thrust faults seaward of the toe of the prism. There is no evidence of diagenetic deposition associated with the frontal thrust that rises from the dècollement. If the dècollement is an important aquifer, apparently the fluids are passed either to the strike-slip faults which intersect the dècollement or to the incipient faults in Cascadia Basin for expulsion. Diagenetic deposits seaward of the prism toe probably consist dominantly of gas hydrates.

Cascadia Basin↗

USGS2014

Ice sheet load cycling and fluid underpressures in the Eastern Michigan Basin, Ontario, Canada

Strong fluid underpressures have been detected in Paleozoic strata in the eastern Michigan Basin, with hydraulic heads reaching ~400 m below land surface (~4 MPa underpressure) and ~200 m below sea level in strata where unusually low permeabilities (~10 −20 –10 −23  m 2 ) were measured in situ. Multiple glaciations, including three with as much as 3 km of ice cover at the site in the last 120 ka, suggest a causal link with the underpressures. We examined this possibility using a one-dimensional groundwater flow model incorporating mechanical loading from both ice weight and lithospheric flexure. Because hydrologic and mechanical changes during glaciation are not well characterized and subsurface properties are imperfectly known, the model was used inversely to estimate flexural loads and loosely constrained permeabilities by matching observed pressures. Acceptable matches were obtained for a surprisingly wide range of scenarios with permeabilities close to measured values and plausible flexural loads. Matches were not obtained when too many parameters were preselected, or when permeabilities were constrained to be significantly larger than measured values. In successful model runs groundwater expulsion under glacial-mechanical loads caused the underpressuring, and flexural loads were important if aquifer and sub-glacial pressures were significantly elevated during glaciation. Simulated fluid pressures in the low-permeability strata fluctuated by 30–40 MPa during glacial cycles but resulted in advective transport of only tens of meters or less. Although other mechanisms cannot be ruled out, we conclude that glacial-mechanical forcing of a water-saturated system can explain the observed underpressures.

Ontario↗

USGS2012

Heat flow in vapor dominated areas of the Yellowstone Plateau volcanic field: implications for the thermal budget of the Yellowstone Caldera

Characterizing the vigor of magmatic activity in Yellowstone requires knowledge of the mechanisms and rates of heat transport between magma and the ground surface. We present results from a heat flow study in two vapor dominated, acid-sulfate thermal areas in the Yellowstone Caldera, the 0.11 km 2 Obsidian Pool Thermal Area (OPTA) and the 0.25 km 2 Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 ± 0.4 MW and 7.5 ± 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 ± 0.3 MW and 1.2 ± 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 ± 1.4 MW. These estimates result in a total heat output of 11.8 ± 1.4 MW and 8.8 ± 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. Extrapolation of the average heat flux from the OPTA (103 ± 2 W·m −2 ) and SPTA (35 ± 3 W·m −2 ) to the ~35 km 2 of vapor dominated areas in Yellowstone yields 3.6 and 1.2 GW, respectively, which is less than the total heat output transported by steam from the Yellowstone Caldera as estimated by the chloride inventory method (4.0 to 8.0 GW).

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 1. Strong motions

We develop a methodology to perform finite fault source inversions from strong motion data using Green's functions (GFs) calculated for a three-dimensional (3-D) velocity structure. The 3-D GFs are calculated numerically by inserting body forces at each of the strong motion sites and then recording the resulting strains along the target fault surface. Using reciprocity, these GFs can be recombined to represent the ground motion at each site for any (heterogeneous) slip distribution on the fault. The reciprocal formulation significantly reduces the required number of 3-D finite difference computations to at most 3NS, where NS is the number of strong motion sites used in the inversion. Using controlled numerical resolution tests, we have examined the relative importance of accurate GFs for finite fault source inversions which rely on near-source ground motions. These experiments use both 1-D and 3-D GFs in inversions for hypothetical rupture models in order (1) to analyze the ability of the 3-D methodology to resolve trade-offs between complex source phenomena and 3-D path effects, (2) to address the sensitivity of the inversion results to uncertainties in the 3-D velocity structure, and (3) to test the adequacy of the 1-D GF method when propagation effects are known to be three-dimensional. We find that given "data" from a prescribed 3-D Earth structure, the use of well-calibrated 3-D GFs in the inversion provides very good resolution of the assumed slip distribution, thus adequately separating source and 3-D propagation effects. In contrast, using a set of inexact 3-D GFs or a set of hybrid 1-D GFs allows only partial recovery of the slip distribution. These findings suggest that in regions of complex geology the use of well-calibrated 3-D GFs has the potential for increased resolution of the rupture process relative to 1-D GFs. However, realizing this full potential requires that the 3-D velocity model and associated GFs should be carefully validated against the true 3-D Earth structure before performing the inverse problem with actual data.

Journal of Geophysical Research B: Solid Earth↗

USGS1998

Paleoseismic investigations in the Santa Cruz mountains, California: Implications for recurrence of large-magnitude earthquakes on the San Andreas fault

Trenching, microgeomorphic mapping, and tree ring analysis provide information on timing of paleoearthquakes and behavior of the San Andreas fault in the Santa Cruz mountains. At the Grizzly Flat site alluvial units dated at 1640–1659 A.D., 1679–1894 A.D., 1668–1893 A.D., and the present ground surface are displaced by a single event. This was the 1906 surface rupture. Combined trench dates and tree ring analysis suggest that the penultimate event occurred in the mid-1600 s, possibly in an interval as narrow as 1632–1659 A.D. There is no direct evidence in the trenches for the 1838 or 1865 earthquakes, which have been proposed as occurring on this part of the fault zone. In a minimum time of about 340 years only one large surface faulting event (1906) occurred at Grizzly Flat, in contrast to previous recurrence estimates of 95–110 years for the Santa Cruz mountains segment. Comparison with dates of the penultimate San Andreas earthquake at sites north of San Francisco suggests that the San Andreas fault between Point Arena and the Santa Cruz mountains may have failed either as a sequence of closely timed earthquakes on adjacent segments or as a single long rupture similar in length to the 1906 rupture around the mid-1600 s. The 1906 coseismic geodetic slip and the late Holocene geologic slip rate on the San Francisco peninsula and southward are about 50–70% and 70% of their values north of San Francisco, respectively. The slip gradient along the 1906 rupture section of the San Andreas reflects partitioning of plate boundary slip onto the San Gregorio, Sargent, and other faults south of the Golden Gate. If a mid-1600 s event ruptured the same section of the fault that failed in 1906, it supports the concept that long strike-slip faults can contain master rupture segments that repeat in both length and slip distribution. Recognition of a persistent slip rate gradient along the northern San Andreas fault and the concept of a master segment remove the requirement that lower slip sections of large events such as 1906 must fill in on a periodic basis with smaller and more frequent earthquakes.

Journal of Geophysical Research B: Solid Earth↗

USGS1999

Formation of natural gas hydrates in marine sediments 1. Conceptual model of gas hydrate growth conditioned by host sediment properties

The stability of submarine gas hydrates is largely dictated by pressure and temperature, gas composition, and pore water salinity. However, the physical properties and surface chemistry of deep marine sediments may also affect the thermodynamic state, growth kinetics, spatial distributions, and growth forms of clathrates. Our conceptual model presumes that gas hydrate behaves in a way analogous to ice in a freezing soil. Hydrate growth is inhibited within fine-grained sediments by a combination of reduced pore water activity in the vicinity of hydrophilic mineral surfaces, and the excess internal energy of small crystals confined in pores. The excess energy can be thought of as a "capillary pressure" in the hydrate crystal, related to the pore size distribution and the state of stress in the sediment framework. The base of gas hydrate stability in a sequence of fine sediments is predicted by our model to occur at a lower temperature (nearer to the seabed) than would be calculated from bulk thermodynamic equilibrium. Capillary effects or a build up of salt in the system can expand the phase boundary between hydrate and free gas into a divariant field extending over a finite depth range dictated by total methane content and pore-size distribution. Hysteresis between the temperatures of crystallization and dissociation of the clathrate is also predicted. Growth forms commonly observed in hydrate samples recovered from marine sediments (nodules, and lenses in muds; cements in sands) can largely be explained by capillary effects, but kinetics of nucleation and growth are also important. The formation of concentrated gas hydrates in a partially closed system with respect to material transport, or where gas can flush through the system, may lead to water depletion in the host sediment. This "freeze-drying" may be detectable through physical changes to the sediment (low water content and overconsolidation) and/or chemical anomalies in the pore waters and metastable presence of free gas within the normal zone of hydrate stability.

Journal of Geophysical Research B: Solid Earth↗

USGS1996

Type of faulting and orientation of stress and strain as a function of space and time in Kilauea's south flank, Hawaii

Earthquake focal mechanisms of events occurring between 1972 and 1992 in the south flank of Kilauea volcano, Hawaii, are used to infer the state of stress and strain as a function of time and space. We have determined 870 fault plane solutions from P wave first motion polarities for events with magnitudes M L ≥ 2.5 and depth ranging between 6 and 12 km. Faulting is characterized by a mixture of decollement, reverse, and normal faults. Most large earthquakes with magnitude M < 7 slip on reverse faults striking NE at 40° and dipping SE between 60° and 70°. In Hawaii, the earthquakes with M > 7 rupture the decollement plane, since it is the only surface large enough to generate magnitude 7 or larger earthquakes. The percentage of reverse faulting events is high compared to the decollement and normal faulting mechanisms for the period 1972–1983. The percentage of decollement type focal mechanisms becomes dominant after 1983. This pattern of faulting activity suggests that pressure was building up within Kilauea's rift zone prior to the 1983 Puu'Oo eruption. Overall, a single stress orientation with the maximum compressive stress oriented SE perpendicular to the rift and dipping at 45° is compatible with the coeval existence of decollement, reverse, and normal faults. However, in a crustal volume east of longitude 155°10′W, we find a change of the orientation of σ 1 from nearly horizontal to plunging 45° SE occurring in 1979. This stress rotation suggests magma movements within the aseismic part of Kilauea's east rift zone. The strain and stress orientations are coaxial in the south flank except within the volume where the stress rotation is observed. We observe a change in the relationship between stress and strain directions caused either by the shifting of seismic activity from reverse faults to decollements, while stress stays constant, or by a rotation of stress, while strain remains constant. Assuming that the model of a noncohesive Coulomb wedge is appropriate for Kilauea's south flank, we find that high pore pressures are prevalent along the decollement and within the wedge for a coefficient of friction equal to 0.85.

Hawaii↗