SEARCH · USGS Science

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 433 records · Page 24Linked to original sources

USGS2012

A seismic refraction study of the Oregon Cascades

A 275‐km‐long reversed refraction profile in the Oregon Cascades, two shallow earthquakes of magnitude 5 in southern Washington, a shallow earthquake of magnitude 4.6 in northern California, and a previously published analysis of the Bouguer gravity field are used to develop a crustal P wave velocity model for the Oregon Cascades. Travel time analysis of the refraction profile indicates a crustal structure characterized by surface layers with P wave velocities that vary from 2.9 to 5.2 km/s and thicknesses that vary from 2.5 to 5.0 km, upper crustal velocities of 6.1 to 6.5 km/s between the depths of 3 and 29 km, lower crustal velocities near 7.0 km/s between the depths of 29 and 44 km, and a mantle reflector at a depth of 44 km. Comparison with synthetic seismograms supports this structure and shows that the lower crustal and Mono transitions can be better modeled by continuous velocity gradients than by first‐order discontinuities. Arrival times from two shallow earthquakes in southern Washington in 1981 across 14 telemetry stations in the Oregon Cascades show an apparent velocity of 7.62 km/s. On the basis of a previously published gravity analysis, the mean north‐south component of the dip to the crust‐mantle boundary is estimated at 1° down dip to the south. The apparent velocity from the earthquakes and the dip estimate from gravity indicate that the true Pn velocity is 7.70 km/s. Arrival times from a shallow earthquake in northern California in 1978 across telemetry stations in the Washington Cascades are consistent with this model and indicate that the upper mantle velocity reaches 8.20 km/s at a depth of approximately 100 km.

Oregon↗

USGS2002

Fault structure and mechanics of the Hayward Fault, California from double-difference earthquake locations

The relationship between small-magnitude seismicity and large-scale crustal faulting along the Hayward Fault, California, is investigated using a double-difference (DD) earthquake location algorithm. We used the DD method to determine high-resolution hypocenter locations of the seismicity that occurred between 1967 and 1998. The DD technique incorporates catalog travel time data and relative P and S wave arrival time measurements from waveform cross correlation to solve for the hypocentral separation between events. The relocated seismicity reveals a narrow, near-vertical fault zone at most locations. This zone follows the Hayward Fault along its northern half and then diverges from it to the east near San Leandro, forming the Mission trend. The relocated seismicity is consistent with the idea that slip from the Calaveras Fault is transferred over the Mission trend onto the northern Hayward Fault. The Mission trend is not clearly associated with any mapped active fault as it continues to the south and joins the Calaveras Fault at Calaveras Reservoir. In some locations, discrete structures adjacent to the main trace are seen, features that were previously hidden in the uncertainty of the network locations. The fine structure of the seismicity suggests that the fault surface on the northern Hayward Fault is curved or that the events occur on several substructures. Near San Leandro, where the more westerly striking trend of the Mission seismicity intersects with the surface trace of the (aseismic) southern Hayward Fault, the seismicity remains diffuse after relocation, with strong variation in focal mechanisms between adjacent events indicating a highly fractured zone of deformation. The seismicity is highly organized in space, especially on the northern Hayward Fault, where it forms horizontal, slip-parallel streaks of hypocenters of only a few tens of meters width, bounded by areas almost absent of seismic activity. During the interval from 1984 to 1998, when digital waveforms are available, we find that fewer than 6.5% of the earthquakes can be classified as repeating earthquakes, events that rupture the same fault patch more than one time. These most commonly are located in the shallow creeping part of the fault, or within the streaks at greater depth. The slow repeat rate of 2–3 times within the 15-year observation period for events with magnitudes around M = 1.5 is indicative of a low slip rate or a high stress drop. The absence of microearthquakes over large, contiguous areas of the northern Hayward Fault plane in the depth interval from ∼5 to 10 km and the concentrations of seismicity at these depths suggest that the aseismic regions are either locked or retarded and are storing strain energy for release in future large-magnitude earthquakes.

California↗

USGS2024

Hyperspectral imaging predicts differences in carbon and nitrogen status among representative biocrust functional groups of the Colorado Plateau

Biological soil crusts (biocrusts) are widespread soil photosynthetic communities covering about 12% of Earth's land surface and play crucial roles in terrestrial carbon (C) and nitrogen (N) cycles, yet scalable quantifications of biocrusts and their biogeochemical contributions are notably lacking. While remote sensing has enormous potential to assess, scale, and contextualize biocrusts and their functions, the applicability of hyperspectral data in predicting C- and N-related biocrust traits remains largely unexplored. We address this issue by evaluating the potential of in situ hyperspectral data to predict C and N across a range of biocrust species and different environmental conditions. We found that in situ hyperspectral reflectance measurements can be used to predict biocrust tissue C/N ratios and N concentrations with relatively high accuracy but to a lesser extent for potential biocrust N2 fixation rates. Critical wavelength domains included the visible region of the spectrum from roughly 490–600 nm, which most effectively captured variations in biocrust tissue C, and the shortwave infrared region from 1,150 to 1,350 nm and 1,550–1,650 nm, which most effectively captured biocrust tissue N and N2 fixation potential. Finally, we provide evidence that multi- and hyperspectral missions with targeted band placement, such as the proposed 26-band Landsat Next, could be effective in predicting biocrust traits. This work provides a critical step in understanding how to apply data from new and upcoming satellite missions to the monitoring of biocrusts.

Utah↗

USGS2002

Crustal structure and relocated earthquakes in the Puget Lowland, Washington, from high-resolution seismic tomography

The availability of regional earthquake data from the Pacific Northwest Seismograph Network (PNSN), together with active source data from the Seismic Hazards Investigation in Puget Sound (SHIPS) seismic experiments, has allowed us to construct a new high-resolution 3-D, P wave velocity model of the crust to a depth of about 30 km in the central Puget Lowland. In our method, earthquake hypocenters and velocity model are jointly coupled in a fully nonlinear tomographic inversion. Active source data constrain the upper 10–15 km of the model, and earthquakes constrain the deepest portion of the model. A number of sedimentary basins are imaged, including the previously unrecognized Muckleshoot basin, and the previously incompletely defined Possession and Sequim basins. Various features of the shallow crust are imaged in detail and their structural transitions to the mid and lower crust are revealed. These include the Tacoma basin and fault zone, the Seattle basin and fault zone, the Seattle and Port Ludlow velocity highs, the Port Townsend basin, the Kingston Arch, and the Crescent basement, which is arched beneath the Lowland from its surface exposure in the eastern Olympics. Strong lateral velocity gradients, consistent with the existence of previously inferred faults, are observed, bounding the southern Port Townsend basin, the western edge of the Seattle basin beneath Dabob Bay, and portions of the Port Ludlow velocity high and the Tacoma basin. Significant velocity gradients are not observed across the southern Whidbey Island fault, the Lofall fault, or along most of the inferred location of the Hood Canal fault. Using improved earthquake locations resulting from our inversion, we determined focal mechanisms for a number of the best recorded earthquakes in the data set, revealing a complex pattern of deformation dominated by general arc-parallel regional tectonic compression. Most earthquakes occur in the basement rocks inferred to be the lower Tertiary Crescent formation. The sedimentary basins and the eastern part of the Olympic subduction complex are largely devoid of earthquakes. Clear association of hypocenters and focal mechanisms with previously mapped or proposed faults is difficult; however, seismicity, structure, and focal mechanisms associated with the Seattle fault zone suggest a possible high-angle mode of deformation with the north side up. We suggest that this deformation may be driven by isostatic readjustment of the Seattle basin.

Washington↗

USGS2020

Geological and thermal control of the hydrothermal system in northern Yellowstone Lake: Inferences from high resolution magnetic surveys

A multiscale magnetic survey of the northern basin of Yellowstone Lake was undertaken in 2016 as part of the Hydrothermal Dynamics of Yellowstone Lake Project (HD‐YLAKE)—a broad research effort to characterize the cause‐and‐effect relationships between geologic and environmental processes and hydrothermal activity on the lake floor. The magnetic survey includes lake surface, regional aeromagnetic, and near‐bottom autonomous underwater vehicle (AUV) data. The study reveals a strong contrast between the northeastern lake basin, characterized by a regional magnetic low punctuated by stronger local magnetic lows, many of which host hydrothermal vent activity, and the northwestern lake basin with higher‐amplitude magnetic anomalies and no obvious hydrothermal activity or punctuated magnetic lows. The boundary between these two regions is marked by a steep gradient in heat flow and magnetic values, likely reflecting a significant structure within the currently active ~20‐km‐long Eagle Bay‐Lake Hotel fault zone that may be related to the ~2.08‐Ma Huckleberry Ridge caldera rim. Modeling suggests that the broad northeastern magnetic low reflects both a shallower Curie isotherm and widespread hydrothermal activity that has demagnetized the rock. Along the western lake shoreline are sinuous‐shaped, high‐amplitude magnetic anomaly highs, interpreted as lava flow fronts of upper units of the West Thumb rhyolite. The AUV magnetic survey shows decreased magnetization at the periphery of the active Deep Hole hydrothermal vent. We postulate that lower magnetization in the outer zone results from enhanced hydrothermal alteration of rhyolite by hydrothermal condensates while the vapor‐dominated center of the vent is less altered.

Wyoming↗

USGS2024

Brittle regime slip partitioned damage and deformation mechanisms along the eastern Denali fault zone in southwestern, Yukon

Rare bedrock exposures of the eastern Denali fault zone in southwestern Yukon allow for the measurement, sampling, and analyses of brittle regime fault slip data and deformation mechanisms to explore relations to far field, oblique plate motions. Host rock lithologies and associated slip surfaces show episodic damage zone‐related deformation and calcite ± hematite ± chlorite related hydrothermal fluid flow. This regional scale network of asymmetric fault damage is spatially and kinematically linked to a discrete and narrow fault core. Fault network observations, orientations, slip data, and strain inversions document a slip partitioned strike‐slip fault system with locally and mutually overprinting strike‐, oblique‐, and dip‐slip components. Microstructural analyses reveal crystal plastic and co‐seismic brittle deformation mechanisms active in a narrow range of upper crustal temperature, pressure, fluid, and chemical conditions. The net damage related slip is not exclusively formed by a single kinematic system, but rather a fully partitioned, time integrated system likely operative for much of the fault's brittle regime evolution temporally constrained by previously published thermochronometric data. Although the fault slip data was collected from outcrop‐scale exposures at sites tens of kilometers apart, results show remarkable correlation between fault kinematics and plate motions along the ∼580 km long eastern Denali fault segment. End member, subhorizontal, northeast directed reverse and north directed dextral strike slip fault strain axes closely reflect relative plate motion interactions over at least the last 30 m.y. and act as a proxy for far‐field stresses compatible with the kinematics of the damage zone network.

Alaska↗

USGS2015

Direct measurement of asperity contact growth in quartz at hydrothermal conditions

Earthquake recurrence requires interseismic fault restrengthening which results from solid state deformation in room-temperature friction and indentation experiments. In contrast exhumed fault zones show solution-transport processes such as pressure solution and contact overgrowths influence fault zone properties . In the absence of fluid flow, overgrowths are driven by gradients in surface curvature where material is dissolved, diffuses, and precipitates at the contact without convergence normal to the contact. To determine the rate of overgrowth for quartz, we conducted single contact experiments in an externally heated pressure vessel. Convergence was continuously monitored using reflected-light interferometry through a long-working-distance microscope. Contact normal force was constant with an initial effective normal stress of 1.7 MPa, temperature was between 350 and 530{degree sign}C, and water pressure was constant at 150 MPa. Two control experiments were conducted: one dry at 425{degree sign}C and one bi-material (sapphire) at 425{degree sign}C and 150 MPa water pressure. No contact growth or convergence was observed in the controls. For wet single-phase contacts, growth was initially rapid and then decreased with time. No convergence was observed. Fluid inclusions indicate that the contact is not uniformly wetted. The contact is bounded by small regions of high aperture, reflecting local free-face dissolution as the source for the overgrowth. The apparent activation energy is ~125 kJ/mol. Extrapolation predicts rates of contact area increase orders of magnitude faster than in dry, room-temperature and hydrothermal friction experiments, suggesting that natural strength recovery near the base of the seismogenic zone could be dominated by contact overgrowth.

Journal of Geophysical Research B: Solid Earth↗

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↗