USGS Science⌕ Search

SEARCH · USGS Science

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 1,657 records · Page 92Linked to original sources
USGS2012

Analysis of ERS 1 synthetic aperture radar data of frozen lakes in northern Montana and implications for climate studies

Lakes that freeze each winter are good indicators of regional climate change if key parameters, such as freeze-up and breakup date and maximum ice thickness, are measured over a decade-scale time frame. Synthetic aperture radar (SAR) satellite data have proven to be especially useful for measurement of climatologically significant parameters characteristic of frozen lakes. In this paper, five lakes in Glacier National Park, Montana, have been studied both in the field and using Earth Remote-Sensing Satellite (ERS) 1 SAR data during the 1992-1993 winter. The lakes are characterized by clear ice, sometimes with tubular or rounded bubbles, and often with a layer of snow ice on top of the clear ice. They are also often snow covered. Freeze-up is detected quite easily using ERS 1 SAR data as soon as a thin layer of ice forms. The effect of snow ice on the backscatter is thought to be significant but is, as yet, undetermined. On the five lakes studied, relative backscatter was found to increase with ice thickness until a maximum was reached in February. Breakup, an often ill-defined occurrence, is difficult to detect because surface water causes the SAR signal to be absorbed, thus masking the ice below. Comparison of the bubble structure of thaw lakes in northern Alaska with lakes in northern Montana has shown that the ice structure is quite different, and this difference may contribute to differential SAR signature evolution in the lakes of the two areas.

Journal of Geophysical Research↗
USGS2012

The M2 tide on the Amazon Shelf

As part of A Multidisciplinary Amazon Shelf Sediment Study (AMASSEDS), moored and shipboard current measurements made over the Amazon shelf during 1990–1991 have been analyzed to determine the dominant semidiurnal tidal constituent, the M 2 . These results have been combined with coastal sea level data from within the Amazon and Para Rivers, the adjacent shelf, and with satellite‐derived tidal elevation data from off the shelf to provide a more complete description of the M 2 tide in this complex river/shelf system. Near the Amazon River mouth the M 2 tide propagates across the shelf and through the mouth as a damped progressive wave, with its amplitude decreasing and phase increasing upriver. Over the adjacent shelf north of Cabo Norte, the M 2 tide approaches a damped standing wave, with large amplitudes (greater than 1.5 m) near the coast due to near resonance within the coastal embayment formed by the Cabo Norte shoal to the south and Cabo Cassipore to the north. The observed M 2 tidal currents are nearly rectilinear and oriented primarily across the local isobaths. Comparisons between tidal observations in both the North Channel and the Cabo Norte‐Cabo Cassipore embayment and a simple variable‐width channel tidal model indicate that (1) most of the M 2 tidal energy dissipation occurs over the mid‐ and inner shelf (in water depths less than 20 m) and (2) fluid muds found there cause a significant reduction (of order 50%) in the effective bottom friction felt by the M 2 tide. The approximate resonant period of the Cabo Norte‐Cabo Cassipore embayment is 11.9 hours, and at resonance the average energy dissipation per forcing period is roughly 2.2 times the average mechanical energy in the embayment. This damping rate is large enough that the tidal amplification is rather insensitive to forcing frequency, so that the response of the embayment to forcing over the semidiurnal band should be essentially the same. The vertical structure of the M 2 tidal current is examined at one outer shelf site located in 65‐m water depth. The observed semimajor axis increases logarithmically with height above bottom within the lowest 1–2 m and reaches a maximum in excess of 0.5 m/s at approximately 11 m above bottom. The mean ellipticity is small (less than 0.1) and positive, indicating clockwise rotation of a nearly rectilinear current, and the semimajor axis is oriented within 10° of the local cross‐isobath direction. The M 2 phase increases with height above bottom, with flood at the bottom leading flood at the surface by about 1 hour. A simple, local homogeneous tidal model with time‐ and space‐dependent eddy viscosity simulates the observed near‐bottom velocity reasonably well, however, the model suggests that stratification above the lowest few meters may significantly affect the tidal boundary layer structure at this site. The M 2 energy flux onto the Amazon shelf and into the Amazon and Para Rivers has been estimated using current and surface elevation data and the best fit variable‐width channel model results. The net M 2 energy flux into the mouths of the Amazon and Para Rivers is 0.47×10 10 W and 0.19×10 10 W, respectively. A net M 2 energy flux of about 3.3×10 10 W occurs onto the shelf between the North Channel of the Amazon River and Cabo Cassipore. This stretch of the Amazon shelf accounts for about 1.3% of the global dissipation of the M 2 tide.

Journal of Geophysical Research C: Oceans↗
USGS2012

Dust deposition in southern Nevada and California, 1984–1989: Relations to climate, source area, and source lithology

Dust samples collected annually for 5 years from 55 sites in southern Nevada and California provide the first regional source of information on modern rates of dust deposition, grain size, and mineralogical and chemical composition relative to climate and to type and lithology of dust source. The average silt and clay flux (rate of deposition) in southern Nevada and southeastern California ranges from 4.3 to 15.7 g/m 2 /yr, but in southwestern California the average silt and clay flux is as high as 30 g/m 2 /yr. The climatic factors that affect dust flux interact with each other and with the factors of source type (playas versus alluvium), source lithology, geographic area, and human disturbance. Average dust flux increases with mean annual temperature but is not correlated to decreases in mean annual precipitation because the regional winds bring dust to relatively wet areas. In contrast, annual dust flux mostly reflects changes in annual precipitation (relative drought) rather than temperature. Although playa and alluvial sources produce about the same amount of dust per unit area, the total volume of dust from the more extensive alluvial sources is much larger. In addition, playa and alluvial sources respond differently to annual changes in precipitation. Most playas produce dust that is richer in soluble salts and carbonate than that from alluvial sources (except carbonate‐rich alluvium). Gypsum dust may be produced by the interaction of carbonate dust and anthropogenic or marine sulfates. The dust flux in an arid urbanizing area may be as much as twice that before disturbance but decreases when construction stops. The mineralogic and major‐oxide composition of the dust samples indicates that sand and some silt is locally derived and deposited, whereas clay and some silt from different sources can be far‐traveled. Dust deposited in the Transverse Ranges of California by the Santa Ana winds appears to be mainly derived from sources to the north and east.

California, Nevada↗
USGS2012

Gravity survey of the Mount Toondina impact structure, South Australia

Gravity and seismic reflection data, together with geologic mapping, indicate that the Mount Toondina feature in South Australia is best interpreted as an eroded 4-km-diameter impact structure consisting of a ring structural depression surrounding a pronounced central uplift. Beds at the center of the structure within the central uplift have been raised as much as 200 m from depth and deformed by convergent flow. Seismic reflection data indicate that deformation extends to depths of only ≈800 m; at greater depths the reflectors are nearly flat lying, indicating little or no deformation. Gravity data show residual anomalies of +1.0 mGal coincident with the central uplift and a −0.5 Mgal low associated with the ring structural depression. Modeling of the gravity data indicates that relatively high-density material occurs within the central uplift, whereas the ring depression is filled with low-density material. The deformation at Mount Toondina is typical of a complex impact crater; the 4-km diameter is consistent with the expected threshold size for complex craters formed in weak to moderate strength sedimentary rocks.

South Australia↗
USGS2012

Hydrous carbonates on Mars?: Evidence from Mariner 6/7 infrared spectrometer and ground‐based telescopic spectra

Absorption features at 2.28 and 5.4 μm identified in Mariner 6/7 infrared spectrometer and terrestrial telescopic spectra are consistent with the spectra of hydrous magnesium carbonates such as hydromagnesite and artinite. Spectral characteristics of these hydrous carbonates are different from those of the anhydrous carbonates, as the former do not have the strong spectral features typically associated with anhydrous carbonates such as calcite and siderite. Theoretical mixing indicates that, depending on the type of hydrous carbonate, 10–20 wt % can be incorporated into the regolith without contradicting the spectral observations or the Viking x ray fluorescence chemical analysis. Hydrous carbonates form as weathering products of mafic minerals in the presence of H 2 O and CO 2 , even in the Antarctic. Their formation as evaporite minerals from either original magmas or hydrothermally altered rocks is consistent with the Martian environment, provided liquid water is or has been at least transiently present. On Earth, formation of hydrous Mg carbonates is associated with the production of amorphous iron oxides, which is consistent with both the environment and the inferred surface mineralogy of Mars. These minerals are about 60 wt % H 2 O, CO 3 , and OH; if they are abundant everywhere at the 10% level, then about 6% of the surface weight could be volatiles bound in this type of mineral. Although the stability of hydrous carbonates in a Martian environment is uncertain, there may be kinetic factors inhibiting the dehydration of these minerals, which may persist metastably in the current environment. Although the spectroscopic evidence for anhydrous carbonates is scant, the possible presence of hydrous carbonates provides an appealing mechanism for the existence of carbonates on Mars.

Journal of Geophysical Research↗
USGS2012

Superposed local and regional paleostresses: fault-slip analysis of Neogene extensional faulting near coeval caldera complexes, Yucca Flat, Nevada

Numerous reduced stress tensors are computed by multiple inversions of 906 temporally and spatially partitioned fault-slip data from the Yucca Flat region in the southwest Nevada volcanic field to constrain the Neogene paleostress and faulting history and to investigate how the regional tectonic stress field was affected by local caldera magmatism. Perturbed, shallow (<400 m), pre-11 Ma paleostress configurations, determined west and northwest of present (post-11 Ma) Yucca Flat basin, existed during mild extensional faulting and are attributed to superposition of transient caldera-magmatic stresses on the regional stress field. Northwest of Yucca Flat a progressive shift in least principal stress (σ 3 ) directions near known calderas located 5–15 km to the west occurred under a normal-slip stress state during caldera development between about 15 and 13 Ma. A brief (∼0.5 m.y.) change to a strike-slip stress state occurred at about 13 Ma and was accompanied by small-offset, quasi-conjugate strike-slip faulting. This stress state was most distinct, relative to a normal-slip state, near calderas where stress solutions and fault relations indicate closer affinities to a reverse-slip state. Inferred 11.6–11.45 Ma paleostress tensors indicate radial tension associated with either initial caldera collapse or local post-collapse topographic modification of the stress field. Post-11 Ma normal-slip stress tensors are associated with normal- and oblique-slip faults that accommodated subsidence and eastward extension of Yucca Flat basin away from the caldera complexes. These tensors do not indicate stress modifications due to residual caldera-related effects and thus were used to infer post-11 Ma regional stress changes. The stress field has rotated as much as 65° clockwise since 11 Ma during extensional development of Yucca Flat basin, with most of the rotation and extension occurring before about 8.5 Ma. Results suggest that shallow magmatism and caldera development can strongly alter extensional tectonic stress fields, fault patterns, and slip directions in the uppermost crust out to distances of roughly two magma chamber radii away from a magma body.

Journal of Geophysical Research↗
USGS2012

Interseismic uplift at the Nankai subduction zone, southwest Japan, 1951–1990

Uplift as a function of time from 1951 through 1990 has been deduced from annual mean sea level measured at 15 tide gages along the Nankai subduction zone, southwest Japan. The recurrence interval for rupture of the Nankai subduction zone is about 100 years, and the most recent rupture was in late 1946. Thus the 1951–1990 uplift record covers most of the first half of the earthquake cycle. The precision of the uplift record is better than could be obtained currently by annual Global Positioning System (GPS) measurements. The pre‐1960 uplift record shows rapid deformation that appears to terminate in 1959. After 1959 the uplift record is remarkably linear in time: Significant curvature in the uplift‐versus‐time plots is detected at only three of the 15 tide gages. The inferred uplift rates are not quantitatively consistent with the predictions of either the viscoelastic coupling or elastic half‐space models of subduction, but the agreement is sufficient to suggest that the causative physical mechanisms have been identified. The immediate postseismic response is controlled by the propagation of slip downward along the downdip extension of the coseismic rupture, and the interseismic response is caused by the accumulation of a slip deficit on the main thrust zone. Asthenosphere relaxation is not required to explain the observations.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Seismic velocity structure and composition of the continental crust: A global view

Seismic techniques provide the highest‐resolution measurements of the structure of the crust and have been conducted on a worldwide basis. We summarize the structure of the continental crust based on the results of seismic refraction profiles and infer crustal composition as a function of depth by comparing these results with high‐pressure laboratory measurements of seismic velocity for a wide range of rocks that are commonly found in the crust. The thickness and velocity structure of the crust are well correlated with tectonic province, with extended crust showing an average thickness of 30.5 km and orogens an average of 46.3 km. Shields and platforms have an average crustal thickness nearly equal to the global average. We have corrected for the nonuniform geographical distribution of seismic refraction profiles by estimating the global area of each major crustal type. The weighted average crustal thickness based on these values is 41.1 km. This value is 10% to 20% greater than previous estimates which underrepresented shields, platforms, and orogens. The average compressional wave velocity of the crust is 6.45 km/s, and the average velocity of the uppermost mantle ( Pn velocity) is 8.09 km/s. We summarize the velocity structure of the crust at 5‐km depth intervals, both in the form of histograms and as an average velocity‐depth curve, and compare these determinations with new measurements of compressional wave velocities and densities of over 3000 igneous and metamorphic rock cores made to confining pressures of 1 GPa. On the basis of petrographic studies and chemical analyses, the rocks have been classified into 29 groups. Average velocities, densities, and standard deviations are presented for each group at 5‐km depth intervals to crustal depths of 50 km along three different geotherms. This allows us to develop a model for the composition of the continental crust. Velocities in the upper continental crust are matched by velocities of a large number of lithologies, including many low‐grade metamorphic rocks and relatively silicic gneisses of amphibolite facies grade. In midcrustal regions, velocity gradients appear to originate from an increase in metamorphic grade, as well as a decrease in silica content. Tonalitic gneiss, granitic gneiss, and amphibolite are abundant midcrustal lithologies. Anisotropy due to preferred mineral orientation is likely to be significant in upper and midcrustal regions. The bulk of the lower continental crust is chemically equivalent to gabbro, with velocities in agreement with laboratory measurements of mafic granulite. Garnet becomes increasingly abundant with depth, and mafic garnet granulite is the dominant rock type immediately above the Mohorovicic discontinuity. Average compressional wave velocities of common crustal rock types show excellent correlations with density. The mean crustal density calculated from our model is 2830 kg/m 3 , and the average SiO 2 content is 61.8%.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Geodetic monitoring of the southern San Andreas Fault, California, 1980-1991

Five geodetic arrays (10 to 40 km aperture) located along the San Andreas fault have been surveyed frequently (several times in most years) over the 1980–1991 interval to detect possible fluctuations in the deformation rate. In each survey of an array the distances between the same four to seven pairs of geodetic monuments were measured. The distances measured (with corresponding standard deviation) range from 8.4 km (3.4 mm) to 38 km (8.2 mm). The data are displayed as plots of measured distance as a function of time. Linear fits in such plots furnish a satisfactory representation of the data. That is, the scatter of the data about the linear fits is within the range expected for the estimated standard deviations in measurement. However, there are coherent low‐amplitude (within the observational error) fluctuations apparent in some of the measured distances. Those fluctuations need not be tectonic effects but rather may be either random patterns or artifacts introduced by systematic drift in instrument calibration or wander of the geodetic monuments. We conclude that the measurements are consistent with steady deformation of the arrays over the 1980–1991 interval, and we find no convincing evidence in the data for fluctuations in the rate of deformation.

California↗
USGS2012

Estimation of rod scale errors in geodetic leveling

Comparisons among repeated geodetic levelings have often been used for detecting and estimating residual rod scale errors in leveled heights. Individual rod-pair scale errors are estimated by a two-step procedure using a model based on either differences in heights, differences in section height differences, or differences in section tilts. It is shown that the estimated rod-pair scale errors derived from each model are identical only when the data are correctly weighted, and the mathematical correlations are accounted for in the model based on heights. Analyses based on simple regressions of changes in height versus height can easily lead to incorrect conclusions. We also show that the statistically estimated scale errors are not a simple function of height, height difference, or tilt. The models are valid only when terrain slope is constant over adjacent pairs of setups (i.e., smoothly varying terrain). In order to discriminate between rod scale errors and vertical displacements due to crustal motion, the individual rod-pairs should be used in more than one leveling, preferably in areas of contrasting tectonic activity. From an analysis of 37 separately calibrated rod-pairs used in 55 levelings in southern California, we found eight statistically significant coefficients that could be reasonably attributed to rod scale errors, only one of which was larger than the expected random error in the applied calibration-based scale correction. However, significant differences with other independent checks indicate that caution should be exercised before accepting these results as evidence of scale error. Further refinements of the technique are clearly needed if the results are to be routinely applied in practice.

Journal of Geophysical Research B: Solid Earth↗
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

Effects of geothermal development on deformation in the Long Valley Caldera, eastern California, 1985-1994

Long Valley caldera in east central California has been the site of crustal unrest in the form of seismicity and ground deformation in response to magmatic inflation since 1980. Uplift of the resurgent dome has totaled ∼0.6 m over the 1975–1992 period. Within this region of uplift, and near the southwestern edge of the resurgent dome, is the Casa Diablo area, which experienced relative subsidence of 0.17 m between 1985 and 1992. Geothermal fluid production began at Casa Diablo in 1985 to supply a binary electric power plant; currently, three such plants generate a total of about 40 MW of electricity. The plants are supplied by 170°C water pumped from depths near 150m; all the produced water is reinjected at cooler temperatures at depths near 600 m. Analyses of data from (1) regional leveling lines and a network of bench marks in the geothermal well field area, (2) tilt observations from an L‐shaped array of bench marks, and (3) calculations relating pressure and temperature changes to subsidence indicate that relative subsidence at Casa Diablo results mainly from reductions in pressure in the shallow production reservoir and overlying formations and reductions in temperature in the underlying injection reservoir. Secondary effects on deformation are caused by offsets along the Taylor‐Bryant fault bounding the west side of the well field.

California↗
USGS1996

Identifying coseismic subsidence in tidal-wetland stratigraphic sequences at the Cascadia subduction zone of western North America

Tidal-wetland stratigraphy reveals that great plate boundary earthquakes have caused hundreds of kilometers of coast to subside at the Cascadia subduction zone. However, determining earthquake recurrence intervals and mapping the coastal extent of past great earthquake ruptures in this region are complicated by the effects of many sedimentologic, hydrographic, and oceanographic processes that occur on the coasts of tectonically passive as well as active continental margins. Tidal-wetland stratigraphy at many Cascadia estuaries differs little from that at similar sites on passive-margin coasts where stratigraphic sequences form through nonseismic processes unrelated to coseismic land level changes. Methods developed through study of similar stratigraphic sequences in Europe provide a framework for investigating the Cascadia estuarine record. Five kinds of criteria must be evaluated when inferring regional coastal subsidence due to great plate boundary earthquakes: the suddenness and amount of submergence, the lateral extent of submerged tidal-wetland soils, the coincidence of submergence with tsunami deposits, and the degree of synchroneity of submergence events at widely spaced sites. Evaluation of such criteria at the Cascadia subduction zone indicates regional coastal subsidence during at least two great earthquakes. Evidence for a coseismic origin remains equivocal, however, for the many peat-mud contacts in Cascadia stratigraphic sequences that lack (1) contrasts in lithology or fossils indicative of more than half a meter of submergence, (2) well-studied tsunami deposits, or (3) precise ages needed for regional correlation. Paleoecologic studies of fossil assemblages are particularly important in estimating the size of sudden sea level changes recorded by abrupt peat-mud contacts and in helping to distinguish erosional and gradually formed contacts from coseismic contacts. Reconstruction of a history of great earthquakes for the Cascadia subduction zone will require rigorous application of the above criteria and many detailed investigations.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Interseismic deformation along the San Andreas Fault in southern California

Eight trilateration networks located along the San Andreas fault in southern California have been surveyed 8 to 19 times within a 14- to 17-year interval between 1971 and 1992. The data, measurements of distances between the same 10 to 32 pairs of geodetic monuments within a network in each of the surveys, have been corrected for coseismic offsets from nearby earthquakes calculated from dislocation models. The corrected data, a total of 2027 measurements, are displayed in the form of plots of measured distance versus time for each of the 167 lines measured. The hypothesis that the interseismic deformation is steady is tested by examining whether deviations from linear fits to the data in the plots are within the range expected for observational error. A significant deviation from steady deformation is found only for the network located neax Palmdale, California. In that network, many of the measurements made in the survey of early 1982 deviate from the trend defined by the measurements in other surveys. The deviations are not of the form (fixed proportional error) one would generally expect from systematic survey error. If the survey of early 1982 is excluded, the remaining data at Palmdale are consistent with steady deformation. Thus the apparent strain event observed near Palmdale in 1982 was transient.

Journal of Geophysical Research↗
USGS2012

Principal component analysis of interseismic deformation in southern California

Eight trilateration networks along the San Andreas fault in southern California have been surveyed repeatedly (8 to 19 times) within an approximately 15-year interval between 1971 and 1992. The data for each network were analyzed by principal component analysis to represent the observed changes in line length as a superposition of individual modes: L ij - L i = ∑ k A ik C k ( t j ), where L ij is the length of the i th line measured at the time t j , L i is the average length of the i th line, and A ik C k ( t j ) is the k th mode. Each mode is described by a common time function Ck(tj) that is evaluated only at the times t j of the surveys, and the participation of each line in that mode is described by a constant amplitude A ik for that line. For all networks the deformation is largely described by the first mode, and the first-mode time function is essentially linear in time. That is, the deformation in each of the networks appears to be steady. Two anomalies in the accumulation of deformation were observed: In the network near Palmdale a transient strain event occurred in the second mode in early 1982, and in the network within Cajon Pass, fluctuations in the accumulation of deformation somewhat greater than would be attributed to measurement error were found. Although both anomalies are formally significant given the estimated measurement error, I am not confident that either anomaly is real. The principal conclusion is that the accumulation of interseismic deformation in the southern California trilateration networks is steady within the precision of measurement.

Journal of Geophysical Research↗
USGS2012

Strain accumulation across the central Nevada seismic zone, 1973–1994

Five trilateration networks extending for 280 km along the central Nevada seismic zone (1915 Pleasant Valley, M = 7.3; 1954 Dixie Valley, M = 6.8; 1954 Stillwater, M = 6.8; 1954 Rainbow Mountain, M = 6.6; 1954 Fairview Peak, M = 7.1; and 1932 Cedar Mountain, M = 7.2) have been surveyed 6 times since 1973 to determine deformation along the zone. Within the precision of measurement the deformation appears uniform along the zone and is described by the principal strain rates 0.036±0.008 μstrain/yr N60°W±3° and −0.031±0.008 μstrain/yr N30°E±3°, extension reckoned positive. The observed strain rates are consistent with simple, right‐lateral, tensor shear at the rate of 0.033 μstrain/yr across a shear zone striking N15°W. This central Nevada shear zone appears to be the northward continuation of the eastern California shear zone. The orientation of the strike‐slip and normal‐slip ruptures within the central Nevada seismic zone are consistent with principal stress axes parallel to the measured principal strain rate axes. Space‐based geodetic measurements (very long baseline interferometry) indicate that the relative motion accommodated across the Basin and Range province west of Ely, Nevada, is about 9.1±1.5 mm/yr N16°W±8° (Dixon et al., 1995.) Notice that the right‐lateral shear zone postulated to explain deformation in the central Nevada seismic zone is properly oriented to accommodate that relative motion. However, a 135‐km effective width of the shear zone would be required to accommodate all of the 9.1 mm/yr relative motion at the strain rates observed in the Nevada seismic zone; only about 3 mm/yr of that relative motion is accommodated within the span of the trilateration networks.

California, Nevada↗
USGS1996

Crustal velocity field near the big bend of California's San Andreas fault

We use geodetic data spanning the 1920–1992 interval to estimate the horizontal velocity field near the big bend segment of California's San Andreas fault (SAF). More specifically, we estimate a horizontal velocity vector for each node of a two-dimensional grid that has a 15-min-by-15-min mesh and that extends between latitudes 34.0°N and 36.0°N and longitudes 117.5°W and 120.5°W. For this estimation process, we apply bilinear interpolation to transfer crustal deformation information from geodetic sites to the grid nodes. The data include over a half century of triangulation measurements, over two decades of repeated electronic distance measurements, a decade of repeated very long baseline interferometry measurements, and several years of Global Positioning System measurements. Magnitudes for our estimated velocity vectors have formal standard errors ranging from 0.7 to 6.8 mm/yr. Our derived velocity field shows that (1) relative motion associated with the SAF exceeds 30 mm/yr and is distributed on the Earth's surface across a band (>100 km wide) that is roughly centered on this fault; (2) when velocities are expressed relative to a fixed North America plate, the motion within our primary study region has a mean orientation of N44°W ± 2° and the surface trace of the SAF is congruent in shape to nearby contours of constant speed yet this trace is oriented between 5° and 10° counterclockwise relative to these contours; and (3) large strain rates (shear rates > 150 nrad/yr and/or areal dilatation rates < −150 nstr/yr) exist near the Garlock fault, near the White Wolf fault, and in the Ventura basin.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Tholeiitic‐alkalic transition at subglacial volcanoes, Tuya region, British Columbia, Canada

Ash Mountain, South Tuya, and Tuya Butte are three small basaltic volcanoes in the Stikine volcanic belt of northern British Columbia. The volcanoes rise 700, 500, and 400 m above their bases and are about 3.2, 1.6, and 2.6 km 3 in volume, respectively. They began eruptive activity under several hundred meters of overlying glacial ice, or water in an ice‐impounded lake, and undegassed pillow lava was erupted and forms the bases of all three. Later, as the vents grew into shallow water, explosive phreatomagmatic activity erupted partly degassed glassy tuffs. Finally, when the volcano emerged through the surface of the ice or water (or the water was drained), degassed subaerial lava flows were erupted and were converted to assemblages of foreset‐bedded pillow breccia and pillow lava when subaerial flows crossed a shoreline and flowed into meltwater lakes. The undegassed subglacial pillow base of Ash Mountain is overlain by partly degassed pillows and hyaloclastite tuff cut by dikes; at South Tuya the pillow base is overlain by hyaloclastite tuffs and lenses of pillow lava; at Tuya Butte the pillow base is overlain by foreset‐bedded pillow lava, pillow breccias, and hyaloclastite tuffs, which in turn are overlain by subaerial lava flows composing a small shield volcano. The undegassed basal subglacial pillow lava of the three volcanoes contain 0.10 ± 0.01 wt % sulfur and ∼0.5 wt % H 2 O. The overlying partly degassed assemblages contain 0.06 ± 0.02% sulfur and ∼0.2% H 2 O at Ash Mountain, 0.07±0.01% sulfur at South Tuya, and 0.03±0.01% sulfur at Tuya Butte. The differences in the degree of degassing can be related to the nature of eruption and quenching and the distance of flow of the subaerial lava. When the volcanoes switched from subglacial to shallow water or subaerial eruptions, as shown by change to more explosive activity and then to subaerial lava flows (and by a marked reduction of sulfur in volcanic glass), the magma shifted from tholeiitic to alkalic composition. This transition occurs at each of the three volcanoes. The tholeiitic and alkalic magmas cannot be related by shallow crystal fractionation and apparently originated by differing degrees of deep melting at a mantle source. Prior to eruption the tholeiitic melts overlay alkalic melts in shallow chambers underlying each of the volcanoes because of their lower density and were, therefore, the first to erupt under subglacial conditions. As the volcano grew through the ice (or ice‐impounded water), the volcanic conduit vented to the atmosphere, producing a partial depressurization of the conduit and the subsurface chamber. This sudden reduction in confining pressure caused enhanced vesiculation of volatile saturated melts, particularly of the more volatile‐rich alkalic melts, causing them to rise to the top of the chamber and erupt.

British Columbia↗