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 289 records · Page 16Linked to original sources
USGS2012

Lead and strontium isotopes and related trace elements as genetic tracers in the Upper Cenozoic rhyolite-basalt association of the Yellowstone Plateau volcanic field

Supported by various field geologic and petrologic data, the contents of Pb, U, Th, Rb, and Sr and the isotopic compositions of Pb and Sr for upper Cenozoic volcanic rocks of the Yellowstone Plateau volcanic field are consistent with the hypothesis of derivation of the basaltic and rhyolitic magmas by partial melting of distinct source regions in the upper mantle and lower crust, respectively. All the basalt samples analyzed but one have systematically lower values of 207 Pb/ 204 Pb and 87 Sr/ 86 Sr than the rhyolites. The values of 206 Pb/ 204 Pb are smaller, and 87 Sr/ 86 Sr are mostly larger than known values in oceanic basalts. In all but one case, the values of 207 Pb/ 204 Pb are higher than expected from an extrapolation of known values in oceanic basalts to less radiogenic values of 206 Pb/ 204 Pb. Because there are no xenoliths, phenocrysts are only moderate to sparse in abundance, REE patterns are low and flat at the radiogenic end of lead isotopic compositions, several values of Rb/Sr are low, and 80% of the basalt samples form a well-developed secondary isochron separate from the rhyolites, we favor an interpretation for basalt genesis wherein isotopic signatures of most mafic magmas were attained in a continental ‘keel’ of mantlelike character about 2.6 b.y. old or somewhat older attached to the crust, and these signatures were unaltered by magma passage through the crust. At the very least, the current data continue to cast serious doubt as to the inevitability of crustal contamination for basaltic magma intruding the continental environment and postulate that much can be learned about the mantle under continents through the study of continental basalts. One basalt unit with an unusually low value of 207 Pb/ 204 Pb and an 87 Ar/ 86 Ar less than 0.704 may represent subcontinental ‘keel’-derived magma that rose unaltered to the surface. Our data also are not consistent with formation of this rhyolite-basalt association primarily by such processes as crystal fractionation, separation of immiscible silicate liquids from a common parental magma, or fractional melting of a homogeneous source. Rather as a conceptual model, we envision large mafic intrusions to have been injected into the lower crust resulting in rhyolite generation through partial anatexis of the adjacent wall rocks which probably had a 206 Pb/ 204 Pb < 17 and 87 Sr/ 86 Sr > 0.709; a model that has much in common with that proposed by Holmes (1931). All the other hypotheses listed have the necessary added complication that either the basalt or the rhyolite or both become contaminated after the two magma types separated, have problems accounting for the lack of igneous rocks of intermediate compositions or production of such large volumes of rhyolitic material (∼5000 km 3 ), and fail to explain why rhyolitic magma is not a more common occurrence in the ocean basin. We appeal to bouyancy of rhyolites to generate a barrier for basalt magma migration and account for the great preponderance of rhyolite relative to basalt at the surface. Furthermore, the complex isotopic picture in the rhyolites indicates that many of these magmas interacted with the upper crustal geologic units that they traversed. The interactions involved diverse processes, probably including reacton with hydrothermal fluids or hydrothermally altered rocks at high levels as well as by contamination with Phanerozoic sedimentary and Precambrian crystalline rocks at deeper levels. At the very least, we feel our study adds a cautionary note to the currently increasingly popular hypothesis that differentiation of basalt or gabbro magmas to rhyolite or granite (as distinct from tonalite or dacite) is a common occurrence and is therefore an important continential building process. Models for formation of rhyolite and granite predominantly by reworking of crust (anatexis) must still be considered. The primitive Archean mantle of the region was characterized by higher Rb/Sr, U/Pb, and Th/U values than are typical of modern suboceanic mantle. The mantle residuum within the continental subcrustal lithosperic ‘keel’ that resulted from the Archean crustal differentiation event probably was depleted in Rb/Sr and U/Pb, and the crust was correspondingly enriched in these ratios. The crust probably was further differentiated by an Archean high-grade metamorphism, during or after the primary event, into a granulitic lower crust depleted in U/Pb and Rb/Sr and a lower-grade upper crust enriched in these ratios.

Journal of Geophysical Research Solid Earth↗
USGS2021

The seismo-acoustics of submarine volcanic eruptions

Many of the world’s volcanoes are hidden beneath the ocean’s surface where eruptions are difficult to observe. However, seismo‐acoustic signals produced by these eruptions provide a useful means of identifying active submarine volcanism. A literature survey revealed reports of 119 seismo‐acoustically recorded submarine eruptions since 1939. Submarine eruptions have been recorded in all major tectonic settings, with a range of geochemistries, and at a variety of water depths, but the reports are dominated by eruptions in the Pacific and at only a few locations. Many of the reports offer little detail, with over half of the observations made from distances >500 km, and only about half were confirmed as eruptions by non‐seismo‐acoustic evidence. The reported seismo‐acoustic signals cover a wide variety of processes, including earthquakes, explosions, various types of tremor, signals related to lava extrusion, and landslides. Recorded signals can sometimes be difficult to classify or confidently associate with an eruption, although there has been progress in this regard. Real‐time monitoring of submarine eruptions has been on‐going for several decades on regional and global scales with growing interest and effort in local networks. Real‐time networks are complemented by short‐term instrument deployments that often give more detailed insights into the dynamics and processes of submarine eruptions. Thorough seismo‐acoustic monitoring and study has increased our understanding of submarine eruptions, especially of deep‐sea volcanoes and spreading centers. Despite this, there are still many outstanding questions that need to be addressed for submarine volcanoes to be as well understood and monitored as their terrestrial counterparts.

Journal of Geophysical Research↗
USGS2024

GRACE and GRACE Follow-On gravity observations of intermediate-depth earthquakes contrasted with those of shallow events

Earthquakes involve mass redistribution within the solid Earth and the ocean, and as a result, perturb the Earth's gravitational field. For most of the shallow (<60 km) earthquakes with M w > 8.0, the GRACE satellite gravity measurements suggest considerable volumetric disturbance of rocks. At a spatial scale of hundreds of km, the effect of volumetric change exceeds gravity change by vertical deformation; for example, negative gravity anomalies associated with volumetric expansion are characteristic patterns after shallow thrust events. In this study, however, we report contrasting observations of gravity change from two intermediate-depth (100–150 km) earthquakes of 2016 & 2017 Mw 8.0 (two combined) Papua New Guinea thrust faulting events and 2019 Mw 8.0 Peru normal faulting and highlight the importance of compressibility in earthquake deformation. The combined 2016/17 thrust events resulted in a positive gravity anomaly of 5–6 microGal around the epicenter, while the 2019 normal faulting produced a negative gravity anomaly of 3–4 microGal. Our modeling found that these gravity changes are manifestation of vertical deformation with limited volumetric change, distinct from gravity changes after the shallow earthquakes. The stronger resistance of rocks to volume change at intermediate-depth results in largely incompressible deformation and thus in a gravity change dominated by vertical deformation. In addition, malleable rocks under high pressure and temperature at depth facilitated substantial afterslip and/or fast viscoelastic relaxation causing additional vertical deformation and gravity change equivalent to the coseismic change. For the Papua New Guinea events, this means that postseismic relaxation enhanced coseismic uplift and relative sea level decrease.

Journal of Geophysical Research: Solid Earth↗
USGS2003

Fault trends on the seaward slope of the Aleutian Trench: Implications for a laterally changing stress field tied to a westward increase in oblique convergence

Normal faults along the seaward trench slope (STS) commonly strike parallel to the trench in response to bending of the oceanic plate into the subduction zone. This is not the circumstance for the Aleutian Trench, where the direction of convergence gradually changes westward, from normal to transform motion. GLORIA side-scan sonar images document that the Aleutian STS is dominated by faults striking oblique to the trench, west of 179??E and east of 172??W. These images also show a pattern of east-west trending seafloor faults that are aligned parallel to the spreading fabric defined by magnetic anomalies. The stress-strain field along the STS is divided into two domains west and east, respectively, of 179??E. Over the western domain, STS faults and nodal planes of earthquakes are oriented oblique (9??-46??) to the trench axis and (69??-90??) to the magnetic fabric. West of 179??E, STS fault strikes change by 36?? from the E-W trend of STS where the trench-parallel slip gets larger than its orthogonal component of convergence. This rotation indicates that horizontal stresses along the western domain of the STS are deflected by the increasing obliquity in convergence. An analytical model supports the idea that strikes of STS faults result from a superposition of stresses associated with the dextral shear couple of the oblique convergence and stresses caused by plate bending. For the eastern domain, most nodal planes of earthquakes strike parallel to the outer rise, indicating bending as the prevailing mechanism causing normal faulting. East of 172??W, STS faults strike parallel to the magnetic fabric but oblique (10??-26??) to the axis of the trench. On the basis of a Coulomb failure criterion the trench-oblique strikes probably result from reactivation of crustal faults generated by spreading. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Propagation and linkage of oceanic ridge segments

We have investigated the propagation of spreading ridges and the development of structures that link ridge segments using an analogy between ridges and cracks in elastic plates. The ridge-propagation force and a path factor that controls propagation direction were calculated for echelon ridge segments propagating toward each other. The ridge-propagation force increases as ridge ends approach but then declines sharply as the ends pass, so ridge segments may overlap somewhat. The sign of the path factor changes as ridge ends approach and pass, so the overlapping ridge ends may diverge and then converge following a hook-shaped path. The magnitudes of shear stresses in the plane of the plate and orientations of maximum shear planes between adjacent ridge segments were calculated to study transform faulting. For different loading conditions simulating ridge push, plate pull, and ridge suction, we identify a zone of intense mechanical interaction between adjacent ridge ends in which stresses are concentrated. For all conditions, the shear stress in the interaction zone increases as ends approach and remains large as the ends overlap; thus crust in this zone may fracture and weaken in preparation for the formation of a through-going transform fault. The calculated shear planes rotate toward an orientation about 90° from the strike of ridges as the ends pass, thus favoring the orthogonal arrangement of ridges and transforms. The magnitudes of mean stresses in the plane of the plate and orientations of principal stress planes were also calculated. The mean stress is tensile in the interaction zone, so basins may form there, except in the case of ridge push loading. The planes across which the maximum tension acts are oblique to ridges, thus favoring obliquely oriented normal faults bounding the transform valley.

Journal of Geophysical Research Solid Earth↗
USGS2001

Asian dust events of April 1998

On April 15 and 19, 1998, two intense dust storms were generated over the Gobi desert by springtime low-pressure systems descending from the northwest. The windblown dust was detected and its evolution followed by its yellow color on SeaWiFS satellite images, routine surface-based monitoring, and through serendipitous observations. The April 15 dust cloud was recirculating, and it was removed by a precipitating weather system over east Asia. The April 19 dust cloud crossed the Pacific Ocean in 5 days, subsided to the surface along the mountain ranges between British Columbia and California, and impacted severely the optical and the concentration environments of the region. In east Asia the dust clouds increased the albedo over the cloudless ocean and land by up to 10–20%, but it reduced the near-UV cloud reflectance, causing a yellow coloration of all surfaces. The yellow colored backscattering by the dust eludes a plausible explanation using simple Mie theory with constant refractive index. Over the West Coast the dust layer has increased the spectrally uniform optical depth to about 0.4, reduced the direct solar radiation by 30–40%, doubled the diffuse radiation, and caused a whitish discoloration of the blue sky. On April 29 the average excess surface-level dust aerosol concentration over the valleys of the West Coast was about 20–50 μg/m 3 with local peaks >100 μg/m 3 . The dust mass mean diameter was 2–3 μm, and the dust chemical fingerprints were evident throughout the West Coast and extended to Minnesota. The April 1998 dust event has impacted the surface aerosol concentration 2–4 times more than any other dust event since 1988. The dust events were observed and interpreted by an ad hoc international web-based virtual community. It would be useful to set up a community-supported web-based infrastructure to monitor the global aerosol pattern for such extreme aerosol events, to alert and to inform the interested communities, and to facilitate collaborative analysis for improved air quality and disaster management.

Journal of Geophysical Research D: Atmospheres↗
USGS2012

Northern East Pacific Rise: Magnetic anomaly and bathymetric framework

The oceanic crust in the eastern Pacific between 7°N and 30°N and east of 127°W contains a fairly complete history of the spreading centers associated with the East Pacific Rise since 25 m.y. B.P. (late Oligocene). In this paper, we have summarized the seafloor spreading magnetic-anomaly data and the bathymetric data that reflect the record of this tectonic history. The well-defined magnetic lineations north of the Clarion fracture zone, in the mouth of the Gulf of California, and on the east flank of the East Pacific Rise (EPR) are carefully examined and used to provide a guide for interpreting the spreading pattern between the Clarion and Clipperton fracture zones, southward of the Rivera fracture zone over the Mathematician Ridge, and over the entire EPR east of the Mathematician Ridge between the Rivera and Siqueiros fracture zones. The bathymetric data provide a trace of the fracture zone pattern in each of the above mentioned areas. The fracture zone bathymetry and the seafloor spreading magnetic lineations on the EPR south of the Rivera fracture zone have a distinctive fanning pattern caused by close poles of rotation and plate boundary reorganizations. All these data provide a good record of the plate reorganizations in the middle Miocene at magnetic anomaly 5 A time (12.5 to 11 m.y. B.P.), in the late Miocene at magnetic anomaly 3′−4 time (6.5 m.y. B.P.), and in the Pliocene at magnetic anomaly 2′−3 time (3.5 m.y. B.P.). Several abandoned spreading centers, including the Mathematician Ridge, were left behind as a result of these reorganizations. The Mathematician Ridge is shown to be a set of ridges and trough whose origin is related to the tectonic activity associated with each of the above mentioned reorganizations since anomaly 5A.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Origin and depositional environment of clastic deposits in the Hilo drill hole, Hawaii

Volcaniclastic units cored at depths of about 87, 164, 178, 226, and 246 m below sea level and carbonate units located between depths of 27 and 53 m below sea level in the Hilo drill core were found to be deposited at or near sea level. Four of these units are hydroclastic deposits, formed when subaerially erupted Mauna Loa lava flows entered the ocean and fragmented to produce quenched, glassy fragments during hydrovolcanic explosions. Ash units 24 and 26, at 178 m depth, accumulated at sea level in a freshwater bog. They contain pyroxenes crystallized from tholeiitic magma that we infer erupted explosively at the summit of Kilauea volcano. Two carbon‐rich layers from these ashes have a weighted average radiocarbon age of 38.6 ± 0.9 ka; the ashes probably correlate with the oldest and thickest part of the Pahala ash. Ash unit 44, at the transition from Mauna Kea to Mauna Loa lava flows, was probably nearly 3.2 m thick and is inferred to be equivalent to the lower thick part of the composite Homelani ash mapped in Hilo and on the flanks of Mauna Kea. The age of this part of Homelani ash is between 128 ± 33 and 200 ± 10 ka; it may have erupted subglacially during the Pohakuloa glacial maxima on Mauna Kea. Beach sand units 12 and 22 were derived from nearby Mauna Loa and Mauna Kea lava flows. The middle of beach sand unit 38 was derived mainly from lava erupted near the distal end of the subaerial east rift zone of Kilauea volcano; these sands were transported about 33 km northwest to Hilo Bay by prevailing longshore currents. Combined age, depth, and sea level markers in the core allow us to determine that lava flow recurrence intervals averaged one flow every 4 kyr during the past 86 kyr and one flow every 16 kyr between 86 and 200 ka at the drill site and that major explosive eruptions that deposit thick ash in Hilo have occurred only twice in the last 400 kyr. These recurrence intervals support the moderate lava flow hazard zonation (zone 3) for coastal Hilo previously determined from surficial mapping.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Structural evidence for northeastward movement on the Chocolate Mountains Thrust, southeasternmost California

The Late Cretaceous Chocolate Mountains thrust of southeastern California and southwestern Arizona places a block of Proterozoic and Mesozoic continental crust over the late Mesozoic continental margin oceanic sedimentary and volcanic rocks of the regionally distinctive Orocopia Schist. The Chocolate Mountains thrust is interpreted as a thrust (burial, subduction) fault rather than a low-angle normal (exhumation, unroofing, uplift) fault. An important parameter required to understand the tectonic significance of the Chocolate Mountains and related thrusts is their sense of movement. The Chocolate Mountains thrust zone contains sparse to locally abundant mesoscopic asymmetric folds. Fabric relations, supported by regional geologic evidence, indicate that these folds are an integral part of and coeval with the thrust zone. On a lower hemisphere equal-area plot representing the orientation and sense of asymmetry of 80 thrust zone folds from 36 localities, spread over an area 60 by 10 km, Z folds plot northwest of and S folds plot southeast of a northeast-southwest striking vertical plane of overall monoclinic symmetry. The only sense of movement consistent with the collective asymmetry of the thrust zone folds is top to the northeast. Asymmetric microstructures studied at several localities also indicate top to the northeast movement. Paleomagnetic data suggest that the original sense of thrusting, prior to Neogene vertical axis tectonic rotation related to the San Andreas fault system, was northward. The essential point is that movement of the upper plate of the Chocolate Mountains thrust evidently was continentward. Continentward thrusting suggests a tectonic scenario in which an insular or peninsular microcontinental fragment collided with mainland southern California. The suture predicted by this model is elusive; but the Chocolate Mountains thrust and underlying Orocopia Schist themselves may represent the suture, at the present level of exposure. Alternative tectonic models involving subduction of the Orocopia Schist eastward beneath continental southern California circumvent the suture problem but are presently not supported by any direct structural evidence.

Journal of Geophysical Research↗
USGS2012

Trace element and isotopic constraints on magmatic evolution at Lassen volcanic center

Magmatic evolution at the Lassen volcanic center (LVC) is characterized by a transition from predominantly andesitic to predominantly silicic volcanism with time. Magmas of the andesitic, or “Brokeoff phase” of volcanism range in composition from basaltic andesite io dacite, whereas those of silicic, or “Lassen phase” range in composition from basaltic andesite to rhyolite. The compositions of magmas from each phase define well organized but distinct variation trends. Compared with Brokeoff-phase magmas of similar SiO 2 content, most Lassen-phase magmas contain lower concentrations of most incompatible minor and trace elements. Based on the behavior of both incompatible and compatible trace elements, the geochemical trends defined by the Brokeoff-phase magmas cannot be ascribed merely to fractional crystallization from a single or multiple mafic parental magmas, Moreover, the Lassen-phase magmas cannot be derived from the Brokeoff-phase magmas by fractional crystallization. Rather, the geochemical trends that characterize each volcanic phase define arrays that primarily indicate mixing between well-homogenized silicic and heterogeneous mafic magmas. The distinctive mixing-dominated arrays for each volcanic phase manifest the generation and evolution of two physically distinct, but genetically related magma systems. The LVC magmas have Sr, Nd, and Pb isotopic characteristics that approximate two-component mixing arrays. One isotopic component is similar in composition to that of NE Pacific Ocean ridge and seamount basalts (“MORB” component), the other to mafic Mesozoic granitoids sampled from the neighboring KSamath and Sierra Nevada provinces (“KSN” component). The isotopic compositions of the most silicic LVC magmas lie within the ranges defined by the more mafic LVC magmas, which in turn lie within broad ranges defined by primitive mafic lavas sampled from the Lassen region. The lack of a correlation between the major element and isotopic compositions of LVC magmas seriously limits any model for magmatic evolution that relies on assimilation of old middle to upper crust by isotopically homogeneous mafic magmas during their ascent through the crust. Alternatively, the isotopic and geochemical uniformity of the most silicic magmas of the Brokeoff and Lassen phases suggests that they are well-homogenized partial melts. The likely source region for these silicic melts is the lower crust, which we envision to consist primarily of mafic igneous rocks that are similar in geochemical and isotopic diversity to the regional mafic lavas. Magmatic evolution at LVC can be viewed in terms of a series of mantle melting events that subsequently stimulated meiting in a progressively increasing volume of the lower crust. In general, the LVC magmas represent slightly fractionated mixtures of the mantle-derived mafic magmas and silicic partial melts of the lower crust, the latter melts increasing in relative proportion over the history of the volcanic center. The voluminous rhyolitic lavas and pyroclastic materials erupted during the early Lassen phase represent lower crustal melts thai pooled into sufficient volumes to avoid significant blending and dilution with intruding mafic magmas. The geochemical and isotopic heterogeneity of both Brokeoff-phase andésites and mafic magmatic inclusions in the silicic Lassen-phase magmas must in part record the corresponding diversity of the mantle-derived contributions to arc magmagenesis in this region.

Journal of Geophysical Research Solid Earth↗
USGS1997

Transform push, oblique subduction resistance, and intraplate stress of the Juan de Fuca plate

The Juan de Fuca plate is a small oceanic plate between the Pacific and North America plates. In the southernmost region, referred to as the Gorda deformation zone, the maximum compressive stress σ 1 constrained by earthquake focal mechanisms is N-S. Off Oregon, and possibly off Washington, NW trending left-lateral faults cutting the Juan de Fuca plate indicate a σ 1 in a NE-SW to E-W direction. The magnitude of differential stress increases from north to south; this is inferred from the plastic yielding and distribution of earthquakes throughout the Gorda deformation zone. To understand how tectonic forces determine the stress field of the Juan de Fuca plate, we have modeled the intraplate stress using both elastic and elastic-perfectly plastic plane-stress finite element models. We conclude that the right-lateral shear motion of the Pacific and North America plates is primarily responsible for the stress pattern of the Juan de Fuca plate. The most important roles are played by a compressional force normal to the Mendocino transform fault, a result of the northward push by the Pacific plate and a horizontal resistance operating against the northward, or margin-parallel, component of oblique subduction. Margin-parallel subduction resistance results in large N-S compression in the Gorda deformation zone because the force is integrated over the full length of the Cascadia subduction zone. The Mendocino transform fault serves as a strong buttress that is very weak in shear but capable of transmitting large strike-normal compressive stresses. Internal failure of the Gorda deformation zone potentially places limits on the magnitude of the fault-normal stresses being transmitted and correspondingly on the magnitude of strike-parallel subduction resistance. Transform faults and oblique subduction zones in other parts of the world can be expected to transmit and create stresses in the same manner.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Analysis of stress-induced oval fractures in a borehole at Deep Sea Drilling Project Site 504, eastern equatorial Pacific

Deep Sea Drilling Project (DSDP) Hole 504B is located in the eastern equatorial Pacific Ocean and extends to a total depth of 2111 m beneath the seafloor (mbsf). Several acoustic televiewer logs have been obtained in this well during successive stages of drilling, and the resulting digital images have revealed numerous oval-shaped fractures seemingly etched into the borehole wall. A theoretical examination of these stress-induced features identifies a unique and ephemeral set of stress distributions and magnitudes that are necessary for their production. Consequently, the ovals provide a basis for quantifying the magnitudes and orientations of the maximum and minimum horizontal principal stresses, S H and S h at this site. Vertical, truncated breakouts and horizontal tensile fractures define the spatial boundaries of the ovals. Explicit criteria for their occurrence are combined with estimates for various physical properties of the rock to yield a range of possible values for the horizontal principal stresses. The conspicuous oval geometry is completed by a curved fracture that joins the vertical and horizontal components. Its degree of curvature is delineated by the modified Griffith failure criterion and is directly related to the principal stress difference ( S H − S h ). Matching a series of type curves corresponding to specific values for ( S H − S h ) with the actual undistorted well bore images allows the magnitude of the stress difference to be further constrained. With a value for ( S H − S h ) of 45 ± 5 MPa the individual magnitudes of S H and S h are determined more precisely. Final estimates for the horizontal principal stresses in DSDP Hole 504B at a depth of 1200 mbsf are 141 MPa ≤ S H ≤ 149 MPa and 91 MPa ≤ S h ≤ 109 MPa. Stress magnitudes derived from this approach rely heavily upon the values of a variety of physical properties, and complementary laboratory measurements performed on relevant rock samples provide critical information. Uncertainties in estimating these properties translate into less precise predictions of principal stresses.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Alteration and mineralization of an oceanic forearc and the ophiolite-ocean crust analogy

Mineralogical, chemical, and isotopic (O, C, S, and Sr) analyses were performed on minerals and bulk rocks from a forearc basement section to understand alteration processes and compare with mid-ocean ridges (MOR) and ophiolites. Ocean Drilling Program Hole 786B in the Izu-Bonin forearc penetrates 103 m of sediment and 725 m into volcanic flows, breccias, and basal dikes. The rocks comprise boninites and andesites to rhyolites. Most of the section was affected by low-temperature (<100°C) seawater alteration, with temperatures increasing downward. The rocks are partly (5–25%) altered to smectite, Fe-oxyhydroxide, calcite, and phillipsite, and exhibit gains of K, Rb, and P, loss of Ca, variable changes in Si, Na, Mg, Fe, Sr, and Y, and elevated δ 18 O and 87 Sr/ 86 Sr. Higher temperatures (∼150°C) in the basal dikes below 750 m led to more intense alteration and formation of chlorite-smectite, corrensite, albite, K-feldspar, and quartz (±chlorite). A 5 m thick hydrothermally altered and pyritized zone at 815 m in the basal dikes reacted with mixtures of seawater and hydrothermal fluids to Mg-chlorite, albite, and pyrite, and gained Mg and S and lost Si and Ca. Focused flow of hydrothermal fluids produced serialization halos (Na-K sericite, quartz, pyrophyllite, K-feldspar, and pyrite) along quartz veins at temperatures of 200°–250°C. High 87 Sr/ 86 Sr ratios of chloritized (∼0.7055) and sericitized (∼0.7065) rocks indicate involvement of seawater via mixing with hydrothermal fluids. Low δ 34 S of sulfide (−2 to −5.5‰) and sulfate (12.5‰) are consistent with input of magmatic SO 2 into hydrothermal fluids and disproportionation to sulfide and sulfate. Alteration processes were generally similar to those at MORs, but the arc section is more intensively altered, in part because of the presence of abundant glassy rocks and mafic phases. The increase in alteration grade below 750 m and the mineralization in the basal dikes are analogous to changes that occur near the base of the volcanic section in MOR and the Troodos ophiolite.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Cyclic stressing and seismicity at strongly coupled subduction zones

We use the finite element method to analyze stress variations in and near a strongly coupled subduction zone during an earthquake cycle. Deformation is assumed to be uniform along strike (plane strain on a cross section normal to the trench axis), and periodic earthquake slip is imposed consistent with the long-term rate of plate convergence and degree of coupling. Simulations of stress and displacement rate fields represent periodic fluctuations in time superimposed on an average field. The oceanic plate, descending slab, and continental lithosphere are assumed here to respond elastically to these fluctuations, and the remaining mantle under and between plates is assumed to respond as Maxwell viscoelastic. In the first part of the analysis we find that computed stress fluctuations in space and time are generally consistent with observed earthquake mechanism variations with time since a great thrust event. In particular, trench-normal extensional earthquakes tend to occur early in the earthquake cycle toward the outer rise but occur more abundantly late in the cycle in the subducting slab downdip of the main thrust zone. Compressional earthquakes, when they occur at all, have the opposite pattern. Our results suggest also that the actual timing of extensional outer rise events is controlled by the rheology of the shallow aseismic portion of the thrust interface. The second part of the analysis shows the effects of mantle relaxation on the rate of ground surface deformation during the earthquake cycle. Models without relaxation predict a strong overall compressional strain rate in the continental plate above the main thrust zone, with the strain rate constant between mainshocks. However with significant relaxation present, a localized region of unusually low compressional, or even slightly extensional, strain rate develops along the surface of the continental plate above and somewhat inland from the downdip edge of the locked main thrust zone. The low strain rate starts in the middle or late part of the cycle, depending on position. This result suggests that the negligible or small contraction measured on the Shumagin Islands, Alaska, during 1980 to 1991, may not invalidate an interpretation of that region as being a moderately coupled subduction zone. In contrast, mantle relaxation causes only modest temporal nonuniformity of uplift rates in the overriding plate and of extensional stress rates in the subducting plate, even when the Maxwell time is an order of magnitude less than the recurrence interval.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

CRUST 5.1: A global crustal model at 5° x 5°

We present a new global model for the Earth's crust based on seismic refraction data published in the period 1948-1995 and a detailed compilation of ice and sediment thickness. An extensive compilation of seismic refraction measurements has been used to determine the crustal structure on continents and their margins. Oceanic crust is modeled with both a standard model for normal oceanic crust, and variants for nonstandard regions, such as oceanic plateaus. Our model (CRUST 5.1) consists of 2592 5° × 5° tiles in which the crust and uppermost mantle are described by eight layers: (1) ice, (2) water, (3) soft sediments, (4) hard sediments, (5) crystalline upper, (6) middle, (7) lower crust, and (8) uppermost mantle. Topography and bathymetry are adopted from a Mandard database (ETOPO-5). Compressional wave velocity in each layer is based on field measurements, and shear wave velocity and density are estimated using recently published empirical V p -V s and V p -density relationships. The crustal model differs from previous models in that (1) the thickness and seismic/density structure of sedimentary basins is accounted for more completely, (2) the velocity structure of unmeasured regions is estimated using statistical averages that are based on a significantly larger database of crustal structure, (3) the compressional wave, shear wave, and density structure have been explicitly specified using newly available constraints from field and laboratory studies. Thus - this global crustal model is based on substantially more data than previous models and differs from them in many important respects. A new map of the thickness of the Earth's crust is presented, and we illustrate the application of this model by using it to provide the crustal correction for surface wave phase velocity maps. Love waves at 40 s are dominantly sensitive to crustal structure, and there is a very close correspondence between observed phase velocities at this period and those predicted by CRUST 5.1. We find that the application of crustal corrections to long-period (167 s) Rayleigh waves significantly increases the variance in the phase velocity maps and strengthens the upper mantle , velocity anomalies beneath stable continental regions. A simple calculation of crustal isostacy indicates significant lateral variations in upper mantle density. The model CRUST 5.1 provides a complete description of the physical properties of the Earth's crust at a scale of 5° × 5° and can be used for a wide range of seismological and nonseismological problems.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

The dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component AM3 of the GFDL global coupled model CM3

The Geophysical Fluid Dynamics Laboratory (GFDL) has developed a coupled general circulation model (CM3) for the atmosphere, oceans, land, and sea ice. The goal of CM3 is to address emerging issues in climate change, including aerosol–cloud interactions, chemistry–climate interactions, and coupling between the troposphere and stratosphere. The model is also designed to serve as the physical system component of earth system models and models for decadal prediction in the near-term future—for example, through improved simulations in tropical land precipitation relative to earlier-generation GFDL models. This paper describes the dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component (AM3) of this model. Relative to GFDL AM2, AM3 includes new treatments of deep and shallow cumulus convection, cloud droplet activation by aerosols, subgrid variability of stratiform vertical velocities for droplet activation, and atmospheric chemistry driven by emissions with advective, convective, and turbulent transport. AM3 employs a cubed-sphere implementation of a finite-volume dynamical core and is coupled to LM3, a new land model with ecosystem dynamics and hydrology. Its horizontal resolution is approximately 200 km, and its vertical resolution ranges approximately from 70 m near the earth’s surface to 1 to 1.5 km near the tropopause and 3 to 4 km in much of the stratosphere. Most basic circulation features in AM3 are simulated as realistically, or more so, as in AM2. In particular, dry biases have been reduced over South America. In coupled mode, the simulation of Arctic sea ice concentration has improved. AM3 aerosol optical depths, scattering properties, and surface clear-sky downward shortwave radiation are more realistic than in AM2. The simulation of marine stratocumulus decks remains problematic, as in AM2. The most intense 0.2% of precipitation rates occur less frequently in AM3 than observed. The last two decades of the twentieth century warm in CM3 by 0.32°C relative to 1881–1920. The Climate Research Unit (CRU) and Goddard Institute for Space Studies analyses of observations show warming of 0.56° and 0.52°C, respectively, over this period. CM3 includes anthropogenic cooling by aerosol–cloud interactions, and its warming by the late twentieth century is somewhat less realistic than in CM2.1, which warmed 0.66°C but did not include aerosol–cloud interactions. The improved simulation of the direct aerosol effect (apparent in surface clear-sky downward radiation) in CM3 evidently acts in concert with its simulation of cloud–aerosol interactions to limit greenhouse gas warming.

Journal of Climate↗
USGS2012

Rift flank uplifts and Hinterland Basins: Comparison of the Transantarctic Mountains with the Great Escarpment of southern Africa

Uplifted rift margins are a common feature of continents and oceans. Two variants of rift flank morphologies have been recognized: One in which the topography warps down from an inland high toward the continental margin, and one where the tropographic peak lies close to the continental margin. The Great Escarpment of southern Africa and the Transantarctic Mountains are examples of the first and the second variants of rift flanks, respectively. Both rift flanks are bordered on their landward side by broad continental basins: the Kalahari and the Wilkes hinterland basins. If these basins are interpreted as flexural “outer lows” that deepen in unison with the uplift of the rift flanks, the lithosphere on the uplifted side is very rigid in both cases (elastic thickness T e of 100 ± 20 km for southern Africa and 110 ± 20 km for East Antarctica). We suggest that the variation in rift flank morphology is caused by the isostatic response to uplift forces of elastic plates sharing different boundary conditions. We model the uplift of the Transantarctic Mountains as an upward deflection of an elastic plate which is broken at the front of the Transantarctic Mountains, and we model the uplift of the Great Escarpment as an upward deflection of a continuous elastic plate that is modified by the downward load of sediments on the continental margin. Although the Transantarctic Mountain uplift is young (60–0 Ma) and the southern African uplift is old (<100 Ma), the different isostatic responses of the two margins are not a function of age, because most loading (sedimentation) and unloading (erosion) took place shortly after rifting. Detailed modeling of topography, gravity, geological markers, and the locations of depocenters suggests that lithospheric rigidity decreases under the Transantarctic Mountains, whereas in southern Africa the decrease occurs not under the Great Escarpment but far seaward under the continental shelf and slope. If the distribution of lithospheric rigidity is indicative of the thermal regime of the lithosphere, then uplifted rift flanks are not always underlain by a thermal anomaly. This and other geological evidence indicate that a single mechanism cannot explain the uplift of both the Antarctic and the African margins.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Tidal triggering of low frequency earthquakes near Parkfield, California: Implications for fault mechanics within the brittle-ductile transition

Studies of nonvolcanic tremor (NVT) have established the significant impact of small stress perturbations on NVT generation. Here we analyze the influence of the solid earth and ocean tides on a catalog of ∼550,000 low frequency earthquakes (LFEs) distributed along a 150 km section of the San Andreas Fault centered at Parkfield. LFE families are identified in the NVT data on the basis of waveform similarity and are thought to represent small, effectively co-located earthquakes occurring on brittle asperities on an otherwise aseismic fault at depths of 16 to 30 km. We calculate the sensitivity of each of these 88 LFE families to the tidally induced right-lateral shear stress (RLSS), fault-normal stress (FNS), and their time derivatives and use the hypocentral locations of each family to map the spatial variability of this sensitivity. LFE occurrence is most strongly modulated by fluctuations in shear stress, with the majority of families demonstrating a correlation with RLSS at the 99% confidence level or above. Producing the observed LFE rate modulation in response to shear stress perturbations requires low effective stress in the LFE source region. There are substantial lateral and vertical variations in tidal shear stress sensitivity, which we interpret to reflect spatial variation in source region properties, such as friction and pore fluid pressure. Additionally, we find that highly episodic, shallow LFE families are generally less correlated with tidal stresses than their deeper, continuously active counterparts. The majority of families have weaker or insignificant correlation with positive (tensile) FNS. Two groups of families demonstrate a stronger correlation with fault-normal tension to the north and with compression to the south of Parkfield. The families that correlate with fault-normal clamping coincide with a releasing right bend in the surface fault trace and the LFE locations, suggesting that the San Andreas remains localized and contiguous down to near the base of the crust. The deep families that have high sensitivity to both shear and tensile normal stress perturbations may be indicative of an increase in effective fault contact area with depth. Synthesizing our observations with those of other LFE-hosting localities will help to develop a comprehensive understanding of transient fault slip below the “seismogenic zone” by providing constraints on parameters in physical models of slow slip and LFEs.

California↗