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.     

Three thousand years of flank and central vent eruptions of the San Salvador volcanic complex (El Salvador) and their effects on El Cambio archeological site: a review based on tephrostratigraphy

The volcanic events of the last 3,000?years at San Salvador volcanic complex are reviewed using detailed stratigraphic records exposed in new excavations between 2005 and 2007, at El Cambio archeological site (Zapotitán Valley, El Salvador), and in other outcrops on the northern and northwestern sectors of the complex. The sequences that overlie Tierra Blanca Joven (cal. 429?±?107?ad), from the Ilopango caldera, comprise the Loma Caldera (cal. 590?±?90?ad) and El Playón (1658?C1671) deposits and the San Andrés Tuff (cal. 1031?±?29?ad), related to El Boquerón Volcano. The surge deposits within the El Playón, San Andrés Tuff and overlying Talpetate II sequences indicate the significance of phreatomagmatic phases in both central vent and flank eruptions during the last 1,600?years. Newly identified volcanic deposits underlying Tierra Blanca Joven at El Cambio extend the stratigraphic record of the area to 3,000?years?bp. Paleosols interstratified with those deposits contain cultural artifacts which could be associated with the Middle Preclassic period (900?C400?bc). If correct, human occupation of the site during the Preclassic period was more intense than previously known and volcanic eruptions must have affected prehistoric settlements. The archeological findings provide information on how prehistoric populations dealt with volcanic hazards, thousands of years ago in the eastern Zapotitán Valley, where several housing projects are currently being developed. The new stratigraphic and volcanological data can be used as a basis for local and regional hazard assessment related to future secondary vent activity in the San Salvador Volcanic Complex.     

Constructing Constitutive Relationships for Seismic and Aseismic Fault Slip

For the purpose of modeling natural fault slip, a useful result from an experimental fault mechanics study would be a physically-based constitutive relation that well characterizes all the relevant observations. This report describes an approach for constructing such equations. Where possible the construction intends to identify or, at least, attribute physical processes and contact scale physics to the observations such that the resulting relations can be extrapolated in conditions and scale between the laboratory and the Earth. The approach is developed as an alternative but is based on Ruina (1983) and is illustrated initially by constructing a couple of relations from that study. In addition, two example constitutive relationships are constructed; these describe laboratory observations not well-modeled by Ruina’s equations: the unexpected shear-induced weakening of silica-rich rocks at high slip speed (Goldsby and Tullis, 2002) and fault strength in the brittle ductile transition zone (Shimamoto, 1986). The examples, provided as illustration, may also be useful for quantitative modeling.     

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.     

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.     

New Seismic Noise Models Obtained Using Very Broadband Stations

It has been two decades since the last comprehensive standard model of ambient earth noise was published Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The PETERSON model was updated by analyzing the absolute quietest conditions for stations within the GSN (Berger et al. in J Geophys Res 109, 2005; Mcnamara and Buland in Bull Seism Soc Am 94:1517–1527, 2004; Ringler et al. in Seismol Res Lett 81(4) doi:10.1785/gssrl.81.4.605, 2010). Unfortunately, both the original model and the updated models did not include any deployed station in North Africa and Middle East, which reflects the noise levels within the desert environment of those regions. In this study, a survey was conducted to create a new seismic noise model from very broadband stations which recently deployed in North Africa. For this purpose, 1 year of continuous recording of seismic noise data of the Egyptian National Seismic Network (ENSN) was analyzed in order to create a new noise model. Seasonal and diurnal variations in noise spectra were recorded at each station. Moreover, we constructed a new noise model for each individual station. Finally, we obtained a new cumulative noise model for all the stations. We compared the new high-noise model (EHNM) and new low-noise model (ELNM) with both the high-noise model (NHNM) and low-noise model (NLNM) of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The obtained noise levels are considerably lower than low-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993) at ultra long period band (ULP band), but they are still below the high-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The results of this study could be considered as a first step to create permanent seismic noise models for North Africa and Middle East regions.     

Gradual Fault Weakening with Seismic Slip: Inferences from the Seismic Sequences of L’Aquila, 2009, and Northridge, 1994

We estimate seismological fracture energies from two subsets of events selected from the seismic sequences of L’Aquila (2009), and Northridge (1994): 57 and 16 selected events, respectively, including the main shocks. Following Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), we postulate that fracture energy (G) represents the post-failure integral of the dynamic weakening curve, which is described by the evolution of shear traction as a function of slip. Following a direct-wave approach, we compute mainshock-/aftershock-source spectral ratios, and analyze them using the approach proposed by Malagnini et al. (Pure Appl. Geophys., this issue, 2014) to infer corner frequencies and seismic moment. Our estimates of source parameters (including fracture energies) are based on best-fit grid-searches performed over empirical source spectral ratios. We quantify the source scaling of spectra from small and large earthquakes by using the MDAC formulation of Walter and Taylor (A revised Magnitude and Distance Amplitude Correction (MDAC2) procedure for regional seismic discriminants, 2001). The source parameters presented in this paper must be considered as point-source estimates representing averages calculated over specific ruptured portions of the fault area. In order to constrain the scaling of fracture energy with coseismic slip, we investigate two different slip-weakening functions to model the shear traction as a function of slip: (i) a power law, as suggested by Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), and (ii) an exponential decay. Our results show that the exponential decay of stress on the fault allows a good fit between measured and predicted fracture energies, both for the main events and for their aftershocks, regardless of the significant differences in the energy budgets between the large (main) and small earthquakes (aftershocks). Using the power-law slip-weakening function would lead us to a very different situation: in our two investigated sequences, if the aftershock scaling is extrapolated to events with large slips, a power law (a la Abercrombie and Rice) would predict unrealistically large stress drops for large, main earthquakes. We conclude that the exponential stress evolution law has the advantage of avoiding unrealistic stress drops and unbounded fracture energies at large slip values, while still describing the abrupt shear-stress degradation observed in high-velocity laboratory experiments (e.g., Di Toro et al., Fault lubrication during earthquakes, Nature 2011).     

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.     

Characterizing and Locating Very Weak (−2.2 ≥ M L ≥ −3.4) Induced Seismicity in Unstable Sandstone Cliffs by Nanoseismic Monitoring

Unloaded natural rock masses are known to generate seismic signals (Green et al., 2006; Hainzl et al., 2006; Husen et al., 2007; Kraft et al., 2006). Following a 1,000 m³ mass failure into the Mediterranean Sea, centimeter-wide tensile cracks were observed to have developed on top of an unstable segment of the coastal cliff. Nanoseismic monitoring techniques (Wust-Bloch and Joswig, 2006; Joswig, 2008), which function as a seismic microscope for extremely weak seismic events, were applied to verify whether brittle failure is still generated within this unconsolidated sandstone mass and to determine whether it can be detected. Sixteen days after the initial mass failure, three small-aperture sparse arrays (Seismic Navigation Systems-SNS) were deployed on top of this 40-m high shoreline cliff. This paper analyzes dozens of spiky nanoseismic (?2.2 ≥ M _L ≥ ?3.4) signals recorded over one night in continuous mode (at 200 Hz) at very short slant distances (3–67 m). Waveform characterization by sonogram analysis (Joswig, 2008) shows that these spiky signals are all short in duration (>0.5 s). Most of their signal energy is concentrated in the 10–75 Hz frequency range and the waveforms display high signal similarity. The detection threshold of the data set reaches M _L ?3.4 at 15 m and M _L ?2.7 at 67 m. The spatial distribution of source signals shows 3-D clustering within 10 m from the cliff edge. The time distribution of M _L magnitude does not display any decay pattern of M _L over time. This corroborates an unusual event decay over time (modified Omori’s law), whereby an initial quiet period is followed by regained activity, which then fades again. The polarization of maximal waveform amplitude was used to estimate spatial stress distribution. The orientation of ellipses displaying maximal signal energy is consistent with that of tensile cracks observed in the field and agrees with rock mechanics predictions. The M _L– surface rupture length relationship displayed by our data fits a constant-slope extrapolation of empirical data collected by Wells and Coppersmith (1994) for normal fault features at much larger scale. Signal characterization and location as well as the absence of direct anthropogenic noise sources near the monitoring site, all indicate that these nanoseismic signals are generated by brittle failure within the top section of the cliff. The atypical event decay over time that was observed suggests that the cliff material is undergoing post-collapse bulk strain accommodation. This feasibility study demonstrates the potential of nanoseismic monitoring in rapidly detecting, locating and analyzing brittle failure generated within unconsolidated material before total collapse occurs.     

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.     

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.     

The interdependence of seismic and volcanic phenomena: Some space — Time relationships in seismicity and volcanism

Research is described aimed at detecting significant relationships between patterns of earthquakes and volcanic eruptions. A large data file of seismic and volcanic events has been assembled on punched cards, and two computer programmes used to relate earthquakes to active or potentially active volcanic centres, and deep earthquakes to succeeding shallower ones and volcanic eruptions after the manner suggested byBlot. It is concluded that recognisable sequences occur in some areas, but that the process is neither universal nor invariable. The presence of random events or events due to other causes almost always obscures patterns of seismicity which seem to be related to volcanic eruptions. The «Blot process» appears to operate in some cases but is probably a second order phenomenon: primarily, correlated sequences of seismic and volcanic events result from periods of tectonic instability and perhaps increased tensional conditions which affect very wide areas of the Earth’s surface for periods of several months to several years at a time.     

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.     

Petrological features of boseti mountains,a complex volcanic system in the axial portion of the main Ethiopian rift

A petrographical and chemical research carried out mainly for the rocks belonging to the 2nd cycle of volcanic activity of Boseti Mts., shows the predominance of the salic over the mafic and the intermediate products. The magmatic evolution, which chiefly took place in a thermodynamic environment with initial low \(P_{H_2 O} \) and \(f_{O_2 } \) , has brought about subaluminiferous mugearitic types, also with trachytic affinity. Among peralkali-rhyolites, an anorthoclase fractionation can explain both their chemical variation and the production of pantelleritic liquids.     

The Stress Accumulation Method and the Pattern Informatics Index: Complementary Approaches to Earthquake Forecasting

In this work, we provide a joint study of the stress accumulation method (SAM) (King and Bowman, 2003) and the Pattern Informatics (PI) index (Tiampo et al., 2002b). We examine the theoretical underpinnings for the similarities between the two techniques, as well as the differences in their application. The SAM technique is employed to determine likely mechanisms for smaller areas of increased probability identified by the PI index, while a modified version of the PI index can be used to locate regions where the smaller magnitude associated with the anomaly is below the resolution of the SAM. Finally, we present three case studies from different regions of the San Andreas fault system to illustrate both their complementary nature, as well as the advantages to combining them in one synthesized analysis.     

Attenuation and Shock Waves in Linear Hereditary Viscoelastic Media; Strick–Mainardi,Jeffreys–Lomnitz–Strick and Andrade Creep Compliances

Dispersion, attenuation and wavefronts in a class of linear viscoelastic media proposed by Strick and Mainardi (Geophys J R Astr Soc 69:415–429, 1982) and a related class of models due to Lomnitz, Jeffreys and Strick are studied by a new method due to the author. Unlike the previously studied explicit models of relaxation modulus or creep compliance, these two classes support propagation of discontinuities. Due to an extension made by Strick, either of these two classes of models comprise both viscoelastic solids and fluids. We also discuss the Andrade viscoelastic media. The Andrade media do not support discontinuity waves and exhibit the pedestal effect.     

An Elliptical Model for Deformation Due to Groundwater Fluctuations

Historically, surface subsidence as a result of subsurface groundwater fluctuations have produced important and, at times, catastrophic effects, whether natural or anthropogenic. Over the past 30?years, numerical and analytical techniques for the modeling of this surface deformation, based upon elastic and poroelastic theory, have been remarkably successful in predicting the magnitude of that deformation (Le Mouélic and Adragna in Geophys Res Lett 29:1853, 2002). In this work we have extended the formula for a circular-shaped aquifer (Geertsma in J Petroleum Tech 25:734–744, 1973) to a more realistic elliptical shape. We have improved the accuracy of the approximation by making use of the cross terms of the expansion for the elliptic coordinates in terms of the eccentricity, e, and the mean anomaly angle, M, widely used in astronomy. Results of a number of simulations, in terms of e and M developed from the transcendental Kepler equation, are encouraging, giving realistic values for the elliptical approximation of the vertical deformation due to groundwater change. Finally, we have applied the algorithm to modeling of groundwater in southern California.     

Neogene calc-alkaline extrusive and intrusive rocks of the karalar-yesiller area

The Karalar-Ye?iller area lies on the estern flank of the Kazda \(\bar g\) i massif in northwestern Anatolia (Asiatic Turkey), and includes about 600 square kilometers, of which approximately 80% is covered by rhyodacite-quartz latite extrusive rocks, and a comagmatic granodiorite-quartz monzonite batholith, and stocks, all of middle Miocene age. Extrusive rocks consist chiefly of rhyodacite lava flows of the Hallaçlar Formation and quartz latite-rhyodacite ash-flow tuff, lava, and mudflow deposits of the overlying Dede Tepe Formation. These volcanic rocks lie on a basement composed, in asceding stratigraphic order, of: 1) pre-Permian Kalabak sequence, 2) Upper Triassic Halilar Formation, 3) post-Upper Triassic Ba \(\bar g\) burun Formation, and 4) allochthons of middle Permian: and Upper Jurassic limestone. Intrusion of a granodiorite-quartz monzonite batholith during the middle Miocene was accompanied or shortly followed by extensive alteration and base metal mineralization of the Hallaçlar and older rocks. Intrusion of six rhyodacitequartz latite stocks followed the main phase of hydrothermal alteration. At least one of these stocks may have contributed material to the Dede Tepe Formation of middle Miocene age. Field relations, petrologic, and geochemical data as well as radiometric age dates suggest that the intrusive and extrusive rocks are comagmatic. Parent magma is probably derived from partial melting of subducted oceanic crust and accompanying oceanic sediments. Geologic relations locally indicate assimilation of sialic crust by the grandiorite-quartz monzonite batholith.     

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.