USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geophysical Journal”

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,693 records · Page 94Linked to original sources
USGS1996

The 40Ar/39Ar and K/Ar dating of lavas from the Hilo 1-km core hole, Hawaii Scientific Drilling Project

Mauna Kea lava flows cored in the Hilo hole range in age from <200 ka to about 400 ka based on 40 Ar/ 39 Ar incremental heating and K‐Ar analyses of 16 groundmass samples and one coexisting plagioclase. The lavas, all subaerially deposited, include a lower section consisting only of tholeiitic basalts and an upper section of interbedded alkalic, transitional tholeiitic, and tholeiitic basalts. The lower section has yielded predominantly complex, discordant 40 Ar/ 39 Ar age spectra that result from mobility of 40 Ar and perhaps K, the presence of excess 40 Ar, and redistribution of 39 Ar by recoil. Comparison of K‐Ar ages with 40 Ar/ 39 Ar integrated ages indicates that some of these samples have also lost 39 Ar. Nevertheless, two plateau ages of 391 ± 40 and 400 ± 26 ka from deep in the hole, combined with data from the upper section, show that the tholeiitic section accumulated at an average rate of about 7 to 8 m/kyr and has a mean recurrence interval of 0.5 kyr/flow unit. Samples from the upper section yield relatively precise 40 Ar/ 39 Ar plateau and isotope correlation ages of 326 ± 23, 241 ± 5, 232 ± 4, and 199 ± 9 ka for depths of −415.7 m to −299.2 m. Within their uncertainty, these ages define a linear relationship with depth, with an average accumulation rate of 0.9 m/kyr and an average recurrence interval of 4.8 kyr/flow unit. The top of the Mauna Kea sequence at −280 m must be older than the plateau age of 132 ± 32 ka, obtained for the basal Mauna Loa flow in the corehole. The upward decrease in lava accumulation rate is a consequence of the decreasing magma supply available to Mauna Kea as it rode the Pacific plate away from its magma source, the Hawaiian mantle plume. The age‐depth relation in the core hole may be used to test and refine models that relate the growth of Mauna Kea to the thermal and compositional structure of the mantle plume.

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↗
USGS1996

Crustal structure of the Colorado Plateau, Arizona: Application of new long-offset seismic data analysis techniques

The Colorado Plateau is a large crustal block in the southwestern United States that has been raised intact nearly 2 km above sea level since Cretaceous marine sediments were deposited on its surface. Controversy exists concerning the thickness of the plateau crust and the source of its buoyancy. Interpretations of seismic data collected on the plateau vary as to whether the crust is closer to 40 or 50 km thick. A thick crust could support the observed topography of the Colorado Plateau isostatically, while a thinner crust would indicate the presence of an underlying low-density mantle. This paper reports results on long-offset seismic data collected during the 1989 segment of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that extended from the Transition Zone into the Colorado Plateau in northwest Arizona. We apply two new methods to analyze long-offset data that employ finite difference travel time calculations: (1) a first-arrival time inverter to find upper crustal velocity structure and (2) a forward-modeling technique that allows the direct use of the inverted upper crustal solution in modeling secondary reflected arrivals. We find that the crustal thickness increases from 30 km beneath the metamorphic core complexes in the southern Basin and Range province to about 42 km beneath the northern Transition Zone and southern Colorado Plateau margin. We observe some crustal thinning (to ∼37 km thick) and slightly higher lower crustal velocities farther inboard; beneath the Kaibab uplift on the north rim of the Grand Canyon the crust thickens to a maximum of 48 km. We observe a nonuniform crustal thickness beneath the Colorado Plateau that varies by ∼15% and corresponds approximately to variations in topography with the thickest crust underlying the highest elevations. Crustal compositions (as inferred from seismic velocities) appear to be the same beneath the Colorado Plateau as those in the Basin and Range province to the southwest, implying that the plateau crust represents an unextended version of the Basin and Range. Some of the variability in crustal structure appears to correspond to preserved lithospheric discontinuities that date back to the Proterozoic Era.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Hydrogeology of the Hawaii Scientific Drilling Project borehole KP-1 2. Groundwater geochemistry and regional flow patterns

A series of downhole and surface water samples were taken from the 1‐km‐deep KP‐1 borehole located on the eastern flank of the island of Hawaii. Early samples from depths of more than 700 m showed salinities nearly equivalent to seawater but having anomalous cation concentrations that are attributed to ion exchange between formation fluids and residual drilling mud clays. Later deep samples found only minor variations from seawater cation chemistry that are consistent with low‐temperature weathering of basalts; δ 18 O values are equivalent to seawater values and are consistent with this interpretation. Carbon 14 activities of dissolved inorganic carbonate indicate a water age ranging from 5890 to 7170 years B.P. and fluid transport rates of 1.8 to 2.2 m/yr. Fluid samples from perforations at 310 m in the borehole demonstrate that a freshwater aquifer is present at the Mauna Kea/Mauna Loa interface; borehole resistivity logs indicate that it is ∼200 m thick. Although it has not yet been possible to obtain samples of the freshwater zone without contamination from the deep saline fluids, the chloride concentrations of the low‐salinity zone are estimated using a mixing enthalpy calculation to be less than 100 mg/L. Light stable isotope data indicate that the fresh water at 320 m is derived from recharge entering the island at an average elevation of 2000 m. Inferred 14 C activities of the dissolved bicarbonate in the freshwater zone indicate an average calibrated age of 2200 years B.P. and an average fluid velocity of at least 14 m/yr. A regional water flow model is proposed that suggests that the fresh water found at the 320‐m depth is derived from rainfall recharge from the middle elevations of Mauna Kea volcano. This rainfall is channeled beneath the Mauna Loa lavas by the thick soil layer separating the two volcanoes. A second shallow fresh‐to‐brackish water zone, derived from Mauna Loa recharge, is also inferred to exist below the carbonate formation that underlies the shallow basal lens. The results of our preliminary study of the groundwater system below the KP‐1 drill site demonstrate that intervolcano and interflow aquicludes can have a substantial impact on water circulation and discharge from young island volcanoes.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Analysis of fractures intersecting Kahi Puka Well 1 and its relation to the growth of the island of Hawaii

As part of the Hawaii Scientific Drilling Project, Kahi Puka Well 1 penetrated about 275 m of Mauna Loa basalts overlying a sequence of Mauna Kea flow units as it was drilled and cored to a total depth of 1053 m below land surface. A borehole televiewer (BHTV) was run in most of the well in successive stages prior to casing in order to obtain magnetically oriented acoustic images of the borehole wall. A total of 283 individual fractures were identified from this log and characterized in terms of strike and dip. These data are divided into three vertical sections based upon age and volcanic source, and lower hemisphere stereographic plots identify two predominant, subparallel fracture subsets common to each section. Assuming that most of the steeply dipping fractures observed in the BHTV log are tensile features generated within basalt flows during deposition and cooling, this fracture information can be combined with models of the evolution of the island of Hawaii to investigate the depositional history of these Mauna Loa and Mauna Kea basalts over the past 400 kyr. The directions of high-angle fractures appear to be generally parallel to topography or to the coastline at the time of deposition, as is supported by surface mapping of modern flows. Consequently, an overall counterclockwise rotation of about 75?? in the strike of these fractures from the bottom to the top of the well represents a systematic change in depositional slope direction over time. We attribute the observed rotation in the orientations of the two predominant fracture subsets over the past 400 kyr to changes in the configurations of volcanic sources during shield building and to the structural interference of adjacent volcanoes that produces shifts in topographic patterns.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

The effects of snowpack grain size on satellite passive microwave observations from the Upper Colorado River Basin

Understanding the passive microwave emissions of a snowpack, as observed by satellite sensors, requires knowledge of the snowpack properties: water equivalent, grain size, density, and stratigraphy. For the snowpack in the Upper Colorado River Basin, measurements of snow depth and water equivalent are routinely available from the U.S. Department of Agriculture, but extremely limited information is available for the other properties. To provide this information, a field program from 1984 to 1995 obtained profiles of snowpack grain size, density, and temperature near the time of maximum snow accumulation, at sites distributed across the basin. A synoptic basin-wide sampling program in 1985 showed that the snowpack exhibits consistent properties across large regions. Typically, the snowpack in the Wyoming region contains large amounts of depth hoar, with grain sizes up to 5 mm, while the snowpack in Colorado and Utah is dominated by rounded snow grains less than 2 mm in diameter. In the Wyoming region, large depth hoar crystals in shallow snowpacks yield the lowest emissivities or coldest brightness temperatures observed across the entire basin. Yearly differences in the average grain sizes result primarily from variations in the relative amount of depth hoar within the snowpack. The average grain size for the Colorado and Utah regions shows much less variation than do the grain sizes from the Wyoming region. Furthermore, the greatest amounts of depth hoar occur in the Wyoming region during 1987 and 1992, years with strong El Niño Southern Oscillation, but the Colorado and Utah regions do not show this behavior.

Journal of Geophysical Research C: Oceans↗
USGS1996

Satellite imagery of the onset of streaming flow of ice streams C and D, West Antarctica

Five overlapping Landsat multispectral scanner satellite images of the interior of the West Antarctic ice sheet were enhanced with principal component analysis, high-pass filtering, and linear contrast stretching and merged into a mosaic by aligning surface features in the overlap areas. The mosaic was registered to geodetic coordinates, to an accuracy of about 1 km, using the five scene centers as control points. The onset of streaming flow of two tributaries of ice stream C and one tributary of ice stream D is visible in the mosaic. The onset appears to occur within a relatively short distance, less than the width of the ice stream, typically at a subglacial topographic feature, such as a step or ridge. The ice streams extend farther up into the interior than previously mapped. Ice stream D starts about 150 km from the ice divide, at an altitude of about 1500 m, approximately halfway up the convex-upward dome shape of the interior ice sheet. Ice stream D is relatively much longer than ice stream C, possibly because ice stream D is currently active, whereas ice stream C is currently inactive. The grounded portion of the West Antarctic ice sheet is perhaps best conceptualized as an ice sheet in which ice streams are embedded over most of its area, with slow-moving ice converging into fast-moving ice streams in a widely distributed pattern, much like that of streams and rivers in a hydrologic basin. A relic margin appears to parallel most of the south margin of the tributary of ice stream D, separated from the active shear margin by about 10 km or less for a distance of over 200 km. This means there is now evidence for recent changes having occurred in three of the five major ice streams which drain most of West Antarctica (B, C, and D), two of which (B and D) are currently active.

Journal of Geophysical Research C: Oceans↗
USGS1996

Southern Ocean monthly wave fields for austral winters 1985-1988 by Geosat radar altimeter

Four years of monthly averaged wave height fields for the austral winters 1985–1988 derived from the Geosat altimeter data show a spatial variability of the scale of 500–1000 km that varies monthly and annually. This variability is superimposed on the zonal patterns surrounding the Antarctic continent and characteristic of the climatology derived from the U.S. Navy [1992] Marine Climatic Atlas of the World. The location and the intensity of these large-scale features, which are not found in the climatological fields, exhibit strong monthly and yearly variations. A global underestimation of the climatological mean wave heights by more than 1 m is also found over large regions of the Southern Ocean. The largest monthly averaged significant wave heights are above 5 m and are found during August of every year in the Indian Ocean, south of 40°S. The monthly wave fields show more variability in the Atlantic and Pacific Oceans than in the Indian Ocean. The Seasat data from 1978 and the Geosat data from 1985 and 1988 show an eastward rotation of the largest wave heights. However, this rotation is absent in 1986 and 1987; the former was a year of unusually low sea states, and the latter was a year of unusually high sea states, which suggests a link to the El Niño-Southern Oscillation event of 1986.

Journal of Geophysical Research C: Oceans↗
USGS1996

Cadmium in the California Current system: Tracer of past and present upwelling

Over 100 samples were collected off the west coast of North America during 1991–1993 to determine the relation between wind‐driven upwelling and nearshore concentrations of dissolved silicate (Si), phosphate (P), and cadmium (Cd). Highly enriched in deep water offshore, these constituents are sensitive indicators of upwelling. Coastal water was sampled from the shore in January and June 1992 at 12 sites distributed between 36° and 48°N latitude. In January the composition of nearshore water along this transect was fairly uniform: 5–15 μmol/kg for Si, 0.5 to 1.0 μmol/kg for P, and 0.1–0.3 nmol/kg for Cd. In June, elevated concentrations of Si (30 μmol/kg), P (2.0 μmol/kg), and Cd (0.6 nmol/kg) revealed a region of intense upwelling between 38° and 40°N. The pattern is broadly consistent with meridional gradients in coastal upwelling calculated from the long‐term mean of alongshore winds compiled from ship reports. Nearshore water was also collected biweekly to monthly at two sites 3 km apart near San Francisco Bay (37.5°N) during 1991–1993. The variability seen in the time series suggests that the composition of nearshore water integrates the effect of alongshore winds over timescales of several weeks. Seasonal variations in Si (5–50 μmol/kg), P (0.5–2.5 μmol/kg), and Cd (0.1–0.8 nmol/kg) concentrations were consistent with upwelling during spring and summer. Maximum Si, P, and Cd concentrations reached in May 1991 were consistent with advection to the very nearshore region from a depth of about 300 m relative to a vertical profile at a distance of 200 km from the coast. Nearshore Si, P, and Cd concentrations were reduced relative to 1991 in 1992, and, to a lesser extent, in 1993 due to weaker upwelling linked to the warm phase of the El Niño‐Southern Oscillation. During periods of weaker upwelling or downwelling, variations in P, Si, and Cd concentrations became uncoupled. There is a good correlation between the coastal Cd time series near San Francisco Bay (37.5°N) and a second order polynomial function of the upwelling index of Bakun [1975] at 36°N, filtered with a 30‐day running mean ( r 2 = 0.71, n = 39). The index is a daily estimate of coastal upwelling calculated from 6‐hourly mean atmospheric pressure distributions at 36°N. From this function and a record of daily upwelling indices, we infer a range of annually averaged coastal Cd concentrations of at least 0.3–0.5 nmol/kg since 1967. Cd/Ca ratios in shells of foraminifera from San Francisco Bay suggest that average coastal Cd concentrations 3500–4500 years ago were at the upper end of this range.

Journal of Geophysical Research C: Oceans↗
USGS1996

Potential climatic impacts of vegetation change: A regional modeling study

The human species has been modifying the landscape long before the development of modern agrarian techniques. Much of the land area of the conterminous United States is currently used for agricultural production. In certain regions this change in vegetative cover from its natural state may have led to local climatic change. A regional climate version of the Colorado State University Regional Atmospheric Modeling System was used to assess the impact of a natural versus current vegetation distribution on the weather and climate of July 1989. The results indicate that coherent regions of substantial changes, of both positive and negative sign, in screen height temperature, humidity, wind speed, and precipitation are a possible consequence of land use change throughout the United States. The simulated changes in the screen height quantities were closely related to changes in the vegetation parameters of albedo, roughness length, leaf area index, and fractional coverage.

Journal of Geophysical Research D: Atmospheres↗
USGS1996

Modeling of fluidized ejecta emplacement over digital topography on Venus

The FLOW computer model of McEwen and Malin (1989) modified for application to the study of Venus fluidized ejecta blankets (FEBs) demonstrates that relatively low viscosities, yield strengths, and initial velocities are required to duplicate the observed flow paths of the outflow materials. The model calculates the velocities and simulated flow paths of gravity flows over Magellan topography. The model is formulated to determine flow movements from initial conditions, gravitational acceleration, and resistance to motion as described by Coulomb, viscous, and turbulent resistance forces. Successful duplication of observed FEB flow paths has been achieved for the FEB craters Addams, Isabella, and Cochran. When used as a simple energy‐line model, the model requires low coefficients of friction to extend FEBs to near their observed termini in the synthetic aperture radar (SAR) imagery, although the resulting straight flow lines do not follow the observed flow paths well. For Bingham flow, the model requires low values of viscosity and yield strength which are more similar to pyroclastic or debris flows than basaltic lavas. Flows of 100‐m depth require 1 to 2 orders of magnitude higher values of both viscosity and yield strength than 10‐m‐deep flows. The complicated nature of the flow lines for the low velocity model suggests that FEBs were probably emplaced under variably laminar and turbulent flow conditions, where underlying topography influenced both the direction and energy of flow materials.

Journal of Geophysical Research E: Planets↗
USGS2012

Compositional variations on the Moon: Recalibration of Galileo solid‐state imaging data for the Orientale region and farside

Updated radiometric calibration and systematic processing procedures for Galileo solid‐state imaging (SSI) data from the first (1990) Earth‐Moon encounter are presented. These procedures were applied to a whole‐disk imaging sequence of the Moon centered near Mare Orientale, called Lunmap 14 (L14). Processing of L14 data included radiometric calibration, subpixel coregistration, scattered light removal, geometric control and reprojection, photometric normalization, and calibration to Earth‐based spectra. Coregistration and scattered‐light removal procedures are improvements over the initial calibration of the SSI mosaics. The effects of scattered light correction are best seen using a whole‐Moon view such as L14; resolution of the debated amounts of light scattering from within or outside the camera field of view is not necessary. Scattered light removal particularly affects the 1‐μm spectral region and has implications for interpretation of mafic mineral signatures in mare deposits of the lunar limb and farside. Recalibrated spectra indicate that mare ponds of the limb show moderately deep 1‐μm absorptions, and thus mafic mineral contents, comparable to those of other nearside basalts. Mafic mineral contents of Schiller‐Shickard cryptomaria are higher than previously thought and are similar to some low‐Ti nearside basalts. Many of the recalibrated spectra from South Pole/Aitken are similar to those of Schiller‐Schickard cryptomaria, suggesting that many of these soils represent a mixed mare/highland lithology. The hypothesis that there is an olivine enrichment in southern South Pole/Aitken basin is not supported strongly by spectra shown here.

Journal of Geophysical Research: Planets↗
USGS1996

Volatile emissions from the crater and flank of Oldoinyo Lengai volcano, Tanzania

As a comparison to airborne infrared (IR) flux measurements, ground-based sampling of fumarole and soil gases was used to characterize the quiescent degassing of CO 2 from Oldoinyo Lengai volcano. Aerial and ground-based measurements are in good agreement: ∼75% of the aerially measured CO 2 flux at Lengai (0.05–0.06 × 10 12 mol yr −1 or 6000–7200 tonnes CO 2 d −1 ) can be attributed to seven large crater vents. In contrast to Etna and Vulcano Island, where 15–50% of the total CO 2 flux emanates diffusely through the volcanic flanks, diffuse emissions were measured only within 500 m of the crater rim at Lengai, contributing <2% of the total flux. The lack of extensive flank emissions may reflect the dimensions of the magma chamber and/or the lack of a shallow fluid flow system. Thermodynamic restoration of fumarole analyses shows that gases are the most CO 2 -rich and H 2 O-poor reported for any volcano, containing 64–74% CO 2 , 24–34% H 2 O, 0.88–1.0% H 2 , 0.1–0.4% CO and <0.1% H 2 S, HCl, HF, and CH 4 . Volatile emissions of S, Cl, and F at Oldoiyno Lengai are estimated as 4.5, 1.5, and 1.0 × 10 7 mol yr −1 , respectively. Accuracy of the airborne technique was also assessed by measuring the C emission rate from a coal-burning power plant. CO 2 fluxes were measured within ±10% near the plant; however, poor resolution at increased distances caused an underestimation of the flux by a factor of 2. The relatively large CO 2 fluxes measured for alkaline volcanoes such as Oldoinyo Lengai or Etna may indicate that midplate volcanoes represent a large, yet relatively unknown, natural source of CO 2 .

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Isotopic evolution of Mauna Kea volcano: Results from the initial phase of the Hawaii Scientific Drilling Project

We have examined the Sr, Nd, and Pb isotopic compositions of Mauna Kea lavas recovered by the first drilling phase of the Hawaii Scientific Drilling Project. These lavas, which range in age from ∼200 to 400 ka, provide a detailed record of chemical and isotopic changes in basalt composition during the shield/postshield transition and extend our record of Mauna Kea volcanism to a late-shield period roughly equivalent to the last ∼100 ka of Mauna Loa activity. Stratigraphic variations in isotopic composition reveal a gradual shift over time toward a more depleted source composition (e.g., higher 143 Nd/ 144 Nd, lower 87 Sr/ 86 Sr, and lower 3 He/ 4 He). This gradual evolution is in sharp contrast with the abrupt appearance of alkalic lavas at ∼240 ka recorded by the upper 50 m of Mauna Kea lavas from the core. Intercalated tholeiitic and alkalic lavas from the uppermost Mauna Kea section are isotopically indistinguishable. Combined with major element evidence (e.g., decreasing SiO 2 and increasing FeO) that the depth of melt segregation increased during the transition from tholeiitic to alkalic volcanism, the isotopic similarity of tholeiitic and alkalic lavas argues against significant lithosphere involvement during melt generation. Instead, the depleted isotopic signatures found in late shield-stage lavas are best explained by increasing the proportion of melt generated from a depleted upper mantle component entrained and heated by the rising central plume. Direct comparison of Mauna Kea and Mauna Loa lavas erupted at equivalent stages in these volcanoes' life cycles reveals persistent chemical and isotopic differences independent of the temporal evolution of each volcano. The oldest lavas recovered from the drillcore are similar to modern Kilauea lavas, but are distinct from Mauna Loa lavas. Mauna Kea lavas have higher 143 Nd/ 144 Nd and 206 Pb/ 204 Pb and lower 87 Sr/ 86 Sr. Higher concentrations of incompatible trace elements in primary magmas, lower SiO 2 , and higher FeO also indicate that Mauna Kea lavas formed through smaller degrees of partial melting at greater depth than Mauna Loa lavas. These chemical and isotopic differences are consistently found between volcanoes along the western “Loa” and eastern “Kea” trends and reflect large-scale variations in source composition and melting environment. We propose a simple model of a radially zoned plume centered beneath the Loa trend. Loa trend lavas generated from the hot plume axis reflect high degrees of partial melting from a source containing a mixture of enriched plume-source material and entrained lower mantle. Kea trend lavas, in contrast, are generated from the cooler, peripheral portions of the plume, record lower degrees of partial melting, and tap a source containing a greater proportion of depleted upper mantle.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

An 84-kyr paleomagnetic record from the sediments of Lake Baikal, Siberia

We have conducted a paleomagnetic study of sediment cores obtained from the Selenga prodelta region of Lake Baikal, Russia. This record, which spans approximately the last 84 kyr, contributes to a better understanding of the nature of geomagnetic field behavior in Siberia and is a useful correlation and dating tool. We demonstrate that the Lake Baikal sediments are recording variations in the geomagnetic field. The directional record displays secular variation behavior with a geomagnetic excursion at 20 ka and additional excursions appearing as large-amplitude secular variation at 41, 61, and 67 ka. Smoothing of the geomagnetic excursion behavior occurs in Lake Baikal sediments owing to the intermediate sedimentation rate (13 cm kyr −1 ). The Lake Baikal relative paleointensity record correlates to absolute paleointensity data for the last 10 kyr and to relative paleointensity records from the Mediterranean Sea and Indian Ocean for the last 84 kyr. This correlation suggests a strong global (i.e., dipole) component to these records and further supports the reliability of sediments as recorders of relative geomagnetic paleointensity. We show that a relative geomagnetic intensity stratigraphy has a potential resolution of 7 kyr by correlating continental and marine records. The geomagnetic intensity stratigraphy helps constrain the age of the difficult to date Lake Baikal sediments.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Frictional constraints on crustal faulting

We consider how variations in fault frictional properties affect the phenomenology of earthquake faulting. In particular, we propose that lateral variations in fault friction produce the marked heterogeneity of slip observed in large earthquakes. We model these variations using a rate- and state-dependent friction law, where we differentiate velocity-weakening behavior into two fields: the strong seismic field is very velocity weakening and the weak seismic field is slightly velocity weakening. Similarly, we differentiate velocity-strengthening behavior into two fields: the compliant field is slightly velocity strengthening and the viscous field is very velocity strengthening. The strong seismic field comprises the seismic slip concentrations, or asperities. The two “intermediate” fields, weak seismic and compliant, have frictional velocity dependences that are close to velocity neutral: these fields modulate both the tectonic loading and the dynamic rupture process. During the interseismic period, the weak seismic and compliant regions slip aseismically, while the strong seismic regions remain locked, evolving into stress concentrations that fail only in main shocks. The weak seismic areas exhibit most of the interseismic activity and aftershocks but can also creep seismically. This “mixed” frictional behavior can be obtained from a sufficiently heterogeneous distribution of the critical slip distance. The model also provides a mechanism for rupture arrest: dynamic rupture fronts decelerate as they penetrate into unloaded complaint or weak seismic areas, producing broad areas of accelerated afterslip. Aftershocks occur on both the weak seismic and compliant areas around a fault, but most of the stress is diffused through aseismic slip. Rapid afterslip on these peripheral areas can also produce aftershocks within the main shock rupture area by reloading weak fault areas that slipped in the main shock and then healed. We test this frictional model by comparing the seismicity and the coseismic slip for the 1966 Parkfield, 1979 Coyote Lake, and 1984 Morgan Hill earthquakes. The interevent seismicity and aftershocks appear to occur on fault areas outside the regions of significant slip: these regions are interpreted as either weak seismic or compliant, depending on whether or not they manifest interevent seismicity.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Frictional behavior of large displacement experimental faults

The coefficient of friction and velocity dependence of friction of initially bare surfaces and 1-mm-thick simulated fault gouges (<90 μm) of Westerly granite were determined as a function of displacement to >400 mm at 25°C and 25 MPa normal stress. Steady state negative friction velocity dependence and a steady state fault zone microstructure are achieved after ∼18 mm displacement, and an approximately constant strength is reached after a few tens of millimeters of sliding on initially bare surfaces. Simulated fault gouges show a large but systematic variation of friction, velocity dependence of friction, dilatancy, and degree of localization with displacement. At short displacement (<10 mm), simulated gouge is strong, velocity strengthening and changes in sliding velocity are accompanied by relatively large changes in dilatancy rate. With continued displacement, simulated gouges become progressively weaker and less velocity strengthening, the velocity dependence of dilatancy rate decreases, and deformation becomes localized into a narrow basal shear which at its most localized is observed to be velocity weakening. With subsequent displacement, the fault restrengthens, returns to velocity strengthening, or to velocity neutral, the velocity dependence of dilatancy rate becomes larger, and deformation becomes distributed. Correlation of friction, velocity dependence of friction and of dilatancy rate, and degree of localization at all displacements in simulated gouge suggest that all quantities are interrelated. The observations do not distinguish the independent variables but suggest that the degree of localization is controlled by the fault strength, not by the friction velocity dependence. The friction velocity dependence and velocity dependence of dilatancy rate can be used as qualitative measures of the degree of localization in simulated gouge, in agreement with previous studies. Theory equating the friction velocity dependence of simulated gouge to the sum of the friction velocity dependence of bare surfaces and the velocity dependence of dilatancy rate of simulated gouge fails to quantitatively account for the experimental observations.

Journal of Geophysical Research B: Solid Earth↗