Coseismic slip in the 1964 Prince William Sound earthquake: A new geodetic inversion

The 1964 Prince William Sound earthquake (March 28, 1964;M _w =9.2) caused crustal deformation over an area of approximately 140,000 km² in south central Alaska. In this study geodetic and geologic measurements of this surface deformation were inverted for the slip distribution on the 1964 rupture surface. Previous seismologic, geologic, and geodetic studies of this region were used to constrain the geometry of the fault surface. In the Kodiak Island region, 28 rectangular planes (50 by 50 km each) oriented 218°N, with a dip varying from 8^o nearest the Aleutian trench to 9^o below Kodiak Island, define the rupture surface. In the Prince William Sound region 39 planes with variable dimensions (40 by 50 km near the trench, 64 by 50 km inland) and orientation (218°N in the west and 270°N in the east) were used to approximate the complex faulting. Prior information was introduced to constrain offshore dip-slip values, the strike-slip component, and slip variation between adjacent planes. Our results suggest a variable dip-slip component with local slip maximums occurring near Montague Island (up to 30 m), further to the east near Kayak Island (up to 14 m), and trenchward of the northeast segment of Kodiak Island (up to 17m). A single fault plane dipping 30°NW, corresponding to the Patton Bay fault, with a slip value of 8 m modeled the localized but large uplift on Montague Island. The moment calculated on the basis of our geodetically derived slip model of 5.0×10²⁹ dyne cm is 30% less than the seismic moment of 7.5×10²⁹ dyne cm calculated from long-period surface waves (Kanamori, 1970) but is close to the seismic moment of 5.9×10²⁹ dyne cm obtained byKikuchi andFukao (1987).     

Cooling history and differentiation of a thick North Mountain Basalt flow (Nova Scotia,Canada)

A thick (<175 m) North Mountain Basalt flow at McKay Head, Nova Scotia (Canada) shows 25-cm-thick differentiated layers separated by 130 cm of basalt in its upper 34m. Upper layers (5 m below the lava top) are highly vesicular whereas lower ones are pegmatitic and contain a thin (2 cm) rhyolite band. The layering of the flow closely resemble that of some Hawaiian lava lakes. The eesicular basalts and mafic pegmatites are inferred to be liquid-rich segregations which drained into horizontal cracks that formed within a crystalline mush. The cracks resulted from a thermal contraction associated with cooling and shrinkage of the mush. Rhyolites were formed by in situ differentiation. Gas overpressures fractured the pegmatites and gas effervescence filter pressing forced silicarich residual liquid from pegmatite interstices into the fractures creating bands. Chemical differences between the pegmatitic layers and early formed, highly differentiated upper vesicular layers may reflect a role for volatiles in the differentiation process along with crystal fractionation.     

Distribution of thermal areas on an active lava flow field: Landsat observations of Kilauea,Hawaii, July 1991

A Landsat Thematic Mapper (TM) image acquired on 23 July 1991 recorded widespread activity associated with the Episode 48 of the Pu'u 'O'o-Kupaianaha eruption of Kilauea Volcano, Hawaii. The scene contains a very large number (>3500) of thermally elevated near infrared (0.8–2.35 m) pixels (each 900 m²), which enable the spatial distribution of volcanic activity to be identified. This activity includes a lava lake within Pu'u 'O'o cone, an active lava tube system (7.9 km in length) with skylights between the Kupaianaha lava shield and several ocean entry points, and extensive active surface flows (total area of 1.3 km²) within a much larger area of cooling flows (total16 km²). The production of an average flux density map from the TM data of the flow field, wherein the average flux density is defined in units of Wm^-2, allows for the chronology of emplacement of active and cooling flows to be determined. The flux density map reveals that there were at least three breakouts (>5000 Wm^-2) feeding active flows, but on the day that the data were collected the TM recorded a waning phase of surface activity in this area, based on the relatively large amount of intermediate power-emitting (cooling) flows compared to high power-emitting (active) flows. The production of a comparable flux density map for future eruptions would aid in the assessment of volcanic hazards if the data were available in near-real time.     

Scattering attenuation and the fractal geometry of fracture systems

Scattering of seismic waves can be shown to have a frequency dependenceQ ^{–1 3–v} if scattering is produced by arrays of inhomogeneities with a 3D power spectrumW _3D(k) k ^–v. In the earth's crust and upper mantle the total attenuation is often dominated by scattering rather than intrinsic absorption, and is found to be frequency dependent according toQ ^–1 , where –1<–0.5. IfD ₁ is the fractal dimension of the surface of the 3D inhomogeneities measured on a 2D section, then this corresponds respectively to 1.5<D ₁1.75, since it can be shown that =2(D ₁–2). Laboratory results show that such a distribution of inhomogeneities, if due to microcracking, can be produced only at low stress intensities and slow crack velocities controlled by stress corrosion reactions. Thus it is likely that the earth's brittle crust is pervaded by tensile microcracks, at least partially filled by a chemically active fluid, and preferentially aligned parallel to the maximum principal compressive stress. The possibility of stress corrosion implies that microcracks may grow under conditions which are very sensitive to pre-existing heterogeneities in material constants, and hence it may be difficult in practice to separate the relative contribution of crack-induced heterogeneity from more permanent geological heterogeneities.By constrast, shear faults formed by dynamic rupture at critical stress intensities produceD ₁=1, consistent with a dynamic rupture criterion for a power law distribution of fault lengths with negative exponentD. The results presented here suggest empirically thatD ₁-1/2(D+1), thereby providing the basis for a possible framework to unify the interpretation of temporal variations in seismicb-value (b-D/2) and the frequency dependence of scattering attenuation ().This is PRIS contribution 046.     

Spatial distribution of upstream magnetospheric ≥50 keV ions

We present for the first time a statistical study of 50 keV ion events of a magnetospheric origin upstream from Earths bow shock. The statistical analysis of the 50–220 keV ion events observed by the IMP-8 spacecraft shows: (1) a dawn-dusk asymmetry in ion distributions, with most events and lower intensities upstream from the quasi-parallel pre-dawn side (4 LT-6 LT) of the bow shock, (2) highest ion fluxes upstream from the nose/dusk side of the bow shock under an almost radial interplanetary magnetic field (IMF) configuration, and (3) a positive correlation of the ion intensities with the solar wind speed and the index of geomagnetic index Kp, with an average solar wind speed as high as 620 km s^–1 and values of the index Kp 2. The statistical results are consistent with (1) preferential leakage of 50 keV magnetospheric ions from the dusk magnetopause, (2) nearly scatter free motion of 50 keV ions within the magnetosheath, and (3) final escape of magnetospheric ions from the quasi-parallel dawn side of the bow shock. An additional statistical analysis of higher energy (290–500 keV) upstream ion events also shows a dawn-dusk asymmetry in the occurrence frequency of these events, with the occurrence frequency ranging between 16%-34% in the upstream region.     

Mid-latitude ionospheric perturbation associated with the Spacelab-2 plasma depletion experiment at Millstone Hill

Elevation scans across geomagnetic mid latitudes by the incoherent scatter radar at Millstone Hill captured the ionospheric response to the firing of the Space Shuttle Challenger OMS thrusters near the peak of the F layer on July 30, 1985. Details of the excitation of airglow and the formation of an ionospheric hole during this event have been reported in an earlier paper by Mendillo et al.. The depletion (factor 2) near the 320 km Shuttle orbital altitude persisted for 35 min and then recovered to near normal levels, while at 265 km the density was reduced by a factor of 6; this significant reduction in the bottomside F-region density persisted for more than 3 hours. Total electron content in the vicinity of the hole was reduced by more than a factor of 2, and an oscillation of the F-region densities with 40-min period ensued and persisted for several hours. Plasma vertical Doppler velocity varied quasi-periodically with a 80-min period, while magnetic field variations observed on the field line through the Shuttle-burn position exhibited a similar 80-min periodicity. An interval of magnetic field variations at hydromagnetic frequencies (95 s period) accompanied the ionospheric perturbations on this field line. Radar observations revealed a downward phase progression of the 40-min period density enhancements of -1.12° km–1, corresponding to a 320-km vertical wavelength. An auroral-latitude geomagnetic disturbance began near the time of the Spacelab-2 experiment and was associated with the imposition of a strong southward IMF Bz across the magnetosphere. This created an additional complication in the interpretation of the active ionospheric experiment. It cannot be determined uniquely whether the ionospheric oscillations, which followed the Spacelab-2 experiment, were related to the active experiment or were the result of a propagating ionospheric disturbance (TID) launched by the enhanced auroral activity. The most reasonable conclusion is that the ionospheric oscillations were a result of the coincident geomagnetic disturbance. The pronounced depletion of the bottomside ionosphere, however, accentuated the oscillatory behavior during the interval following the Shuttle OMS burn.     

Far-field Gravity and Tilt Signals by Large Earthquakes: Real or Instrumental Effects?

A wide set of dynamics phenomena (i.e., Geodynamics, Post Glacial Rebound, seismicity and volcanic activity) can produce time gravity changes, which spectrum varies from short to long (more than 1 year) periods. The amplitude of the gravity variations is generally in the order of consequently their detection requires instruments with high sensitivity and stability: then, high quality experimental data. Spring and superconducting gravimeters are intensively used with this target and they are frequently jointed with tiltmeters recording stations in order to measure the elasto-gravitational perturbation of the Earth. The far-field effects produced by large earthquakes on records collected by spring gravimeters and tiltmeters are investigated here. Gravity and tilt records were analyzed on time windows spanning the occurrence of large worldwide earthquakes; the gravity records have been collected on two stations approximately 600 km distant. The background noise level at the stations was characterized, in each season, in order to detect a possible seasonal dependence and the presence of spectral components which could hide or mask other geophysical signals, such as, for instance, the highest mode of the Seismic Free Oscillation (SFO) of the Earth. Some spectral components (6.5; 8; 9; 14, 20, 51) have been detected in gravity and tilt records on the occasion of large earthquakes and the effect of the SFO has been hypothesized. A quite different spectral content of the EW and NS tiltmeter components has been detected and interpreted as a consequence of the radiation pattern of the disturbances due to the earthquakes. Through the analysis of the instrumental sensitivity, instrumental effects have been detected for gravity meters at very low frequency.F.S.E. (Fondo Sociale Europeo -European Community -)     

Stratospheric ozone variability near the IMF sector boundary

Summary We search for the effects of the interplanetary magnetic field (IMF) sector boundary crossing (SBC) in upper stratospheric ozone. The SBUV data (Nimbus-7) at the 10, 3 and 1 hPa levels are analysed for latitudes 45° N and 55° N for winters of the period December 1979 to December 1982. An effect of the IMF SBC wos only found at the 10 hPa level. These first results concerning the IMF SBC effect in upper stratospheric ozone are rather preliminary.

---

¶rt;m uu nu mau () nam aum n () a mam. SBUV ¶rt;a (u-7) a nm ¶rt;au 10, 3 u 1a aauum ¶rt; um 45° u 55° . . ua nu¶rt;a ¶rt;a 1979 – ¶rt;a 1982. m uu a¶rt; m a 10a. mu n mam n uuu a mam m n¶rt;aumu.

     

Is there any relation between the sun's motion and global seismic activity?

Summary The seismic energy released by global earthquake activity with time was correlated with F=||, where || is the absolute value of the change of the Sun's acceleration with time. For deep earthquakes and probability P=0.95, the coefficient of correlation was found to be around 0.4. For shallow earthquakes, the dependence was not proved.

---

¶rt;a u uu (1900–1982.) ua ¶rt; F=||,¶rt; || am au uu u a. mu (h>60 ) u ¶rt; mmu =0.95 uum uu ¶rt;muam0.45. nm mu auum a n¶rt;m¶rt;a.

     

Spherical harmonic analysis and Fourier transformation

Summary The paper deals with the computation of spherical harmonic coefficients from surface measurements of the magnetic or gravity field of the Earth when the measurements are distributed regularly. The Fourier representation of associated Legendre functions which this procedure makes use of, then enables the harmonic analysis to be transformed to Fourier analysis which has better numerical properties.

---

ama na ama uu uu uum n nm uu aum uuaumau n u a, ¶rt;a mu uu an ammu ma. ¶rt;mau n a¶rt; , m unm mm n¶rt;¶rt;, nm nmauu aau naam aau , m a¶rt;am uu uumu mau.

     

Microseismicity, tectonics and seismic potential in southern Caribbean and northern Venezuela

Approximately one thousand microearthquakes with body-wave magnitude mb have been located in northern Venezuela and the southern Caribbean region (9–12° N; 64–70° W) since the installation in 1980 of the Venezuelan Seismological Array, together with forty events of mb 4, one of them with surface-wave magnitude Ms 6. Focal depths are in the range of 0 to <15 km. This geologically complex region is part of the boundary between the Caribbean and the South American Plates. Epicentral locations indicate that this E–W oriented portion of the boundary is formed by two 400 km long subparallel fault zones: San Sebastián fault zone (SSF), 20 km north of Caracas along the coast; and La Victoria fault zone (LVF), 25 km south of the city. They are clearly delineated by the microseismicity. New composite focal mechanism solutions (CFMS) along these faults show right-lateral strike-slip (RLSS) motion on nearly E–W oriented fault planes. NW-striking subsidiary active faults occur in the region and intercept the two main E–W fault zones. These interceptions show high levels of microearthquake activity and seismic moment release when compared to other portions of both, the main and subsidiary faults. New CFMS at those fault crossing sites show NW-striking RLSS motion and normal faulting, in an en-echelon-like structural behavior. Geological data and quantitative comparisons with other transcurrent plate boundaries in the world suggest that the rate of plate motion in this area is on the order of 20 mm/y. Several moderate and large shocks have occurred along the SSF and LVF since 1640, including an Ms 7.6 event in 1900 on SSF. Although the region may be relatively far from a repeat of this earthquake, seismicity data indicate that strong shocks could take place along segments of the seismically active faults identified in this study.     

Predictions of magnetic fields for Uranus and Neptune

Summary The magnetic moments of Uranus and Neptune have been predicted using different scaling laws of planetary magnetism. The predictions for Uranus cover a broad band of values from very weak magnetic fields (tidal relations) to moderate fields (thermal convection hypothesis). Therefore, the direct measurements of this field by Voyager 2 (January 1986) will be very important for testing the individual hapotheses.

---

a m ama a ¶rt; amua nam ¶rt;a n aum m nam a a nm. aa n¶rt;nmu nuu mu, ¶rt;m u¶rt;am a aum n a nmu aa u u a nmu nma.

---

---

Presented at the Fifth Scientific Assembly of IAGA in Prague 1985.     

Topographic effects on tidal strains and tilts

Summary Topographic effects on tidal strains and tilts are studied using a homogeneous elastic spherical model. Expressions for local perturbing strains and tilts are derive das functions of the physical and geometrical parameters of the model. It is demonstrated that tidal tilts are affected more by the topography than tidal strains.

---

n ¶rt;¶rt; n u ¶rt;u u¶rt;a uu a mmu a nuu ¶rt;auu u a. ¶rt; au, nuau a uau nuu ¶rt;au u a auumu m uuu umuu naam ¶rt;u. aa, m uau nuu a uau nuu ¶rt;au.

     

Linear modelling in vertical ionospheric sounding

Summary Applying the methods of computing N(h) profiles to scalar product spaces provides a more general view of the differences between the individual ionospheric models, which enables a better selection of the optimum model.

---

u n¶rt; m¶rt;uu ama N(h) nu nmama a nu¶rt;u anauam u ¶rt; a au ¶rt; m¶rt;u uu ¶rt;u, m nm nmm uam nmua ¶rt;.

     

Effects of the imf sector structure in the midlatitude troposphere and stratosphere

Summary The data on geopotential heights and temperatures at 7 pressure levels between 1000-10 hPa above Berlin(52.5 °N, 13.4 °E) are analysed for the winters of 1963–1973. No demonstrable effect of the interplanetary magnetic field sector boundary crossing (IMF SBC) is found in the lower and middle stratosphere, but there is a demonstrable effect in the middle troposphere at the 500 hPa level. This effect is less important than the IMF SBC effect in the tropospheric vorticity area index and seems to be of a different type.

---

auum ¶rt;a nnmua m u mnam a 7 nm ¶rt;au ¶rt; 1000-10 a a¶rt; u(52,5 °.., 13,4 °.¶rt;.) ¶rt; u 1963–1973. ua ¶rt;aam m nu mau nam aum n( ) ¶rt;a amu u u ¶rt; mam, ma m a¶rt; ¶rt; mn a 500 a. mm m a, m u¶rt; na¶rt;u aumu am, u am m ¶rt; muna.

     

Travel times of Friuli earthquakes propagating to the Bohemian Massif

17 mmu u uma uu 1976. anua 5 ¶rt;numu mauumu u¶rt;a mu u m na Pn, Pg, Sn u Sg. u¶rt; numm muam u mum ma¶rt;am¶rt;a ¶rt; uu mmu n¶rt;naam nu m m uamm aumm. ¶rt;am nu m mmu maua.     

Normal density earth models

Summary Models of the Earth's density, close to thePREM model, have been derived, they reproduce the external normal gravitational field of the Earth and its dynamic flattening, and are referred to as normal density models. The Earth's surface is approximated by an ellipsoid of the order of the flattening, or of its square. Of the group of normal models sgtisfying the solution of the inverse problem, the normal density modelHME2 is recommended. The spherically symmetric density modelPREM, which was corrected in the course of solving the inverse problem, thus creating the modifiedPREM-E2 model, was used as the a priori information.

---

¶rt; ¶rt;u an¶rt;u nmmu uu ¶rt;uPREM (m. a. a ¶rt;u nmmu), aumau n m u¶rt;mu na¶rt;am auaumau n u. m u annuum am unu¶rt; au. uau amu a ¶rt; mam H==0.003 273 994. ma ¶rt; a ¶rt; ¶rt;m ¶rt;HME2. am anu u a ¶rt; nmmu a unaa ¶rt; a¶rt;ua umua ¶rt;PREM. ¶rt;aam ¶rt;uuau m ¶rt;u n¶rt; aauPREM-E2.

     

Iterative geophysical tomography

Summary The algorithm of iterative geophysical tomography is presented. The medium is approximated smoothly by means of B-splines. The tww-point problem of ray computation is solved with the aid of paraxial approximation. The parameters of the medium are obtained from the iterative algorithm of minimizing the quadratic form. Two numerical 2-D examples are given.

---

u¶rt; au umamuuu mauu. ¶rt;a annuuaa n nu nu -na. ma na aa a nu nu naaua annuauu. aam ¶rt; n a umamu aua uuauauu a¶rt;amu . am nu¶rt; ¶rt;a 2-D u nua.

     

Deep structure of the bohemian massif from phase velocities of Rayleigh and Love waves

Summary Phase velocities of Rayleigh and Love waves have been measured between the broadband seismic stations KHC (Kaperské Hory, South Bohemia) and KSP (Ksi, Lower Silesia), a profile that nearly coincides with the Interactional DSS Profile VII. The data for both wave types were separately inverted into models of shear-wave velocity versus depth. Novotný's model KHKS 82[1] for the DSS Profile VII was used as a start model. While the crustal section of Novotný's model is compatible with both of our data sets, our Rayleigh-wave data require smaller shear-wave velocities, on the average by 0.24 km/s, in the top 180 km of the mantle. The average difference between Novotný's model and our Love-wave model in that depth range is only 0.06 km/s. If our identification of the observed Love waves as the fundamental mode is correct, this result indicates the presence of polarization anisotropy in the uppermost mantle.

---

u a mu u a ¶rt; unu uuu mauu (an , a u) u S (, ua uu). u - S nuuum mmmm ¶rt;a¶rt; nu VII . a ¶rt; u mun m¶rt; umua ¶rt;u auumu mu nn mu. am mam¶rt;u a unaa ¶rt; S 82 m[1] ¶rt; nu VII . a am ¶rt;u m aua uu aau au ¶rt;a, ¶rt;a au ¶rt;a ¶rt; mm uu m nn — ¶rt; a 0,24 / 180 uma amuu. ¶rt;a aua ¶rt; ¶rt; m u a ¶rt; ¶rt; a m ¶rt;uanau n¶rt;mam 0,06 /. u aa u¶rt;muuau a¶rt;a a a ¶rt;ama ¶rt;a m nau, m mm mam u¶rt;ummm numuu nuau aumnuu amuu.

     

ТЕОРИЯ ВОЗНИКНОВЕНИЯ КВАЗИДВУХЛЕТНЕГО ЦИКЛА В АТМОСФЕРЕ

a aamaa ¶rt;uaua mu ¶rt;n nauuu a mu mu u ¶rt;uau uma am. m muu mam, m au¶rt;mu u au uuu am a mnuuu nau ¶rt;au. a a mu na an¶rt; aumuu mum ama, m am um m nu¶rt;u au u uu au¶rt;m ua.