Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysical Journal”

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 901 records · Page 50Linked to original sources
USGS2002

Fault structure and kinematics of the Long Valley Caldera region, California, revealed by high-accuracy earthquake hypocenters and focal mechanism stress inversions

We have determined high-resolution hypocenters for 45,000+ earthquakes that occurred between 1980 and 2000 in the Long Valley caldera area using a double-difference earthquake location algorithm and routinely determined arrival times. The locations reveal numerous discrete fault planes in the southern caldera and adjacent Sierra Nevada block (SNB). Intracaldera faults include a series of east/west-striking right-lateral strike-slip faults beneath the caldera's south moat and a series of more northerly striking strike-slip/normal faults beneath the caldera's resurgent dome. Seismicity in the SNB south of the caldera is confined to a crustal block bounded on the west by an east-dipping oblique normal fault and on the east by the Hilton Creek fault. Two NE-striking left-lateral strike-slip faults are responsible for most seismicity within this block. To understand better the stresses driving seismicity, we performed stress inversions using focal mechanisms with 50 or more first motions. This analysis reveals that the least principal stress direction systematically rotates across the studied region, from NE to SW in the caldera's south moat to WNW-ESE in Round Valley, 25 km to the SE. Because WNW-ESE extension is characteristic of the western boundary of the Basin and Range province, caldera area stresses appear to be locally perturbed. This stress perturbation does not seem to result from magma chamber inflation but may be related to the significant (???20 km) left step in the locus of extension along the Sierra Nevada/Basin and Range province boundary. This implies that regional-scale tectonic processes are driving seismic deformation in the Long Valley caldera.

California↗
USGS2003

A revised dislocation model of interseismic deformation of the Cascadia subduction zone

CAS3D‐2, a new three‐dimensional (3‐D) dislocation model, is developed to model interseismic deformation rates at the Cascadia subduction zone. The model is considered a snapshot description of the deformation field that changes with time. The effect of northward secular motion of the central and southern Cascadia forearc sliver is subtracted to obtain the effective convergence between the subducting plate and the forearc. Horizontal deformation data, including strain rates and surface velocities from Global Positioning System (GPS) measurements, provide primary geodetic constraints, but uplift rate data from tide gauges and leveling also provide important validations for the model. A locked zone, based on the results of previous thermal models constrained by heat flow observations, is located entirely offshore beneath the continental slope. Similar to previous dislocation models, an effective zone of downdip transition from locking to full slip is used, but the slip deficit rate is assumed to decrease exponentially with downdip distance. The exponential function resolves the problem of overpredicting coastal GPS velocities and underpredicting inland velocities by previous models that used a linear downdip transition. A wide effective transition zone (ETZ) partially accounts for stress relaxation in the mantle wedge that cannot be simulated by the elastic model. The pattern of coseismic deformation is expected to be different from that of interseismic deformation at present, 300 years after the last great subduction earthquake. The downdip transition from full rupture to no slip should take place over a much narrower zone.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Undrained poroelastic response of sandstones to deviatoric stress change

Deformation of porous crustal rock through diagenesis, tectonic loading, or other processes can change pore volume and affect fluid pressure. The largest stress-induced pore pressure changes occur when fluid is trapped in pores in an “undrained” condition. We have measured the undrained poroelastic response of two sandstones to changes in mean and deviatoric stress. Pore pressure was found to respond to mean stress σ m in the usual manner: Δ p = B Δσ m ( B ranging from 0.4 to 0.7), nearly independent of the ambient deviatoric stress state. However, variations in deviatoric stress (σ d = (σ 1 − σ 3 )/2) at constant mean stress were also found to induce a reversible (elastic) pore pressure response to stress levels up to and exceeding 80% failure stress (i.e., Δ p = = const.). The coefficient η became more negative with increasing deviatoric stress level in sandstone and Ottawa sand samples. That is, η represents a dilatant response where increased deviatoric stress causes a decrease in pore pressure. The poroelastic response to deviatoric stress is explained in terms of anisotropic matrix stiffening due to closure of crack-like pore space or flattening of grain contacts at high ambient stress levels and can be important in calculations of earthquake stress transfer.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Concordant paleolatitudes for Neoproterozoic ophiolitic rocks of the Trinity Complex, Klamath Mountains, California

New paleomagnetic results from the eastern Klamath Mountains of northern California show that Neoproterozoic rocks of the Trinity ophiolitic complex and overlying Middle Devonian volcanic rocks are latitudinally concordant with cratonal North America. Combining paleomagnetic data with regional geologic and faunal evidence suggests that the Trinity Complex and related terranes of the eastern Klamath plate were linked in some fashion to the North American craton throughout that time, but that distance between them may have varied considerably. A possible model that is consistent with our paleomagnetic results and the geologic evidence is that the Trinity Complex formed and migrated parallel to paleolatitude in the basin between Laurasia and Australia–East Antarctica as the Rodinian supercontinent began to break up. It then continued to move parallel to paleolatitude at least through Middle Devonian time. Although the eastern Klamath plate served as a nucleus against which more western components of the Klamath Mountains province amalgamated, the Klamath superterrane was not accreted to North America until Early Cretaceous time.

California↗
USGS2002

Episodic thermal perturbations associated with groundwater flow: An example from Kilauea Volcano, Hawaii

Temperature measurements in deep drill holes on volcano summits or upper flanks allow a quantitative analysis of groundwater induced heat transport within the edifice. We present a new temperature-depth profile from a deep well on the summit of Kilauea Volcano, Hawaii, and analyze it in conjunction with a temperature profile measured 26 years earlier. We propose two groundwater flow models to interpret the complex temperature profiles. The first is a modified confined lateral flow model (CLFM) with a continuous flux of hydrothermal fluid. In the second, transient flow model (TFM), slow conductive cooling follows a brief, advective heating event. We carry out numerical simulations to examine the timescales associated with each of the models. Results for both models are sensitive to the initial conditions, and with realistic initial conditions it takes between 750 and 1000 simulation years for either model to match the measured temperature profiles. With somewhat hotter initial conditions, results are consistent with onset of a hydrothermal plume ∼550 years ago, coincident with initiation of caldera subsidence. We show that the TFM is consistent with other data from hydrothermal systems and laboratory experiments and perhaps is more appropriate for this highly dynamic environment. The TFM implies that volcano-hydrothermal systems may be dominated by episodic events and that thermal perturbations may persist for several thousand years after hydrothermal flow has ceased.

Hawaii↗
USGS2003

Importance of a stochastic distribution of floods and erosion thresholds in the bedrock river incision problem

Fluvial erosion of bedrock occurs during occasional flood events when boundary shear stress exceeds a critical threshold to initiate incision. Therefore efforts to model the evolution of topography over long timescales should include an erosion threshold and should be driven by a stochastic distribution of erosive events. However, most bedrock incision models ignore the threshold as a second‐order detail. In addition, climate is poorly represented in most landscape evolution models, so the quantitative relationship between erosion rate and measurable climatic variables has been elusive. Here we show that the presence of an erosion threshold, when combined with a well‐constrained, probabilistic model of storm and flood occurrence, has first‐order implications for the dynamics of river incision in tectonically active areas. First, we make a direct calculation of the critical shear stress required to pluck bedrock blocks for a field site in New York. Second, we apply a recently proposed stochastic, threshold, bedrock incision model to a series of streams in California, with known tectonic and climatic forcing. Previous work in the area has identified a weak relationship between channel gradient or relief and rock uplift rate that is not easily explained by simpler detachment‐limited models. The results with the stochastic threshold model show that even low erosion thresholds, which are exceeded in steep channels during high‐frequency flood events, fundamentally affect the predicted relationship between gradient and uplift rate in steady state rivers, in a manner consistent with the observed topography. This correspondence between theory and data is, however, nonunique; models in which a thin alluvial cover may act to inhibit channel incision in the low uplift rate zone also provide plausible explanations for the observed topography. Third, we explore the broader implications of the stochastic threshold model to the development of fluvial topography in active tectonic settings. We suggest that continued field applications of geomorphic models, including physically meaningful thresholds and stochastic climate distributions, are required to advance our knowledge of interactions among surficial, climatic, and crustal processes.

California↗
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / n , n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Fire dynamics and implications for nitrogen cycling in boreal forests

We used a dynamic, long-term mass balance approach to track cumulative carbon (C) and nitrogen (N) losses to fire in boreal Manitoba over the 6500 years since deglaciation. Estimated C losses to decomposition and fire, combined with measurements of N pools in mature and burned forest floors, suggest that loss of N by combustion has likely resulted in a long-term loss that exceeds the amount of N stored in soil today by 2 to 3 times. These estimates imply that biological N fixation rates could be as high as 5 to 10 times atmospheric deposition rates in boreal regions. At the site scale, the amount of N lost is due to N content of fuels, which varies by stand type and fire severity, which in turn vary with climate and fire dynamics. The interplay of fire frequency, fire severity, and N partitioning during regrowth are important for understanding rates and sustainability of nutrient and carbon cycling over millenia and over broad regions.

Manitoba↗
USGS2002

Stable isotope composition of waters in the Great Basin, United States 1. Air-mass trajectories

Isentropic trajectories, calculated using the NOAA/Climate Monitoring and Diagnostics Laboratory's isentropic transport model, were used to determine air-parcel origins and the influence of air mass trajectories on the isotopic composition of precipitation events that occurred between October 1991 and September 1993 at Cedar City, Utah, and Winnemucca, Nevada. Examination of trajectories that trace the position of air parcels backward in time for 10 days indicated five distinct regions of water vapor origin: (1) Gulf of Alaska and North Pacific, (2) central Pacific, (3) tropical Pacific, (4) Gulf of Mexico, and (5) continental land mass. Deuterium (??D) and oxygen-18 (??18O) analyses were made of precipitation representing 99% of all Cedar City events. Similar analyses were made on precipitation representing 66% of the precipitation falling at Winnemucca during the same period. The average isotopic composition of precipitation derived from each water vapor source was determined. More than half of the precipitation that fell at both sites during the study period originated in the tropical Pacific and traveled northeast to the Great Basin; only a small proportion traversed the Sierra Nevada. The isotopic composition of precipitation is determined by air-mass origin and its track to the collection station, mechanism of droplet formation, reequilibration within clouds, and evaporation during its passage from cloud to ground. The Rayleigh distillation model can explain the changes in isotopic composition of precipitation as an air mass is cooled pseudo-adiabatically during uplift. However, the complicated processes that take place in the rapidly convecting environment of cumulonimbus and other clouds that are common in the Great Basin, especially in summer, require modification of this model because raindrops that form in the lower portion of those clouds undergo isotopic change as they are elevated to upper levels of the clouds from where they eventually drop to the ground.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Stable isotope compositions of waters in the Great Basin, United States 2. Modern precipitation

Precipitation was collected between 1991 and 1997 at 41 locations within and adjacent to parts of the Great Basin lying in California, Oregon, Nevada, and Utah. These samples were analyzed for their deuterium (??D) and oxygen-18 (??18O) contents. Separate collections were made of summer and winter season precipitation at stations ranging in elevation from -65 m to 3246 m. The ??D per mil values of stations that were closely spaced but at different elevations showed an average ??D decrease of approximately 10???/km rise in elevation. Data for all samples representing winter precipitation, when plotted on a ??D versus ??18O plot, fall close to the Meteoric Water Line (??D = 8 ??18O + 10); samples representing summer precipitation define a line of slightly lower slope due to evaporation of the raindrops during their passage from cloud to ground. Comparison of our 1991-1997 ??D data with those from the same three stations reported by an earlier study in the southeastern California shows seasonal differences ranging from 0 per mil to 19??? (average: 15) and annual differences ranging from 0 to 13 per mil (average: 2), illustrating the degree of annual and seasonal variability in this region. When contoured, the ??D values display gradients indicating a north to northwest decrease in deuterium, with values ranging from -60 to -125??? in winter precipitation and from -40 to -110??? in summer precipitation. These gradient trends can be explained by the predominance of air mass trajectories originating in the tropical Pacific, the Gulf of California, and (in summer) the Gulf of Mexico.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Stable isotope compositions of waters in the Great Basin, United States 3. Comparison of groundwaters with modern precipitation

Groundwater samples from wells and springs, scattered over most of the Great Basin province, were collected and analyzed for their isotopic makeup. They were augmented by previously published isotopic data on groundwaters from southeast California and by several hundred unpublished isotopic analyses. The ratio of 2H (deuterium, D) to 1H, in water samples from valleys in parts of California, Idaho, Nevada, Oregon, and Utah, are here compared with the winter, summer, and annual isotopic compositions of precipitation falling in or near the sampled areas. The main goal of this study was to identify basins where the groundwaters have isotopic compositions that are "lighter" (depleted in the heavier isotope, D) relative to modern winter precipitation. Where these basins do not adjoin substantially higher terrain, we consider those light groundwaters to be of Pleistocene age and thus more than 10,000 years old. Where the groundwater is 10 to 19??? lighter than local winter precipitation, we consider it to be possibly an indication of Pleistocene water; where the ??D makeup is >20??? lighter, we consider it to be probably Pleistocene water. More than 80 sites underlain by waters of possible or probable Pleistocene age were identified.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Modeling soil thermal and carbon dynamics of a fire chronosequence in interior Alaska

In this study, the dynamics of soil thermal, hydrologic, and ecosystem processes were coupled to project how the carbon budgets of boreal forests will respond to changes in atmospheric CO 2 , climate, and fire disturbance. The ability of the model to simulate gross primary production and ecosystem respiration was verified for a mature black spruce ecosystem in Canada, the age-dependent pattern of the simulated vegetation carbon was verified with inventory data on aboveground growth of Alaskan black spruce forests, and the model was applied to a postfire chronosequence in interior Alaska. The comparison between the simulated soil temperature and field-based estimates during the growing season (May to September) of 1997 revealed that the model was able to accurately simulate monthly temperatures at 10 cm ( R > 0.93) for control and burned stands of the fire chronosequence. Similarly, the simulated and field-based estimates of soil respiration for control and burned stands were correlated ( R = 0.84 and 0.74 for control and burned stands, respectively). The simulated and observed decadal to century-scale dynamics of soil temperature and carbon dynamics, which are represented by mean monthly values of these variables during the growing season, were correlated among stands ( R = 0.93 and 0.71 for soil temperature at 20- and 10-cm depths, R = 0.95 and 0.91 for soil respiration and soil carbon, respectively). Sensitivity analyses indicate that along with differences in fire and climate history a number of other factors influence the response of carbon dynamics to fire disturbance. These factors include nitrogen fixation, the growth of moss, changes in the depth of the organic layer, soil drainage, and fire severity.

Alaska↗
USGS2001

Dating the Vostok ice core record by importing the Devils Hole chronology

The development of an accurate chronology for the Vostok record continues to be an open research question because these invaluable ice cores cannot be dated directly. Depth-to-age relationships have been developed using many different approaches, but published age estimates are inconsistent, even for major paleoclimatic events. We have developed a chronology for the Vostok deuterium paleotemperature record using a simple and objective algorithm to transfer ages of major paleoclimatic events from the radiometrically dated 500,000-year δ 18 O-paleotemperature record from Devils Hole, Nevada. The method is based only on a strong inference that major shifts in paleotemperature recorded at both locations occurred synchronously, consistent with an atmospheric teleconnection. The derived depth-to-age relationship conforms with the physics of ice compaction, and internally produces ages for climatic events 5.4 and 11.24 which are consistent with the externally assigned ages that the Vostok team needed to assume in order to derive their most recent chronology, GT4. Indeed, the resulting V-DH chronology is highly correlated with GT4 because of the unexpected correspondence even in the timing of second-order climatic events that were not constrained by the algorithm. Furthermore, the algorithm developed herein is not specific to this problem; rather, the procedure can be used whenever two paleoclimate records are proxies for the same physical phenomenon, and paleoclimatic conditions forcing the two records can be considered to have occurred contemporaneously. The ability of the algorithm to date the East Antarctic Dome Fuji core is also demonstrated.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Thermal infrared spectroscopy of experimentally shocked anorthosite and pyroxenite: Implications for remote sensing of Mars

The feldspar and pyroxene mineralogies on Mars revealed by the Thermal Emission Spectrometer (TES) on Mars Global Surveyor likely record a variety of shock effects, as suggested by petrologic analyses of the Martian meteorites and the abundance of impact craters on the planet's surface. To study the effects of shock pressures on thermal infrared spectra of these minerals, we performed shock recovery experiments on orthopyroxenite and anorthosite samples from the Stillwater Complex (Montana) over peak pressures from 17 to 63 GPa. We acquired emissivity and hemispherical reflectance spectra (350–1400 cm −1 ; ∼7–29 μm) of both coherent chips and fine-grained powders of shocked and unshocked samples. These spectra are more directly comparable to remotely sensed data of Mars (e.g., TES) than previously acquired absorption or transmission spectra of shocked minerals. The spectra of experimentally shocked feldspar show systematic changes with increasing pressure due to depolymerization of the silica tetrahedra. For the spectra of chips, this includes the disappearance of small bands in the 500–650 cm −1 region and a strong band at 1115 cm −1 , and changes in positions of a strong band near 940 cm −1 and the Christiansen feature near 1250 cm −1 . Spectra of the shocked powders show the gradual disappearance of a transparency feature near 830 cm −1 . Fewer changes are observed in the pyroxene spectra at pressures as high as 63 GPa. Spectra of experimentally shocked minerals will help identify more precisely the mineralogy of rocks and soils not only from TES but also from Mars instruments such as miniTES and THEMIS.

Journal of Geophysical Research E: Planets↗
USGS2003

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Stress-induced, time-dependent fracture closure at hydrothermal conditions

Time-dependent closure of fractures in quartz was measured in situ at 22–530°C temperature and 0.1-150 MPa water pressure. Unlike previous crack healing and rock permeability studies, in this study, fracture aperture is monitored directly and continuously using a windowed pressure vessel, a long-working-distance microscope, and reflected-light interferometry. Thus the fracture volume and geometry can be measured as a function of time, temperature, and water pressure. Relatively uniform closure occurs rapidly at temperatures and pressures where quartz becomes significantly soluble in water. During closure the aperture is reduced by as much as 80% in a few hours. We infer that this closure results from the dissolution of small particles or asperities that prop the fracture open. The driving force for closure via dissolution of the prop is the sum of three chemical potential terms: (1) the dissolution potential, proportional to the logarithm of the degree of undersaturation of the solution; (2) the coarsening potential, proportional to the radius of curvature of the prop; and (3) the pressure solution potential, proportional to the effective normal stress at the contact between propping particles and the fracture wall. Our observations suggest that closure is controlled by a pressure solution-like process. The aperture of dilatant fractures and microcracks in the Earth that are similar to those in our experiments, such as ones generated from thermal stressing or brittle failure during earthquake rupture and slip, will decrease rapidly with time, especially if the macroscopic stress is nonhydrostatic.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

A fault constitutive relation accounting for thermal pressurization of pore fluid

The heat generated in a slip zone during an earthquake can raise fluid pressure and thereby reduce frictional resistance to slip. The amount of fluid pressure rise depends on the associated fluid flow. The heat generated at a given time produces fluid pressure that decreases inversely with the square root of hydraulic diffusivity times the elapsed time. If the slip velocity function is crack-like, there is a prompt fluid pressure rise at the onset of slip, followed by a slower increase. The stress drop associated with the prompt fluid pressure rise increases with rupture propagation distance. The threshold propagation distance at which thermally induced stress drop starts to dominate over frictionally induced stress drop is proportional to hydraulic diffusivity. If hydraulic diffusivity is 0.02 m 2 /s, estimated from borehole samples of fault zone material, the threshold propagation distance is 300 m. The stress wave in an earthquake will induce an unknown amount of dilatancy and will increase hydraulic diffusivity, both of which will lessen the fluid pressure effect. Nevertheless, if hydraulic diffusivity is no more than two orders of magnitude larger than the laboratory value, then stress drop is complete in large earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Magma supply dynamics at Westdahl volcano, Alaska, modeled from satellite radar interferometry

A group of satellite radar interferograms that span the time period from 1991 to 2000 shows that Westdahl volcano, Alaska, deflated during its 1991-1992 eruption and is reinflating at a rate that could produce another eruption within the next several years. The rates of inflation and deflation are approximated by exponential decay functions having time constants of about 6 years and a few days, respectively. This behavior is consistent with a deep, constant-pressure magma source connected to a shallow reservoir by a magma-filled conduit. An elastic deformation model indicates that the reservoir is located about 6 km below sea level and beneath Westdahl Peak. We propose that the magma flow rate through the conduit is governed by the pressure gradient between the deep source and the reservoir. The pressure gradient, and hence the flow rate, are greatest immediately after eruptions. Pressurization of the reservoir decreases both the pressure gradient and the flow rate, but eventually the reservoir ruptures and an eruption or intrusion ensues. The eruption rate is controlled partly by the pressure gradient between the reservoir and surface, and therefore it, too, decreases with time. When the supply of eruptible magma is exhausted, the eruption stops, the reservoir begins to repressurize at a high rate, and the cycle repeats. This model might also be appropriate for other frequently active volcanoes with stable magma sources and relatively simple magma storage systems.

Journal of Geophysical Research B: Solid Earth↗