A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

Experimental investigation of the elastic modulus of a fractal system—A model of fractured rocks

The elastic properties of a physical model representing a damaged rock matrix were studied using a square lattice deformed under tensile stress. The elastic modulusM of such a system varies in agreement with percolation theory as|x–x _c | ^f , wherex is the damage parameter andx _c the threshold value of the damage parameter,f3.6. Atxx _c the scale dependence ofM can be expressed asML ^–f/v , whereL is the size of the sample andv the correlation exponent in percolation theory.The experimental results are of interest in assessing elastic properties in earthquake focal zones and fault zones in general.     

Footprint prediction of scalar fluxes from analytical solutions of the diffusion equation

The use of analytical solutions of the diffusion equation for footprint prediction is explored. Quantitative information about the footprint, i.e., the upwind area most likely to affect a downwind flux measurement at a given height z, is essential when flux measurements from different platforms, particularly airborne ones, are compared. Analytical predictions are evaluated against numerical Lagrangian trajectory simulations which are detailed in a companion paper (Leclerc and Thurtell, 1990). For neutral stability, the structurally simple solutions proposed by Gash (1986) are shown to be capable of satisfactory approximation to numerical simulations over a wide range of heights, zero displacements and roughness lengths. Until more sophisticated practical solutions become available, it is suggested that apparent limitations in the validity of some assumptions underlying the Gash solutions for the case of very large surface roughness (forests) and tentative application of the solutions to cases of small thermal instability be dealt with by semi-empirical adjustment of the ratio of horizontal wind to friction velocity. An upper limit of validity of these solutions for z has yet to be established.     

Petrogenesis and paleotectonic history of the Wild Bight Group,an Ordovician rifted island arc in central Newfoundland

The Wild Bight Group (WBG) is a sequence of early and middle Ordovician volcanic, subvolcanic and epiclastic rocks, part of the Dunnage Tectonostratigraphic Zone of the Newfoundland Appalachians. A detailed geochemical and Nd-isotopic study of the volcanic and subvolcanic rocks has been carried out to determine the geochemical characteristics of the rocks, interpret their palcotectonic environments and constrain their petrogenetic history. The lower and central stratigraphic levels of the WBG contain mafic volcanic rocks with island-arc geochemical signatures, including LREE-enriched are tholeiites with _Nd(t) =-0.1 to +2.2 (type A-I), LREE-depleted arc tholeiites with _Nd(t) =+5.6 to +7.1 (type A-II) and an unusual suite of strongly incompatible-element depleted tholeiites in which _Nd(t) ranges from-0.9 to +4.6 and is negatively correlated with¹⁴⁷Sm/¹⁴⁴Nd (type A-III). High-silica, low-K rhyolites occur locally in the central part of the stratigraphy, associated with mafic rocks of arc affinity, and have _Nd(t) =+4.7 to +5.4. The upper stratigraphic levels of the WBG dominantly contain rocks with non-arc geochemical signatures, including alkalic basalts with _Nd(t) =+4.6 to +5.5 (type N-I), strongly LREE- and incompatible element-enriched tholeiites that are transitional between alkalic and non-alkalic rocks with _Nd(t) =+4.4 to +7.0 (type N-II) and rocks with flat to slightly LREE-enriched patterns and _Nd(t) =+5.1 to +7.4 (type N-III). Rocks with non-arc and arc signatures are locally interbedded near the stratigraphic type of the WBG. Nd-isotopic data in the type A-I and A-II rocks are generally compatible with mixing/partial melting models involving depleted mantle, variably contaminated by a subducted crustally-derived sediment. The petrogenesis of type A-III rocks must involve source mixing and multi-stage partial melting, but the details are not clear. The geochemistry and Nd isotope data for types N-I, N-II and N-III rocks are compatible with petrogenetic models involving variable partial melting of a source similar to that postulated for modern oceanic island basalts. Comparison of the WBG with modern analogues suggests a 3-stage developmental model: stage 1) island-arc volcanism (eruption of type mafic volcancs); stage 2) arc-rifting (continued eruption of type A-I, A-I, eruption of types A-II and A-III mafic volcanics and high-silica, low-K rhyolites); and stage 3) back-arc basin volcanism (continued minor eruption of type A-I basalts, eruption of types N-I, N-II, N-III basalts). Stages 1 and 2 volcanism involved partial melting of subduction contaminated mantle, while stage 3 volcanism utilized depleted-mantle sources not affected by the subducting slab. This model provides a basis for interpreting coeval sequences in central Newfoundland and a comparative framework for some early Paleozoic oceanic volcanic sequences elsewhere in the Appalachian orogen.     

Petrology,isotope characteristics,and K-Ar ages of the Maranhão,northern Brazil,Mesozoic basalt province

Northern Brazil contains remnants of Mesozoic flood basalts and hypabyssal rocks that were apparently emplaced during tectonism related to opening of the Atlantic Ocean. Analyses and new K-Ar ages reveal that this 700x250 km Maranhão province (5°–8°S) has low-Ti basalts (1.1 wt% TiO₂) in the western part that range about 160 to 190 Ma, and high-Ti basalts (3.4–4.4 wt% TiO₂) in the eastern part about 115–122 Ma. Low-Ti basalt compositions are less evolved and have a smaller range, Mg# 62-56, than the high-Ti basalts, Mg# 44–33. General characteristics of the least evolved members of low- and high-Ti groups include, respectively, Zr 100 and 250 ppm, Sr 225 and 475 ppm, Ba 200 and 500 ppm, Nb 10 and 26 ppm, Y 29 and 36 ppm, La/Yb_(n) 4.2 and 8.8, where La_(n) is 30 and 90. Overall compositions resemble the low- and high-Ti basaltic rocks of the Mesozoic Serra Geral (Paraná) province in southern Brazil. The Maranhão low-Ti basalts have more radiogenic Sr and Pb and higher ¹⁸O than the high-Ti basalts. Respectively, low- vs high-Ti: _Sr26–54 vs 15–18; ²⁰⁶Pb/²⁰⁴Pb=18.25–.78 vs 18.22–.24; and ¹⁸O 8.9–12.6 vs 6.5–8.6. Nd isotopes overlap: _Nd–1.6 to –3.8 vs –2.1 to –3. Ages, compositions, and isotopes indicate that the low- and high-Ti groups had independent parentages from enriched subcontinental mantle. However, both groups can be modeled from one source composition if low-Ti basalt isotopes reflect crustal contamination, and if the parentages for each group were picritic liquids that represent either higher (for low-Ti) or lower (for high-Ti) percentages of melting of that single source. When comparing Pb isotopes of Maranhão and Serra Geral high-Ti basalts (uncontaminated) to evaluate the DUPAL anomaly, Maranhão has Pb 7/4=4.6–11, and Pb 8/4=72–87; Serra Geral has Pb 7/4=10–13, and Pb 8/4=95–125. The small difference is not enough to conform to DUPAL contours, and is inconsistent with large-scale isotopic heterogeneity of mantle beneath Brazil prior to rifting of South America from Africa. Maranhão low-Ti magmas probably relate to the opening of central North Atlantic, and high-Ti magmas to the opening of equatorial Atlantic. The proposed greater percentage of source melting for low-Ti basalts may reflect a Triassic-Jurassic hotspot, while lesser melting for high-Ti magmas may relate to Cretaceous decompressional (rifting) melting.     

Carbon isotopic fractionation between CO2 vapour,silicate and carbonate melts: an experimental study to 30 kbar

The carbon isotopic fractionation between CO₂ vapour and sodamelilite (NaCaAlSi₂O₇) melt over a range of pressures and temperatures has been investigated using solid-media piston-cylinder high pressure apparatus. Ag₂C₂O₄ was the source of CO₂ and experimental oxygen fugacity was buffered at hematite-magnetite by the double capsule technique. The abundance and isotopic composition of carbon dissolved in sodamelilite (SM) glass were determined by stepped heating and the ¹³C of coexisting vapour was determined directly by capsule piercing. CO₂ solubility in SM displays a complex behavior with temperature. At pressures up to 10 kbars CO₂ dissolves in SM to form carbonate ion complexes and the solubility data suggest slight negative temperature dependence. Above 20 kbars CO₂ reacts with SM to form immiscible Na-rich silicate and Ca-rich carbonate melts and CO₂ solubility in Na-enriched silicate melt rises with increasing temperature above the liquidus. Measured values for carbon isotopic fractionation between CO₂ vapour and carbonate ions dissoived in sodamelilite melt at 1200°–1400° C and 5–30 kbars average 2.4±0.2, favouring¹³C enrichment in CO₂ vapour. The results are maxima and are independent of pressure and temperature. Similar values of 2 are obtained for the carbon isotopic fractionation between CO₂ vapour and carbonate melts at 1300°–1400° C and 20–30 kbars.     

Spectrum of cosmic fireballs

A progress report on cosmic fireballs is presented. The main new results are: (a) the phenomenon should be almost universal, and most explosive -ray sources should show the characteristic fireball spectrum; (b) even if the radiation density is insufficient, pair production in electron-proton or electron-electron scattering might start the fireball; (c) some computed fireball spectra are shown. They all have in common a 1/E low-energy behaviour, a 100 keV flattening, and a 0.5 MeV cut-off.Paper presented at the Symposium on Cosmic Gamma-Ray Bursts, held at Toulouse, France, 26–29 November, 1979.     

A Monte Carlo investigation of Jovian perturbations on short-period comet orbits

We present statistical distributions of Jovian perturbations on short-period comet orbits resulting from accurate numerical integrations. Our sample of 60, 000 cometary orbits with low inclinations and random orientations is characterized by perihelia between 0 and 7 AU and aphelia between 4 and 13 AU. The perturbations considered are those experienced because of Jupiter's gravitation per orbital revolution by the comets. Regularization and accurate step-length control in the numerical integration gives statistical results appreciably different from those computed by Rickman and Vaghi (1978). Their use of a crude method of integration led to erroneous results for close encounters. Strong asymmetries of the $\text{δ(}\text{1}\text{a}\text{)}$ distributions, in particular for the extreme tails, are observed for perihelion- or aphelion-tangent orbits. These orbits are also shown to experience the strongest energy perturbations on the average. Some results concerning the perturbations of Tisserand parameters are indicated. The perturbation distributions for the angular elements are described and discussed. The role of the minimum distance from Jupiter as an indicator of perturbations is investigated.     

Radius of convergence of Lie series for some elliptic elements

For equatorial orbits about an oblate body, we show that the Lie series for the elliptic elementse,f,l and diverge when the oblateness exceeds a critical multiple of the transformed eccentricity constant. The use of similar truncated series expansions for such elliptic elements by Brouwer accounts for the first-order errors at low eccentricity in his derived coordinates for an artificial satellite.