USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research- Oceans”

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 145 records · Page 8Linked to original sources

Lu-Hf constraints on the evolution of lunar basalts

Very low Ti basalts and green glass samples from the moon show high Lu/Hf ratios and low Hf concentrations. Low-Ti lunar basalts show high and variable Lu/Hf ratios and higher Hf concentrations, whereas high-Ti lunar basalts show low Lu/Hf ratios and high Hf concentrations. KREEP basalts have constant Lu/Hf ratios and high but variable Hf concentrations. Using the Lu-Hf behavior as a constraint, we propose a model for the mare basalts evolution. This constraint requires extensive crystallization of the primary lunar magma ocean prior to formation of the lunar mare basalt sources and the KREEP basalts. Mare basalts are produced by the melting of the cumulate rocks, and KREEP basalts represent the residual liquid of the magma ocean. Lu and Hf concentrations and the Hf isotopic data of lunar rocks suggest that assimilation cannot be accepted as a major process to explain the diversity of the lunar mare basalts. The urKREEP hypothesis is also unnecessary. Both high- and low-Ti basalts show enough iron enrichment to be regarded as melting products of the last stage cumulate rocks from the lunar magma ocean. The KREEP basalts are also rich in iron and may be regarded as the final, residual liquid left after the crystallization of the major portion of the primary lunar magma ocean.

Journal of Geophysical Research Solid Earth↗

Three‐dimensional model of Hellenic Arc deformation and origin of the Cretan uplift

[1] The Hellenic Arc of Greece is the most seismically active part of Europe, but little is know about its mechanics. We modeled deformation along the arc using a finite element model. The model was intended to capture large‐scale 3‐D structure of Nubian plate subduction beneath the Aegean block and its deformational consequences. The shape of the interface was developed using mapped traces at the surface and earthquake hypocenters at depth. Model block motions were constrained by recent compilations of GPS velocity vectors. We simulated a 10 ka period of convergence between Nubia and the Aegean and calculated the strain field in the overriding plate as well as the spatial distribution and orientation of differential stress (∣ σ 1 − σ 3 ∣). From these calculations we derived testable quantities such as the expected seismic moment rate on the interplate contact, uplift pattern, and distribution of strain modes. Our relatively simple model broadly reproduced observed uplift patterns, earthquake activity, and loci of extension and contraction. The model showed a localization of uplift near the island of Crete, where the fastest Aegean uplift rates are well documented. Comparison of calculated expected seismic moment and observed earthquake catalogs implies a nearly fully coupled interplate contact. On the basis of our modeling results, we suggest that south Aegean deformation is driven primarily by the fast moving (∼33 mm a −1 ) Aegean upper plate overriding a nearly stalled (∼5 mm a −1 ) Nubian lower plate. This tectonic setting thus more closely resembles a continental thrust than it does a typical oceanic subduction zone.

Aegean Sea, Hellenic Arc↗

Early Cretaceous paleolatitude of the Yukon-Koyukuk province, Alaska

We report results from a paleomagnetic study of the Yukon-Koyukuk province, a key region for reconstructing the Cretaceous paleogeography of northern Alaska and northeastern Asia. The province lies between the displaced continental fragment of Arctic Alaska and the accreted terranes of southern Alaska. Although Lower Cretaceous volcanogenic rocks of the Yukon-Koyukuk province preserve primary remanent magnetizations, the primary component is masked by a pervasive overprint. The over-printing magnetization, which has dual polarity, was acquired after a Late Cretaceous orogeny left the region highly deformed. By fitting planes to tilt-corrected thermal demagnetization data and solving for the best intersection of the corresponding remagnetization circles, we separated the primary component, presumably a vestige of the original thermoremanent magnetization, from the overprinting component of the natural remanence. The primary magnetization corresponds to a mean paleopole of Early Cretaceous age at 60.2°N, 159°E. The mean paleolatitude is 69°N with the 95% confidence ellipse ranging from 60° to 75°N. This analysis was successful because the tilt corrections of the sampled strata had sufficient variation to ensure a strong intersection of the remagnetization circles. We infer that Lower Cretaceous rocks of the Brooks Range and the Yukon-Koyukuk province were once separated by a wide ocean basin that was closed by subduction, probably during the Cretaceous. The northern limb of the Angayucham terrane, which lies between the Brooks Range and the Yukon-Koyukuk province, may contain remnants of the now-closed ocean basin.

Journal of Geophysical Research Solid Earth↗

Reflected and mode-converted seismic waves within the shallow Aleutian subduction zone, southern Kenai Peninsula, Alaska

Pronounced secondary phases observed in local recordings of quarry shots and earthquakes on the southern Kenai Peninsula, Alaska, are identified as reflected P and S and converted S-to-P phases originating within four depth ranges: in the upper few kilometers of the Cook Inlet Tertiary basin, at midcrustal depths within the overthrust North American plate, at about 35 km depth near the top of the Wadati-Benioff seismic zone in proximity to the inferred interplate megathrust, and at about 5–10 km below the megathrust in the subducted Pacific plate. Beneath the landward margin of the accreted Chugach terrane, the mid-upper plate (MUP) discontinuity dips as steeply as 20°–30° to the west-northwest over the depth range 12–18 km. At shallower depths it flattens and possibly arches with the crest at about 10 km depth. Similar midcrustal reflectors have been imaged about 125 km to the southwest and about 350 km to the northeast along the structural trend of the Chugach terrane. The extensive reflectors may have a common origin, possible caused by the presence of underplated rocks. Relative amplitudes and polarities of the secondary phases originating at MUP depths provide few constraints on the nature of the discontinuity. It appears that the MUP discontinuity is seismically inactive and does not represent a brittle-ductile transition zone within the upper plate. The two converted S -to- P phases generated near the top of the subducted plate could indicate a low-velocity zone associated with subducted oceanic crust.

Journal of Geophysical Research Solid Earth↗

Deuterium in interstitial water from deep-sea cores

As part of the Joint Oceanographic Institutions Deep Earth Sampling project, the interstitial waters of cores from 69 holes were sampled for deuterium analysis. Sixteen of the cores penetrated sediments as old as Eocene, and several sampled Cretaceous sediments, which allowed us to examine changes in the deuterium content of the oceans with time. Deuterium is shown to be a conservative constituent of the pore water. Its abundance in the pore fluids can be changed by diffusion, but the rate is slow, and corrections for this effect have been made. Changes in the abundance of deuterium can be related to changes in the amount of ice stored in continental glaciers, inasmuch as precipitation in the form of snow is highly depleted in deuterium compared with the oceans. Many of the cores show a change in isotopic composition of samples from early to late Miocene that can be ascribed to the buildup of the Antarctic ice sheets. After correcting for the role of diffusion in reducing the isotopic contrast between samples from a single core, we estimate an increase of 10 per mil (‰) δD (corresponding to a° δ 18 O change of about 1.2‰) between the early and late Miocene. A similiar analysis of Pleistocene to Holocene changes indicates a δD rise of 8‰ during the time of maximum continental ice, which corresponds to a δ 18 O increase of about 1.0‰. On the basis of limited data, we find no δD change in the oceans from Cretaceous to Miocene.

Journal of Geophysical Research Solid Earth↗

Finite element model predictions of static deformation from dislocation sources in a subduction zone: Sensitivities to homogeneous, isotropic, Poisson-solid, and half-space assumptions

Dislocation models can simulate static deformation caused by slip along a fault. These models usually take the form of a dislocation embedded in a homogeneous, isotropic, Poisson-solid half-space (HIPSHS). However, the widely accepted HIPSHS assumptions poorly approximate subduction zone systems of converging oceanic and continental crust. This study uses three-dimensional finite element models (FEMs) that allow for any combination (including none) of the HIPSHS assumptions to compute synthetic Green's functions for displacement. Using the 1995 Mw = 8.0 Jalisco-Colima, Mexico, subduction zone earthquake and associated measurements from a nearby GPS array as an example, FEM-generated synthetic Green's functions are combined with standard linear inverse methods to estimate dislocation distributions along the subduction interface. Loading a forward HIPSHS model with dislocation distributions, estimated from FEMs that sequentially relax the HIPSHS assumptions, yields the sensitivity of predicted displacements to each of the HIPSHS assumptions. For the subduction zone models tested and the specific field situation considered, sensitivities to the individual Poisson-solid, isotropy, and homogeneity assumptions can be substantially greater than GPS. measurement uncertainties. Forward modeling quantifies stress coupling between the Mw = 8.0 earthquake and a nearby Mw = 6.3 earthquake that occurred 63 days later. Coulomb stress changes predicted from static HIPSHS models cannot account for the 63-day lag time between events. Alternatively, an FEM that includes a poroelastic oceanic crust, which allows for postseismic pore fluid pressure recovery, can account for the lag time. The pore fluid pressure recovery rate puts an upper limit of 10-17 m2 on the bulk permeability of the oceanic crust. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Influence of welded boundaries in anelastic media on energy flow, and characteristics of P, S-I, and S-II waves: Observational evidence for inhomogeneous body waves in low-loss solids

A general computer code, developed to calculate anelastic reflection-refraction coefficients, energy flow, and the physical characteristics for general P , S -I, and S -II waves, quantitatively describes physical characteristics for wave fields in anelastic media that do not exist in elastic media. Consideration of wave fields incident on boundaries between anelastic media shows that scattered wave fields experience reductions in phase and energy speeds, increases in maximum attenuation and Q −1 , and directions of maximum energy flow distinct from phase propagation. Each of these changes in physical characteristics are shown to vary with angle of incidence. Finite relaxation times for anelastic media result in energy flow due to interaction of superimposed radiation fields and contribute to energy flow across anelastic boundaries for all angles of incidence. Agreement of theoretical and numerical results with laboratory measurements argues for the validity of the theoretical and numerical formulations incorporating inhomogeneous wave fields. The agreement attests to the applicability of the model and helps confirm the existence of inhomogeneous body waves and their associated set of distinct physical characteristics in the earth. The existence of such body waves in layered, low-loss anelastic solids implies the need to reformulate some seismological models of the earth. The exact anelastic formulation for a liquid-solid interface with no low-loss approximations predicts the existence of a range of angles of incidence or an anelastic Rayleigh window, through which significant amounts of energy are transmitted across the boundary. The window accounts for the discrepancy apparent between measured reflection data presented in early textbooks and predictions based on classical elasticity theory. Characteristics of the anelastic Rayleigh window are expected to be evident in certain sets of wide-angle, ocean-bottom reflection data and to be useful in estimating Q −1 for some ocean bottom reflectors.

Journal of Geophysical Research B: Solid Earth↗

Hydrothermal circulation at the Cleft-Vance overlapping spreading center: Results of a magnetometric resistivity survey

We report on a magnetometric resistivity sounding carried out in the overlapping spreading center between the Cleft and Vance segments of the Juan de Fuca Ridge. The data collected reveal a strong three dimensionality in the crustal electrical resistivity structure on wavelengths of a few kilometers. Areas of reduced crustal electrical resistivities, with values approaching that of seawater, are seen beneath the neovolcanic zones of both active spreading centers. We interpret these reduced resistivities as evidence of active hydrothermal circulation within the uppermost 1 km of hot, young oceanic crust.

Journal of Geophysical Research B: Solid Earth↗

Tectonic setting of the Yukon-Koyukuk basin and its borderlands, western Alaska

The Yukon-Koyukuk basin of western Alaska is composed of an arcuate belt of Jurassic and Early Cretaceous subduction-related volcanic and plutonic rocks (Koyukuk terrane) flanked by deep subbasins filled with mid-Cretaceous terrigenous sedimentary rocks. The basin is bordered on three sides by metamorphosed Proterozoic and Paleozoic continental rocks (Seward, Arctic Alaska, and Ruby terranes) and is separated from the metamorphic borderlands by a narrow, highly tectonized belt of oceanic crust and mantle rocks (composite Angayucham-Tozitna terrane). The oceanic and mantle rocks, which dip inward beneath the basin and are thrust outward onto the borderlands, are divided into three separate thrust panels: (1) a structurally lowest panel (Slate Creek) composed of phyllite and metagraywacke of probable Devonian age, (2) a middle panel (Narvak) composed of imbricated basalt, chert, and gabbro of Devonian to Early Jurassic age, and (3) a structurally highest panel (Kanuti) composed of gabbro and peridotite of probable Middle and Late Jurassic age. The three thrust panels appear to represent a reversely stacked sequence that progresses from continental slope deposits in the lower panel to cumulus and mantle peridotites in the upper. Metamorphic mineral K-Ar ages from garnet amphibolite on the sole of the upper panel suggest that the upper panel was emplaced on the middle panel in the Middle to Late Jurassic. Subsequent accretion of all three panels to the continental rocks of the borderlands occurred in the latest Jurassic and Early Cretaceous, synchronous with arc volcanism within the basin. Arc volcanism waned and the accretionary phase ended in the middle of Early Cretaceous time. Uplift and erosion of the metamorphic borderlands and the obducted oceanic rocks began in late Early Cretaceous and was accompanied by the rapid filling of two flanking subbasins with turbiditic sediments. In the latest Early Cretaceous and early Late Cretaceous, shallow marine and nonmarine conglomerates were deposited around the margins of the basin, and a prograding delta was built out from the southeast margin of the basin across the turbiditic subbasins and the remnant volcanic arc. In the Late Cretaceous, western Alaska was subjected to strong east-west compression which severely deformed both the Yukon-Koyukuk basin and the borderlands. Several widespread magmatic episodes in the mid- and Late Cretaceous and in early Tertiary transgress the basin boundaries and stitch together the accreted arc and oceanic terranes and the continental borderlands.

Journal of Geophysical Research Solid Earth↗

Pb, Sr, and Nd isotopes in seamount basalts from the Juan de Fuca Ridge and Kodiak-Bowie seamount chain, northeast Pacific

Pb, Sr, and Nd isotopic ratios and their parent/daughter element concentrations for 28 basalts from 10 hotspot and nonhotspot seamounts are reported. Nd and Sr isotopic compositions ( 143 Nd/ 144 Nd = 0.51325–0.51304; 87 Sr/ 86 Sr = 0.70237–0.70275) plot in the envelope for Juan de Fuca-Gorda ridge basalts with tholeiitic basalts showing more depleted sources and a better negative correlation than transitional to alkalic basalts. Pb isotopic ratios in tholeiitic and alkalic basalts overlap ( 206 Pb/ 204 Pb = 18.29–19.44) and display a trend toward more radiogenic Pb in alkalic basalts. The isotopic data for hotspot and nonhotspot basalts are indistinguishable and correlate broadly with rock composition, implying that they are controlled by partial melting. The isotopic variation in the seamount basalts is about 60% (Nd-Sr) to 100% (Pb) of that in East Pacific Rise basalts and is interpreted as a lower limit for the magnitude of mantle heterogeneity in the northeast Pacific. The data indicate absence of a chemically distinct plume component in the linear seamount chains and strongly suggest an origin from mid-ocean ridge basalt-like east Pacific mantle.

Journal of Geophysical Research Solid Earth↗

Density structure of the lithosphere in the southwestern United States and its tectonic significance

We calculate a density model of the lithosphere of the southwestern United States through an integrated analysis of gravity, seismic refraction, drill hole, and geological data. Deviations from the average upper mantle density are as much as ?? 3%. A comparison with tomographic images of seismic velocities indicates that a substantial part (>50%) of these density variations is due to changes in composition rather than temperature. Pronounced mass deficits are found in the upper mantle under the Basin and Range Province and the northern part of the California Coast Ranges and adjacent ocean. The density structure of the northern and central/southern Sierra Nevada is remarkably different. The central/southern part is anomalous and is characterized by a relatively light crust underlain by a higher-density upper mantle that may be associated with a cold, stalled subducted plate. High densities are also determined within the uppermost mantle beneath the central Transverse Ranges and adjoining continental slope. The average density of the crystalline crust under the Great Valley and western Sierra Nevada is estimated to be up to 200 kg m~3 higher than the regional average, consistent with tectonic models for the obduction of oceanic crust and uppermost mantle in this region.

Journal of Geophysical Research B: Solid Earth↗

The crustal structure of the Wrangellia Terrane along the East Glenn Highway, eastern‐southern Alaska

Recently acquired seismic refraction data from eastern‐southern Alaska provide new information on the structure and composition of the Wrangellia and adjacent terranes. The data comprise a 160‐km‐long refraction profile along the East Glenn (Tok‐Cutoff) Highway that was collected as part of the U.S. Geological Survey's multidisciplinary Trans‐Alaska Crustal Transect program. The upper 3 km of the Wrangellia terrane and associated rocks is characterized by low compressional wave velocities (V p = 1.9, 3.3, 4.6, 5.6 km s −1 ) and high‐velocity gradients common to most onshore seismic refraction profiles. There is also clear seismic expression of the West Fork fault system as a steep, down‐to‐the‐southwest fault that separates the Peninsular terrane to the southwest and the metamorphic complex of Gulkana River to the northeast. In contrast, no seismic expression occurs for the Paxson Lake fault, which separates the Wrangellia terrane from the metamorphic complex of Gulkana River. Adjacent to the Denali fault, within the Wrangellia terrane, two high‐velocity bodies (V p = 6.6 km s −1 ) occur in the upper crust. One of these extends to ∼10‐km depth and correlates with a late Paleozoic dioritic complex, suggesting that the Wrangellia terrane is at least 10 km thick in this part of Alaska. From 5 to 23 km depth, the crust appears seismically homogeneous, with velocity increasing from V p = 6.2 to V p = 6.6 km s −1 . Beneath this level, the crust is less well resolved, although evidence exists for a low‐velocity zone between 23 and 26 km and a possible southwest dipping interface at 35 km. No identifiable mantle refraction or reflection is observed, possibly indicating a crust as thick as 55 km. The relatively low seismic velocities in the upper 23 km of the crust compare favorably with laboratory‐measured velocities on pelitic schists and intermediate‐composition plutonic rocks (granites and granodiorites), both of which are recognized in Wrangellia. We interpret the seismic velocities to indicate that silicic‐to‐intermediate‐composition rocks are important constituents of the basement of this part of Wrangellia. Geologic evidence indicates that the Alaskan part of the Wrangellia terrane is a Paleozoic and Mesozoic island arc: our seismic evidence indicates it may have been built mostly on continental crust as opposed to the fragment of Wrangellia from Vancouver Island which was probably built on oceanic crust.

Alaska↗

Olivine friction at the base of oceanic seismogenic zones

We investigate the strength and frictional behavior of olivine aggregates at temperatures and effective confining pressures similar to those at the base of the seismogenic zone on a typical ridge transform fault. Triaxial compression tests were conducted on dry olivine powder (grain size ???60 ??m) at effective confining pressures between 50 and 300 MPa (using Argon as a pore fluid), temperatures between 600??C and 1000??C, and axial displacement rates from 0.06 to 60 ??m/s (axial strain rates from 3 ?? 10-6 to 3 ?? 10-3 s-1). Yielding shows a negative pressure dependence, consistent with predictions for shear enhanced compaction and with the observation that samples exhibit compaction during the initial stages of the experiments. A combination of mechanical data and microstructural observations demonstrate that deformation was accommodated by frictional processes. Sample strengths were pressure-dependent and nearly independent of temperature. Localized shear zones formed in initially homogeneous aggregates early in the experiments. The frictional response to changes in loading rate is well described by rate and state constitutive laws, with a transition from velocity-weakening to velocity-strengthening at 1000??C. Microstructural observations and physical models indicate that plastic yielding of asperities at high temperatures and low axial strain rates stabilizes frictional sliding. Extrapolation of our experimental data to geologic strain rates indicates that a transition from velocity weakening to velocity strengthening occurs at approximately 600??C, consistent with the focal depths of earthquakes in the oceanic lithosphere. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Tectonic stresses in the lithosphere: constraints provided by the experimental deformation of rocks.

The strengths of rocks clearly place an upper limit on the stress that can be sustained by the upper half of the lithosphere. Laboratory data on rock rheology are generally lacking at intermediate temperatures and pressures on the important rock types expected in the lithosphere, so a definitive accounting of the strength distribution with depth in the upper lithosphere is presently unattainable. Analogies are drawn between the fragmentary strength data on silicates at intermediate temperature and the more extensive experimental data on marble and limestone, and several tentative conclusions are drawn: First, brittle processes, such as faulting and cataclasis, are expected to control rock strength at low pressures and temperatures. The strengths associated with these brittle mechanisms increase rapidly with increasing effective pressure and are relatively insensitive to temperature and strain rate. Second, the transitions between brittle and ductile processes occur at critical values of the least principal stress σ 3 . I suggest that the concept of the deformation mechanism map of Ashby (1972) be extended to brittle-ductile transitions by normalizing the applied differential stress σ by σ 3 , i.e., the transitions occur at critical values of σ/σ 3 . Third, the high temperature flow law of olivine, widely applied to problems involving flow in the asthenosphere, is inappropriate to the conditions of temperature found in the lithosphere, because different dislocation flow mechanisms dominate at low to intermediate temperatures. The fragmentary rheological data suggest the following rheological structure of the lithosphere where it is inelastically deforming: A rapid pressure-driven increase in rock strength with depth culminates with a shear strength maximum of up to 8 kbar at depths that depend on the state of stress and on the temperature distribution. A review of the mechanisms of weakening associated with water suggests that water weakening effects are probably not important in the oceanic lithosphere but are likely to be controlling in the continental crust.

Journal of Geophysical Research Solid Earth↗

Association between atmospheric circulation patterns and firn-ice core records from the Inilchek glacierized area, central Tien Shan, Asia

Glacioclimatological research in the central Tien Shan was performed in the summers of 1998 and 1999 on the South Inilchek Glacier at 5100-5460 m. A 14.36 m firn-ice core and snow samples were collected and used for stratigraphic, isotopic, and chemical analyses. The firn-ice core and snow records were related to snow pit measurements at an event scale and to meteorological data and synoptic indices of atmospheric circulation at annual and seasonal scales. Linear relationships between the seasonal air temperature and seasonal isotopic composition in accumulated precipitation were established. Changes in the ??18O air temperature relationship, in major ion concentration and in the ratios between chemical species, were used to identify different sources of moisture and investigate changes in atmospheric circulation patterns. Precipitation over the central Tien Shan is characterized by the lowest ionic content among the Tien Shan glaciers and indicates its mainly marine origin. In seasons of minimum precipitation, autumn and winter, water vapor was derived from the and and semiarid regions in central Eurasia and contributed annual maximal solute content to snow accumulation in Tien Shan. The lowest content of major ions was observed in spring and summer layers, which represent maximum seasonal accumulation when moisture originates over the Atlantic Ocean and Mediterranean and Black Seas. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Airborne microorganisms in the African desert dust corridor over the mid-Atlantic ridge, Ocean Drilling Program, Leg 209

The objective of this study was to enhance our understanding of the fate and trans-Atlantic transport of dustborne microorganisms from Northern Africa to the Caribbean and Americas, and more specifically to determine if culturable populations could be detected at a mid-ocean site, closer to the source of dust relative to land-based Caribbean sites, during the early summer months of May and June. Between the dates of 22 May and 30 June 2003, daily air samples were collected and evaluated for the presence of culturable bacterial and fungal colony-forming units (CFU). Here we report a statistically significant correlation between daily atmospheric CFU counts at a mid-ocean research site (???15??N, 45??W) and daily desert dust concentrations as determined by the U.S. Navy's Naval Aerosol Analysis and Prediction System (NAAPS) Global Aerosol Model (Honrath et al. (2004). Journal of Geophysical Research, 109; Johnson et al. (2003). Global Biogeochemical Cycles, 17, 1063; Reid et al. (2004). Geophysical Research Letters, 31; Schollaert, Yoder, Westphal, & O'Reilly (2003). Journal of Geophysical Research, 108, 3191). ?? Springer Science+Business Media B.V. 2006.

Aerobiologia↗

Variations in the reflectivity of the moho transition zone beneath the Midcontinent Rift System of North America: results from true amplitude analysis of GLIMPCE data

True amplitude processing of The Great Lakes International Multidisciplinary Program on Crustal Evolution seismic reflection data from the Midcontinent Rift System of North America shows large differences in the reflectivity of the Moho transition zone beneath the axial rift, beneath the rift flanks, and outside of the rift. The Moho reflection from the axial rift has a discontinuous, diffractive character marginally stronger (several decibels) than an otherwise transparent lower crust and upper mantle. Beneath the axial rift, Moho is interpreted to be a synrift igneous feature. Beneath the rift flanks, the reflectivity of the Moho transition is generally well developed with two identifiable boundaries, although in places it is weakly reflective to nonreflective, similar to Moho outside the rift. The two boundaries are interpreted as the base of essentially intact, although stretched, prerift Archean crust (upper boundary) and new synrift Moho 1-2 s (6-7 km) deeper (lower boundary). Beneath the rift flanks, the layered reflection Moho transition results from the preexisting crustal composition and fabric modified by synrift igneous processes and extensional tectonic/metamorphic processes. The geologic evidence for extensive basaltic magmatism in the rift is the basis for interpreting the Moho signature as a Keweenawan structure that has been preserved for 1.1 b.y. Extension and magmatism appear to enhance reflectivity in the lower crust and Moho transition zone only where stretching factors are moderate (rift flanks) and not where they are extreme (axial rift). This leads to the prediction that the reflectivity across analogous volcanic passive continental margins should be greatest beneath the moderately stretched continental shelves and should decrease towards the ocean-continent boundary.

Journal of Geophysical Research B: Solid Earth↗

Annual modulation of non-volcanic tremor in northern Cascadia

Two catalogs of episodic tremor events in northern Cascadia, one from 2006 to 2012 and the other from 1997 to 2011, reveal two systematic patterns of tremor occurrence in southern Vancouver Island: (1) most individual events tend to occur in the third quarter of the year; (2) the number of events in prolonged episodes (i.e., episodic tremor and slip events), which generally propagate to Vancouver Island from elsewhere along the Cascadia subduction zone, is inversely correlated with the amount of precipitation that occurred in the preceding 2 months. We rationalize these patterns as the product of hydrologic loading of the crust of southern Vancouver Island and the surrounding continental region, superimposed with annual variations from oceanic tidal loading. Loading of the Vancouver Island crust in the winter (when the land surface receives ample precipitation) and unloading in the summer tends to inhibit and enhance downdip shear stress, respectively. Quantitatively, for an annually variable surface load, the predicted stress perturbation depends on mantle viscoelastic rheology. A mechanical model of downdip shear stress on the transition zone beneath Vancouver Island—driven predominantly by the annual hydrologic cycle—is consistent with the 1997–2012 tremor observations, with peak-to-peak downdip shear stress of about 0.4 kPa. This seasonal dependence of tremor occurrence appears to be restricted to southern Vancouver Island because of its unique situation as an elongated narrow-width land mass surrounded by ocean, which permits seasonal perturbations in shear stress at depth.

Journal of Geophysical Research B: Solid Earth↗