Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysics Journal International”

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

Limitations of quantitative analysis of deep crustal seismic reflection data: Examples from GLIMPCE

Amplitude preservation in seismic reflection data can be obtained by a relative true amplitude (RTA) processing technique in which the relative strength of reflection amplitudes is preserved vertically as well as horizontally, after compensating for amplitude distortion by near-surface effects and propagation effects. Quantitative analysis of relative true amplitudes of the Great Lakes International Multidisciplinary Program on Crustal Evolution seismic data is hampered by large uncertainties in estimates of the water bottom reflection coefficient and the vertical amplitude correction and by inadequate noise suppression. Processing techniques such as deconvolution, F-K filtering, and migration significantly change the overall shape of amplitude curves and hence calculation of reflection coefficients and average reflectance. Thus lithological interpretation of deep crustal seismic data based on the absolute value of estimated reflection strength alone is meaningless. The relative strength of individual events, however, is preserved on curves generated at different stages in the processing. We suggest that qualitative comparisons of relative strength, if used carefully, provide a meaningful measure of variations in reflectivity. Simple theoretical models indicate that peg-leg multiples rather than water bottom multiples are the most severe source of noise contamination. These multiples are extremely difficult to remove when the water bottom reflection coefficient is large (>0.6), a condition that exists beneath parts of Lake Superior and most of Lake Huron.

Great Lakes↗

Detection, attribution, and sensitivity of trends toward earlier streamflow in the Sierra Nevada

Observed changes in the timing of snowmelt dominated streamflow in the western United States are often linked to anthropogenic or other external causes. We assess whether observed streamflow timing changes can be statistically attributed to external forcing, or whether they still lie within the bounds of natural (internal) variability for four large Sierra Nevada (CA) basins, at inflow points to major reservoirs. Streamflow timing is measured by "center timing" (CT), the day when half the annual flow has passed a given point. We use a physically based hydrology model driven by meteorological input from a global climate model to quantify the natural variability in CT trends. Estimated 50-year trends in CT due to natural climate variability often exceed estimated actual CT trends from 1950 to 1999. Thus, although observed trends in CT to date may be statistically significant, they cannot yet be statistically attributed to external influences on climate. We estimate that projected CT changes at the four major reservoir inflows will, with 90% confidence, exceed those from natural variability within 1-4 decades or 4-8 decades, depending on rates of future greenhouse gas emissions. To identify areas most likely to exhibit CT changes in response to rising temperatures, we calculate changes in CT under temperature increases from 1 to 5??. We find that areas with average winter temperatures between -2??C and -4??C are most likely to respond with significant CT shifts. Correspondingly, elevations from 2000 to 2800 in are most sensitive to temperature increases, with CT changes exceeding 45 days (earlier) relative to 1961-1990. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Generation and propagation of nonlinear internal waves in Massachusetts Bay

During the summer, nonlinear internal waves (NLIWs) are commonly observed propagating in Massachusetts Bay. The topography of the area is unique in the sense that the generation area (over Stellwagen Bank) is only 25 km away from the shoaling area, and thus it represents an excellent natural laboratory to study the life cycle of NLIWs. To assist in the interpretation of the data collected during the 1998 Massachusetts Bay Internal Wave Experiment (MBIWE98), a fully nonlinear and nonhydrostatic model covering the generation/shoaling region was developed, to investigate the response of the system to the range of background and driving conditions observed. Simplified models were also used to elucidate the role of nonlinearity and dispersion in shaping the NLIW field. This paper concentrates on the generation process and the subsequent evolution in the basin. The model was found to reproduce well the range of propagation characteristics observed (arrival time, propagation speed, amplitude), and provided a coherent framework to interpret the observations. Comparison with a fully nonlinear hydrostatic model shows that during the generation and initial evolution of the waves as they move away from Stellwagen Bank, dispersive effects play a negligible role. Thus the problem can be well understood considering the geometry of the characteristics along which the Riemann invariants of the hydrostatic problem propagate. Dispersion plays a role only during the evolution of the undular bore in the middle of Stellwagen Basin. The consequences for modeling NLIWs within hydrostatic models are briefly discussed at the end.

Massachusetts↗

Modeling effects of bank friction and woody bank vegetation on channel flow and boundary shear stress in the Rio Puerco, New Mexico

[1] We have applied a physically based model for steady, horizontally uniform flow to calculate reach-averaged velocity and boundary shear-stress distributions in a natural stream with woody vegetation on the channel banks. The model calculates explicitly the form drag on woody plant stems and includes the effects of vegetation on the boundary shear stress, velocity, and turbulence fields. Average channel shapes, bed gradients, and shrub characteristics were determined for four long, internally similar reaches covering an 81 -km segment of the lower Rio Puerco, in north-central New Mexico. Hydraulic geometries for each of three flow events were determined from reach-average channel shapes, slopes, and measured silt-line heights above the bed. Model results show that friction on the lateral boundaries reduced the boundary shear stress in the center of the Rio Puerco channel by as much as 20%. In reaches with moderate to dense bank shrubs, perimeter-averaged boundary shear stresses were reduced by almost 40% during near bankfull flows. Model-calculated discharges for all three flow events indicate there was a loss of about 40% of the upstream discharge through the 81-km river segment. Infiltration rates estimated from discharge losses during the quasi-steady peak flows and areas of the bed and banks suggest rates of flow loss into the silty sand of the upper banks are about 8 times greater than rates of loss into the clayey silt-covered bed and lower banks.

Journal of Geophysical Research F: Earth Surface↗

Evidence for fluctuating wind in shaping an ancient Martian dune field: The Stimson formation at the Greenheugh pediment, Gale crater

Temporal fluctuations of wind strength and direction can influence aeolian bedform morphology and orientation, which can be encoded into the architecture of aeolian deposits. These strata represent a direct record of atmospheric processes and can be used to understand ancient Martian atmospheric processes as well as those on Earth. The strata can: give insight to ancient atmospheric circulation, how the atmosphere evolved in response to global changes in habitability, and how ancient processes differ from modern processes. The Stimson formation at the Greenheugh pediment (Gale crater) records evidence of fluctuating wind across multiple temporal scales. The strata can be subdivided into three intervals–Gleann Beag, Ladder, and Edinburgh intervals. Internally, the intervals record changes of dune morphology and orientation, correlatable to wind fluctuations at multiple temporal scales. The basal Gleann Beag interval comprises compound cross-strata, deposited by oblique compound dunes. These dunes record a bimodal wind regime, resulting in net sediment transport toward the north. The Ladder interval records a reversal of sediment transport to the south, where straight-crested simple-dunes shaped by a seasonally variable winds formed. Finally, the Edinburgh interval records sediment transport to the west, where a unimodal wind formed sinuous-crested simple dunes. These observations demonstrate active and variable atmospheric circulation in Gale crater during the accumulation of the Stimson dune field, at multiple temporal scales from seasonally driven winds to much longer time-frames, during the Hesperian. These observations can be used to further understand ancient atmospheric conditions and processes, at a high temporal resolution on Mars.

Journal of Geophysical Research - Planets↗

Block and shear-zone architecture of the Minnesota River Valley subprovince: Implications for late Archean accretionary tectonics

The Minnesota River Valley subprovince of the Superior Province is an Archean gneiss terrane composed internally of four crustal blocks bounded by three zones of east-northeast-trending linear geophysical anomalies. Two of the block-bounding zones are verified regional-scale shears. The geological nature of the third boundary has not been established. Potential-field geophysical models portray the boundary zones as moderately north-dipping surfaces or thin slabs similar in strike and dip to the Morris fault segment of the Great Lakes tectonic zone at the north margin of the subprovince. The central two blocks of the subprovince (Morton and Montevideo) are predominantly high-grade quartzofeldspathic gneiss, some as old as 3.6 Ga, and late-tectonic granite. The northern and southern blocks (Benson and Jeffers, respectively) are judged to contain less gneiss than the central blocks and a larger diversity of syntectonic and late-tectonic plutons. A belt of moderately metamorphosed mafic and ultramafic rocks having some attributes of a dismembered ophiolite is partly within the boundary zone between the Morton and Montevideo blocks. This and the other block boundaries are interpreted as late Archean structures that were reactivated in the Early Proterozoic. The Minnesota River Valley subprovince is interpreted as a late accretionary addition to the Superior Province. Because it was continental crust, it was not subductible when it impinged on the convergent southern margin of the Superior Craton in late Archean time, and it may have accommodated to convergent-margin stresses by dividing into blocks and shear zones capable of independent movement.

Minnesota↗

Prodigious submarine landslides on the Hawaiian Ridge

The extensive area covered by major submarine mass wasting deposits on or near the Hawaiian Ridge has been delimited by systematic mapping of the Hawaiian exclusive economic zone using the side‐looking sonar system GLORIA. These surveys show that slumps and debris avalanche deposits are exposed over about 100,000 km 2 of the ridge and adjacent seafloor from Kauai to Hawaii, covering an area more than 5 times the land area of the islands. Some of the individual debris avalanches are more than 200 km long and about 5000 km 3 in volume, ranking them among the largest on Earth. The slope failures that produce these deposits begin early in the history of individual volcanoes when they are small submarine seamounts, culminate near the end of subaerial shield building, and apparently continue long after dormancy. Consequently, landslide debris is an important element in the internal structure of the volcanoes. The dynamic behavior of the volcanoes can be modulated by slope failure, and the structural features of the landslides are related to elements of the volcanoes including rift zones and fault systems. The landslides are of two general types, slumps and debris avalanches. The slumps are slow moving, wide (up to 110 km), and thick (about 10 km) with transverse blocky ridges and steep toes. The debris avalanches are fast moving, long (up to 230 km) compared to width, and thinner (0.05–2 km); they commonly have a well‐defined amphitheater at their head and hummocky terrain in the lower part. Oceanic disturbance caused by rapid emplacement of debris avalanches may have produced high‐level wave deposits (such as the 365‐m elevation Hulopoe Gravel on Lanai) that are found on several islands. Most present‐day submarine canyons were originally carved subaerially in the upper parts of debris avalanches. Subaerial canyon cutting was apparently promoted by the recently steepened and stripped slopes of the landslide amphitheaters.

Journal of Geophysical Research B: Solid Earth↗

Imaging the Ferron Member of the Mancos Shale formation using reprocessed high-resolution 2-D seismic reflection data: Emery County, Utah

Late in 1982 and early in 1983, Arco Exploration contracted with Rocky Mountain Geophysical to acquired four high-resolution 2-D multichannel seismic reflection lines in Emery County, Utah. The primary goal in acquiring this data was an attempt to image the Ferron Member of the Upper Cretaceous Mancos Shale. Design of the high-resolution 2-D seismic reflection data acquisition used both a short geophone group interval and a short sample interval. An explosive energy source was used which provided an input pulse with broad frequency content and higher frequencies than typical non-explosive Vibroseis?? sources. Reflections produced by using this high-frequency energy source when sampled at a short interval are usually able to resolve shallow horizons that are relatively thin compared to those that can be resolved using more typical oil and gas exploration seismic reflection methods.The U.S. Geological Survey-Energy Resources Program, Geophysical Processing Group used the processing sequence originally applied by Arco in 1984 as a guide and experimented with processing steps applied in a different order using slightly different parameters in an effort to improve imaging the Ferron Member horizon. As with the Arco processed data there are sections along all four seismic lines where the data quality cannot be improved upon, and in fact the data quality is so poor that the Ferron horizon cannot be imaged at all.Interpretation of the seismic and core hole data indicates that the Ferron Member in the study area represent a deltaic sequence including delta front, lower delta plain, and upper delta plain environments. Correlating the depositional environments for the Ferron Member as indicated in the core holes with the thickness of Ferron Member suggests the presence of a delta lobe running from the northwest to the southeast through the study area. The presence of a deltaic channel system within the delta lobe complex might prove to be an interesting conventional exploration target along with the coal-bed methane production already proven in the area. ?? 2003 Elsevier B.V. All rights reserved.

International Journal of Coal Geology↗

Are the stress drops of small earthquakes good predictors of the stress drops of moderate-to-large earthquakes?

The stress drops of small earthquakes often exhibit spatial patterns of variability. If moderate and large earthquakes follow the same spatial patterns, the stress drops of possible future damaging earthquakes could be better predicted by considering the stress drops of nearby small events. Better stress drop predictability could reduce ground-motion uncertainty in Probabilistic Seismic Hazard Assessment (PSHA) and Earthquake Early Warning (EEW). I find that for an internally consistent stress drop catalog of M1.8-3.1 events in southern California, the stress drops of the bigger earthquakes are predictable from the nearby smaller events. However, this catalog only weakly spatially correlates with another catalog of M3.0-5.8 earthquakes, and is spatially uncorrelated with five other stress drop catalogs of M≥3.4 earthquakes. For southern California events M5.5-7.5, stress drops compiled from the literature are weakly spatially correlated with the stress drops of the M1.8-3.1 events, although the correlations are not statistically significant. The lack of strong spatial correlation may be due to actual differences in the controlling factors of stress drop, for example dynamic weakening in moderate-to-large earthquakes. Alternatively, a stronger spatial correlation may exist that is obscured by methodological heterogeneity and large errors in the stress drop estimates. Either way, the stress drops of small earthquakes do not appear to be good predictors of the stress drops of nearby moderate-to-large earthquakes, at least for current techniques of stress drop estimation. If these results are representative, small-earthquake stress drops are not currently useful for substantially reducing uncertainty in PSHA and EEW.

California↗

San Andreas fault zone drilling project: scientific objectives and technological challenges

We are leading a new international initiative to conduct scientific drilling within the San Andreas fault zone at depths of up to 10 km. This project is motivated by the need to understand the physical and chemical processes operating within the fault zone and to answer fundamental questions about earthquake generation along major plate-boundary faults. Through a comprehensive program of coring, fluid sampling, downhole measurements, laboratory experimentation, and long-term monitoring, we hope to obtain critical information on the structure, composition, mechanical behavior and physical state of the San Andreas fault system at depths comparable to the nucleation zones of great earthquakes. The drilling, sampling and observational requirements needed to ensure the success of this project are stringent. These include: 1) drilling stable vertical holes to depths of about 9 km in fractured rock at temperatures of up to 300°C; 2) continuous coring and completion of inclined holes branched off these vertical boreholes to intersect the fault at depths of 3, 6, and 9 km; 3) conducting sophisticated borehole geophysical measurements and fluid/rock sampling at high temperatures and pressures; and 4) instrumenting some or all of these inclined core holes for continuous monitoring of earthquake activity, fluid pressure, deformation and other parameters for periods of up to several decades. For all of these tasks, because of the overpressured clay-rich formations anticipated within the fault zone at depth, we expect to encounter difficult drilling, coring and hole-completion conditions in the region of greatest scientific interest.

California↗

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↗

The strength and rheology of methane clathrate hydrate

Methane clathrate hydrate (structure I) is found to be very strong, based on laboratory triaxial deformation experiments we have carried out on samples of synthetic, high-purity, polycrystalline material. Samples were deformed in compressional creep tests (i.e., constant applied stress, ??), at conditions of confining pressure P = 50 and 100 MPa, strain rate 4.5 ?? 10-8 ??? ?? ??? 4.3 ?? 10-4 s-1, temperature 260 ??? T ??? 287 K, and internal methane pressure 10 ??? PCH4 ??? 15 MPa. At steady state, typically reached in a few percent strain, methane hydrate exhibited strength that was far higher than expected on the basis of published work. In terms of the standard high-temperature creep law, ?? = A??ne-(E*+PV*)/RT the rheology is described by the constants A = 108.55 MPa-n s-1, n = 2.2, E* = 90,000 J mol-1, and V* = 19 cm3 mol-1. For comparison at temperatures just below the ice point, methane hydrate at a given strain rate is over 20 times stronger than ice, and the contrast increases at lower temperatures. The possible occurrence of syntectonic dissociation of methane hydrate to methane plus free water in these experiments suggests that the high strength measured here may be only a lower bound. On Earth, high strength in hydrate-bearing formations implies higher energy release upon decomposition and subsequent failure. In the outer solar system, if Titan has a 100-km-thick near-surface layer of high-strength, low-thermal conductivity methane hydrate as has been suggested, its interior is likely to be considerably warmer than previously expected.

Journal of Geophysical Research B: Solid Earth↗

Antarctica: A Keystone in a Changing World--Online Proceedings for the Tenth International Symposium on Antarctic Earth Sciences. Santa Barbara, California, U.S.A.--August 26 to September 1, 2007

Overview: The International Symposium on Antarctic Earth Sciences (ISAES) is held once every four years to provide an international forum for presenting research results and new ideas and for planning future Antarctic geoscience research projects. This Tenth ISAES coincides with the International Polar Year (IPY; 50th Anniversary of the International Geophysical Year) and has been structured to showcase the great breadth of geoscience research being done in Antarctic regions by more than more than 100 institutions located in over 30 countries. The science program of the Symposium encompasses six broad themes that cover key topics on evolution and interactions of the geosphere, cryosphere and biosphere and their cross-linkages with past and historic paleoclimates. Emphasis is also on deciphering the climate records in ice cores, geologic cores, rock outcrops and those inferred from climate models. New technologies for the coming decades of geoscience data collection are also highlighted. Ten keynote presentations at the symposium outline the foundation for the research sessions of the symposium and the structure of the Online Proceedings and Proceedings Book for the Tenth ISAES. The ISAES is traditionally a cornerstone meeting for the Scientific Committee on Antarctic Research (SCAR). In recognition of the Tenth ISAES being held in the U.S. for the first time in 30 years and during IPY, the publication of the symposium proceedings is being handled as a special collaborative effort of the U.S. National Science Foundation, the U.S. Geological Survey, The National Academies Polar Research Board and The National Academies Press. The National Academies Polar Research Board oversees the activities of SCAR in the U.S. Special attention has been directed at publication formats for the symposium, to expedite the open and wide sharing of mature and preliminary research results presented in talks and posters at the Tenth ISAES. All symposium presentations are documented by a short Summary printed in the Symposium Program Booklet and either a short research paper or an extended abstract. The research papers and extended abstracts are compiled in this Online Proceedings and are replicated on a DVD-ROM that is placed in the back of the Tenth ISAES Proceedings Book. The Proceedings Book has printed versions of the keynote talks and an overview paper of the symposium. The short research papers and extended abstracts have been handled differently. Research papers present mature research results and syntheses. They have been peer-reviewed using standard journal procedures. Following revisions and acceptances by co-editors, the papers have been formatted for publication and proofread by authors. Each paper has been assigned a Digital Object Identification (DOI) number and separate html link, and posted online (as part of this USGS Open-File Report series) to ensure open and wide access to the research results. Extended abstracts focus on preliminary research results and have not been peer reviewed. They have had only minimal editorial review and revision. Authors have formatted and proofread their papers. Extended abstracts were not given DOI numbers and are included together in a separate chapter of this online Proceedings. USGS publications staff and Stanford-student editorial assistants indexed and compiled the PDF versions of the short research papers and extended abstracts for inclusion in this online Proceedings. USGS staff created the master DVD-ROM that contains a replica of the Online Proceedings for the Tenth ISAES, and provided the DVD-ROM copies that are included in the Tenth ISAES Proceedings Book published by The National Academies Press in the U.S. A team of more than 25 co-editors coordinated with the numerous authors and peer reviewers in handling the many research papers and extended abstracts that are included herein. Handling the large volume of short papers and extended abstracts, getting most of them onli

Open-File Report↗

Earthquake-derived seismic velocity changes during the 2018 caldera collapse of Kīlauea volcano

The 2018 Kīlauea caldera collapse produced extraordinary sequences of seismicity and deformation, with 62 episodic collapse events which significantly altered the landscape of the summit region. Despite decades of focused scientific studies at Kīlauea, detailed information about the internal structure of the volcano is limited. Recently developed techniques in seismic interferometry can be used to monitor the internal structure of an active volcano more directly by detecting subtle spatiotemporal changes in seismic wave velocity, but their utility relies on accurate interpretations of the underlying phenomena causing those velocity changes. Here, we retrospectively apply repeating-earthquake-based seismic interferometry to the 2018 Kīlauea eruption sequence. We find that seismic velocities changed over two distinct time scales: a sudden increase followed by a slower decrease in velocity in the hours following each collapse event, and a gradual, long-term decrease in velocity over several weeks that ceased approximately 1 month prior to the end of the eruption. Modeling suggests that short-term changes can be explained by magma reservoir pressurization which specifically closed vertical ring fractures. Long-term changes are related to subsidence of the caldera and likely include the influence of inelastic strain from the formation of new fractures. These observations provide new insights into the evolution of Kīlauea during its progressive collapse and will inform future interpretations for near-real-time monitoring at hazardous volcanoes around the world using similar techniques, especially where a dominant fracture orientation is present.

Hawaii↗

General circulation model simulations of winter and summer sea-level pressures over North America

Temporal variability of synoptic-scale circulation is a dominant factor in regional temporal climate variations. In the midlatitudes, temporal variability of synoptic-scale circulation has been found to be related to temporal variations in regional air temperature and precipitation. General circulation model (GCM) simulations of long-term changes in synoptic-scale circulation may be useful in assessing the effects of climate change on regional air temperature and precipitation. Atmospheric pressure is an important component of synoptic-scale circulation. In this paper, observed sea-level pressures were used to evaluate winter and summer sea-level pressures over North America simulated by the Goddard Institute for Space Studies (GISS) and the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation models. The objective of the study is to determine how similar the spatial and temporal distributions of GCM-simulated daily sea-level pressures over North America are to observed distributions. Overall, both models are better at reproducing observed within-season variance of winter and summer sea-level pressures than they are at simulating the magnitude of mean winter and summer sea-level pressures. Results indicate that winter sea-level pressures are better simulated by the GISS model than by the GFDL model. The GFDL model simulates anomalously high winter sea-level pressures over north-eastern North America. In summer, the GISS model underestimates sea-level pressures, whereas the GFDL model overestimates sea-level pressures over most of North America. Kolmogorov-Smirnov tests indicate that much of the error in model simulations can be attributed to differences between simulated and observed means. Comparisons of spatial correlation matrices of observed and simulated sea-level pressures indicate that GCMs simulate spatial variability of sea-level pressures that is similar to observed. These results suggest that GCMs are able to simulate major components of synoptic-scale circulation over North America, especially for winter.

International Journal of Climatology↗

Borehole techniques identifying subsurface chimney heights in loose ground-some experiences above underground nuclear explosions

The location of the subsurface top of the chimney formed by the collapse of the cavity resulting from an underground nuclear explosion is examined at five sites at the Nevada Test Site. The chimneys were investigated by drilling, coring, geophysical logging (density, gamma-ray, caliper), and seismic velocity surveys. The identification of the top of the chimney can be complicated by chimney termination in friable volcanic rock of relatively high porosity. The presence of an apical void in three of the five cases is confirmed as the chimney horizon by coincidence with anomalies observed in coring, caliper and gamma-ray logging (two cases), seismic velocity, and drilling. In the two cases where an apical void is not present, several of these techniques yield anomalies at identical horizons, however, the exact depth of chimney penetration is subject to some degree of uncertainty. This is due chiefly to the extent to which core recovery and seismic velocity may be affected by perturbations in the tuff above the chimney due to the explosion and collapse. The data suggest, however, that the depth uncertainty may be only of the order of 10 m if several indicators are available. Of all indicators, core recovery and seismic velocity indicate anomalous horizons in every case. Because radiation products associated with the explosion are contained within the immediate vicinity of the cavity, gamma-ray logs are generally not diagnostic of chimney penetration. In no case is the denisty log indicative of the presence of the chimney. ?? 1993.

International Journal of Rock Mechanics and Mining↗

Structural control on megathrust rupture and slip behavior: Insights from the 2016 Mw 7.8 Pedernales Ecuador earthquake

The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 Mw 7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906 Mw 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a Mw 7.7 earthquake.

Journal of Geophysical Research: Earth Surface↗

Travel times of P and S from the global digital seismic networks: Implications for the relative variation of P and S velocity in the mantle

We present new data sets of P and S arrival times which have been handpicked from long-period vertical and transverse component recordings of the various global seismic networks. Using events which occurred from 1976 to 1994 results in ∼38,000 globally well-distributed measurements of teleseismic P and ∼41,000 measurements of S . These data are particularly useful for looking at the relative variation of S and P velocities in the lower mantle. We describe both the measurement techniques and the gross characteristics of the data sets. The size of our data sets allows us to exploit the internal consistency of the data to identify outliers using a summary ray analysis. Since the polarity of each arrival is also known, we can construct fault plane solutions and/or compare with polarities predicted by the Harvard centroid moment tensor solutions to further diagnose phase misidentification. This analysis results in ∼5% of the data being identified as outliers. An analysis of variance indicates that the S residual travel times are dominated by the effects of three-dimensional structure but the P data have comparable contributions from noise and source mislocation effects. The summary ray analysis reveals the basic character of lower mantle structure, and there are large-scale patterns in both the S and P data sets that correlate quite well with each other. This analysis suggests that on average, d ln v s / d ln v p is an increasing function of depth in the mantle going from a value of ∼1.7 at the top of the lower mantle to an apparent value of 4 near the base of the mantle. This latter extreme value of R seems to result mainly from data which sample one region in the lowermost mantle under the central Pacific, where large positive S residuals are associated with very small P residuals. Such an anomaly cannot be thermal in origin.

Journal of Geophysical Research B: Solid Earth↗