Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Planets”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72Linked to original sources

Nonlinear strain buildup and the earthquake cycle on the San Andreas fault

Two contrasting models of the earthquake deformation cycle on strike slip faults predict significant temporal declines in shear strain rate near the fault, accompanied by a progressive broadening of the zone of deformation adjacent to it. In the thin lithosphere model, transient deformation results from flow in the asthenosphere due to stress relaxation following faulting through most or all of the lithosphere. For an earth model with a thick elastic lithosphere (plate thickness » depth of seismic slip), transient motions are due to postearthquake aseisrnic slip below the coseismic fault plane. Data from the San Andreas fault indicate a long-term temporal decrease in strain rate that persists for at least 30 years and may extend through the entire earthquake cycle. Observations support a cycle-long rate decrease and a temporal spreading of the deformation profile only if movement cycles on the northern and southern locked sections of the fault are basically similar. If so, the usually lower strain rates and broader deformation zone currently observed on the southern San Andreas represent a later evolutionary stage of the northern locked section, where a great earthquake is a more recent occurrence. Although the data allow some extreme models to be discarded, no sufficiently strong constraints exist to decide between the thin and thick lithosphere models. Regardless of the appropriate model the geodetic observations themselves indicate that strain buildup is sufficiently nonlinear to cause significant departures from recurrence estimates based on linear strain accumulation and the time-predictable model.

Journal of Geophysical Research Solid Earth↗

The use of MAGSAT data to determine secular variation

A combined spatial and secular variation model of the geomagnetic field, labeled M061581, is derived from a selection of Magsat data. Secular variation (SV) data computed from linear fits to midnight hourly values from 19 magnetic observatories were also included in the analysis but were seen to have little effect on the model. The SV patterns from this new model are compared with those from the 1980 IGRF and with those for 1970 computed by the DGRF and with the 1960 patterns computed using the GSFC(12/66) model. Most of the features of the M061581 are identical in location and level with those of the 1980 IGRF. Together they confirm that the reversals in sign of field change seen over Asia and North America between 1965 and 1975 are reverting to the pre-1965 states. The M061581 model gives −32 nT/yr for the dipole decay rate, larger than the 70% increase already reported since 1965. This abnormally high value is interpreted as being a defect of the model because it appears to result from a much larger (−100 nT/yr) drop in field over the polar regions not indicated by the 1980 IGRF. This north polar decrease is shown to be of external origin as the result of a combination of the seasonal effect of the north polar ionospheric (counterclockwise) afternoon Sq p cell increasing in intensity from the beginning (November 1979) to the end (June 1980) of the Magsat data collection period, coupled with an enhancement of its effect as the orbit lowered from the 350– 550 km initial altitudes to near 200 km just prior to burnup. This experiment indicates that secular variation can be obtained from satellite data for intervals of less than a full year if corrections can be made for seasonal effects and that ‘annual snapshots’ of the field by a satellite would allow easy and accurate models of secular change without the use of any surface data.

Journal of Geophysical Research Solid Earth↗

The 1979 Homestead Valley earthquake sequence, California: Control of aftershocks and postseismic deformation

The coseismic slip and geometry of the March 15, 1979, Homestead Valley, California, earthquake sequence are well constrained by precise horizontal and vertical geodetic observations and by data from a dense local seismic network. These observations indicate 0.52±0.10 m of right-lateral slip and 0.17±0.04 m of reverse slip on a buried vertical 6-km-long and 5-km-deep fault and yield a mean static stress drop of 7.2±1.3 MPa. The largest shock had M S = 5.6. Observations of the ground rupture revealed up to 0.1 m of right-lateral slip on two mapped faults that are subparallel to the modeled seismic slip plane. In the 1.9 years since the earthquakes, geodetic network displacements indicate that an additional 60±10 mm of postseismic creep took place. The rate of postseismic shear strain (0.53±0.13 μrad/yr) measured within a 30×30-km network centered on the principal events was anomalously high compared to its preearthquake value and the postseismic rate in the adjacent network. This transient cannot be explained by postseismic slip on the seismic fault but rather indicates that broadscale release of strain followed the earthquake sequence. We have calculated the postearthquake stress field caused by the modeled coseismic slip. We assume that failure is promoted when the sum of the shear stress plus 0.75 times the faultopening stress increases. Most aftershocks concentrate at points where the stresses are enhanced by 0.3 MPa (3 bars) or more; aftershocks are nearly absent where postearthquake stresses decrease by 0.3–0.5 MPa. Isolated off-fault clusters of aftershocks that locate at one fault length from the rupture plane are explainable by this hypothesis. We find that ground rupture and postseismic creep take place where near-surface stresses are calculated to increase within the preexisting fault zones. Two patches that extend 4 km from both ends of the seismic fault exhibited neither aftershocks nor measurable postseismic creep. The sensitivity of aftershocks and ground rupture to changes in stress that are less than 5% of the earthquake stress drop demonstrates that the region around the earthquakes was within a few percent of its failure threshold before the main shocks. The preearthquake stress field and the stress required for failure must also have been nearly uniform.

Journal of Geophysical Research Solid Earth↗

Earthquakes in the Orozco transform zone: seismicity, source mechanisms, and tectonics

As part of the Rivera Ocean Seismic Experiment, a network of ocean bottom seismometers and hydrophones was deployed in order to determine the seismic characteristics of the Orozco transform fault in the central eastern Pacific. We present hypocentral locations and source mechanisms for 70 earthquakes recorded by this network. All epicenters are within the transform region of the Orozco Fracture Zone and clearly delineate the active plate boundary. About half of the epicenters define a narrow line of activity parallel to the spreading direction and situated along a deep topographic trough that forms the northern boundary of the transform zone (region 1). Most focal depths for these events are very shallow, within 4 km of the seafloor; several well-determined focal depths, however, are as great as 7 km. No shallowing of seismic activity is observed as the rise-transform intersection is approached; to the contrary, the deepest events are within 10 km of the intersection. First motion polarities for most of the earthquakes in region 1 are compatible with right-lateral strike slip faulting along a nearly vertical plane, striking parallel to the spreading direction. Another zone of activity is observed in the central part of the transform (region 2). The apparent horizontal and vertical distribution of activity in this region is more scattered than in the first, and the first motion radiation patterns of these events do not appear to be compatible with any known fault mechanism. Pronounced lateral variations in crustal velocity structure are indicated for the transform region from refraction data and measurements of wave propagation directions. The effect of this lateral heterogeneity on hypocenters and fault plane solutions is evaluated by tracing rays through a three-dimensional velocity grid. While findings for events in region 1 are not significantly affected, in region 2, epicentral mislocations of up to 10 km and azimuthal deflections of up to 45° may result from assuming a laterally homogeneous velocity structure. When corrected for the effects of lateral heterogeneity, the epicenters and fault plane solutions for earthquakes in region 2 are compatible with predominantly normal faulting along a topographic trough trending NW–SE; the focal depths, however, are poorly constrained. These results suggest an en echelon spreading center or leaky transform regime in the central transform region.

Journal of Geophysical Research B: Solid Earth↗

Precision of geodolite surveys: A reply to Jackson and Cheng

Jackson and Cheng have suggested that the changes in areal dilatation measured in Geodolite surveys by the U.S. Geological Survey may be simply an artifact of the measuring system. Although systematic error could conceivably account for the observed excursions in dilatation, we maintain that the specific criticisms by Jackson and Cheng are incorrect: the excursions in dilatation cannot be attributed to the offset correction nor to proportional error associated with temperature. The absence of both errors is demonstrated by using data that are particularly sensitive to the two effects.

Journal of Geophysical Research B: Solid Earth↗

Experimental deformation of polycrystalline H2O ice at high pressure and low temperature: Preliminary results

Interest in the mechanical properties of water ice under the conditions in which it exists in the outer solar system has motivated the development and use of a new high‐pressure, low‐temperature triaxial deformation apparatus. Constant displacement rate tests on 70 samples of pure polycrystalline water ice have been performed at temperatures 77≤≤258 K, confining pressures 0.1≤≤350 MPa, and strain rates 3.5×10≤ ≤3.5×10 s. In most cases, the ice polymorph tested was ice . Both brittle and ductile behavior have been observed. Brittle behavior of ice, promoted by lower pressure, lower temperature, and higher strain rate, is analogous to that in rocks, with the important exception that brittle fracture strength becomes independent of confining pressure above 50 MPa pressure and the fracture angle is approximately 45° to the loading direction (i.e., the coefficient of internal friction is approximately zero). Ductile flow, the predominant behavior in our tests at ≥195 K, follows a law of form = σ exp (−*/) (σ is stress; is the gas constant; , * are material constants). Three sets of material constants are required to fit the data, with changes in sets (or mechanisms) occurring near 243 K and 195 K. The value of remains near 4 throughout the measured ductile field, but * drops from 91 to 61 to 31 /mole as temperature decreases. The maximum brittle strength measured was 171 MPa; the maximum ductile strength measured was 91 MPa. At confining pressures near the phase transition pressure of ice → ice II, the ductile strength is observed to drop dramatically. Some overlap with previous work occurs at higher temperatures and lower pressures. Agreement with present work is generally good, both quantitatively in the values of and *, and qualitatively in the mechanism of deformation. Although the ductile strengths measured here are somewhat higher than expected on the basis of extrapolations of previous work, the low value of * at <195 K indicates that the ice layer on icy bodies in the solar system is much weaker than has generally been predicted.

Journal of Geophysical Research B: Solid Earth↗

Chryse Basin channels: low-gradients and ponded flows.

Gradients on the floors of the Martian outflow channels that are derived from radar-elevation profiles across Lunae Planum and Chryse Basin have much lower values than those obtained from the U.S. Geological Survey's topographic map. Whereas the gradients of Maja and Ares Valles are similar to those of the catastrophic flood channels in the Scablands of Washington State, the gradients of Simud and Tiu Valles are essentially level, and the movement of fluids to the N poses problems. It is proposed that ponding may have formed lakes in depressions associated with the Valles Marineris grabens, ancient craters in the chaotic terrain area, and possibly even the regional low where most chaotic terrains occur. It is envisaged that lakes eventually overflowed, forming the present channels. When dams broke, floods were released catastrophically, with a final gigantic flood from the Valles Marineris system of troughs, which would have had sufficient head to move fluids across nearly level gradients through the Simud and Tiu channels.

Journal of Geophysical Research↗

Petrology and comparative thermal and mechanical histories of clasts in breccia 62236

Lunar breccia 62236 contains large lithic fragments of troctolite, norite, and anorthosite. The mafic phases, olivine, inverted pigeonite, and augite, fill interstitial areas between larger plagioclases and appear to be cumulate phases with extensive adcumulus growth. Pyroxene compositional homogeneity indicates that cation exchange during cooling was limited to an area of about 1 mm. Primary augite and pigeonite both contain 30 μm-wide lamellae of the other along ‘001’. Pigeonite inverted to orthopyroxene without retaining any crystallographic orientation and subsequently exsolved fine lamellae of augite on (100). Profiles across orthopyroxene-augite interfaces obtained in the analytical transmission electron microscope show an increase of ∼5% wollastonite in augite within 0.5 μm of the interface, suggesting that short-range cation exchange continued to temperatures below 500°C. The entire sample has undergone heterogeneous shock deformation. Shock melting of the troctolitic clast suggests pressures of 200–300 kb and well-developed basal twinning in augite from the norite clasts suggests pressures of 50–300 kbars. The present evidence indicates that 62236 contains parts of a slowly cooled microlayered adcumulate that has been heterogeneously shocked several times and combined into the present breccia.

Journal of Geophysical Research Solid Earth↗

Strangways Crater, Northern Territory, Australia: Siderophile element enrichment and lithophile element fractionation

The Strangways Crater, Northern Territory, Australia (15°12′S, 133°35′ E ), has a central core, about 10 km in diameter, of shocked granitic gneiss and amphibolite, and some remnants of a melt rock sheet, surrounded by outer rings of quartzite and siltstone to a diameter of 20–25 km. Seven samples of melt rock (six granitic melts, one shale melt clast) and four samples of country rock (granitic gneiss, amphibolite, shale, quartzite) were analyzed by neutron-activation analysis: for Sc, Cr, Fe, Co, Zn, Rb, Zr, Sb, Cs, Ba, rare earth elements, Hf, Ta, Th, and U, the samples were analyzed instrumentally; and for Ni, Se, Pd, Ag, Cd, Re, Os, Ir, and Au, they were analyzed radiochemically. Siderophile elements are significantly enriched in the granitic melt rocks relative to country rocks; for example, the Ir enrichments range from 0.6 to 2.8 ppb. The low Ir/Ni ratio (∼0.16 relative to C1 chondrites) excludes a chondritic impacting body, and Cr enrichment argues against impact by an iron meteorite. The Strangways Crater may have been formed by the impact of an olivine-rich achondrite and melt rocks appear to contain about 3 wt.% of projectile material. The composition of the granitic melt rocks cannot be reproduced by any simple mixture of analyzed country rock types and chemical fractionation by selective shock melting appears to have taken place.

Journal of Geophysical Research Solid Earth↗

Small-scale features in the Earth's magnetic field observed by Magsat

A spherical harmonic expansion to degree and order 29 is derived using a selected magnetically quiet sample of Magsat data. Global maps representing the contribution due to terms of the expansion above n = 13 at 400 km altitude are compared with previously published residual anomaly maps and shown to be similar, even in polar regions. An expansion with such a high degree and order displays all but the sharpest features seen by the satellite and gives a more consistent picture of the high-order field structure at a constant altitude than do component maps derived independently.

Journal of Geophysical Research Solid Earth↗

A seismic refraction survey of the Imperial Valley Region, California

The U.S. Geological Survey conducted an extensive seismic refraction survey in the Imperial Valley region of California in 1979. The Imperial Valley is located in the Salton Trough, an active rift between the Pacific and North American plates. Forty shots fired at seven shot points were recorded by 100 portable seismic instruments at typical spacing of 0.5–1 km. More than 1300 recording locations were occupied, and more than 3000 usable seismograms were obtained. We analyzed five profiles using a standard ray‐tracing program, constructed a contour map of reduced travel times from our most widely recorded shot point, and modeled an existing gravity profile across the Salton Trough. Results are itemized: (1) All models have in common a sedimentary layer ( V p = 1.8–5.0 km/s), a “transition zone” ( V p = 5.0–5.65 km/s), a basement ( V p = 5.65 km/s in the Imperial Valley, 5.9 km/s on the bordering mesas), and subbasement ( V p = 7.2 km/s). (2) The sedimentary layer ranges in thickness along the axis of the Salton Trough from 3.7 km (Salton Sea) to 4.8 km (U.S.‐Mexican border). On the bordering mesas it is quite variable in thickness. (3) The “transition” zone is about 1 km thick in most places. In the Imperial Valley there are no marked velocity discontinuities in this zone between the sedimentary layer and basement. On the bordering mesas, however, there is a discontinuity at the top of this zone. (4) There are apparently two types of basement. On the bordering mesas, basement is crystalline igneous and metamorphic rocks. In the Imperial Valley, basement is mostly lower‐greenshist‐facies sedimentary rocks, based primarily on the smooth transition in character from sediment to basement arrivals, the low value of basement velocity, and the fact that deep (4 km) wells in the valley penetrate only the upper part of the known Cenozoic stratigraphic column for the Salton Trough. (5) The subbasement, or intermediate crustal layer, ranges in depth along the axis of the Salton Trough from 16 km (Salton Sea) to 10 km (U.S.‐Mexican border). Gravity modeling requires that this layer deepen and/or pinch out beneath the bordering mesas and mountain ranges. Based on its high velocity and the presence of intrusive basaltic rocks in the sedimentary section in the Imperial Valley, the subbasement is thought to be a mafic intrusive complex similar to oceanic middle crust. (6) Several structures are seen that affect basement, transition zone, and deeper parts of the sedimentary layer. They include a scarp along the Imperial fault, as much as 1 km down to the northeast, and a scarp passing roughly along the topographic boundary between the Imperial Valley and the bordering mesa to the west, as much as 3½ km down to the east. We interpret the latter scarp to be the suture, or rift boundary, between the older crystalline basement on the mesa and the younger metasedimentary basement in the Imperial Valley. (7) On a contour map of reduced travel time from our most widely recorded shot point, subtle patches of early arrivals among otherwise late arrivals in the central Imperial Valley correlate well with known geothermal resource areas having reservoir temperatures of more than 150°C. Apparently the Salton Trough is a location where new crust is being generated. As the rift opens, mafic intrusive rocks fill it from below as sedimentary rocks fill it from above. Rifting and intrusion produce high heat flow that metamorphoses the sedimentary rocks to shallow depth (metasedimentary basement in the Imperial Valley) and thus consolidates the new crust.

California↗

Local gravity anomalies produced by dislocation sources

Rundle (1978) and Walsh and Rice (1979) have shown that the change in the vertical component of gravity is proportional to uplift for a spherical source of dilatation and for slip on an infinitely long dip-slip fault. In the first case, no free air gravity anomaly is produced and in the second case no Bouguer gravity anomaly. Gravity anomalies due to other dislocation sources in three dimensions are as follows: For strike-slip faulting the ratio of the gravity change to uplift depends upon position; however, the gravity change contours are roughly similar to those corresponding to a zero free air gravity anomaly. Nor is the ratio constant for dip-slip faulting except for the two special cases of dip slip on a vertical fault and horizontal slip on a horizontal fault, neither of which produce a Bouguer anomaly. The Bouguer anomaly produced by an open horizontal crack is the same as would be produced had the material within the crack been mined out without deforming the solid. If the horizontal crack were filled with material of density equal to the host rock (a good approximation to sill formation), no Bouguer anomaly is produced. For cracks of other inclinations the ratio of gravity change to uplift is not constant. Thus, dilatancy, in general, does not correspond to the absence of a free air anomaly, as might be suggested by the special case of a spherical source of dilatation. For two-dimensional models a cylindrical source of dilatation produces no free air gravity anomaly, dip-slip faulting produces no Bouguer anomaly, and open cracks produce a Bouguer anomaly equal to that which would be produced had the material within the crack been mined out without deforming the solid. A two-dimensional crack filled with material of density equal to that of the host rock would produce no Bouguer anomaly. Jachens et al. (1983) have reported temporal changes in gravity, elevation, and area strain along the San Andreas fault in southern California such that the Bouguer anomaly apparently remains unchanged and the uplift-to-strain ratio is about −100 km. Several dislocation mechanisms are proposed that fulfill these constraints, but these mechanisms appear to be rather contrived and are not regarded as satisfactory explanations.

Journal of Geophysical Research Solid Earth↗

Regional deformation near Palmdale, California, 1973-1983 (USA)

The Tehachapi trilateration network spans the intersection of the San Andreas and Garlock faults in southern California in the “Big Bend” region of the San Andreas fault. Analysis of data from 1973–1983 shows strain differences between the northwest and southeast regions of the network and slip at depth on both faults. The Palmdale network, spanning the San Andreas fault entirely within the Tehachapi network, showed increases of about 1 μstrain in both east‐west and north‐south extension in late 1979. The Tehachapi strains also jumped at this time, but the magnitude of the increase was only about one third that of Palmdale. The principal strain rates for Tehachapi over the time interval 1973–1983 are μstrain/yr and μstrain/yr, with the 1 axis directed N76°E. Strains were also computed for two subregions. The principal strain rates for the southeast Tehachapi region are μstrain/yr and μStrain/yr, with the 1 axis directed N73°E. This result differs significantly from the principal strain rates at Palmdale ( μstrain/yr, μstrain/yr, with the 1 axis directed N71°E). For the northwest Tehachapi region, μstrain/yr and μstrain/yr, with the 1 axis directed N93°E.

California↗

A seismic refraction study of the Oregon Cascades

A 275‐km‐long reversed refraction profile in the Oregon Cascades, two shallow earthquakes of magnitude 5 in southern Washington, a shallow earthquake of magnitude 4.6 in northern California, and a previously published analysis of the Bouguer gravity field are used to develop a crustal P wave velocity model for the Oregon Cascades. Travel time analysis of the refraction profile indicates a crustal structure characterized by surface layers with P wave velocities that vary from 2.9 to 5.2 km/s and thicknesses that vary from 2.5 to 5.0 km, upper crustal velocities of 6.1 to 6.5 km/s between the depths of 3 and 29 km, lower crustal velocities near 7.0 km/s between the depths of 29 and 44 km, and a mantle reflector at a depth of 44 km. Comparison with synthetic seismograms supports this structure and shows that the lower crustal and Mono transitions can be better modeled by continuous velocity gradients than by first‐order discontinuities. Arrival times from two shallow earthquakes in southern Washington in 1981 across 14 telemetry stations in the Oregon Cascades show an apparent velocity of 7.62 km/s. On the basis of a previously published gravity analysis, the mean north‐south component of the dip to the crust‐mantle boundary is estimated at 1° down dip to the south. The apparent velocity from the earthquakes and the dip estimate from gravity indicate that the true Pn velocity is 7.70 km/s. Arrival times from a shallow earthquake in northern California in 1978 across telemetry stations in the Washington Cascades are consistent with this model and indicate that the upper mantle velocity reaches 8.20 km/s at a depth of approximately 100 km.

Oregon↗

Deformation of clinopyroxenite: Evidence for a transition in flow mechanisms and semibrittle behavior

A systematic suite of constant strain rate experiments was performed on a vacuum-dried, high-purity, fine-grained clinopyroxenite using NaCl and NaF as confining media in a Griggs-type piston-cylinder apparatus. The experiments were carried out over a range of temperatures from 400° to 1100°C, strain rates from 10 −3 to 10 −7 s −1 , and confining pressures from 170 to 1990 MPa. At T = 600°C and = 1.1 × 10 −5 s −1 , three modes of deformation occur with increasing confining pressure: (1) Macroscopic faulting associated with low strength and stress drops, (2) stable microfracturing and plastic deformation associated with pressure-dependent strength, and (3) plastic deformation (mechanical twinning and 〈001〉 slip) with high strengths which are insensitive to pressure variations. In experiments at P = 1500 MPa, within this high-pressure plastic mode, two regimes of flow are clearly defined. At low to intermediate temperatures and high strain rates, flow strengths are insensitive to changes in strain rate and temperature. Optical and transmission electron microscope observations indicate that plastic strain is accomplished by mechanical twinning on (100) and (001) and by {100}〈001〉 slip. In contrast, at high temperatures and low strain rates the flow stress is strongly dependent on temperature and strain rate. Specimens deformed in this regime show evidence of recovery, multiple slip, and recrystallization; and plastic strain is much more homogeneous. The flow data within each regime can be satisfactorily fit to thermally activated power laws. In the low-temperature regime n (the stress exponent) = 83 ± 16 and E * (the activation energy for flow) = 220 ± 40 kJ/mol. We believe that these parameters reflect flow dominated by the kinetics of dislocation glide associated with mechanical twinning and (100)〈001〉 slip. In the high-temperature regime, n = 5.3 ± 1.1 and E * = 380 ± 30 kJ/mol. These parameters describe creep by multiple slip accompanied by increased rates of diffusion and recovery.

Journal of Geophysical Research Solid Earth↗

Permeability of fault gouge under confining pressure and shear stress

The permeability of both clay-rich and non-clay gouges, as well as several pure clays, was studied as a function of confining pressures from 5 to 200 MPa and shear strain to 10. Permeability ranged over four orders of magnitude, from around 10 −22 to 10 −18 m 2 (1 darcy = 0.987 × 10 −12 m 2 ). The lowest values were characteristic of the montmorillonite-rich and finer grained non-clay gouges. Illite, kaolinite, and chlorite had intermediate permeabilities, while the highest values were typical of the serpentine and coarser grained non-clay gouges. Grain size was an important factor in determining permeability, particularly for the clay-rich samples. The coarse grained gouges were the most permeable and decreased in permeability after shearing. Conversely, the fine grained gouges had characteristically lower permeabilities that did not vary significantly after various amounts of shearing. The permeabilities of the non-clay samples were not significantly different than those of the clays. Therefore, comminuted rock flours can be equally as effective in reducing the flow of water as the characteristically low permeability clay gouges. The strengths of the samples were quite variable. The non-clay gouges were consistently the strongest, with yield points (beginning of nonelastic behavior) around 850 MPa, while montmorillonite had an anomalously low strength in relation to all the other gouges at 250 MPa. Strength of the saturated samples under drained (low pore pressure) conditions did not correlate with high or low permeability. However, the low permeabilities of these gouges could be a factor in the measured low shear stresses along fault regions if excess pore pressures were created as a result of shearing or compaction, and this pressure was unable to dissipate through a thick section of the material.

Journal of Geophysical Research Solid Earth↗

Origin of Hawaiian tholeiite: A metasomatic model

Two voluminous magma types generated in the mantle underlying the Pacific plate are mid‐ocean ridge tholeiite (MORB) erupted at the East Pacific Rise spreading center and Hawaiian tholeiite (HT) erupted above the Hawaiian hot spot or melting anomaly. MORB has low initial 87 Sr/ 86 Sr ratios and low amounts of all incompatible trace elements including rare earths; chondrite‐normalized patterns are depleted in light rare earths. HT, by contrast, has higher initial 87 Sr/ 86 Sr and higher amounts of incompatible trace elements; chondrite‐nor‐malized patterns are enriched in the middle and light rare earths. HT is generally poorer in CaO and Al 2 O 3 and much richer in total iron and TiO 2 compared with MORB having the same MgO content. Primary magma compositions for the two volcanic systems are calculated in Fe‐Mg equilibrium with residual olivine (Fo 92 ). MORB is generated by partial melting of a trace element depleted Iherzolite source leaving a residual assemblage dominated by olivine and orthopyroxene. The percentage of partial melting for a primary magma containing 15% MgO is calculated to be 35–42% in a source mantle having a heavy rare earth content of 3×chondrite and 33–35% MgO. HT, represented by Kilauea tholeiite, is generated by partial melting of a mixture of unmelted and residual mantle for MORB which has been modified by metasomatic addition of a nephelinitic fluid, amphibole, and minor amounts of apatite and Fe‐bearing phases such as sulfide and magnetite/ilmenite. This model yields a picritic magma in equilbrium with magnesian dunite at high (>40%) degrees of partial melting. The source also has 35% MgO before partial melting. Melting in both systems in polyvariant and not controlled by lower‐temperature invariant equilibria. The low‐velocity zone is considered to be the source of metasomatic fluids that are driven upward into the lowermost lithosphere in response to a thermal plume. Picritic primary magmas are produced by shear melting, localized in the zone of thinned and metasomatized lithosphere beneath Hawaii. Melt extraction is rapid and episodic at intervals of months to decades; magma is not stored in the mantle but passes upward to a plexus of storage reservoirs located 2–6 km beneath the surface of Kilauea. Kilauea primary magmas fractionate olivine during upward transport to reach bulk compositions of 13–14% MgO in storage. Different magma batches erupted to the surface, distinguished by different major and minor element compositons compared at similar MgO content, represent combinations of differing degrees of metasomatic enrichment, differing degrees of partial melting, and some effects of premelting mantle heterogeneity.

Hawaii↗

Introduction and digest to the Special Issue on Chemical Effects of Water on the Deformation and Strengths of Rocks

The important role of pore pressure in promoting such brittle processes as cataclasis, hydraulic fracturing, large‐scale faulting, and earthquakes within the crust is widely accepted in geology and geophysics [, 1957; , 1959; , 1958; ., 1963; , 1968; ., 1968; ., 1976; , 1973, 1980; , 1981]. Provided that fluid pressure is fully communicated with rock pore space, the effective normal stresses that control crack growth, macroscopic fracture, and friction are reduced by the magnitude of the fluid pressui'e. Beyond this physical effect of pore fluids, there are chemical effects of water on the strength of rocks that are also important in governing differential stresseg and flow in the continental crust. Some of these chemical effects of water on rock deformation have long been recognized.

Journal of Geophysical Research B: Solid Earth↗