S-wave travel times for a spherically averaged earth

Summary. Shear-wave travel times in a spherically averaged earth are estimated using 'differential' S minus P ( S – P ) travel-time measurements and detailed statistical procedures. Fourteen earthquakes and 48 stations are specially selected, yielding 302 S - P times for 6° < Δ < 111°. Analysis of variance techniques are used to estimate simultaneously azimuthally varying source and station adjustments while constructing an S – P travel-time model. A method of weighting the equations of condition based on the distribution of stations and epicentres is developed to reduce the effects of systematic errors due to non-random sampling of the Earth. The resulting S - P travel times are added to the 1939 Jeffreys–Bullen and the 1968 Herrin P travel times as a function of distance to obtain shear-wave travel-time models. Confidence intervals for the models are estimated from the variance of the observed S – P travel times.
The standard error for a single observed S – P travel time (6° < Δ < 111°) is 2.1 s and the residual distribution is not significantly different from a normal distribution at the 95 per cent confidence level. For 30° < Δ < 80° the mean S travel time is 1.3 s later than the corresponding mean for Jeffreys–Bullen tables, which is significant at the 95 per cent confidence level.     

Seismic travel times for central Asian epicentres

b
Observations of P and S from deep earthquakes in the Hindu Kush region are analysed. They are reduced to a surface focus by a method of S. Mohorovičić. Cubic forms in Δ are found to fit the times of P satisfactorily up to 19|Mo and of S to 21|Mo. The coefficient of the cube term is less for S than for P , indicating that the S rays are much straighter than the P ones.
A cubic from 21|Mo to 37|Mo and a quadratic from 37|Mo to 95|Mo fitted the times of P , with t and dt/dΔ continuous. Beyond 92|Mo it seems that the time curve becomes nearly straight, as had been suspected previously. For S cubics from 22.5|Mo to 37|Mo and from 37|Mo to 94|Mo fitted; there are few data beyond 95|Mo and a linear extrapolation was made to 105|Mo.
For both P and S there are considerable drops in dt/dΔ about distance 20|Mo. Analysis for P at 1|Mo intervals showed that the times beyond 20|Mo vary smoothly with distance, and that the 20|Mo discontinuity cannot be smoothed over more than 2|Mo. The large amplitudes in this range strongly suggest that the variation of dt/dΔ is continuous and not due to the entry of an overtaking branch. The data for S are less clear but would suggest a similar strong d²dΔ² about 22|Mo.
There is evidence that the times of P , extrapolated to distance 0|Mo, are about 0.6 s more than those found by Willmore in the Heligoland explosion. The intervals between P and its near reflexion agree with this, and could be interpreted as an effect of the height of the region. If so it may be desirable at some stage to adapt the results to sea level.     

Detection of the 11 yr sunspot cycle signal in Earth rotation

Summary. A two-channel MESA or maximum entropy spectrum analysis (Morf et al. 1978) between Δ T °= ET – UT (Morrison 1973) and sunspot numbers spanning 1832–1975 yields the following results over a bandwidth 2–20yr: (1) The spectra of Δ T ° and sunspot numbers are both dominated by a narrow band signal at 11.0 yr; (2) On average the coherency over the continuum is 0.14 while at 11.0yr it peaks at 0.83; (3) The 11-yr sunspot cycle signal in length of day (lod) has an amplitude of 0.16 ms, and in time lags that in sunspot numbers by 3.4yr. Estimates obtained from segmenting the series yield extremal values which grossly bound the above estimates: The narrow band signal has period range (10.5–11.4yr) coherency (0.8–1.0), amplitude (0.06–0.31 ms), and time lag (3.0–3.Syr). In addition, a two channel analysis of sunspot numbers with a new Δ T ⁿ series from 1861–1978 (Morrison 1979) and an earlier segment of Δ T °, as well as a single-channel analysis of Δ T ° and Δ T ⁿ individually, further support the conclusion that the solar sunspot cycle in Earth rotation has been detected. These experimental results have implications in astronomy, solar physics, meteorology and climatology.     

S times from deep-focus Indian earthquakes

Summary. One hundred and eight deep-focus earthquakes occurring in the Indian region (6–37° N, 70–100° E) were originally selected for the study of P residuals (Paper I). Of these, 61 earthquakes were found suitable for analysing S residuals.
S times for this region are generally greater than J—B S times. Analysis of residuals reveals that anomalous regions are present both in the upper mantle as well as the lower mantle.     

Seismic anisotropy in the core–mantle transition zone

Split S waves observed at Hockley, Texas from events in the Tonga–Fiji region of the southwest Pacific show predominantly vertically polarized shear-wave ( SV  ) energy arriving earlier than horizontally polarized ( SH ) energy for rays propagating horizontally through D" . After corrections are made for the effects of upper-mantle anisotropy beneath Hockley, a time lag of 1.5 to 2.0  s remains for the furthest events (93.9°–100.6° ), while the time lags of the nearer observations (90.5°–92.9° ) nearly disappear. At closer distances, the S waves from these same events do not penetrate as deeply into the lower mantle, and are not split. These observations suggest that a patch of D" beneath the central Pacific is anisotropic, while the mantle immediately above the patch is isotropic. The thickness of the anisotropic zone appears to be of the order of 100–200  km.
  Observations of shear-wave splitting have previously been made for paths that traverse D" under the Caribbean and under Alaska. SH leads SV , the reverse of the Hockley observations, but in these areas the fact that SV  leads SH in the HKT data shown here suggests a different sort of anisotropy under the central Pacific from that under Alaska and the Caribbean. The case of SH travelling faster than SV  is consistent with transverse isotropy with a vertical axis of symmetry (VTI) and does not require variations with azimuth. The case of SV  leading SH is consistent with transverse isotropy with a horizontal axis of symmetry (HTI), an azimuthally anisotropic medium, and with a VTI medium formed by a hexagonal crystal. Given that (Mg,Fe)SiO₃ perovskite appears unlikely to form anisotropic fabrics on a large scale, the presence of anisotropy may point to chemical heterogeneity in the lowermost mantle, possibly due to mantle–core interactions.     

Geophysical surveys on the East Pacific Rise — Galapagos Rise system

Summary. The East Pacific Rise at 12–15° S is topographically smooth with a crestal horst or linear volcanic peak marking the present axis of spreading. The Galapagos Rise at 14–17° S is topographically rough with a possible central graben marking the extinct spreading axis. The seafloor spreading magnetic anomalies on the East Pacific Rise are of low amplitude, but fracture-zone anomalies at 13–14° S have amplitudes of up to 1250 nT. Anomalies of this amplitude at the magnetic equator must be formed within the fracture zone by some combination of block reversal boundaries, anomalously-high magnetic intensities, and/or anomalously-large thicknesses of the magnetic layers within the fracture zone. Magnetization and major-element chemical analyses of basalts dredged from four locales along the fracture zone indicate that the large magnetic-anomaly amplitudes are caused by the high iron and titanium content of these ferrobasalts. The magnetic-anomaly profiles from the Galapagos Rise and its fracture-zone system are of normal amplitude and are extremely difficult to correlate internally or with the geomagnetic timescale.
Eighty-one heat-flow measurements indicate that the values measured are controlled by sediment thickness. Where the thickness of the sediment blanket is greater than 100 m, high heat flow is measured and possibly is representative of the total heat transfer at the seafloor. Where the sediment thickness is less than 100 m, seawater circulation in the oceanic crust is thought to remove most of the heat convectively; thus causing low conductive heat-flow values to be measured by the usual heat-flow apparatus. The heat loss by convective processes is probably a function also of topographic roughness and sediment permeability.     

The GH branch of PKP from deep-focus earthquakes

summary . Approximately 1000 PKP observations in the range 110° < Δ < 170° have been statistically studied to examine the existence of Bolt's GH branch of the PKP travel-time curve. Data give firm evidence for the existence of this branch from 130 to 153° but the possibility of its extending further cannot be ruled out. Observations for this branch fit the form:
Where a =(1150.3254±0.5328)s, b = (16.4907±0.3803)s, c = (−0.1316±0.0630)s.
Times for this branch are always greater by 2 s than the corresponding Gogna times from Pacific earthquakes.     

A Maastrichtian palaeomagnetic pole for the Pacific plate from a skewness analysis of marine magnetic anomaly 32

The asymmetry (skewness) of marine magnetic anomaly 32 (72.1–73.3  Ma) on the Pacific plate has been analysed in order to estimate a new palaeomagnetic pole. Apparent effective remanent inclinations of the seafloor magnetization were calculated from skewness estimates of 108 crossings of anomaly 32 distributed over the entire Pacific plate and spanning a great-circle distance of ~12  000  km. The data were inverted to obtain a palaeomagnetic pole at 72.1°N, 26.8°E with a 95 per cent confidence ellipse having a 4.0° major semi-axis oriented 98° clockwise of north and a 1.8° minor semi-axis; the anomalous skewness is 14.2° ± 3.7°. The possible dependence of anomalous skewness on spreading rate was investigated with two empirical models and found to have a negligible effect on our palaeopole analysis over the range of relevant spreading half-rates, ~25 to ~90  mm  yr⁻¹ . The new pole is consistent with the northward motion for the Pacific plate indicated by coeval palaeocolatitude and palaeoequatorial data, but differs significantly from, and lies to the northeast of, coeval seamount poles. We attribute the difference to unmodelled errors in the seamount poles, mainly in the declinations. Comparison with the northward motion inferred from dated volcanoes along the Hawaiian–Emperor seamount chain indicates 13° of southward motion of the Hawaiian hotspot since 73  Ma. When the pole is reconstructed with the Pacific plate relative to the Pacific hotspots, it differs by 14°–18° from the position of the pole relative to the Indo–Atlantic hotspots. This has several possible explanations including bias in one or more of the palaeomagnetic poles, motion between the Pacific and Indo–Atlantic hotspots, and errors in plate reconstructions relative to the hotspots.     

Travel times for Pacific earthquakes

Summary. Pacific earthquakes studied by Gogna, also three important explosions in the Tuamotu archipelago, are rediscussed. The results are very consistent, but those from Tuamotu are later than Gogna's by about 1 s in the times of P about 60°. Both sets of data give PKP residuals about -5 s about 140° - 142°, indicating that the observations there referred to the neighbourhood of the cusp of the travel-time curve but the ISS had compared them with the DEF branch. The corresponding difference in the 1940 tables is about 2 s.
Analysis at intervals of 1° indicated that the cusp of PKP is about 141° instead of 143° as in the 1940 tables and the difference between it and the DEF branch at these distances is about - 5 s.
Travel times of S under the Pacific were found but need more data, especially at distances under 10°.
Times of PcP reported by Kogan and Carder were compared with those calculated from P in Gogan's explosions, and indicated a radius of the core of 3479.8 ± 1.8 km.     

Lower-crustal rifting in the Rukwa Graben, East Africa

An M_w 5.9 earthquake occurred in the Lake Rukwa rift, Tanzania, on 1994 August 18, and was well recorded by 20 broad-band seismic stations at distances of 160 to 800 km and 21 broad-band stations at teleseismic distances. The regional and teleseismic waveforms have been used to investigate the source characteristics of the main shock, and also to locate aftershocks that occurred within three weeks of the main shock. Teleseismic body-wave modelling yields the following source parameters for the main shock: source depth of 25 ± 2 km, a normal fault orientation, with a horizontal tension axis striking NE-SW and an almost vertical pressure axis (Nodal Plane I: strike 126°–142°, dip 63°–66°, and rake 280°–290°; Nodal Plane II: strike 273°–289°, dip 28°–31°, and rake 235°–245°), a scalar moment of 4.1 times 10¹⁷ N m, and a 2 s impulsive source time function. Four of the largest aftershocks also nucleated at depths of 25 km, as deduced from regional sPmp–Pmp times. The nodal planes are broadly consistent with the orientation of both the Lupa and Ufipa faults, which bound the Rukwa rift to the northeast and southwest, respectively. The rupture radius of the main shock, assuming a circular fault, is estimated to be 4 km with a corresponding stress drop of 6.5 MPa. Published estimates of crustal thickness beneath the Rukwa rift indicate that the foci of the main shock and aftershocks lie at least 10 km above the Moho. The presence of lower-crustal seismicity beneath the Rukwa rift suggests that the pre-rift thermal structure of the rifted crust has not been strongly modified by the rifting, at least to depths of 25 km.     

Experimental study of the anisotropy of longitudinal and transverse waves from local earthquake records

Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ t_s1-_s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s^-1 and the velocities of S -waves vary from 3.15 to 3.5 km s^-1.
The delay times Δ t_S1 - _S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ t_S1 - _S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.     

DEGLACIATION OF THE EASTERN FLANK OF THE NORTH PATAGONIAN ICEFIELD AND ASSOCIATED CONTINENTAL-SCALE LAKE DIVERSIONS

ABSTRACT. We examine the deglaciation of the eastern flank of the North Patagonian Icefield between latitudes 46° and 48°S in an attempt to link the chronology of the Last Glacial Maximum moraines and those close to present-day outlet glaciers. The main features of the area are three shorelines created by ice-dammed lakes that drained eastwards to the Atlantic. On the basis of 16 ¹⁴C and exposure age dates we conclude that there was rapid glacier retreat at 15–16 ka (calendar ages) that saw glaciers retreat 90–125 km to within 20 km of their present margins. There followed a phase of glacier and lake stability at 13.6–12.8 ka. The final stage of deglaciation occurred at c. 12.8 ka, a time when the lake suddenly drained, discharging nearly 2000 km³ to the Pacific Ocean. This latter event marks the final separation of the North and South Patagonian Icefields. The timing of the onset of deglaciation and its stepped nature are similar to elsewhere in Patagonia and the northern hemisphere. However, the phase of lake stability, coinciding with the Antarctic Cold Reversal and ending during the Younger Dryas interval, mirrors climatic trends as recorded in Antarctic ice cores. The implication is that late-glacial changes in southern Patagonia were under the influence of the Antarctic realm and out of phase with those of the northern hemisphere.     

Localized velocity anomalies in the lower mantle

Summary. Two localized regions of velocity heterogeneity in the lower mantle with scale lengths of 1000–2000 km and 2 per cent velocity contrasts are detected and isolated through comparison of S, ScS, P and PcP travel times and amplitudes from deep earthquakes in Peru, Bolivia, Argentina and the Sea of Okhotsk. Comparison of the relative patterns of ScS-S differential travel times and S travel-time residuals across North American WWSSN and CSN stations for the different source regions provides baselines for interpreting which phases have anomalous times. A region of low S and P velocities is located beneath Northern Brazil and Venezuela at depths of 1700–2700 km. This region produces S -wave delays of up to 4 s for signals from deep Argentine events recorded at eastern North American stations. The localized nature of the anomaly is indicated by the narrow bounds in azimuth (15°) and take-off angle (13°) of the arrivals affected by it. The long period S -waves encountering this anomaly generally show 30–100 per cent amplitude enhancement, while the short-period amplitudes show no obvious effect. The second anomaly is a high-velocity region beneath the Caribbean originally detected by Jordan and Lynn, who used travel times from deep Peruvian events. The data from Argentine and Bolivian events presented here constrain the location of the anomaly quite well, and indicate a possible short- and long-period S -wave amplitude diminution associated with it. When the travel-time data are corrected for the estimated effects of these two anomalies, a systematic regional variation in ScS-S station residuals is apparent between stations east of and west of the Rocky Mountains. One possible explanation of this is a long wavelength lateral variation in the shear velocity structure of the lower mantle at depths greater than 2000 km beneath North America.     

Upper mantle shear structure of North America

Summary. The waveforms and travel times of S and SS phases in the range 10°–60° have been used to derive upper mantle shear velocity structures for two distinct tectonic provinces in North America. Data from earthquakes on the East Pacific Rise recorded at stations in western North America were used to derive a tectonic upper mantle model. Events on the north-west coast of North America and earthquakes off the coast of Greenland provided the data to investigate the upper mantle under the Canadian shield. All branches from the triplications due to velocity jumps near 400 and 660 km were observed in both areas. Using synthetic seismograms to model these observations placed tight constraints on heterogeneity in the upper mantle and on the details of its structure. SS–S travel-time differences of 30 s along with consistent differences in waveforms between the two data sets require substantial heterogeneity to at least 350 km depth. Velocities in the upper 170 km of the shield are about 10 per cent higher than in the tectonic area. At 250 km depth the shield velocities are still greater by about 4.5 per cent and they gradually merge near 400 km. Below 400 km no evidence for heterogeneity was found. The two models both have first-order discontinuities of 4.5 per cent at 405 km and 7.5 per cent at 695 km. Both models also have lids with lower velocities beneath. In the western model the lid is very thin and of relatively low velocity. In the shield the lid is 170 km thick with very high elocity (4.78 km s^-1); below it the velocity decreases to about 4.65 km s^-1. Aside from these features the models are relatively smooth, the major difference between them being a larger gradient in the tectonic region from 200 to 400 km.     

Gravimetric geoid in the Northwest Pacific Ocean

Summary. A total of 3708 1 × 1° free-air gravity anomaly averages have been used to construct a new 1 × 1° gravimetric geoid of the Northwest Pacific Ocean. The 1 × 1° averages are based on a compilation of 147000 surface ship and pendulum gravity measurements. The gravimetric geoid reveals information in the geoid of the Northwest Pacific not present in currently used satellite derived models. The RMS difference between the 1 × 1° geoid and satellite derived models is about ±6 m. Difference geoid undulations range from a maximum of +19 m over the Hawaiian ridge to a minimum of −31 m over the junction of the Kuril and Aleutian trenches. The Hawaiian swell is associated with a geoidal high of up to +15 m with wavelengths of about 2200 km and the topographic rises seaward of deep-sea trenches are associated with geoidal highs of up to 4m with wavelengths of about 220–900 km. The main difference between the gravimetric geoid and the satellite derived models occurs over the Pacific basin where discrepancies reach +10 m with wavelengths of 4000 km. The agreement between the gravi-metric geoid and Skylab-4 and Geos-3 altimeter data is close for wavelengths greater than about 300 km but poor for shorter wavelengths.     

Palaeomagnetism of the continental sector of the Cameroon Volcanic Line, West Africa

Measurement of samples from 154 sites in the continental sector of the Cameroon Volcanic Line yielded six palaeomagnetic poles, at 243.6°E, 84.6°N, α ₉₅ = 6.8°; 224.3°E, 81.2°N, α ₉₅ = 8.4°; 176.1°E, 82.0°N, α ₉₅ = 8.5°; 164.3°E, 86.4°N, α ₉₅ = 3.4°; 169.4°E, 82.6°N, α ₉₅ = 4.6° and 174.7°E, 72.8°N, α ₉₅ = 9.5°, belonging to rocks which have been dated by the K–Ar method at 0.4–0.9  Ma, 2.6  Ma, 6.5–11  Ma, 12–17  Ma, 20–24  Ma and 28–31  Ma, respectively. The results are in general agreement with other palaeomagnetic poles from Oligocene to Recent formations in Africa.
  The first three poles for rocks formed between 0.4 and 11  Ma are not significantly different from the present geographical pole. Together with other African poles for the same period, this suggests that the African continent has moved very little relative to the pole since 11  Ma. The other three poles for rocks dated between 12 and 31  Ma are significantly different from the present geographical pole, showing a 5° polar deviation from the present pole in the Miocene and 13° in the Middle Oligocene.     

New palaeomagnetic results from the Kerguelen Islands

New palaeomagnetic analyses have been carried out in the Kerguelen Islands on 32 lava flows of well-established age (20–22 Ma). Combined with previous studies, they yield a reliable pole for the Lower Miocene for the Antarctic plate: N = 59 flows. 349.9°E, 83.5°S, A⁹⁵= 6.1°. A reversal sequence reversed-normal has been identified in the Port Jeanne d'Arc section.     

Rayleigh-wave amplitudes from earthquakes in the range 0–150°

Summary . Vertical component Rayleigh-wave amplitudes from 1461 shallow earthquakes recorded in the distance range 0–150° are analysed to separate the effects of earthquake size, epicentral distance (Δ) and recording station.
The estimated decay of amplitude with distance has the form of a theoretical curve for the decay of Rayleigh waves with distance if the assumption is made that the decay due to dispersion for the data analysed is that of an Airy phase. Writing the decay due to anelastic attenuation as exp (- k Δ), k is estimated to be 0.676/rad over the whole range of distance. If the distance effects are represented by a straight line of the form h log Δ+ constant, h is estimated to be 1.15. The calibration function for computing M _s derived from the estimated distance effects is very similar to that of Marshall & Basham.
Station effects on Rayleigh-wave amplitudes though statistically significant are small, and can probably be ignored in the computation of M _s.
Comparing the estimated surface-wave magnitudes (earthquake size) obtained in this study with the long and short period body-wave magnitudes ( m ^LP_b and m ^SP_b respectively) obtained by Booth, Marshall & Young for the same earthquake shows that m ^LP_b is about equal to M _s over the magnitude range of interest (˜4.0–7.0). The m ^LP_b and M_s relationship shows that the greater the long-period energy radiated by an earthquake the smaller proportionately is the short-period energy.     

The S receiver functions: synthetics and data example

Recently, the S receiver function method has been successfully developed to identify upper mantle interfaces. S receiver functions have the advantage of being free of S -wave multiple reflections and can be more suitable than P receiver functions for studying mantle lithosphere. However, because of specific ray geometry and interference of diverse phases, the S receiver function method has some technical difficulties and limitations. We use synthetic seismograms to demonstrate the feasibility and limitations of S receiver functions for studying mantle structures. Full-wavefield seismograms were calculated using the reflectivity method and processed to generate synthetic S receiver functions for S , SKS and ScS waves. Results show that S receiver functions can be obtained from waveforms of S , SKS and ScS waves. The synthetic S receiver functions for these incident waves show S -to- P converted phases at all discontinuities in the crust and upper mantle. Useful ranges of epicentral distances for calculation of S receiver functions are: 55°–85° for S , >85° for SKS and 50°–75° for ScS waves. We apply both the S and P receiver function methods to data recorded at broadband station YKW3 in Northwest Canada. The study shows that there is significant agreement among different receiver function methods, and demonstrates the usefulness of S receiver functions for imaging the mantle lithosphere.     

Palaeomagnetism of Upper Cretaceous—Lower Tertiary igneous rocks from Central Argentina

Summary. From nine Upper Cretaceous—Lower Tertiary (85 ± 5–66 ± 5 Ma) volcanic hills in Central Argentina (33°S, 65°W), 26 hand samples were collected yielding a palaeomagnetic pole at 45°E 70°s ( A ₉₅ = 12.1°; k = 13.6; N = 12) after AC cleaning. Three sites show normal and nine reversed polarity. This pole is close to the pole for the late Cretaceous (69 Ma) Andacolo Series.