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 199 records · Page 11Linked to original sources

On plate tectonics and the geologic evolution of southwestern North America

Very rapid subduction of the Farallon plate under southwestern North America between 60 and 40 Ma was accompanied by a relatively low volume of magmatism throughout the southwestern United States and northern Mexico. Between 40 and 20 Ma, when subduction slowed significantly and in one area may have even stopped, magmatism became widespread and voluminous from Nevada and Utah to central Mexico. This correlation of rapid subduction with a relatively low volume of magmatism can be explained by the observation that subduction-related andesitic arc volcanism, often formed in a Laramide-style compressional regime, is relatively low volume compared to continental volcanism. The shallow roots of arc volcanic systems are clearly exposed in the porphyry copper deposits found in currently active arcs and common throughout southwestern North America between 60 and 50 Ma. By 43 Ma, worldwide plate motions changed, the Pacific plate began moving away from North America, and subduction of the Farallon plate slowed. By around 36 Ma, the easternmost part of the East Pacific Rise, which was located between the Pioneer and Murray fracture zones, approached the trench and the young, hot, buoyant lithosphere appears to have clogged part of the subduction zone. Uplift on land became widespread. Voluminous continental magmatism formed the Sierra Madre Occidental (SMO) of Mexico, one of the largest batholiths in the world, as well as volcanic centers now exposed in the San Juan Mountains of Colorado and the Rio Grande Rift of New Mexico. Vectors of motion of the Pacific plate relative to the North American plate determined by Stock and Molnar (1988) are consistent with formation of a transtensional environment along the plate boundary sufficient to create a 100- to 200-km-wide void just landward of the old volcanic arc. While the SMO batholith was forming within this void, the Monterey and Arguello microplates just offshore to the west were broken off from the Farallon plate and rotated so that the East Pacific Rise in this immediate area became nearly perpendicular to the trench and perpendicular to the vector of motion of the Pacific plate relative to North America. Formation of the SMO batholith was followed between 24 and 20 Ma by a major increase in the rate of subduction of the Guadalupe plate, a fragment of the former Farallon plate, and by increasing mylonitization, extension, and uplift in the metamorphic core complexes that extend northwestward through southern Arizona from the northern end of the SMO batholith. The plate margin underwent another major change between 12.5 and 10 Ma when subduction again stopped, strike slip faulting became dominant along the coast, the Basin and Range Province opened, and numerous tectonostratigraphic terranes in southern California underwent large rotations. By 3 Ma a large, new terrane had been severed from North America immediately west of the SMO batholith as the Gulf of California opened. These observations can be explained by a model for the weakening and ultimate falling apart of the uppermost part of the subducted oceanic plate in the 20–30 m.y. after the end of rapid subduction. As the plate falls apart, not only is compressional stress relieved, but significant backslip along the old subduction zone is also possible, perhaps bringing blueschists rapidly upward from 20- to 30-km depths.

Journal of Geophysical Research↗

Dynamics of marsh-derived sediments in lagoon-type estuaries

Salt marshes are valuable ecosystems that must trap sediments and accrete in order to counteract the deleterious effect of sea‐level rise. Previous studies have shown that the capacity of marshes to build up vertically depends on both autogenous and exogenous processes including eco‐geomorphic feedbacks and sediment supply from in‐land and coastal ocean. There have been numerous efforts to quantify the role played by the sediments coming from marsh edge erosion on the resistance of salt marshes to sea‐level rise. However, the majority of existing studies investigating the interplay between lateral and vertical dynamics use simplified modelling approaches and they do not consider that marsh retreat can affect the regional scale hydrodynamics and sediment retention in back‐barrier basins. In this study, we evaluated the fate of the sediments originating from marsh lateral loss by using high‐resolution numerical model simulations of Jamaica Bay, a small lagoonal estuary located in New York City. Our findings show that up to 42% of the sediments released during marsh edge erosion deposits on the shallow areas of the basin and over the vegetated marsh platforms, contributing positively to the sediment budget of the remaining salt marshes. Furthermore, we demonstrate that with the present‐day sediment supply from the ocean the system cannot keep pace with sea‐level rise even accounting for the sediment liberated in the bay through marsh degradation. Our study highlights the relevance of multiple sediment sources for the maintenance of the marsh complex.

New York↗

Sedimentary deposits in the northern lowland plains, Mars

The lowland plains on Mars have surfaces marked by large polygonal fracture patterns. It was recently proposed that the fracture patterns were developed on sedimentary deposits from outflow channels. We support this hypothesis because of the following observations. (1) Polygonal fracture patterns tend to occur in low areas on Mars that apparently received influx of sediments; the area of northern plains, in which the fractures occur, also forms a regional low. (2) The fracture patterns are concentrated near reentrants of the plains projecting into the southern highlands. Not only are these reentrants associated with outflow channels, but channels and polygonally fractured deposits have similar ages. (3) Sinuous ridges at the mouths of Martian outflow channels resemble ridges in Antarctic ice streams and ice shelves; the similarity suggests that the Martian material was transported from the southern highlands toward the northern plains, as the Antarctic material is transported from the continent toward the adjacent ocean. Overall, it appears that a massive transfer of material took place midway in Martian history, either through the outflow channels or elsewhere along the northern highland scarp, and that this material probably formed the deposits now characterized by polygonal fracture patterns.

Journal of Geophysical Research B: Solid Earth↗

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗

Widespread gas hydrate instability on the upper U.S. Beaufort margin

The most climate-sensitive methane hydrate deposits occur on upper continental slopes at depths close to the minimum pressure and maximum temperature for gas hydrate stability. At these water depths, small perturbations in intermediate ocean water temperatures can lead to gas hydrate dissociation. The Arctic Ocean has experienced more dramatic warming than lower latitudes, but observational data have not been used to study the interplay between upper slope gas hydrates and warming ocean waters. Here we use (a) legacy seismic data that constrain upper slope gas hydrate distributions on the U.S. Beaufort Sea margin, (b) Alaskan North Slope borehole data and offshore thermal gradients determined from gas hydrate stability zone thickness to infer regional heat flow, and (c) 1088 direct measurements to characterize multidecadal intermediate ocean warming in the U.S. Beaufort Sea. Combining these data with a three-dimensional thermal model shows that the observed gas hydrate stability zone is too deep by 100 to 250 m. The disparity can be partially attributed to several processes, but the most important is the reequilibration (thinning) of gas hydrates in response to significant (~0.5°C at 2σ certainty) warming of intermediate ocean temperatures over 39 years in a depth range that brackets the upper slope extent of the gas hydrate stability zone. Even in the absence of additional ocean warming, 0.44 to 2.2 Gt of methane could be released from reequilibrating gas hydrates into the sediments underlying an area of ~5–7.5 × 10 3 km 2 on the U.S. Beaufort Sea upper slope during the next century.

Alaska↗

Hydrated salt minerals on Europa's surface from the Galileo near-infrared mapping spectrometer (NIMS) investigation

We reported evidence of heavily hydrated salt minerals present over large areas of Europa's surface from analysis of reflectance spectra returned by the Galileo mission near infrared mapping spectrometer (NIMS) [ McCord et al. , 1997a, b, 1998a, b]. Here we elaborate on this earlier evidence, present spatial distributions of these minerals, examine alternate water-ice interpretations, expand on our hydrated-salts interpretation, consider salt mineral stability on Europa, and discuss the implications. Extensive well-defined areas on Europa show distinct, asymmetric water-related absorption bands in the 1 to 2.5-μm region. Radiative transfer modeling of water ice involving different particle sizes and layers at Europa temperatures does not reproduce the distinctive Europa water bands. However, ice near its melting temperature, such as in terrestrial environments, does have some characteristics of the Europa spectrum. Alternatively, some classes of heavily hydrated minerals do exhibit such water bands. Among plausible materials, heavily hydrated salt minerals, such as magnesium and sodium sulfates, sodium carbonate and their mixtures, are preferred. All Europa spectral features are present in some salt minerals and a very good match to the Europa spectrum can be achieved by mixing several salt spectra. However, no single or mix of salt mineral spectra from the limited library available has so far been found to perfectly match the Europa spectrum in every detail. The material is concentrated at the lineaments and in chaotic terrain, which are technically disrupted areas on the trailing side. Since the spectrum of the material on Europa is nearly the same everywhere so-far studied, the salt or salt-mixture composition may be nearly uniform. This suggests similar sources and processes over at least a near-hemispheric scale. This would suggest that an extensive subsurface ocean containing dissolved salts is the source, and several possible mechanisms for deposit emplacement are considered. The hydrogen bonds associated with hydration of these salts are similar or greater in strength and energy to those in pure water ice. Thus, once on the surface, the salt minerals should be as stable to disruption as water ice at the Europa temperatures, and mechanisms are suggested to enhance the stability of both materials. Spectra obtained of MgSO 4 ·6H 2 O at 77 K show only small differences from room temperature spectra. The main difference is the appearance of the individual absorptions composing the broad, composite water features and associated with the several different H 2 O sites in the salt hydrate molecule. This suggests that the Europa absorption bands are also composites. Thus higher spectral resolution may reveal these diagnostic features in Europa's spectrum. The specific salts present and their relative abundances would be indicators of the chemistry and conditions of an ocean environment, and areas of fresh, heavy concentration of these minerals should make ideal lander mission sampling sites.

Journal of Geophysical Research E: Planets↗

Crustal structure of a transform plate boundary: San Francisco Bay and the central California continental margin

Wide-angle seismic data collected during the Bay Area Seismic Imaging Experiment provide new glimpses of the deep structure of the San Francisco Bay Area Block and across the offshore continental margin. San Francisco Bay is underlain by a veneer (<300 m) of sediments, beneath which P wave velocities increase rapidly from 5.2 km/s to 6.0 km/s at 7 km depth, consistent with rocks of the Franciscan subduction assemblage. The base of the Franciscan at-15-18 km depth is marked by a strong wide-angle reflector, beneath which lies an 8- to 10-km-thick lower crust with an average velocity of 6.75??0.15 km/s. The lower crust of the Bay Area Block may be oceanic in origin, but its structure and reflectivity indicate that it has been modified by shearing and/or magmatic intrusion. Wide-angle reflections define two layers within the lower crust, with velocities of 6.4-6.6 km/s and 6.9-7.3 km/s. Prominent subhorizontal reflectivity observed at near-vertical incidence resides principally in the lowermost layer, the top of which corresponds to the "6-s reflector" of Brocher et al. [1994]. Rheological modeling suggests that the lower crust beneath the 6-s reflector is the weakest part of the lithosphere; the horizontal shear zone suggested by Furlong et al. [1989] to link the San Andreas and Hayward/Calaveras fault systems may actually be a broad zone of shear deformation occupying the lowermost crust. A transect across the continental margin from the paleotrench to the Hayward fault shows a deep crustal structure that is more complex than previously realized. Strong lateral variability in seismic velocity and wide-angle reflectivity suggests that crustal composition changes across major transcurrent fault systems. Pacific oceanic crust extends 40-50 km landward of the paleotrench but, contrary to prior models, probably does not continue beneath the Salinian Block, a Cretaceous arc complex that lies west of the San Andreas fault in the Bay Area. The thickness (10 km) and high lower-crustal velocity of Pacific oceanic crust suggest that it was underplated by magmatism associated with the nearby Pioneer seamount. The Salinian Block consists of a 15-km-thick layer of velocity 6.0-6.2 km/s overlying a 5-km-thick, high-velocity (7.0 km/s) lower crust that may be oceanic crust, Cretaceous arc-derived lower crust, or a magmatically underplated layer. The strong structural variability across the margin attests to the activity of strike-slip faulting prior to and during development of the transcurrent Pacific/North American plate boundary around 29 Ma. Copyright 1996 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Time-averaged near-bed suspended sediment concentrations under waves and currents: Comparison of measured and model estimates

Profiles of suspended sediment concentration and velocity were measured over a 15-day period at a near-shore site off Santa Cruz, CA in Monterey Bay. The concentration and velocity data were collected with an Acoustic Backscattering System (ABS) and Acoustic Current Profiler (ACP) that were mounted on a bottom tripod. High-resolution bottom scanning sonar was also attached to the tripod to provide images of bed features during the experiment. Hourly time-averaged near-bed concentrations of suspended sediment were calculated from three models and compared with the measurements. Surface waves and currents that were generated by a storm of moderate intensity caused bed stresses that exceeded threshold stress for D 50 =0.02 cm, the median size of the moderately well-sorted bottom sediment, over a period of about 7 days. Estimates of the concentration at 1 cm above the bottom, C a1 , were obtained using the ABS measurements. These observations have been compared with predictions for the concentration at 1 cm above the bottom, C 1 . Nielsen's models for reference concentration C o [Nielsen, P., 1986. Suspended sediment concentrations under waves. Coastal Engineering 10, 32-31; Nielsen, P., 1992. Coastal Bottom Boundary Layers and Sediment Transport, Advanced Series on Ocean Engineering. World Scientific, Hackensack, NJ.] are purely wave-based and do not include effects of bottom currents on bed stress and bedform scales. C 1 calculated from this model compared well with measured C a1 when currents were weak and small oscillatory ripples were observed in the sonar images. However, during the 3-day period of highest bottom stresses modeled C 1 did not compare well to C a1 . The other two models for C 1 , Glenn and Grant [Glenn, S.M., Grant, W.D., 1987. A suspended sediment stratification correction for combined wave and current flows. Journal of Geophysical Research 92(C8), 8244-8264.] and van Rijn and Walstra [Van Rijn, L.C., Walstra, D.J.R., 2004. Description of TRANSPOR2004 and implementation in Delft3D-ONLINE. Interim Report prepared for DG Rijkswaterstaat, Rijksinstituut voor Kust en Zee. Delft Hydraulics Institute, The Netherlands.], accounted for combined wave-current stresses and included different formulations for predicting bedform scales. C 1 calculated from these models compared well with measurements throughout the early and middle period of the storm, but underpredicted measured values during the latter part of the experiment. An alternative method based on dimensional analysis provided a straightforward way to relate C 1 to bed stress and sediment parameters. This new relationship based on dimensional analysis explained about 77% of the variability in C a1 derived from the ABS data.

California↗

Separating sea and slow slip signals on the seafloor

Seafloor pressure measurements hold promise for estimating vertical displacements from transient slow slip events on submarine faults. We assess the accuracy of pressure offset estimates that evolve over days to weeks and the confidence with which they may be attributed to tectonic deformation or to the ocean water column. One common approach to resolve this ambiguity assumes water column pressures vary insignificantly over the study region and are represented by stable reference site pressures. Assessing the validity of this assumption requires independent evidence. Correlations between pressures and colocated temperatures collected during the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment suggest temperatures might provide a useful independent proxy for water column pressures. We compared offsets estimated using several methods, with temperature and other proxies. The use of a temperature proxy was unsuccessful, because seafloor temperatures did not track the seasonal signal that contributes significantly to seafloor pressure changes over the slow slip event period. Regardless of the estimation method, offsets varied within a few cm around some uncertain reference level. Commonly used statistical measures are shown not to be reliable indicators of offset accuracy since offsets contribute minimally to the total variance. Offsets estimated using identical methods but with seafloor pressures simulated using a regional ocean model were larger than those derived from the data but had a similar pattern. Since the model simulates only water column processes, this suggests a significant fraction of the estimated pressure offsets are due to seasonal water column signal and are not of tectonic origin.

Journal of Geophysical Research- Solid Earth↗

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↗

Mg-spinel lithology: A new rock type on the lunar farside

High-resolution compositional data from Moon Mineralogy Mapper (M 3) for the Moscoviense region on the lunar farside reveal three unusual, but distinctive, rock types along the inner basin ring. These are designated "OOS" since they are dominated by high concentrations of orthopyroxene, olivine, and Mg-rich spinel, respectively. The OOS occur as small areas, each a few kilometers in size, that are widely separated within the highly feldspathic setting of the basin rim. Although the abundance of plagioclase is not well constrained within the OOS, the mafic mineral content is exceptionally high, and two of the rock types could approach pyroxenite and harzburgite in composition. The third is a new rock type identified on the Moon that is dominated by Mg-rich spinel with no other mafic minerals detectable (<5% pyroxene, olivine). All OOS surfaces are old and undisturbed since basin formation. They are effectively invisible in image data and are only recognized by their distinctive composition identified spectroscopically. The origin of these unusual lithologies appears to be linked to one or more magmatic intrusions into the lower crust, perhaps near the crust-mantle interface. Processes such as fractional crystallization and gravity settling within such intrusions may provide a mechanism for concentrating the mafic components within zones several kilometers in dimension. The OOS are embedded within highly anorthositic material from the lunar crust; they may thus be near contemporaneous with crustal products from the cooling magma ocean. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Late Cretaceous paleomagnetism of the Tucson Mountains: Implications for vertical axis rotations in south central Arizona

The Tucson Mountains of southern Arizona are the site of an Upper Cretaceous caldera from which the rhyolitic Cat Mountain Tuff was erupted at about 72 Ma. Two magnetic units within the Cat Mountain Tuff are distinguished by paleomagnetic data in both the northern and southern Tucson Mountains. The resurgent Amole pluton (≃72 Ma) in the northern Tucson Mountains was emplaced soon after eruption of the Cat Mountain Tuff but cooled and was magnetized after northeastward tilting (50°–85°) of the adjacent caldera-fill sequence. Petrologic and paleomagnetic data indicate that the lower magnetic unit of the Cat Mountain Tuff caps the Silver Bell Mountains to the northwest. A previous paleomagnetic investigation ( N = 34) indicates that the Silver Bell Mountains have been rotated clockwise 30°±16° (95% confidence level) about a vertical axis relative to cratonic North America. A similar paleomagnetic study of Upper Cretaceous volcanic, volcaniclastic, and intrusive units in the Tucson Mountains ( N = 26) indicates that these rocks have been rotated 7°±14° clockwise relative to stable North America. A direct comparison of paleomagnetic directions for the lower unit of the Cat Mountain Tuff shows a 17°±10° clockwise rotation between the Silver Bell Mountains and the Tucson Mountains which supports the relative accuracy of the absolute rotations determined for these two mountain ranges. Preliminary paleomagnetic directions for middle Tertiary units from the Silver Bell and Tucson Mountains are consistent with clockwise rotation having occurred prior to deposition of these rocks. Clockwise rotation of crustal blocks in southern Arizona likely was associated with strike-slip movement on major northwest trending faults in the region, and this movement may have been related to oblique subduction of oceanic plates along the western continental margin in Late Cretaceous and early Tertiary time. The available paleomagnetic data indicate that rocks in southern Arizona have not remained unrotated with respect to North America since Late Cretaceous time and that vertical axis rotations may have played an important role in the region during Laramide deformation.

Journal of Geophysical Research↗

Fast rates of subduction erosion along the Costa Rica Pacific margin: Implications for nonsteady rates of crustal recycling at subduction zones

At least since the middle Miocene (???16 Ma), subduction erosion has been the dominant process controlling the tectonic evolution of the Pacific margin of Costa Rica. Ocean Drilling Program Site 1042 recovered 16.5 Ma nearshore sediment at ???3.9 km depth, ???7 km landward of the trench axis. The overlying Miocene to Quaternary sediment contains benthic foraminifera documenting margin subsidence from upper bathyal (???200 m) to abyssal (???2000 m) depth. The rate of subsidence was low during the early to middle Miocene but increased sharply in the late Miocene-early Pliocene (5-6.5 Ma) and at the Pliocene-Pleistocene boundary (2.4 Ma). Foraminifera data, bedding dip, and the geometry of slope sediment indicate that tilting of the forearc occurred coincident with the onset of rapid late Miocene subsidence. Seismic images show that normal faulting is widespread across the continental slope; however, extension by faulting only accounts for a minor amount of the post-6.5 Ma subsidence. Basal tectonic erosion is invoked to explain the subsidence. The short-term rate of removal of rock from the forearc is about 107-123 km3 Myr-1 km-1. Mass removal is a nonsteady state process affecting the chemical balance of the arc: the ocean sediment input, with the short-term erosion rate, is a factor of 10 smaller than the eroded mass input. The low 10Be concentration in the volcanic arc of Costa Rica could be explained by dilution with eroded material. The late Miocene onset of rapid subsidence is coeval with the arrival of the Cocos Ridge at the subduction zone. The underthrusting of thick and thermally younger ocean crust decreased the subduction angle of the slab along a large segment of the margin and changed the dynamic equilibrium of the margin taper. This process may have induced the increase in the rate of subduction erosion and thus the recycling of crustal material to the mantle. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Heat flow in the Western Arctic Ocean (Amerasian Basin)

From 1963 to 1973 the U.S. Geological Survey (USGS) measured heat flow at 356 sites in the Amerasian Basin (Western Arctic Ocean) from a drifting ice island (T-3). The resulting measurements, which are unevenly distributed on Alpha-Mendeleev Ridge (AMR) and in Canada and Nautilus basins, greatly expand available heat flow data for the Arctic Ocean. Average T-3 heat flow is ~54.7 ± 11.3 mW m-2, and Nautilus Basin, including Mendeleev Plain, is the only well-surveyed area (~13% of data) with significantly higher average heat flow (63.8 mW m-2). Heat flow and bathymetry are not correlated at a large scale, and turbiditic surficial sediments (Canada and Nautilus basins) have higher heat flow than the sediments that blanket the AMR. Thermal gradients are mostly near-linear, implying that conductive heat transport dominates and that near-seafloor sediments are in thermal equilibrium with overlying bottom waters. Combining the heat flow data with modern seismic imagery suggests that some of the observed heat flow variability may be explained by local changes in sediment thickness or lithology or the presence of basement faults that channel circulating seawater. A thermal model that incorporates thermal conductivity variations along a profile from Canada Basin (thick sediment on mostly oceanic crust) to Alpha Ridge (thin sediment over thick magmatic units associated with the High Arctic Large Igneous Province) predicts heat flow lower than that observed on Alpha Ridge. This, along with other observations, implies that circulating fluids modulate conductive heat flow and contribute to high variability in the T-3 dataset. .

Journal of Geophysical Research B: Solid Earth↗

Nonlinear refraction and reflection travel time tomography

We develop a rapid nonlinear travel time tomography method that simultaneously inverts refraction and reflection travel times on a regular velocity grid. For travel time and ray path calculations, we apply a wave front method employing graph theory. The first-arrival refraction travel times are calculated on the basis of cell velocities, and the later refraction and reflection travel times are computed using both cell velocities and given interfaces. We solve a regularized nonlinear inverse problem. A Laplacian operator is applied to regularize the model parameters (cell slownesses and reflector geometry) so that the inverse problem is valid for a continuum. The travel times are also regularized such that we invert travel time curves rather than travel time points. A conjugate gradient method is applied to minimize the nonlinear objective function. After obtaining a solution, we perform nonlinear Monte Carlo inversions for uncertainty analysis and compute the posterior model covariance. In numerical experiments, we demonstrate that combining the first arrival refraction travel times with later reflection travel times can better reconstruct the velocity field as well as the reflector geometry. This combination is particularly important for modeling crustal structures where large velocity variations occur in the upper crust. We apply this approach to model the crustal structure of the California Borderland using ocean bottom seismometer and land data collected during the Los Angeles Region Seismic Experiment along two marine survey lines. Details of our image include a high-velocity zone under the Catalina Ridge, but a smooth gradient zone between. Catalina Ridge and San Clemente Ridge. The Moho depth is about 22 km with lateral variations. Copyright 1998 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

The North American upper mantle: Density, composition, and evolution

The upper mantle of North America has been well studied using various seismic methods. Here we investigate the density structure of the North American (NA) upper mantle based on the integrative use of the gravity field and seismic data. The basis of our study is the removal of the gravitational effect of the crust to determine the mantle gravity anomalies. The effect of the crust is removed in three steps by subtracting the gravitational contributions of (1) topography and bathymetry, (2) low-density sedimentary accumulations, and (3) the three-dimensional density structure of the crystalline crust as determined by seismic observations. Information regarding sedimentary accumulations, including thickness and density, are taken from published maps and summaries of borehole measurements of densities; the seismic structure of the crust is based on a recent compilation, with layer densities estimated from P-wave velocities. The resultant mantle gravity anomaly map shows a pronounced negative anomaly (−50 to −400 mGal) beneath western North America and the adjacent oceanic region and positive anomalies (+50 to +350 mGal) east of the NA Cordillera. This pattern reflects the well-known division of North America into the stable eastern region and the tectonically active western region. The close correlation of large-scale features of the mantle anomaly map with those of the topographic map indicates that a significant amount of the topographic uplift in western NA is due to buoyancy in the hot upper mantle, a conclusion supported by previous investigations. To separate the contributions of mantle temperature anomalies from mantle compositional anomalies, we apply an additional correction to the mantle anomaly map for the thermal structure of the uppermost mantle. The thermal model is based on the conversion of seismic shear-wave velocities to temperature and is consistent with mantle temperatures that are independently estimated from heat flow and heat production data. The thermally corrected mantle density map reveals density anomalies that are chiefly due to compositional variations. These compositional density anomalies cause gravitational anomalies that reach ~250 mGal. A pronounced negative anomaly (−50 to −200 mGal) is found over the Canadian shield, which is consistent with chemical depletion and a corresponding low density of the lithospheric mantle, also referred to as the mantle tectosphere. The strongest positive anomaly is coincident with the Gulf of Mexico and indicates a positive density anomaly in the upper mantle, possibly an eclogite layer that has caused subsidence in the Gulf. Two linear positive anomalies are also seen south of 40°N: one with a NE-SW trend in the eastern United States, roughly coincident with the Grenville-Appalachians, and a second with a NW-SE trend beneath the states of Texas, New Mexico, and Colorado. These anomalies are interpreted as being due to (1) the presence of remnants of an oceanic slab in the upper mantle beneath the Grenville-Appalachian suture and (2) mantle thickening caused by a period of shallow, flat subduction during the Laramie orogeny, respectively. Based on these geophysical results, the evolution of the NA upper mantle is depicted in a series of maps and cartoons that display the primary processes that have formed and modified the NA crust and lithospheric upper mantle.

Journal of Geophysical Research B: Solid Earth↗

Seismic structure of the southern Cascadia subduction zone and accretionary prism north of the Mendocino triple junction

Four multichannel-seismic reflection profiles, collected as part of the Mendocino triple junction seismic experiment, image the toe of the southern Cascadia accretionary prism. Today, 250–600 m of sediment is subducting with the Gorda plate, and 1500–3200 m is accreting to the northern California margin. Faults imaged west and east of the deformation front show mixed structural vergence. A north-south trending, 20 km long portion of the central margin is landward vergent for the outer 6–8 km of the toe of the prism. This region of landward vergence exhibits no frontal thrust, is unusually steep and narrow, and is likely caused by a seaward-dipping backstop close to the deformation front. The lack of margin-wide preferred seaward vergence and wedge-taper analysis suggests the prism has low basal shear stress. The three southern lines image wedge-shaped fragments of oceanic crust 1.1–7.3 km in width and 250–700 m thick near the deformation front. These wedges suggest shortening and thickening of the upper oceanic crust. Discontinuities in the seafloor west of the prism provide evidence for mass wasting in the form of slump blocks and debris fans. The southernmost profile extends 75 km west of the prism imaging numerous faults that offset both the Gorda basin oceanic crust and overlying sediments. These high-angle faults, bounding basement highs, are interpreted as strike-slip faults reactivating structures originally formed at the spreading ridge. Northeast or northwest trending strike-slip faults within the basin are consistent with published focal mechanism solutions and are likely caused by north-south Gorda-Pacific plate convergence.

Journal of Geophysical Research B: Solid Earth↗

Geologic implications of topographic, gravity, and aeromagnetic data in the northern Yukon-Koyukuk province and its borderlands, Alaska

The northern Yukon-Koyukuk province is characterized by low elevation and high Bouguer gravity and aeromagnetic anomalies in contrast to the adjacent Brooks Range and Ruby geanticline. Using newly compiled digital topographic, gravity, and aeromagnetic maps, I have divided the province into three geophysical domains. The Koyukuk domain, which is nearly equivalent to the Koyukuk lithotectonic terrane, is a horseshoe-shaped area, open to the south, of low topography, high gravity, and high-amplitude magnetic anomalies caused by an intraoceanic magmatic arc. The Angayucham and Kanuti domains are geophysical subdivisions of the Angayucham lithotectonic terrane that occur along the northern and southeastern margins of the Yukon-Koyukuk province, where oceanic rocks have been thrust over continental rocks of the Brooks Range and Ruby geanticline. Basalt of the Angayucham domain causes strong gravity highs and weak magnetic highs. The Kanuti domain is distinguished from the Angayucham domain by intense magnetic highs caused by cumulus mafic and ultramafic plutonic rocks, abundant ultramafic mantle tectonites, and magnetic syenite and monzonite. Long-wavelength, low-intensity magnetic highs and undulating gravity anomalies indicate an undulating basement surface of varied lithology beneath the Kobuk-Koyukuk and Lower Yukon basins. Modeling of gravity and magnetic anomalies shows that oceanic rocks of the Angayucham and Kanuti domains dip inward beneath the Kobuk-Koyukuk basin. The modeling supports, but does not prove, the hypothesis that the crust of the Kobuk-Koyukuk basin is 32–35 km thick, consisting of a tectonically thickened section of Cretaceous volcanic and sedimentary rocks and older oceanic crust. Plutons of the Brooks Range and the southern Ruby geanticline are nonmagnetic, ilmenite series, S-type granites that cause magnetic lows. Plutons of the northern Ruby geanticline are variable in their magnetic properties and cause both highs and lows. Plutons of both the eastern and western Yukon-Koyukuk province are variable in their magnetic expression but commonly cause magnetic lows in contrast to andesite.

Journal of Geophysical Research Solid Earth↗