USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Planets”

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 487 records · Page 27Linked to original sources
USGS2003

Oxygen-18 concentrations in recent precipitation and ice cores on the Tibetan Plateau

A detailed study of the climatic significance of δ 18 O in precipitation was completed on a 1500 km southwest-northeast transect of the Tibetan Plateau in central Asia. Precipitation samples were collected at four meteorological stations for up to 9 years. This study shows that the gradual impact of monsoon precipitation affects the spatial variation of δ 18 O-T relationship along the transect. Strong monsoon activity in the southern Tibetan Plateau results in high precipitation rates and more depleted heavy isotopes. This depletion mechanism is described as a precipitation “amount effect” and results in a poor δ 18 O-T relationship at both seasonal and annual scales. In the middle of the Tibetan Plateau, the effects of the monsoon are diminished but continue to cause a reduced correlation of δ 18 O and temperature at the annual scale. At the monthly scale, however, a significant δ 18 O-T relationship does exist. To the north of the Tibetan Plateau beyond the extent of the effects of monsoon precipitation, δ 18 O in precipitation shows a strong temperature dependence. δ 18 O records from two shallow ice cores and historic air temperature data were compared to verify the modern δ 18 O-T relationship. δ 18 O in Dunde ice core was positively correlated with air temperature from a nearby meteorological station in the north of the plateau. The δ 18 O variation in an ice core from the southern Plateau, however, was inversely correlated with precipitation amount at a nearby meteorological station and also the accumulation record in the ice core. The long-term variation of δ 18 O in the ice core record in the monsoon regions of the southern Tibetan Plateau suggest past monsoon seasons were probably more expansive. It is still unclear, however, how changes in large-scale atmosphere circulation might influence summer monsoon precipitation on the Tibetan Plateau.

Journal of Geophysical Research D: Atmospheres↗
USGS2003

Climate change and Arctic ecosystems: 2. Modeling, paleodata-model comparisons, and future projections

Large variations in the composition, structure, and function of Arctic ecosystems are determined by climatic gradients, especially of growing-season warmth, soil moisture, and snow cover. A unified circumpolar classification recognizing five types of tundra was developed. The geographic distributions of vegetation types north of 55°N, including the position of the forest limit and the distributions of the tundra types, could be predicted from climatology using a small set of plant functional types embedded in the biogeochemistry-biogeography model BIOME4. Several palaeoclimate simulations for the last glacial maximum (LGM) and mid-Holocene were used to explore the possibility of simulating past vegetation patterns, which are independently known based on pollen data. The broad outlines of observed changes in vegetation were captured. LGM simulations showed the major reduction of forest, the great extension of graminoid and forb tundra, and the restriction of low- and high-shrub tundra (although not all models produced sufficiently dry conditions to mimic the full observed change). Mid-Holocene simulations reproduced the contrast between northward forest extension in western and central Siberia and stability of the forest limit in Beringia. Projection of the effect of a continued exponential increase in atmospheric CO 2 concentration, based on a transient ocean-atmosphere simulation including sulfate aerosol effects, suggests a potential for larger changes in Arctic ecosystems during the 21st century than have occurred between mid-Holocene and present. Simulated physiological effects of the CO 2 increase (to >700 ppm) at high latitudes were slight compared with the effects of the change in climate.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Satellite-based detection of global urban heat-island temperature influence

This study utilizes a satellite-based methodology to assess the urban heat-island influence during warm season months for over 4400 stations included in the Global Historical Climatology Network of climate stations. The methodology includes local and regional satellite retrievals of an indicator of the presence green photosynthetically active vegetation at and around the stations. The difference in local and regional samples of the normalized difference vegetation index (NDVI) is used to estimate differences in mean air temperature. Stations classified as urban averaged 0.90??C (N. Hemisphere) and 0.92??C (S. Hemisphere) warmer than the surrounding environment on the basis of the NDVI-derived temperature estimates. Additionally, stations classified as rural averaged 0.19??C (N. Hemisphere) and 0.16??C (S. Hemisphere) warmer than the surrounding environment. The NDVI-derived temperature estimates were found to be in reasonable agreement with temperature differences observed between climate stations. The results suggest that satellite-derived data sets can be used to estimate the urban heat-island temperature influence on a global basis and that a more detailed analysis of rural stations and their surrounding environment may be necessary to assure that temperature trends derived from assumed rural environments are not influenced by changes in land use/land cover. Copyright 2002 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2003

Multiple oxygen and sulfur isotopic analyses on water-soluble sulfate in bulk atmospheric deposition from the southwestern United States

Sulfate is a major component of bulk atmospheric deposition (including dust, aerosol, fog, and rain). We analyzed sulfur and oxygen isotopic compositions of water-soluble sulfate from 40 sites where year-round dust traps collect bulk atmospheric deposition in the southwestern United States. Average sulfur and oxygen isotopic compositions (δ 34 S and δ 18 O) are 5.8 ± 1.4 (CDT) and 11.2 ± 1.9 (SMOW) (n = 47), respectively. Samples have an oxygen 17 anomaly (Δ 17 O), with an average value of 1.0 ± 0.6‰. Except for a weak positive correlation between δ 18 O and Δ 17 O values (r 2 ≈ 0.4), no correlation exists for δ 18 O versus δ 34 S, Δ 17 O versus δ 34 S, or any of the three isotopic compositions versus elevation of the sample site. Exceptional positive Δ 17 O values (up to 4.23‰) are found in samples from sites in the vicinity of large cities or major highways, and near-zero Δ 17 O values are found in samples close to dry lakes. Comparison of isotopic values of dust trap sulfate and desert varnish sulfate from the region reveals that varnish sulfate has average isotopic values that are ∼4.8‰ lower for δ 18 O, ∼2.1‰ higher for δ 34 S, and ∼0.3‰ lower for Δ 17 O than those of the present-day bulk deposition sulfate. Although other factors could cause the disparity, this observation suggests a possibility that varnish sulfate may have recorded a long-term atmospheric sulfate deposition during the Holocene or Pleistocene, as well as the differences between sulfur and oxygen isotopic compositions of the preindustrial bulk deposition sulfate and those of the industrial era.

Journal of Geophysical Research D: Atmospheres↗
USGS2004

Posteruption suspended sediment transport at Mount St. Helens: Decadal‐scale relationships with landscape adjustments and river discharges

Widespread landscape disturbance by the cataclysmic 1980 eruption of Mount St. Helens abruptly increased sediment supply in surrounding watersheds. The magnitude and duration of the redistribution of sediment deposited by the eruption as well as decades‐ to centuries‐old sediment remobilized from storage have varied chiefly with the style of disturbance. Posteruption suspended sediment transport has been greater and more persistent from zones of channel disturbance than from zones of hillslope disturbance. Despite the severe landscape disturbances caused by the eruption, relationships between discharge magnitudes and frequencies and suspended sediment transport have been remarkably consistent. Discharges smaller than mean annual flows generally have transported <5%, but locally ∼15%, of the annual suspended sediment loads, and infrequent ( p < 0.01), large floods have transported as much as 50% of the annual suspended sediment loads in a single day. However, moderate‐magnitude discharges (those greater than mean annual flows but less than 2‐year floods) have transported the greatest amounts of sediment from all disturbance zones. Such discharges have transported, on average, 60% to ∼95% of the annual suspended sediment loads, usually within cumulative periods of 1–3 weeks each year. Although small‐magnitude and large‐magnitude discharges have locally and episodically transported considerable amounts of suspended sediment, there has not been any notable change in the overall nature of the effective discharges; moderate‐magnitude flows have been the predominant discharges responsible for transporting the majority of suspended sediment during 20 years of posteruption landscape adjustment.

Washington↗
USGS2004

Paleogeodetic records of seismic and aseismic subduction from central Sumatran microatolls, Indonesia

We utilize coral microatolls in western Sumatra to document vertical deformation associated with subduction. Microatolls are very sensitive to fluctuations in sea level and thus act as natural tide gauges. They record not only the magnitude of vertical deformation associated with earthquakes (paleoseismic data), but also continuously track the long-term aseismic deformation that occurs during the intervals between earthquakes (paleogeodetic data). This paper focuses on the twentieth century paleogeodetic history of the equatorial region. Our coral paleogeodetic record of the 1935 event reveals a classical example of deformations produced by seismic rupture of a shallow subduction interface. The site closest to the trench rose 90 cm, whereas sites further east sank by as much as 35 cm. Our model reproduces these paleogeodetic data with a 2.3 m slip event on the interface 88 to 125 km from the trench axis. Our coral paleogeodetic data reveal slow submergence during the decades before and after the event in the areas of coseismic emergence. Likewise, interseismic emergence occurred before and after the 1935 event in areas of coseismic submergence. Among the interesting phenomenon we have discovered in the coral record is evidence of a large aseismic slip or "silent even" in 1962, 27 years after the 1935 event. Paleogeodetic deformation rates in the decades before, after, and between the 1935 and 1962 events have varied both temporally and spatially. During the 25 years following the 1935 event, submergence rates were dramatically greater than in prior decades. During the past four decades, however, rates have been lower than in the preceding decades, but are still higher than they were prior to 1935. These paleogeodetic records enable us to model the kinematics of the subduction interface throughout the twentieth century. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Structural model of the San Bernardino basin, California, from analysis of gravity, aeromagnetic, and seismicity data

The San Bernardino basin is an area of Quaternary extension between the San Jacinto and San Andreas Fault zones in southern California. New gravity data are combined with aeromagnetic data to produce two- and three-dimensional models of the basin floor. These models are used to identify specific faults that have normal displacements. In addition, aeromagnetic maps of the basin constrain strike-slip offset on many faults. Relocated seismicity, focal mechanisms, and a seismic reflection profile for the basin area support interpretations of the gravity and magnetic anomalies. The shape of the basin revealed by our interpretations is different from past interpretations, broadening its areal extent while confining the deepest parts to an area along the modern San Jacinto fault, west of the city of San Bernardino. Through these geophysical observations and related geologic information, we propose a model for the development of the basin. The San Jacinto fault-related strike-slip displacements started on fault strands in the basin having a stepping geometry thus forming a pull-apart graben, and finally cut through the graben in a simpler, bending geometry. In this model, the San Bernardino strand of the San Andreas Fault has little influence on the formation of the basin. The deep, central part of the basin resembles classic pull-apart structures and our model describes a high level of detail for this structure that can be compared to other pull-apart structures as well as analog and numerical models in order to better understand timing and kinematics of pull-apart basin formation. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Stress field variations in the Swiss Alps and the northern Alpine foreland derived from inversion of fault plane solutions

This study is devoted to a systematic analysis of the state of stress of the central European Alps and northern Alpine foreland in Switzerland based on focal mechanisms of 138 earthquakes with magnitudes between 1 and 5. The most robust feature of the results is that the azimuth of the minimum compressive stress, S 3 , is generally well constrained for all data subsets and always lies in the NE quadrant. However, within this quadrant, the orientation of S 3 changes systematically both along the structural strike of the Alpine chain and across it. The variation in stress along the mountain belt from NE to SW involves a progressive, counterclockwise rotation of S 3 and is most clear in the foreland, where it amounts to 45°–50°. This pattern of rotation is compatible with the disturbance to the stress field expected from the indentation of the Adriatic Block into the central European Plate, possibly together with buoyancy forces arising from the strongly arcuate structure of the Moho to the immediate west of our study area. Across the Alps, the variation in azimuth of S 3 is defined by a progressive, counterclockwise rotation of about 45° from the foreland in the north across the Helvetic domain to the Penninic nappes in the south and is accompanied by a change from a slight predominance of strike-slip mechanisms in the foreland to a strong predominance of normal faulting in the high parts of the Alps. The observed rotation can be explained by the perturbation of the large-scale regional stress by a local uniaxial deviatoric tension with a magnitude similar to that of the regional differential stress and with an orientation perpendicular to the strike of the Alpine belt. The tensile nature and orientation of this stress is consistent with the “spreading” stress expected from lateral density changes due to a crustal root beneath the Alps.

European Alps↗
USGS2004

Transient volcano deformation sources imaged with interferometric synthetic aperture radar: Application to Seguam Island, Alaska

Thirty interferometric synthetic aperture radar (InSAR) images, spanning various intervals during 1992&ndash;2000, document coeruptive and posteruptive deformation of the 1992&ndash;1993 eruption on Seguam Island, Alaska. A procedure that combines standard damped least squares inverse methods and collective surfaces, identifies three dominant amorphous clusters of deformation point sources. Predictions generated from these three point source clusters account for both the spatial and temporal complexity of the deformation patterns of the InSAR data. Regularized time series of source strength attribute a distinctive transient behavior to each of the three source clusters. A model that combines magma influx, thermoelastic relaxation, poroelastic effects, and petrologic data accounts for the transient, interrelated behavior of the source clusters and the observed deformation. Basaltic magma pulses, which flow into a storage chamber residing in the lower crust, drive this deformational system. A portion of a magma pulse is injected into the upper crust and remains in storage during both coeruption and posteruption intervals. This injected magma degasses and the volatile products accumulate in a shallow poroelastic storage chamber. During the eruption, another portion of the magma pulse is transported directly to the surface via a conduit roughly centered beneath Pyre Peak on the west side of the island. A small amount of this magma remains in storage during the eruption, and posteruption thermoelastic contraction ensues. This model, made possible by the excellent spatial and temporal coverage of the InSAR data, reveals a relatively simple system of interrelated predictable processes driven by magma dynamics.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Crystallographic controls on the frictional behavior of dry and water-saturated sheet structure minerals

We compare the frictional strengths of 17 sheet structure mineral powders, measured under dry and water-saturated conditions, to identify the factors that cause many of them to be relatively weak. The dry coefficient of friction μ ranges upward from 0.2 for graphite, leveling off at 0.8 for margarite, clintonite, gibbsite, kaolinite, and lizardite. The values of μ (dry) correlate directly with calculated (001) interlayer bond strengths of the minerals. This correlation occurs because shear becomes localized along boundary and Riedel shears and the platy minerals in them rotate into alignment with the shear planes. For those gouges with μ (dry) < 0.8, shear occurs by breaking the interlayer bonds to form new cleavage surfaces. Where μ (dry) = 0.8, consistent with Byerlee's law, the interlayer bonds are sufficiently strong that other frictional processes dominate. The transition in dry friction mechanisms corresponds to calculated surface energies of 2–3 J/m 2 . Adding water causes μ to decrease for every mineral tested except graphite. If the minerals are separated into groups with similar crystal structures, μ (wet) increases with increasing interlayer bond strength within each group. This relationship also holds for the swelling clay montmorillonite, whose water-saturated strength is consistent with the strengths of nonswelling clays of similar crystal structure. Water in the saturated gouges forms thin, structured films between the plate surfaces. The polar water molecules are bonded to the plate surfaces in proportion to the mineral's surface energy, and μ (wet) reflects the stresses required to shear through the water films.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults

We argue that key features of thrust earthquake triggering, inhibition, and clustering can be explained by Coulomb stress changes, which we illustrate by a suite of representative models and by detailed examples. Whereas slip on surface-cutting thrust faults drops the stress in most of the adjacent crust, slip on blind thrust faults increases the stress on some nearby zones, particularly above the source fault. Blind thrusts can thus trigger slip on secondary faults at shallow depth and typically produce broadly distributed aftershocks. Short thrust ruptures are particularly efficient at triggering earthquakes of similar size on adjacent thrust faults. We calculate that during a progressive thrust sequence in central California the 1983 M w = 6.7 Coalinga earthquake brought the subsequent 1983 M w = 6.0 Nuñez and 1985 M w = 6.0 Kettleman Hills ruptures 10 bars and 1 bar closer to Coulomb failure. The idealized stress change calculations also reconcile the distribution of seismicity accompanying large subduction events, in agreement with findings of prior investigations. Subduction zone ruptures are calculated to promote normal faulting events in the outer rise and to promote thrust-faulting events on the periphery of the seismic rupture and its downdip extension. These features are evident in aftershocks of the 1957 M w = 9.1 Aleutian and other large subduction earthquakes. We further examine stress changes on the rupture surface imparted by the 1960 M w = 9.5 and 1995 M w = 8.1 Chile earthquakes, for which detailed slip models are available. Calculated Coulomb stress increases of 2–20 bars correspond closely to sites of aftershocks and postseismic slip, whereas aftershocks are absent where the stress drops by more than 10 bars. We also argue that slip on major strike-slip systems modulates the stress acting on nearby thrust and strike-slip faults. We calculate that the 1857 M w = 7.9 Fort Tejon earthquake on the San Andreas fault and subsequent interseismic slip brought the Coalinga fault ∼1 bar closer to failure but inhibited failure elsewhere on the Coast Ranges thrust faults. The 1857 earthquake also promoted failure on the White Wolf reverse fault by 8 bars, which ruptured in the 1952 M w = 7.3 Kern County shock but inhibited slip on the left-lateral Garlock fault, which has not ruptured since 1857. We thus contend that stress transfer exerts a control on the seismicity of thrust faults across a broad spectrum of spatial and temporal scales.

California↗
USGS2004

Strain accumulation across the Coast Ranges at the latitude of San Francisco, 1994-2000

A 66-monument geodetic array spanning the Coast Ranges near San Francisco has been surveyed more than eight times by GIPS between late 1993 and early 2001. The measured horizontal velocities of the monuments are well represented by uniform, right-lateral, simple shear parallel to N29°W. (The local strike of the San Andreas Fault is ∼N34°W.) The observed areal dilatation rate of 6.9 ± 10.0 nstrain yr −1 (quoted uncertainty is one standard deviation and extension is reckoned positive) is not significantly different from zero, which implies that the observed strain accumulation could be released by strike-slip faulting alone. Our results are consistent with the slip rates assigned by the Working Group on California Earthquake Probabilities [2003] to the principal faults (San Gregorio, San Andreas, Hayward-Rodgers Creek, Calaveras-Concord-Green Valley, and Greenville Faults) cutting across the GPS array. The vector sum of those slip rates is is 39.8 ± 2.6 mm yr −1 N29.8°W ± 2.8°, whereas the motion across the GPS array (breadth 120 km) inferred from the uniform strain rate approximation is 38.7 ± 1.2 mm yr −1 N29.0°W ± 0.9° right-lateral shear and 0.4 ± 0.9 mm yr −1 N61°E ± 0.9° extension. We interpret the near coincidence of these rates and the absence of significant accumulation of areal dilatation to imply that right-lateral slip on the principal faults can release the accumulating strain; major strain release on reverse faults subparallel to the San Andreas Fault within the Coast Ranges is not required.

California↗
USGS2004

Recalculated probability of M ≥ 7 earthquakes beneath the Sea of Marmara, Turkey

New earthquake probability calculations are made for the Sea of Marmara region and the city of Istanbul, providing a revised forecast and an evaluation of time-dependent interaction techniques. Calculations incorporate newly obtained bathymetric images of the North Anatolian fault beneath the Sea of Marmara [Le Pichon et al., 2001; Armijo et al., 2002]. Newly interpreted fault segmentation enables an improved regional A.D. 1500-2000 earthquake catalog and interevent model, which form the basis for time-dependent probability estimates. Calculations presented here also employ detailed models of coseismic and postseismic slip associated with the 17 August 1999 M = 7.4 Izmit earthquake to investigate effects of stress transfer on seismic hazard. Probability changes caused by the 1999 shock depend on Marmara Sea fault-stressing rates, which are calculated with a new finite element model. The combined 2004-2034 regional Poisson probability of M≥7 earthquakes is ~38%, the regional time-dependent probability is 44 ± 18%, and incorporation of stress transfer raises it to 53 ± 18%. The most important effect of adding time dependence and stress transfer to the calculations is an increase in the 30 year probability of a M ??? 7 earthquake affecting Istanbul. The 30 year Poisson probability at Istanbul is 21%, and the addition of time dependence and stress transfer raises it to 41 ± 14%. The ranges given on probability values are sensitivities of the calculations to input parameters determined by Monte Carlo analysis; 1000 calculations are made using parameters drawn at random from distributions. Sensitivities are large relative to mean probability values and enhancements caused by stress transfer, reflecting a poor understanding of large-earthquake aperiodicity.

Sea Of Marmara↗
USGS2004

Interseismic strain and rotation rates in the northeast Mojave domain, eastern California

The northeast Mojave domain, a type locality for bookshelf faulting, is a region of east striking, left-lateral faults in the northeast corner of the Mojave block, a block otherwise dominated by ∼N40°W striking, right-lateral faults. Paleomagnetic evidence suggests that blocks within the domain have rotated clockwise about a vertical axis as much as 60° since 12.8 Ma [Schermer et al., 1996] . In 1994, and again in 2002, the U.S. Geological Survey surveyed an array of 14 geodetic monuments distributed across the northeast Mojave domain. The 2002 survey results were adjusted to remove the coseismic offsets imposed by the nearby Hector Mine earthquake (16 October 1999, M w = 7.1). The adjusted deformation across the array appears to be uniform and can be approximated by the principal strain rates ε 1 = 28.9 ± 9.1 N77.2°W ± 4.8° and ε 2 = −48.2 ± 8.9 N12.8°E ± 4.8° nstrain yr −1 ; extension reckoned positive, and quoted uncertainties are standard deviations. That strain accumulation could be released by slip on faults striking N32°W but not by bookshelf faulting on the east striking faults alone. The vertical axis rotation rate of the northeast Mojave domain as a whole relative to fixed North America is 71.0 ± 6.4 nrad yr −1 (4.07° ± 0.37° Myr −1 ) clockwise, about twice the maximum tensor shear strain rate. The observed rotation rate acting over 12.8 Myr would produce a clockwise rotation of 52.1° ± 4.7°, exclusive of possible coseismic rotations. That rotation is in rough agreement with the paleomagnetic rotation accumulated in the individual fault blocks within the northeast Mojave domain since 12.8 Ma.

California↗
USGS2004

Contemporary tectonic deformation of the Basin and Range province, western United States: 10 years of observation with the Global Positioning System

[1] We have estimated patterns and rates of crustal movement across 800 km of the Basin and Range at ∼39° north latitude with Global Positioning System surveys in 1992, 1996, 1998, and 2002. The total rate of motion tangent to the small circle around the Pacific‐North America pole of rotation is 10.4 ± 1.0 mm/yr, and motion normal to this small circle is 3.9 ± 0.9 mm/yr compared to the east end of our network. On the Colorado Plateau the east end of our network moves by ∼1–2 mm/yr westerly with respect to North America. Transitions in strain rates delimit six major tectonic domains within the province. These deformation zones coincide with areas of modern seismicity and are, from east to west, (1) east‐west extension in the Wasatch Fault zone, (2) low rate east‐west extension centered near the Nevada‐Utah border, (3) low rate east‐west contraction between 114.7°W and 117.9°W, (4) extension normal to and strike‐slip motion across the N10°E striking Central Nevada Seismic Zone, (5) right lateral simple shear oriented N13°W inside the Walker Lane Belt, and (6) shear plus extension near the Sierra Nevada frontal faults. Concentration of shear and dilatational deformation across the three westernmost zones suggests that the Walker Lane Belt lithosphere is rheologically weak. However, we show that linear gradients in viscosity and gravitational potential energy can also effectively concentrate deformation. In the Basin and Range, gradients in gravitational potential are spatially anticorrelated with dilatational strain rates, consistent with the presence of horizontal variations in viscosity of the lithosphere.

Arizona, California, Nevada, Utah↗
USGS2005

Magmatic effects of the Cobb hot spot on the Juan de Fuca Ridge

The interaction of the Juan de Fuca Ridge with the Cobb hot spot has had a considerable influence on the magmatism of the Axial Segment of the ridge, the second-order segment that overlies the hot spot. In addition to the construction of the large volcanic edifice of Axial Seamount, the Axial Segment has shallow bathymetry and a prevalence of constructional volcanic features along its 100-km length, suggesting that hot spot-derived magmas supplement and oversupply the ridge. Lavas are generally more primitive at Axial Seamount and more evolved in the Axial Segment rift zones, suggesting that fractional crystallization is enhanced with increasing distance from the hot spot because of a reduced magma supply and more rapid cooling. Although the Cobb hot spot is not an isotopically enriched plume, it produces lavas with some distinct geochemical characteristics relative to normal mid-ocean ridge basalt, such as enrichments in alkalis and highly incompatible trace elements, that can be used as tracers to identify the presence and prevalence of the hot spot influence along the ridge. These characteristics are most prominent at Axial Seamount and decline in gradients along the Axial Segment. The physical model that can best explain the geochemical observations is a scenario in which hot spot and mid-ocean ridge basalt (MORB) magmas mix to varying degrees, with the proportions controlled by the depth to the MORB source. Modeling of two-component mixing suggests that MORB is the dominant component in most Axial Segment basalts. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Noise in two-color electronic distance meter measurements revisited

Frequent, high-precision geodetic data have temporally correlated errors. Temporal correlations directly affect both the estimate of rate and its standard error; the rate of deformation is a key product from geodetic measurements made in tectonically active areas. Various models of temporally correlated errors are developed and these provide relations between the power spectral density and the data covariance matrix. These relations are applied to two-color electronic distance meter (EDM) measurements made frequently in California over the past 15-20 years. Previous analysis indicated that these data have significant random walk error. Analysis using the noise models developed here indicates that the random walk model is valid for about 30% of the data. A second 30% of the data can be better modeled with power law noise with a spectral index between 1 and 2, while another 30% of the data can be modeled with a combination of band-pass-filtered plus random walk noise. The remaining 10% of the data can be best modeled as a combination of band-pass-filtered plus power law noise. This band-pass-filtered noise is a product of an annual cycle that leaks into adjacent frequency bands. For time spans of more than 1 year these more complex noise models indicate that the precision in rate estimates is better than that inferred by just the simpler, random walk model of noise.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Upper crustal structure of southwestern British Columbia from the 1998 Seismic Hazards Investigation in Puget Sound

This paper applies nonlinear three-dimensional travel time tomography to refraction data recorded during the 1998 Seismic Hazards Investigation in Puget Sound (SHIPS) to derive the first large-scale, high-resolution upper crustal velocity model for southwestern British Columbia. A minimum structure P wave velocity model is constructed using 175,000 first arrival travel times picked from data recorded by 58 temporary onshore stations. The model details forearc crustal structures related to terrane accretion and subsequent basin formation to a depth of about 10 km. The Metchosin igneous complex (correlative with the Eocene Crescent-Siletz Terrane in Washington) is imaged as a laterally extensive WNW trending high-velocity anomaly underlying southernmost Vancouver Island and much of the Strait of Juan du Fuca. Northeast of the Strait of Georgia, the southwesterly dip of the contact between the Wrangellia terrane rocks of Vancouver Island and the Coast Plutonic Complex suggests Wrangellia rocks are down-faulted against the plutonic complex. At the southwestern end of the Strait of Juan de Fuca, the 50 km long WNW trending Clallam basin has a maximum thickness of 5-6 km. Near the eastern end of the Strait of Juan de Fuca, Port Townsend basin has an inferred thickness of approximately 4-5 km. The southern end of the 9 km thick Georgia basin is bounded by a high-velocity basement ridge. Beneath the Strait of Georgia, clusters of well-located earthquakes have a prominent NW trend and coincide spatially with rapid lateral velocity changes. Clusters of microearthquakes there are associated with the intersection of several east trending structural highs within this NW trend. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗