Dike swarms,volcanism and tectonics of the western Afar margin along the Kombolcha-Eloa traverse (Ethiopia)

The rocks outcropping along the Kombolcha-Eloa traverse (western Afar, Ethiopia) can be attributed chiefly to the so-called Oligocenic « Alaji Rhyolitic Series » (32–28 m.y. B.P.) and to the overlying « Termaber Basalts » (28–18 m.y. B.P.), these latter being emitted from a central volcano — Mt. Garfa — rising on the Afar margin and tilted during Miocene. The age and the genesis of the dikes crossing the investigated volcanic succession are discussed in the light of the radiometric ages quoted inMegrue et al. (1972), distinguishing the pre-, the syn- and the post-tilt dikes. The dikes connected with a central volcanism are fed by shallow basins, whereas the others, connected with a fissural volcanism, are to be related to strong tensional movements.     

Karroo basalts of the Tuli Syncline,Rhodesia

Karroo (Triassic-Jurassic) sediments and lavas occupy a N.N.E.-trending trough known as the Tuli Syncline 230 km in length by 45 km wide along the southern margin of Rhodesia. The structure is fault bounded, especially on the north, and comprises a sequence of gently inward dipping lava flows cut by numerous E-W and N.N.E. striking dolerite dykes. The igneous rocks have been subdivided on the basis of mineralogy and chemistry into a stratigraphic sequence of limburgites near the base, olivine basalts and upper augite-plagio-clase-phyric basalts. In addition, rare shoshonite and absarokite types occur. Twentyfour new silicate analyses are comparable to those from rocks in the neighbouring Nuanetsi Province of Rhodesia, and show strong affinities to the northern province Karroo basalts of southern Africa as defined by Cox,McDonald andHornung (1967).     

Some aspects of volcanic rocks of the valley of Mexico

Three periods of volcanic activity connected with tectonic events form the geological history of the Valley of Mexico (Mooser 1963, 1969). An igneous rock suite from rhyodacites to andesites (but lacking rhyolites and basalts) can be observed in each period. During the Tertiary epochs — in the Oligo-Miocene and Upper Miocene-Pliocene — we have a more dacitic volcanism, in the Quaternary epoch a more andesitic volcanism. This result was verified by calculating the average of all available and stratigraphically datable chemical analyses byGunn &Mooser (1971) andNegendank (1972). Using the average chemical composition of the Oligo-Miocene, Upper Miocene-Pliocene and Quaternary products the equivalent igneous rocks were computed using theRittmann-norms in theStreckeisen-Q-A-P-F double triangle with the following result (names in parenthesis are those using the classification ofMiddlemost (1973): Quaternary : quartz-latite-andesite (andesite) Upper Miocene-Pliocene : leuco-quartz-latite-andesite (high lime dacite) Oligo-Miocene : leuco-quartz-latite-andesite (high lime dacite) The equal average composition of the two groups of Tertiary volcanic rocks seems to support the theory of a uniform primary andesite magma apart from which of the two possible theories of petrogenesis one favors. The calculated average trace element abundances show high Cr- and Ni-values which suggests that mantle material was involved if we consider the Tertiary products as partial melting products of the lower crust. A more elegant hypothesis seems to be the model ofGunn &Mooser (1971), who consider these volcanic rocks as partial melting products of oceanic tholeiites or their high pressure derivatives in the sense ofRaleigh &Lee (1969).     

Scaling Transition in Earthquake Sources: A Possible Link Between Seismic and Laboratory Measurements

We estimate the corner frequencies of 20 crustal seismic events from mainshock–aftershock sequences in different tectonic environments (mainshocks 5.7 < M _W < 7.6) using the well-established seismic coda ratio technique (Mayeda et al. in Geophys Res Lett 34:L11303, 2007; Mayeda and Malagnini in Geophys Res Lett, 2010), which provides optimal stability and does not require path or site corrections. For each sequence, we assumed the Brune source model and estimated all the events’ corner frequencies and associated apparent stresses following the MDAC spectral formulation of Walter and Taylor (A revised magnitude and distance amplitude correction (MDAC2) procedure for regional seismic discriminants, 2001), which allows for the possibility of non-self-similar source scaling. Within each sequence, we observe a systematic deviation from the self-similar \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - 3} \) line, all data being rather compatible with \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - (3 + \varepsilon )} \) , where ε > 0 (Kanamori and Rivera in Bull Seismol Soc Am 94:314–319, 2004). The deviation from a strict self-similar behavior within each earthquake sequence of our collection is indicated by a systematic increase in the estimated average static stress drop and apparent stress with increasing seismic moment (moment magnitude). Our favored physical interpretation for the increased apparent stress with earthquake size is a progressive frictional weakening for increasing seismic slip, in agreement with recent results obtained in laboratory experiments performed on state-of-the-art apparatuses at slip rates of the order of 1 m/s or larger. At smaller magnitudes (M _W < 5.5), the overall data set is characterized by a variability in apparent stress of almost three orders of magnitude, mostly from the scatter observed in strike-slip sequences. Larger events (M _W > 5.5) show much less variability: about one order of magnitude. It appears that the apparent stress (and static stress drop) does not grow indefinitely at larger magnitudes: for example, in the case of the Chi–Chi sequence (the best sampled sequence between M _W 5 and 6.5), some roughly constant stress parameters characterize earthquakes larger than M _W ~ 5.5. A representative fault slip for M _W 5.5 is a few tens of centimeters (e.g., Ide and Takeo in J Geophys Res 102:27379–27391, 1997), which corresponds to the slip amount at which effective lubrication is observed, according to recent laboratory friction experiments performed at seismic slip velocities (V ~ 1 m/s) and normal stresses representative of crustal depths (Di Toro et al. in Nature in press, 2011, and references therein). If the observed deviation from self-similar scaling is explained in terms of an asymptotic increase in apparent stress (Malagnini et al. in Pure Appl Geophys, 2014, this volume), which is directly related to dynamic stress drop on the fault, one interpretation is that for a seismic slip of a few tens of centimeters (M _W ~ 5.5) or larger, a fully lubricated frictional state may be asymptotically approached.     

Age-dependent Large-scale Fabric of the Mantle Lithosphere as Derived from Surface-wave Velocity Anisotropy

—Systematic variations of the seismic radial anisotropy ξ to depths of 200–250 km in North America and Eurasia and their surroundings are related to the age of continental provinces, and typical depth dependences of ξ _R are determined. The relative radial anisotropy ξ _R in the mantle lithosphere of Phanerozoic orogenic belts is characterized by ν _SH > ν _SV , with its maximum depth of about 70 km, on the average, while beneath old shields and platforms, it exhibits a maximum deviation from ACY400 model (Montagner and Anderson, 1989) at depths of about 100 km with ν _SV ≥ν _SH signature. An interpretation of the observed seismic anisotropy by the preferred orientation of olivine crystals results in a model of the mantle lithosphere characterized by anisotropic structures plunging steeply beneath old shields and platforms, compared to less inclined anisotropies beneath Phanerozoic regions. This observation supports the idea derived from petrological and geochemical observations that a mode of continental lithosphere generation may have changed throughout earth's history.     

Effects of Location Errors in Pattern Informatics

The effect of location errors in the performance of seismicity-based forecasting methods was studied here using one particular binary forecast technique, the Pattern Informatics (PI) technique (Rundle et al., Proc Nat Acad Sci USA 99, 2514–2521, 2002; Tiampo et al., Pure Appl Geophys 159, 2429–2467, 2002). The Southern Californian dataset was used to generate a series of perturbed catalogs by adding different levels of noise to epicenter locations. The PI technique was applied to these perturbed datasets to perform retrospective forecasts that were evaluated by means of skill scores, commonly used in atmospheric sciences. These results were then compared to the effectiveness obtained from the original dataset. Isolated instances of decline of the PI performance were observed due to the nature of the skill scores themselves, but no clear trend of degradation was identified. Dependence on the total number of events in a catalog also was studied, with no systematic degradation in the performance of the PI for catalogs with events in the cases studied. These results suggest that the stability of the PI method is due to the invariance of the clustering patterns identified by the TM metric (Thirumalai and Mountain, Phys Rev A 39, 3563–3573, 1989) when applied to seismicity.     

Modeling Seismic Attributes of Pn Waves using the Spectral-Element Method

To investigate the nature of Pn propagation, we have implemented the spectral-element method (SEM) for vertically and laterally varying media with and without attenuation. As a practical measure, essential features of the Pn waves are distilled into seismic attributes including arrival times, amplitudes and pulse frequencies. To validate the SEM simulations, we first compare the SEM results with reflectivity calculations of Braile and Smith (Geophys. J.R. Astr. Soc. 40, 145–176, 1975) and then to the asymptotic results of ?erveny and Ravindra (Theory of Seismic Headwaves, University of Toronto Press, pp. 235–250, 1971). Models with random, laterally varying Moho structures are then simulated, where the amplitude and pulse frequency characteristics are found to be stable to small Moho interface perturbations. SEM calculations for models with different upper-mantle velocity gradients are next performed where it is found that interference effects can strongly influence the Pn amplitudes and pulses frequencies. For larger-scale, laterally varying structures, SEM models similar to that found along the Hi-CLIMB array in Tibet are then performed. It is observed that large-scale structures, along with small-scale structures, upper-mantle velocity gradients and attenuation, can all significantly affect the Pn attributes. Ambiguities between upper-mantle velocity gradients and attenuation are also found when using Pn amplitudes and pulse frequency attributes. These ambiguities may be resolved, to some degree, by using the curvature of the travel times at longer regional distance, however, this would also be complicated by lateral variability.     

Flood basalts from dharwars of Mysore,India

The term flood basalt is redefined emphasizing the importance of the subaerial environment. Using the well established physical criteria of aerial extent, internal structures, time of extrusion and associations, flood basalt activity is distinguished in the Archeans (Dharwars) of Mysore from the geosynclinal volcanics. Study of chemical composition of the Dharwar and other Archean volcanics in the light ofSugimura’s (1968) SWS index, and plotting of the chemical analyses on theMacdonald andKatsura’s (1964) alkali-silica diagram,Kuno’s (1968) alkali-alumina-silica diagram and Scheynamann’s silica-Niggli qz diagram shows both geosynclinal and subaerial volcanics are mainly tholeiitic. Therefore in deciphering the environment of volcanism, it is suggested that the physical criteria take precedence over chemical composition.     

Location of Moderate-Sized Earthquakes Recorded by the NARS�CBaja Array in the Gulf of California Region Between 2002 and 2006

We relocated the hypocentral coordinates of small to moderate-sized earthquakes reported by the National Earthquake Information Center (NEIC) between April 2002 and August 2006 in the Gulf of California region and recorded by the broadband stations of the network of autonomously recording seismographs (NARS?CBaja array). The NARS?CBaja array consists of 19 stations installed in the Baja California peninsula, Sonora and Sinaloa, Mexico. The events reported by the preliminary determinations of epicenters (PDE) catalog within the period of interest have moment magnitudes (M _w) ranging between 1.1 and 6.7. We estimated the hypocentral location of these events using P and S wave arrivals recorded by the regional broadband stations of the NARS?CBaja and the RESBAN (Red Sismológica de Banda Ancha) arrays and using a standard location procedure with the HYPOCENTER code (Lienert and Havskov in Seism Res Lett 66:26?C36, 1995) as a preliminary step. To refine the location of the initial hypocenters, we used the shrinking box source-specific station term method of Lin and Shearer (J Geophys Res 110, B04304, 2005). We found that most of the seismicity is distributed in the NW?CSE direction along the axis of the Gulf of California, following a linear trend that, from north to south, steps southward near the main basins (Wagner, Delfin, Guaymas, Carmen, Farallon, Pescadero and Alarcon) and spreading centers. We compared the epicentral locations reported in the PDE with the locations obtained using regional arrival times, and we found that earthquakes with magnitudes in the range 3.2?C5.0?mb differ on the average by as much as 43?km. For the M _w magnitude range between 5 and 6.7 the discrepancy is less, differing on the average by about 25?km. We found that the relocated epicenters correlate well with the main bathymetric features of the Gulf.     

Ergodicity and Earthquake Catalogs: Forecast Testing and Resulting Implications

Recently the equilibrium property of ergodicity was identified in an earthquake fault system (Tiampo et al., Phys. Rev. Lett. 91, 238501, 2003; Phys. Rev. E 75, 066107, 2007). Ergodicity in this context not only requires that the system is stationary for these networks at the applicable spatial and temporal scales, but also implies that they are in a state of metastable equilibrium, one in which the ensemble averages can be substituted for temporal averages when studying their behavior in space and time. In this work we show that this property can be used to identify those regions of parameter space which are stationary when applied to the seismicity of two naturally-occurring earthquake fault networks. We apply this measure to one particular seismicity-based forecasting tool, the Pattern Informatics index (Tiampo et al., Europhys. Lett. 60, 481–487, 2002; Rundle et al., Proc. National Acad. Sci., U.S.A., Suppl. 1, 99, 2463, 2002), in order to test the hypothesis that the identification of ergodic regions can be used to improve and optimize forecasts that rely on historic seismicity catalogs. We also apply the same measure to synthetic catalogs in order to better understand the physical process that affects this accuracy. We show that, in particular, ergodic regions defined by magnitude and time period provide more reliable forecasts of future events in both natural and synthetic catalogs, and that these improvements can be directly related to specific features or properties of the catalogs that impact the behavior of their spatial and temporal statistics.     

Attenuation Tomography Based on Strong Motion Data: Case Study of Central Honshu Region, Japan

Three-dimensional frequency dependent S-wave quality factor (Q_β(f)) value for the central Honshu region of Japan has been determined in this paper using an algorithm based on inversion of strong motion data. The method of inversion for determination of three-dimensional attenuation coefficients is proposed by Hashida and Shimazaki (J Phys Earth. 32, 299–316, 1984) and has been used and modified by Joshi (Curr Sci. 90, 581–585, 2006; Nat Hazards. 43, 129–146, 2007) and Joshi et al. (J. Seismol. 14, 247–272, 2010). Twenty-one earthquakes digitally recorded on strong motion stations of Kik-net network have been used in this work. The magnitude of these earthquake ranges from 3.1 to 4.2 and depth ranging from 5 to 20 km, respectively. The borehole data having high signal to noise ratio and minimum site effect is used in the present work. The attenuation structure is determined by dividing the entire area into twenty-five three-dimensional blocks of uniform thickness having different frequency-dependent shear wave quality factor. Shear wave quality factor values have been determined at frequencies of 2.5, 7.0 and 10 Hz from record in a rectangular grid defined by 35.4°N to 36.4°N and 137.2°E to 138.2°E. The obtained attenuation structure is compared with the available geological features in the region and comparison shows that the obtained structure is capable of resolving important tectonic features present in the area. The proposed attenuation structure is compared with the probabilistic seismic hazard map of the region and shows that it bears some remarkable similarity in the patterns seen in seismic hazard map.     

The Frequency Dependence of Harmonic Fluid Flow Through Networks of Cracks and Pores

—The "dynamic" permeability k(ω) of heterogeneous networks of cracks, tubes and spheres, was determined by numerically simulating the harmonic flow of an interstitial fluid for a wide range of frequencies. For comparison with previous works, this procedure was applied to the 100 network realizations used in Bernabé (1995). In most cases, the calculated frequency dependence of the real and imaginary parts of k(ω) was consistent with the JKD model (Johnson et al., 1987), showing a transition from "viscous", macroscopic flow at low frequencies to "inertial" flow at high frequencies. The viscous skin depth δ _c at the transition was found to be proportional to the critical capillary radius r _c from a capillary invasion (Katz and Thompson, 1986). A simple explanation is that these two length scales arise from the same percolation problem. On the other hand, δ _c was not well correlated with the JKD parameter Λ. The conclusion is that Λ and δ _c (or r _c ?) are two independent parameters, derived from two unrelated approaches (i.e., weighted averaging and percolation theory). Finally, an attempt was made to relax the initial assumptions of a rigid solid matrix and an incompressible fluid. It was observed that the effect of the fluid compressibility could occasionally be very large, especially when networks with large amounts of storage pore space were considered.     

New K-Ar ages,chemical analyses and magnetic data of rocks from the islands of santa maria (azores), porto santo and madeira (madeira archipelago) and gran canaria (Canary Islands)

K/Ar-determinations, major and trace element chemical analyses and magnetic data are reported for rocks from Santa Maria (Azores), Madeira and Porto Santo, and Gran Canaria. Based on these data, the age of the basement of Santa Maria is believed to have formed between about 5.2. and 4.6 m.y.; the unconformably overlying pillow complex interbedded with fossiliferous calcarenites about 3.8 to 3.3 m.y. with the capping subaerial basalt being part of the same magmatic phase. The major erosional phase levelling the basement is thus approximately synchronous with the major Pliocene regression (R2) on Gran Canaria (Lietz andSchmincke, 1975), possibly indicating a widespread eustatic event. The upper part of the submarine, partly fossiliferous series of Porto Santo was dated as ca. 12 to 13 m.y., and the Quaternary age for the major basalt formation in Eastern Madeira (Watkins andAbdel Monem, 1971) is confirmed. The ages of several formations on Gran Canaria were slightly revised. Chemical differences between basement (shield) and later posterosional series on Gran Canaria are re-emphasized by the new data, while such differences are much less pronounced between the basement and younger series on Santa Maria.     

Submarine basaltic volcanism: Morphometric parameters for discriminating hyaloclastites from hyalotuffs

Hyaloclastites are volcaniclastic rocks generated by non-explosive granulation of volcanic glass which takes place when basaltic magmas are quenched by contact with water; hyaloclastites are common products of deep submarine basaltic central volcanoes (seamounts and guyots). We suggested (Honnorlz, 1966) calling hyalotuffs the actual pyroclastic rocks which are generated by phreatomagmatic and phreatic explosions taking place when basaltic volcanoes crupt in shallow waters; hyalotuffs are restricted to shallow subaqueous conditions since no volcanic explosion can occur in deep seas. The distinction between hyaloclastites and hyalotuffs is therefore a useful tool when reconstructing the paleogeography of ancient submarine volcanic edifices and the mechanism by which their lavas were emplaced. We propose using two sets of morphometric parameters to discriminate the hyaloclastites from the hyalotuffs. The granules making up these two volcaniclastic rock types plot in different areas of either 1) a ternary diagram the apexes of which represent the grain planarity (P), convexity (V) and concavity (C) % («roundness» according to Szadeczky-Kardoss); or 2) a binary diagram relating the number of grain corners (N) to their planarity (P) %.     

Phytoplankton and zooplankton of some lakes in northeastern Norway

During the summer of 1966 some limnological aspects of thirteen lakes and one pond in the most northeasterly part of Norway were studied. A total of 418 taxa of algae including 104 genera, and 70 taxa of zooplankton including 46 genera are identified. The composition of the plankton in all the waters is shown to be typical for fundamentally oligotrophic lakes. A few of the lakes show some humic influence. This is in agreement with the hydrographic results from these waters (Bøyum [1]).     

An Independent Assessment of the Load/Unload Response Ratio (LURR) Proposed Method of Earthquake Prediction

The Load/Unload Response Ratio (LURR) method is a proposed technique to predict earthquakes that was first put forward by Yin (1987). LURR is based on the idea that when an area enters the damage regime, the rate of seismic activity during loading of the tidal cycle increases relative to the rate of seismic activity during unloading in the months to one year preceding a large earthquake. Since earth tides generally contribute the largest temporal variations in crustal stress, it seems plausible that earth tides would trigger earthquakes in areas that are close to failure (e.g., Vidale et al., 1998). However, the vast majority of studies have shown that earth tides do not trigger earthquakes (e.g., Vidale et al., 1998; Heaton, 1982; Rydelek et al., 1992). In this study, we conduct an independent test of the LURR method, since there would be important scientific and social implications if it were proven to be a robust method of earthquake prediction. Smith and Sammis (2004) undertook a similar study and found no evidence that there was predictive significance to the LURR method. We have repeated calculations of LURR for the Northridge earthquake in California, following both the parameters of X.C. Yin (personal communication) and the somewhat different ones of Smith and Sammis (2004). Though we have followed both sets of parameters closely, we have been unable to reproduce either set of results. Our examinations have shown that the LURR method is very sensitive to certain parameters. Thus it seems likely that the discrepancies between our results and those of previous studies are due to unaccounted for differences in the calculation parameters. A general agreement was made at the 2004 ACES Workshop in China between research groups studying LURR to work cooperatively to resolve the differences in methods and results, and thus permit more definitive conclusions on the potential usefulness of the LURR method in earthquake prediction.     

Origin of rhyolites by anatectic melting of granitic crustal rocks

Rhyolitic pumice deposits on some East Aegean Islands are the remnants of a tuff sheet which covered formerly more than 2500 km². The rhyolites originated by melting of granitic crustal rocks. The petrological properties leading to this conclusion are as follows:

-mineralogical composition of the pumice being 60 wt % of glass and 40 % of relictic granitic minerals (quartz, plagioclase, potash feldspar, biotite) in «phenocrysts» up to 1 cm Ø.

-phenomena of strong corrosion, resorption and melting of all light colored «phenocrysts» in the rhyolitic pumice.

-mineralogical and chemical identity between pumice and granitic xenolites found in great quantity in the pumice tuff.

-structurally transitional types between pumice and the granites in different state of melting.

-eutectic composition of the pumice glass for a pressure of 2 kb.

     

A Godunov-type Finite Volume Scheme for Meso- and Micro-scale Flows in Three Dimensions

This short note reports the extension of the f-waves approximate Riemann solver (Ahmad and Lindeman, 2007; LeVeque, 2002; Bale et al., 2002) for three-dimensional meso- and micro-scale atmospheric flows. The Riemann solver employs flux-based wave decomposition for the calculation of Godunov fluxes and does not require the explicit definition of the Roe matrix to enforce conservation. The other important feature of the Riemann solver is its ability to incorporate source term due to gravity without introducing discretization errors. The resulting finite volume scheme is second-order accurate in space and time. The finite-difference schemes currently used in atmospheric flow models are neither conservative nor able to resolve regions of sharp gradients. The finite volume scheme described in this paper is fully conservative and has the ability to resolve regions of sharp gradients without introducing spurious oscillations in the solution. The scheme shows promise in accurately resolving flows on the meso- and micro-scales and should be considered for implementation in the dynamical cores of next generation meso- and micro-scale atmospheric flow models.     

Carbonatite-alkalic complex of Panwad-Kawant,Gujarat, and its bearing on the structural characteristics of the area

Carbonatite-alkalic rocks occur in the form of dykes and small volcanic plugs in the area, with major central type volcanic activity restricted to Amba Dongar. The trappean flows of Blanford and Bose are identified as plagioclasecalcite rocks. An attempt is made to explain the origin of these rocks which are extensively cut by dykes of alkaline rocks, carbonatites and dolerites. By far the dominant lavas are «fissure phonolites» (Wright, 1963) and tinguaites. The chemical analyses of these rocks show that the magma is mainly of continental sodic alkaline suite, probably turning sodi-potassic, a suggestion drawn from the occurrence of lamprophyres and pseudoleucite tinguaites, and the higher potassium contents of some rocks. Bagh sediments are mainly represented by sandstones which show mild contact effects with carbonatite, especially in the south.