inverse kinematic and forward dynamic models of the 2002 denali fault earthquake alaska we perform inverse kinematic and forward dynamic models of the m 7 9 2002 denali fault alaska earthquake to shed light on the rupture process and dynamics of this event which took place on a geometrically complex fault system in central alaska we use a combination of local seismic and global positioning system gps data for our kinematic inversion and find that the slip distribution of this event is characterized by three major asperities on the denali fault the rupture nucleated on the susitna glacier thrust fault and after a pause propagated onto the strike slip denali fault approximately 216 km to the east the rupture abandoned the denali fault in favor of the more southwesterly directed totschunda fault three dimensional dynamic models of this event indicate that the abandonment of the denali fault for the totschunda fault can be explained by the totschunda fault s more favorable orientation with respect to the local stress field however a uniform tectonic stress field cannot explain the complex slip pattern in this event we also find that our dynamic models predict discontinuous rupture from the denali to totschunda fault segments such discontinuous rupture helps to qualitatively improve our kinematic inverse models two principal implications of our study are 1 a combination of inverse and forward modeling can bring insight into earthquake processes that are not possible with either technique alone and 2 the stress field on geometrically complex fault systems is most likely not due to a uniform tectonic stress field that is resolved onto fault segments of different orientations rather other forms of stress heterogeneity must be invoked to explain the observed slip patterns

flow directions in dikes from anisotropy of magnetic susceptibility data the bootstrap way one of the first applications of anisotropy of magnetic susceptibility ams was an attempt to determine flow directions from mafic dikes khan 1962 since the seminal work of knight and walker 1988 defining the expected behavior of ams in response to magma flow there has been increasing interest in using ams for this purpose here we present a quantitative method for interpretation of ams data from dikes using a parametric bootstrap first dikes must be sampled with at least five and preferrably more samples from within 10 cm of the dike margin the distributions of the eigenvalues and eigenvectors of the ams tensor are delineated by calculating eigenparameters of many bootstrapped paradata sets we generate paradata sets by first selecting a sample at random then calculating a replacement set of data by drawing tensor elements from normal distributions with the mean and standard deviation of the entire site the bounds containing 95 of the eigenparameters of the bootstrapped data serve as confidence limits for the parameter of interest classification of dikes proceeds as follows sites whose maximum and intermediate eigenvalues could not be distinguished are deemed uninterpretable in addition sites with principal eigenvectors with angles 45 away from the dike margin inverse or with markedly different directions on either side of the dike scissored are excluded the remaining dikes are classified as having unique flow direction information if the principal eigenvectors from at least one side are distinct from the dike plane based on the distribution of the bootstrapped principal eigenvectors if neither side has principal eigenvectors distinct from the dike plane the dikes are classified as having lineation information only a study comprising 251 dikes from the troodos ophiolite has 151 sites with directional data 38 sites with lineations only 7 inverse sites 5 scissored sites and 55 sites not fitting into any other category the flow directions interpreted from the data were generally southerly toward a fossil transform zone

comment on magnetic studies of magma supply and sea floor metamorphism troodos ophiolite dikes by g j borradaile and d gauthier anisotropy of magnetic susceptibility is used as a proxy for the determination of magmatic flow direction in mafic dykes here we take advantage of the dataset from g j borradaile and d gauthier to comment three points 1 the sampling strategy 2 the geometric relationship between magnetic axis dyke and 3 an alternative interpretation to obtain a flow direction the magnetic lineations published by borradaile and gauthier correspond to the zone axis of the dyke and magnetic foliation poles questioning the reliability of the magnetic lineation as a flow estimate an alternative interpretation is based on the use of the tiling of the magnetic foliation plane against the dyke wall

the effect of bandwidth limitations on the inference of earthquake slip weakening distance from seismograms numerous researchers have obtained estimates of slip weakening distance dc and fracture energy for recent earthquakes dc is often observed to be a significant fraction of the total slip and tends to correlate with total slip although these observations may well be true of real earthquakes we show that low pass filtering of strong motion seismograms can also produce some of these effects in inverted rupture models we test the accuracy of dc estimates by calculating them in low pass filtered versions of models a and b of guatteri and spudich 2000 models a and b are two different rupture models for a hypothetical m 6 5 earthquake and they have nearly identical rupture time slip and stress drop distributions and nearly identical predicted seismograms but dc for model b is about twice that for model a by low pass filtering slip models a and b at 1 0 hz we simulate the blurring effects of band limited waveform inversions on these slip models at each point on a fault graphic 1 is defined to be the slip at the time of the peak slip speed at that point low pass filtering the slip models causes an upward bias in dc inferred from stress slip curves and it causes an artificial correlation between graphic 2 and the total slip low pass filtering might also bias fracture energy high and radiated energy low these biases should be considered when interpreting dc derived from band limited slip models of real earthquakes 1 embed inline graphic 1 gif 2 embed inline graphic 2 gif

three dimensional dynamic simulation of the 1992 landers earthquake the 1992 landers california earthquake magnitude 7 3 was modeled as the propagation of a spontaneous rupture controlled by a realistic prestress distribution with the use of a three dimensional finite difference method the dynamic rupture reproduces the general slip pattern used to compute the initial stress distribution and generates near fault ground motions at the surface similar to observations the simulated rupture propagates on the fault along a complex path with highly variable speed and rise time changing the magnitude and pattern of the stress significantly the method provides the framework to estimate earthquake rupture parameters from recorded seismic and geodetic data

potbellies wasp waists and superparamagnetism in magnetic hysteresis because the response of a magnetic substance to an applied field depends strongly on the physical properties of the material much can be learned by monitoring that response through what is known as a magnetic hysteresis loop the measurements are rapid and quickly becoming part of the standard set of tools supporting paleomagnetic research yet the interpretation of hysteresis loops is not simple it has become apparent that although classic single domain pseudo single domain and multidomain loops described in textbooks occur in natural samples loops are frequently distorted having constricted middles wasp waisted loops or spreading middles and slouching shoulders potbellies such complicated loops are often interpreted in oversimplified ways leading to erroneous conclusions the physics of the problem have been understood for nearly half a century yet numerical simulations appropriate to geological materials are almost unavailable in this paper we discuss results of numerical simulations using the simplest of systems the single domain superparamagnetic sd sp system examination of the synthetic hysteresis loops leads to the following observations 1 wasp waisting and potbellies can easily be generated from populations of sd and sp grains 2 wasp waisting requires an sp contribution that saturates quickly resulting in a steep initial slope and potbellies require low initial slopes the sp contribution approaching saturation at higher fields the approach to saturation is dependent on volume hence the cube of grain diameter therefore there is a very strong dependence of hysteresis loop shape on the assumed threshold size 3 we were unable to generate potbellies using an sp sd threshold size as large as 30 nm and wasp waists cannot be generated using a threshold size as small as 8 nm the occurrence of both potbellies and wasp waists in natural samples is consistent with a room temperature threshold size of some 15 nm 5 4 simulations using a threshold size of 15 20 nm with populations dominated by sp grain sizes that is with a small number of sd grains produce synthetic hysteresis loops consistent with measured hysteresis loops and transmission electron microscopic observations from submarine basaltic glass 5 simulations and measurements using two populations with distinct coercivity spectra can also generate wasp waisted loops a relatively straightforward analysis of the resulting loops can distinguish the latter case from wasp waisting resulting from sp sd behavior

the dependence of traction evolution on the earthquake source time function adopted in kinematic rupture models 1 we compute the temporal evolution of traction by solving the elasto dynamic equation and by using the slip velocity history as a boundary condition on the fault plane we use different source time functions to derive a suite of kinematic source models to image the spatial distribution of dynamic and breakdown stress drop strength excess and critical slip weakening distance d c our results show that the source time functions adopted in kinematic source models affect the inferred dynamic parameters the critical slip weakening distance characterizing the constitutive relation ranges between 30 and 80 of the total slip the ratio between d c and total slip depends on the adopted source time functions and in these applications is nearly constant over the fault we propose that source time functions compatible with earthquake dynamics should be used to inter the traction time history

dynamic source inversion of the m6 5 intermediate depth zumpango earthquake in central mexico a parallel genetic algorithm we introduce a method for imaging the earthquake source dynamics from the inversion of ground motion records based on a parallel genetic algorithm the source model follows an elliptical patch approach and uses the staggered grid split node method to simulate the earthquake dynamics a statistical analysis is used to estimate errors in both inverted and derived source parameters synthetic inversion tests reveal that the average rupture speed vr the rupture area and the stress drop may be determined with formal errors of 30 12 and 10 respectively in contrast derived parameters such as the radiated energy er the radiation efficiency r and the fracture energy g have larger errors around 70 40 and 25 respectively we applied the method to the mw 6 5 intermediate depth 62 km normal faulting earthquake of 11 december 2011 in guerrero mexico inferred values of 29 2 6 2 mpa and r 0 26 0 1 are significantly higher and lower respectively than those of typical subduction thrust events fracture energy is large so that more than 73 of the available potential energy for the dynamic process of faulting was deposited in the focal region i e g 14 4 3 5 1014j producing a slow rupture process vr vs 0 47 0 09 despite the relatively high energy radiation er 0 54 0 31 1015 j and energy moment ratio er m0 5 7 10 5 it is interesting to point out that such a slow and inefficient rupture along with the large stress drop in a small focal region are features also observed in both the 1994 deep bolivian earthquake and the seismicity of the intermediate depth bucaramanga nest

laboratory earthquakes along inhomogeneous faults directionality and supershear we report on the experimental observation of spontaneously nucleated ruptures occurring on frictionally held bimaterial interfaces with small amounts of wave speed mismatch rupture is always found to be asymmetric bilateral in one direction rupture always propagates at the generalized rayleigh wave speed whereas in the opposite direction it is subshear or it transitions to supershear the lack of a preferred rupture direction and the conditions leading to supershear are discussed in relation to existing theory and to the earthquake sequence in parkfield california and in north anatolia

the 2002 denali fault earthquake alaska a large magnitude slip partitioned event the mw moment magnitude 7 9 denali fault earthquake on 3 november 2002 was associated with 340 kilometers of surface rupture and was the largest strike slip earthquake in north america in almost 150 years it illuminates earthquake mechanics and hazards of large strike slip faults it began with thrusting on the previously unrecognized susitna glacier fault continued with right slip on the denali fault then took a right step and continued with right slip on the totschunda fault there is good correlation between geologically observed and geophysically inferred moment release the earthquake produced unusually strong distal effects in the rupture propagation direction including triggered seismicity

fabric analysis using the intercept method abstract the intercept method is a convenient yet powerful method for the numerical analysis of fabrics for geological applications it can be used when a set of objects e g pores grains or aggregates of a specific mineral lineaments etc with unique characteristics can be identified in an image the intercept method analyses boundaries of objects as a population of lines how intercepts are definined is therefore important in order to establish the significance of the intercept counts when the method is applied to digital images the subdivision of the information into square pixels must also be considered we apply a linear filter which minimizes the effects of grid anisotropy on the counting of intercepts analysis of the resulting rose of intercept counts by a fourier series considerably facilitates the interpretation of the data a rose of directions can readily be derived for the phase boundaries if a population of objects can be assumed to have started with an originally isotropic orientation and to have been deformed passively the fourier components of the rose of intercepts also give the sectional strain the fourier components also permit an objective test of the validity of describing a fabric as resulting from a homogeneous strain on an initially isotropic section two computer programs for the counting and analysis of intercepts on digital images are presented and applied in four examples which illustrate the usefulness of the intercept method for magmatic fabric measurement strain analysis and lineament analysis we also illustrate the importance of fourier series decomposition of rose of intercept counts for fabric interpretations intercept results are also shown to be compatible with those derived by other methods eigenvector surfor inertia tensor autocorrelation methods

magma flow directions in dikes of the koolau complex oahu determined from magnetic fabric studies we compared the magnetic fabric with the macroscopic surface lineation direction previously assumed to be parallel with magma flow direction for 35 dikes of the koolau complex oahu anisotropy of magnetic susceptibility ams was determined for 530 samples from 59 dikes including a sill and a pluglike intrusive body and two filled lava tubes and statistically significant ams clusters were found in all but five dikes for the 35 lineated dikes the mean maximum ams direction 1 coincides with the macroscopic lineations to within 25 in 18 and to within 15 in eight the ams ellipsoid shape in about half of the samples is prolate cigar shaped and in about half is oblate pancake shaped we interpret the ams to be determined by shape anisotropy of ferromagnetic grains the mean total anisotropy is 2 2 and total anisotropy h ranges from 0 24 to as high as 14 12 the principal ams axes for samples from a dike are tightly grouped on equal area plots and for most of the dikes the minimum ams axes plot near to the pole of the dike our ams measurements enable us to infer absolute magma flow directions for half the dikes studied the maximum ams axes of samples from one side of the dike form a tight cluster on an equal area plot and those from the other side of the dike form another tight cluster the two clusters are close together separated by 10 30 for narrow dikes 100 cm wide and are farther apart 30 60 for wide dikes and are symmetrically disposed on either side of the plane of the dike we infer that these paired clusters are caused by imbrication of nonequant ferromagnetic crystals against each margin of the dike during deposition in a velocity gradient and hence indicate the flow azimuth i e the absolute magma flow direction nearly half of the dikes measured 25 gave paired ams plots in 12 the magma flowed up at low oblique angles 5 30 from the horizontal in five it flowed down at 5 20 in four it flowed horizontally 5 in three it flowed up at a steep angle 50 75 and in one it flowed down at a steep angle the magma flow in general was rarely vertical or horizontal

dynamic faulting studied by a finite difference method we have developed a finite difference method that is especially adapted to the study of dynamic shear cracks we studied a number of simple earthquake source models in two and three dimensions with special emphasis on the modeling of the stress field we compared our numerical results for semi infinite and self similar shear cracks with the few exact solutions that are available in the literature we then studied spontaneous rupture propagation with the help of a maximum stress criterion from dimensional arguments and a few simple examples we showed that the maximum stress criterion depended on the physical dimensions of the fault for a given maximum stress intensity the finer the numerical mesh the higher the maximum stress that had to be adopted a study of in plane cracks showed that at high rupture velocities the numerical results did not resolve the stress concentration due to the rupture front from the stress peak associated with the shear wave propagating in front of the crack we suggest that this is the reason why transonic rupture velocities are found in the numerical solutions of in plane faulting when the rupture resistance is rather low finally we studied the spontaneous propagation of an initially circular rupture two distinct modes of nucleation of the rupture were studied in the first a plane circular shear crack was formed instantaneously in a uniformly prestressed medium after a while once stress concentrations had developed around the crack edge the rupture started to grow in the second type of nucleation a preexisting circular crack became unstable at time t 0 and started to grow the latter model appeared to us as a more realistic simulation of earthquake triggering in this case the initial stress was nonuniform and was the static field of the preexisting fault

shakeout d ground motion estimates using an ensemble of large earthquakes on the southern san andreas fault with spontaneous rupture propagation 1 we simulate ground motion in southern california from an ensemble of 7 spontaneous rupture models of large mw7 8 northwest propagating earthquakes on the southern san andreas fault shakeout d compared to long period spectral accelerations from the next generation attenuation nga empirical relations shakeout d predicts similar average rock site values i e within roughly their epistemic uncertainty but significantly larger values in los angeles and ventura basins due to wave guide focusing effects the shakeout d ground motion predictions differ from those of a kinematically parameterized geometrically similar scenario rupture 1 the kinematic rock site predictions depart significantly from the common distance attenuation trend of the nga and shakeout d results and 2 shakeout d predictions of long period spectral acceleration within the basins of the greater los angeles area are lower by factors of 2 3 than the corresponding kinematic predictions we attribute these differences to a less coherent wavefield excited by the complex rupture paths of the shakeout d sources

crustal structure and magmatism of north atlantic continental margins enormous variations exist along rifted margins of the north atlantic in the magmatic activity that accompanied continental breakup the north atlantic margins exhibit a full range of magmatic behaviour one end of the spectrum commonly termed volcanic margins is found in the northern north atlantic off greenland norway and northwest britain where huge volumes of igneous rock were added to the continental crust as it rifted and the adjacent oceanic crust is considerably thicker than normal the other end of the spectrum so called non volcanic margins occurs further south off france and spain where only minor volcanism accompanied rifting and the oceanic crust immediately adjacent to the continental margin is thinner than normal oceanic crust results from seismic profiles are used to determine the volume and distribution of igneous crust accumulated as the continents broke up and seafloor spreading commenced the volume of melt produced during rifting is controlled mainly by the temperature of the underlying asthenospheric mantle and by the amount and rate of decompression as it rises beneath stretching and thinning lithosphere large quantities of igneous rock were produced as the continents broke up to form the north atlantic where the newly initiated iceland plume brought abnormally hot mantle beneath the rift much smaller quantities of melt were produced further south where mantle temperatures were normal during rifting the distribution of igneous rocks across the continental margin and the adjacent oceanic crust is controlled by the rate at which melt is produced and by the ease with which it can intrude laterally in the crust where the rate of stretching is very slow and the rifting is spread over many millions or even tens of millions of years as is common during the initial stages of rifting then the rising asthenospheric mantle has time to cool by conduction heat is lost from the asthenospheric mantle both vertically into the overlying water layer and laterally into the colder wedges of continental lithosphere on either side of the rifted region as a consequence of the reduced mantle temperatures considerably less partial melting occurs in the slowly rising asthenospheric mantle than in mantle which rises rapidly beneath rifts this explains the abnormally thin oceanic crust found adjacent to the non volcanic rifted continental margins in the north atlantic as the new ocean widens the spreading centre becomes isolated from the continental margins on either side by newly formed warm oceanic lithosphere and the rate of rifting generally increases the oceanic crust then attains its equilibrium thickness as conductive heat loss from the rising asthenospheric mantle reduces and all the melt formed by decompression of the mantle solidifies close to the spreading axis basaltic melts can intrude continental crust far more easily than oceanic crust basaltic dykes may extend many hundreds of kilometres away from the magma source beneath continental rifts this redistributes melt laterally towards the continent in both volcanic and to a lesser extent in non volcanic margins and may be a contributory factor in leaving unusually thin crust in the first formed ocean immediately adjacent to some rifted margins it also allows melt formed under volcanic rifts to migrate large distances laterally as is seen in the rockall area northwest of britain

photo acoustic study of subshear and supershear ruptures in the laboratory abstract we report supershear ruptures observed in spontaneously nucleating laboratory microearthquakes and describe the signature of the associated mach wavefront radiation transducers detect the wavefield both close and at a distance from the fault the rupture velocities are inferred from either photoelastic high speed imaging or using the acoustic recordings both methods yield compatible estimates of sufficient accuracy to discriminate between sub shear and supershear ruptures the acoustic records allow to characterize the mach wavefront radiated from the supershear rupture front in particular its amplitude and decay velocity functions recovered by integrating high frequency accelerometer signals in the case of supershear ruptures consist in a double pulse function a first pulse traveling at supershear velocity followed by a second pulse traveling at the rayleigh wave velocity conversely the sub shear event is characterized by a single pulse finally we perform numerical simulations of our experiment using a prescribed supershear rupture velocity the synthetic waveforms obtained from these simulations including the mach wave amplitude and phase yield a satisfactory fit to the experimental results

evidence for a supershear transient during the 2002 denali fault earthquake elastodynamic considerations suggest that the acceleration of ruptures to supershear velocities is accompanied by the release of rayleigh waves along the fault from the stress breakdown zone these waves generate a secondary slip pulse trailing the rupture front but manifest almost entirely in ground motion perpendicular to the fault in the near source region we construct a spontaneously propagating rupture model exhibiting these features and use it to explain ground motions recorded during the 2002 denali fault earthquake at pump station 10 located 3 km from the fault we show that the initial pulses on both the fault normal and fault parallel components are due to the supershear stress release on the fault whereas the later arriving fault normal pulses result from the trailing subshear slip pulse on the fault online material mpeg movies of rupture history and ground motion

seaward dipping reflectors and the continent ocean boundary at passive continental margins abstract the definition of the continent ocean boundary at passive continental margins has proved to be an elusive task even the relatively direct method of seismic refraction experiments has yielded results that cannot always be interpreted unequivocally multichannel seismic reflection profiles on many passive margins have revealed the presence of remarkable suites of arcuate reflectors dipping seaward to form a wedge shaped structure their general characteristics and velocity structure suggest that they may be largely volcanic in nature but there is no agreed upon model for their origin nevertheless it is generally thought that they lie at or close to the boundary between continent and ocean and as such they are extremely important structural markers that may yield important evidence concerning the structure and evolution of passive margins

constitutive relations between dynamic physical parameters near a tip of the propagating slip zone during stick slip shear failure abstract constitutive relations between physical parameters in the cohesive zone during stick slip shear failure are experimentally investigated stick slip was generated along a 40 cm long precut fault in tsukuba granite samples using a servocontrolled biaxial loading apparatus dynamic behavior during local breakdown processes near a tip of the slipping zone is revealed the slip velocity and acceleration are given as a function of the slip displacement and the cohesive or breakdown shear stress as a function of the slip velocity a cycle of the breakdown and restrengthening process of stick slip is composed of five phases characterized in terms of the cohesive strength and the slip velocity the cohesive strength can degrade regardless of the slip velocity during slip instabilities the maximum slip acceleration u max and the maximum slip velocity u max are obtained experimentally as u max 2 u c u max 2 and u max b g v where u c is the critical displacement b the breakdown stress drop g the rigidity and v the rupture velocity these relations are consistent with ida s theoretical estimation based on the cohesive zone model the above formula gives good estimates for the maximum slip acceleration of actual earthquakes the cutoff frequency max of the power spectral density of the slip acceleration increases with increasing normal stress in particular max is found to be directly proportional to the normal stress n within the normal stress range less than 17 mpa as max khz 4 0 n mpa n max increase with an increase in u max or u max all these results lead to the conclusion that u max u max and max increase with increasing normal stress this is consistent with a previous observation that b increases with increasing normal stress the above empirical linear relation between max and n can be explained by a linear dependence of b on n the size scale dependence of physical parameters is discussed and such size scale dependent parameters as u max and apparent fracture energy are scaled to the breakdown zone size regarded as a characteristic length to extend the results obtained in the laboratory to an earthquake failure in the earth

estimation of dynamic rupture parameters from the radiated seismic energy and apparent stress we propose a new procedure for estimating the critical slip weakening distance dc by evaluating the spatio temporal distribution of the apparent stress calculated from the spatio temporal distribution of the slip velocity function on the fault plane obtained from a kinematic inversion of the earthquake source the idea is based on the fact that the apparent stress can be related to the difference between the earthquake average stress and the frictional dynamical stress on the fault plane during the rupture process from the cumulative slip apparent stress relationship we estimate the critical slip dc of the 1992 landers earthquake we find that the critical slip weakening plays an important role in controlling the rupture velocity in the case of the landers earthquake the very low rupture velocity in the central segment could be explained by the large dc value obtained for the main asperity on that segment

magma flow directions of shallow dykes from the east greenland volcanic margin inferred from magnetic fabric studies abstract the role played by plume generated crustal magmatic complexes in the segmentation of volcanic margins is highlighted by a preliminary study of magma flow directions in shallow intrusives from the east greenland volcanic margin we investigate the magmatic texture using anisotropy of magnetic susceptibility for eight dykes of tholeitic affinity belonging to a dyke swarm associated with the tertiary opening of the north atlantic ocean the thickness of the sampled dykes ranges from 3 to 37 m the dykes are of doleritic texture and contain up to 12 of opaque minerals and 35 of plagioclase laths dykes showing a magnetic foliation plane within the dyke plane i e a magnetic fabric usually attributed to magmatic processes represent 40 of the samples from the studied swarm these dykes have a low degree of anisotropy and their ellipsoid of magnetic susceptibility is strongly oblate inverse magnetic fabrics where the maximum principal susceptibility axis lies near the pole to the dyke account for 60 of the data we interpret these inverse fabrics either as alteration minerals with highly prolate ellipsoids of magnetic susceptibility or primary titanomagnetites with oblate ellipsoid of low anisotropy the flow direction is inferred for normal fabric dykes using the mirror imbrication of magnetic lineation analysis of thin sections shows a good agreement between magnetic fabric directions and phenocryst preferred orientations the inferred flow directions are predominantly horizontal throwing a new light on volcanic margin development

dynamic rupture modeling of the transition from thrust to strike slip motion in the 2002 denali fault earthquake alaska we use three dimensional dynamic spontaneous rupture models to investigate the nearly simultaneous ruptures of the susitna glacier thrust fault and the denali strike slip fault with the 1957 m w 8 3 gobi altay mongolia earthquake as the only other well documented case of significant nearly simultaneous rupture of both thrust and strike slip faults this feature of the 2002 denali fault earthquake provides a unique opportunity to investigate the mechanisms responsible for development of these large complex events we find that the geometry of the faults and the orientation of the regional stress field caused slip on the susitna glacier fault to load the denali fault several different stress orientations with oblique right lateral motion on the susitna glacier fault replicate the triggering of rupture on the denali fault about 10 sec after the rupture nucleates on the susitna glacier fault however generating slip directions compatible with measured surface offsets and kinematic source inversions requires perturbing the stress orientation from that determined with focal mechanisms of regional events adjusting the vertical component of the principal stress tensor for the regional stress field so that it is more consistent with a mixture of strike slip and reverse faulting significantly improves the fit of the slip rake angles to the data rotating the maximum horizontal compressive stress direction westward appears to improve the fit even further

southeast baffin volcanic margin and the north american greenland plate separation plate breakup over plumes is characterized by the development of margins showing extensive magma production both underplated at moho level and extruded as thick piles of seaward dipping lava formations the disko svartenhuk area in west greenland is one of the few places in the world with exposed seaward dipping basalts forming a prism whose thickness increases seaward we present a quantitative tectonic study of this margin which we tentatively restored in its geodynamic position during the different stages of plate separation between greenland and north america our structural data are constrained with recently published 39ar 40ar and new and coherent 40k 40ar geochronology in dikes of different orientations the first order structure is that of a tectonically driven seaward crustal flexure linked to definitive plate breakup between greenland and baffin island during the eocene coeval with the formation of the upper part of the exposed seaward dipping volcanic prism this flexing is associated with a significant crustal stretching associated with arrays of continentward dipping normal faults this across strike structure is correlated to a fundamental along strike segmentation with the three segments adopting a zigzag strike there is a clear increase of extensional strain at the extremities of the segments eocene extension trended n060 on average in the northern and southern segments but was nw trending in the central ne trending nugssuaq segment we discuss the interpretation of such an extension perpendicular to the different margin segments from a regional point of view the n060 extension is distinct in orientation from the approximately n s trend of plate separation between north america and greenland in the baffin bay during the eocene however the extension recorded in the margin is nearly perpendicular to the obliquely spreading eocene accretion axis in the baffin bay this latter point suggests a mantle or ridge control of the development of the margin over regional far field lithospheric stress models

how fast is rupture during an earthquake new insights from the 1999 turkey earthquakes we report that during the two devastating 1999 earthquakes in turkey rupture propagated over a large part of the nearly 200km long fault zone at supershear speed approaching 5km s we present observations and modeling which confirm the original inference of supershear rupture during the izmit earthquake and we show that supershear rupture also occurred during the duzce earthquake we show that the rupture velocity measured about 2 times the shear wave velocity is the value predicted by theoretical studies in fracture dynamics we look for clues to explain these observations

slip weakening distance estimated at near fault stations we estimated the slip weakening distance from the seismograms recorded at near fault stations considering the effect of spatio temporal filtering for a continuous propagation of rupture this effect is usually observed at some distance away from the fault we used a green s function for 2 d anti plane rupture propagating with constant rupture velocity and instantaneous slip at the crack tip as a source time function we used a regularized yoffe function which is a kinematically good approximation of dynamic slip weakening behavior by convolving these functions near fault waveforms are computed as a proxy of slip weakening distance d c we used d c which is defined as a slip at the time of peak slip velocity then d c values can be computed as a function of the distance from the fault this procedure is applied to the observation of the 2000 western tottori japan earthquake m w 6 6 as well as the 2002 denali alaska earthquake m w 7 9 we estimated d at about 0 3 m at gsh station 100 m away from the fault for the tottori earthquake and that at about 2 5 m at ps10 station 3 km away from the fault for the denali earthquake we confirmed that both estimates are not affected by the spatio temporal smoothing

dike surface lineations as magma flow indicators within the sheeted dike complex of the troodos ophiolite cyprus mesoscopic flow lineations and anisotropy of magnetic susceptibility ams have been measured for dikes within the cretaceous age troodos ophiolite with the goal of comparing the direction of initial magma flow through dike conduits immediately following crack propagation with that of flow of subsequent magma emplaced during later stages of dike growth dike margin indicators of flow include cusp axes and elongate vesicles found high in the ophiolite pseudostratigraphy and ridge and groove structures termed hot slickenlines found throughout the complex a unique flow direction is determined where elongate vesicles near dike margins display imbrication with respect to the margin significant changes in vesicle elongation directions across dikes likely indicate either changes in magma flow direction after dike propagation or back flow of magma during the waning stages of intrusion surface lineations generally lie subparallel to the direction of flow inferred from ams determinations on cores within 5 cm of dike margins surface lineations also lie subparallel to the long axis e1 of the orientation ellipsoid defined by long axes of groundmass plagioclase phenocrysts measured in sections from ams cores correlation of surface lineations with interior indicators of flow ams plagioclase trachytic texture indicate that the surface features are good proxies for grain scale magma flow directions during dike propagation in troodos dikes orientations of surface flow features in the dikes of the troodos ophiolite indicate an approximately equal mix of subhorizontal to near vertical magma flow contradicting the paradigm of primarily vertical flow of magma beneath continuous axial magma chambers at oceanic spreading centers our data are consistent with a model of magma emplacement both vertically and horizontally away from isolated magma chambers beneath axial volcanoes spaced along a ridge crest

shallow intrusive directions of sheeted dikes in the troodos ophiolite anisotropy of magnetic susceptibility and structural data sheeted dikes play a central role in the formation of oceanic crust it is commonly assumed that sheeted dikes intrude vertically upward from elongated mid ocean ridge mor magma chambers but there are no direct observational data bearing on this hypothesis this assumption contrasts with the intrusive behavior of subaerial volcanoes where magmas rise into shallow central magma chambers that laterally feed vertically oriented dikes the authors have studied intrusive directions of sheeted dikes in a structural analogue to oceanic crust the troodos ophiolite structural and magnetic fabric data of 65 dikes provide consistent results and suggest a broad distribution of shallow 20 degree to nearly vertical upward magma transport directions these data suggest that horizontal emplacement has to be considered for sheeted dikes at mors implying more centralized mor plumbing systems than previously thought such plumbing systems provide ample opportunity for complex mixing fractionation and contamination of mor lavas in magma chambers and tabular magma storage volumes whether the mor magma supply is linear or centralized also has a fundamental effect on crustal accretion processes and the geometry of hydrothermal convection systems

supershear mach waves expose the fault breakdown slip abstract mikumo et al 2003 showed that it is possible to estimate the breakdown slip dc as the slip at the time of the peak slip rate for rupture propagation with subshear speeds cruz atienza et al 2009 later attempted to extend this method to extract information about dc as the displacement at the time of the peak particle velocity from seismic strong motion records however a reasonably accurate estimate of dc was only possible in a narrow zone adjacent to the fault typically on the order of hundreds of meters due to the fast decay with distance from the fault of the seismic energy related to the stress breakdown process when the rupture propagates with supershear speeds on the other hand this energy is carried much farther away from the fault by mach waves in particular rayleigh mach waves when rupture reaches the earth s surface dunham and bhat 2008 here we present a new approach to estimate dc from strong motion records containing mach waves first we show that the method by mikumo et al is valid for supershear rupture propagation this method is then used to estimate dc via an asymptotic approximation of the slip and slip rate time histories from the mach waves using spontaneous rupture simulations we demonstrate that for a visco elastic half space model dc can be estimated within an error of 40 from mach waves that have propagated a distance of at least 3 km from the fault the method is applied to estimate dc for the 2002 mw7 9 denali alaska earthquake 1 5 m and for the 1999 mw7 6 izmit turkey earthquake 1 7 m

attenuation of radiated ground motion and stresses from three dimensional supershear ruptures 1 radiating shear and rayleigh waves from supershear ruptures form mach waves that transmit large amplitude ground motion and stresses to locations far from the fault we simulate bilateral ruptures on a finite width vertical strike slip fault of width w and half length l with l w breaking the surface of an elastic half space and focus on the wavefield out to distances comparable to l at distances much smaller than w two dimensional plane strain slip pulse models i e models in which the lateral extent of the slip zone is unbounded accurately predict the subsurface wavefield amplitudes in the shear mach wedge of those models are undiminished with distance from the fault when viewed from distances far greater than w rupture is accurately modeled as a moving point source that produces a shear mach cone and on the free surface rayleigh wave mach fronts geometrical spreading of the shear mach cone occurs radially and amplitudes there decrease with the inverse square root of distance the transition between these two asymptotic limits occurs at distances comparable to w similar considerations suggest that rayleigh mach waves suffer no attenuation in the ideally elastic medium studied here the rate at which fault strength weakens at the rupture front exerts a strong influence on the off fault fields only in the immediate vicinity of the fault for both sub rayleigh and supershear ruptures and at the mach fronts of supershear ruptures more rapid weakening generates larger amplitudes at the mach fronts

slip dependent friction law and nucleation processes in earthquake rupture abstract matsu ura m kataoka h and shibazaki b 1992 slip dependent friction law and nucleation processes in earthquake rupture in t mikumo k aki m ohnaka l j ruff and p k p spudich editors earthquake source physics and earthquake precursors tectonophysics 211 135 148 we developed a physical model that interprets slip strengthening and slip weakening processes in brittle shear fracture by considering the microscopic interaction between statistically self similar fault surfaces our model successfully describes general features of the constitutive behaviour observed in laboratory experiments of stick slip with the progress of fault slip frictional stress first increases rapidly up to a peak value and then decreases gradually to a constant level the first slip strengthening process can be explained by the increase of frictional resistance in transition from the in phase contact to the random phase contact of fault surfaces the latter slip weakening process can be explained by the decrease of frictional resistance due to the abrasion of surface asperities the trade off between these two processes determines the peak stress and the critical displacement over which the slip weakening proceeds we examined the nucleation process proceeding in a seismogenic zone with the increase of external shear stress through numerical simulations in our model the seismogenic zone is represented by a broad weak zone with a locally strong part asperity on a fault plane the basic equations governing the nucleation process are the equation of equilibrium and the constitutive relation between frictional stress and fault slip we obtained the following results the nucleation process first proceeds quasi statically at the weak portions with the increase of external shear stress then if the constitutive relation at the asperity has a high and narrow peak of stress brittle dynamic rupture occurs in a limited area around the asperity when the constitutive relation has a low and broad peak the dynamic rupture of the asperity does not occur the process is completely aseismic when the peak of stress is high and broad the dynamic rupture of the asperity leads to the catastrophic rupture which extends over the whole region from these results it is concluded that diversity in the nucleation process can be ascribed to difference in the variation of constitutive behaviour along the fault

terashake2 spontaneous rupture simulations of mw 7 7 earthquakes on the southern san andreas fault abstract previous numerical simulations terashake1 of large m w 7 7 southern san andreas fault earthquakes predicted localized areas of strong amplification in the los angeles area associated with directivity and wave guide effects from northwestward propagating rupture scenarios the terashake1 source was derived from inversions of the 2002 m w 7 9 denali alaska earthquake that source was relatively smooth in its slip distribution and rupture characteristics owing both to resolution limits of the inversions and simplifications imposed by the kinematic parameterization new simulations terashake2 with a more complex source derived from spontaneous rupture modeling with small scale stress drop heterogeneity predict a similar spatial pattern of peak ground velocity pgv but with the pgv extremes decreased by factors of 2 3 relative to terashake1 the terashake2 source excites a less coherent wave field with reduced along strike directivity accompanied by streaks of elevated ground motion extending away from the fault trace the source complexity entails abrupt changes in the direction and speed of rupture correlated to changes in slip velocity amplitude and waveform features that might prove challenging to capture in a purely kinematic parameterization despite the reduced pgv extremes northwest rupturing terashake2 simulations still predict entrainment by basin structure of a strong directivity pulse with pgvs in los angeles and san gabriel basins that are much higher than predicted by empirical methods significant areas of those basins have predicted pgv above the 2 probability of exceedance poe level relative to current attenuation relationships even when the latter includes a site term to account for local sediment depth and wave guide focusing produces localized areas with pgv at roughly 0 1 0 2 poe about a factor of 4 5 above the median in contrast at rock sites in the 0 100 km distance range the median terashake2 pgvs are in very close agreement with the median empirical prediction and extremes nowhere reach the 2 poe level the rock site agreement lends credibility to some of our source modeling assumptions including overall stress drop level and the manner in which we assigned dynamic parameters to represent the mechanical weakness of near surface material future efforts should focus on validating and refining these findings assessing their probabilities of occurrence relative to alternative rupture scenarios for the southern san andreas fault and incorporating them into seismic hazard estimation for southern california

magmatic tectonic evolution of a volcanic rifted margin abstract many north atlantic continental margins are underlain by huge volcanic edifices near the continent ocean boundary a deep crustal hole drilled on the voring plateau odp leg 104 has provided important constraints on the initial breakup history and the subsequent margin evolution by penetrating more than 900 m of igneous rocks and interbedded sediments below a post early eocene cover the rocks constitute two distinct different volcanic series the upper series comprising a seaward dipping reflector wedge consists of transitional mid oceanic tholeiitic lava flows and thin volcaniclastic sediments dacitic flows some dikes and thicker interbedded sediments constitute the lower series the margin evolved by paleocene crustal extension accompanied by uplift and pervasive intrusion of the rift zone just prior to breakup magma from shallow crustal melts produced the lower series the upper series was constructed during an intense rapidly waning subaerial volcanic surge following breakup near the paleocene eocene transition it covers both the new oceanic crust and large areas of the continental crust the dipping wedge was formed by subsidence due to loading and thermal contraction possibly amplified by tectonism when the surge abated the injection center rapidly subsided and normal oceanic crust was formed there is a direct temporal and compositional relationship between the onshore north atlantic volcanic province and the margin volcanism compared to the transverse part of the province including the iceland hotspot which has been active over 60 m y the volcanic margins reflect a 2500 km long transient segment related to continental breakup lasting only about 3 m y the volcanic surge and uplift at breakup are related to a modest regional increased temperature at the base of the lithosphere this mode of margin formation has had important paleoenvironmental implications

coherence of mach fronts during heterogeneous supershear earthquake rupture propagation simulations and comparison with observations 1 we study how heterogeneous rupture propagation affects the coherence of shear and rayleigh mach wavefronts radiated by supershear earthquakes we address this question using numerical simulations of ruptures on a planar vertical strike slip fault embedded in a three dimensional homogeneous linear elastic half space ruptures propagate spontaneously in accordance with a linear slip weakening friction law through both homogeneous and heterogeneous initial shear stress fields in the 3 d homogeneous case rupture fronts are curved owing to interactions with the free surface and the finite fault width however this curvature does not greatly diminish the coherence of mach fronts relative to cases in which the rupture front is constrained to be straight as studied by dunham and bhat 2008a introducing heterogeneity in the initial shear stress distribution causes ruptures to propagate at speeds that locally fluctuate above and below the shear wave speed calculations of the fourier amplitude spectra fas of ground velocity time histories corroborate the kinematic results of bizzarri and spudich 2008a 1 the ground motion of a supershear rupture is richer in high frequency with respect to a subshear one 2 when a mach pulse is present its high frequency content overwhelms that arising from stress heterogeneity present numerical experiments indicate that a mach pulse causes approximately an 1 7 high frequency falloff in the fas of ground displacement moreover within the context of the employed representation of heterogeneities and over the range of parameter space that is accessible with current computational resources our simulations suggest that while heterogeneities reduce peak ground velocity and diminish the coherence of the mach fronts ground motion at stations experiencing mach pulses should be richer in high frequencies compared to stations without mach pulses in contrast to the foregoing theoretical results we find no average elevation of 5 damped absolute response spectral accelerations sa in the period band 0 05 0 4 s observed at stations that presumably experienced mach pulses during the 1979 imperial valley 1999 kocaeli and 2002 denali fault earthquakes compared to sa observed at non mach pulse stations in the same earthquakes a 20 amplification of short period sa is seen only at a few of the imperial valley stations closest to the fault this lack of elevated sa suggests that either mach pulses in real earthquakes are even more incoherent that in our simulations or that mach pulses are vulnerable to attenuation through nonlinear soil response in any case this result might imply that current engineering models of high frequency earthquake ground motions do not need to be modified by more than 20 close to the fault to account for mach pulses provided that the existing data are adequately representative of ground motions from supershear earthquakes

the 14 november 2001 kokoxili tibet earthquake high frequency seismic radiation originating from the transitions between sub rayleigh and supershear rupture velocity regimes seismic array based analysis of the major kokoxili earthquake tibet 14 november 2001 yields an unambiguous reconstruction of the seismic rupture history and relates it to the generated seismic radiation we demonstrate that after a classical sub rayleigh velocity stage the rupture speed has jumped to supershear values close to compressional wave velocity over a 175 km long fault segment before abruptly slowing down in the late part of the earthquake the transition locations between these three phases are correlated with the fault geometry and are associated with the most energetic radiation this observation proves that the rupture velocity changes as theoretically predicted are a primary source of high frequency seismic radiation this result requires reconsidering the origins of seismic damage generally attributed to slip variations

estimation of the breakdown slip from strong motion seismograms insights from numerical experiments fukuyama and mikumo 2007 presented a method to extract informa tion about the friction controlling earthquake rupture directly from near field seismo grams recorded out to several kilometers from a fault this method computes a parameter d 00 approximating the breakdown slip dc from the rake parallel dis placement at the time of the peak particle velocity the validation of the method was based on simple 2d steady state green s functions an approach lacking sufficient physics of the rupture process to support their conclusions here we use 3d simula tions of subshear strike slip spontaneous rupture propagation to demonstrate that d 00 is almost always controlled by rupture and wave propagation effects rather than the fault friction process only if rupture reaches the earth s surface and within a short distance from the fault rc is it possible to extract information about dc from strong motion data due to the fast decay with distance from the fault of the seismic energy related to the stress breakdown process beyond rc d 00 is controlled by the dynamic stress drop from the earthquake rupture without information about dc rc is comparable to the length of the fault cohesive zone where the breakdown process takes place during rupture and approximately equal to 80 of the wavelength associated with the break down frequency defined as the reciprocal of the time span required by the stress to drop to the dynamic level our findings suggest that the dc estimate by fukuyama and mikumo 2007 for the 2000 mw 6 6 tottori earthquake is unrelated to the breakdown slip moreover the very narrow widths for rc obtained from our simulations are in agreement with previous studies challenging the resolution of bandlimited strong motion data to estimate dc from kinematic rupture models guatteri and spudich 2000 spudich and guatteri 2004

north atlantic volcanic margins dimensions and production rates early tertiary lithospheric breakup between eurasia and greenland was accompanied by a transient 3 m y igneous event emplacing both the onshore flood basalts of the north atlantic volcanic province navp and huge extrusive complexes along the continent ocean transition on the rifted continental margins seismic data show that volcanic margins extend 2600 km along the early eocene plate boundary in places underlain by high velocity 7 2 7 7 km s lower crustal bodies quantitative calculations of navp dimensions considered minimum estimates reveal an areal extent of 1 3 106 km2 and a volume of flood basalts of 1 8 106 km3 yielding a mean eruption rate of 0 6 km3 yr or 2 4 km3 yr if two thirds of the basalts were emplaced within 0 5 m y the total crustal volume is 6 6 106 km3 resulting in a mean crustal accretion rate of 2 2 km3 yr thus navp ranks among the world s larger igneous provinces if the volcanic margins are considered the velocity structure of the expanded crust seaward of the continent ocean boundary differs from standard oceanic and continental crustal models based on seismic velocities this volcanic margin crust can be divided into three units of which the upper unit corresponds to basaltic extrusives the regionally consistent velocity structure and geometry of the crustal units suggest that the expanded crust including the high velocity lower crust which extends some distance landward of the continent ocean boundary was emplaced during and subsequent to breakup the volcanic margin crust was formed by excess melting within a wide zone of asthenospheric upwelling probably reflecting the interaction of a mantle plume and a lithosphere already extending

determination of flow directions by combining ams and thin section analyses implications for oligocene volcanism in the kerguelen archipelago southern indian ocean we report an anisotropy of magnetic susceptibility ams study carried out on seven oligocene volcanic sections widely distributed on the surface of the kerguelen archipelago the statistical results at flow scale and at section scale were checked with ams density diagrams and thin section analysis because of the axis permutations frequently observed for lava flows two possible flow directions were nicely defined for the majority of the sections we then used two alternative approaches imbrication of the minimum axis and symmetry of the ams distribution to infer a flow direction for each section these two methods and thin section analysis not being decisive for choosing between the two solutions geological and geomorphological observations were finally used to obtain confident interpretations the results of this study enable the supposed eruptive centres of the kerguelen archipelago to be located more precisely

is this magnetic fabric normal a review and case studies in volcanic formations abstract when studying weakly anisotropic rocks non standard relationships between magnetic and expected petrofabric axes are sometimes observed in particular the so called inverse or intermediate fabrics a review of possible explanations for these effects is presented and leads to the conclusion that a majority of abnormal magnetic fabrics often encountered in magnetite rich volcanic rocks are due to abnormal preferred orientations or distributions of magnetite grains through case studies in dikes from skye island yemen and the iberian peninsula and magma and ash flows from the ethiopian traps and southern corsica the various possible origins of such fabrics are examined it is possible to retrieve reliable magma flow directions from these sites provided proper data selection is performed leading to non unique interpretations

distribution anisotropy the cause of ams in igneous rocks understanding of ams in pristine basic igneous rocks in terms of flow orientation of ellipsoidal magnetite grains is thwarted by the fact that these grains 1 typically are irregularly equant or skeletal rather than distinctly ellipsoidal and 2 usually crystallize interstitially and late after magma flow has ceased however the distribution of such grains will be relatively anisotropic if they grew in residual liquid volumes within a preferredly oriented by flow silicate template this constraint may impose a slight preferred shape orientation to the irregular magnetite grains during their growth but in addition we suggest it is the anisotropic magnetic interaction between the grains which is the basic cause of ams experiments involving the casting of magnetite epoxy mixtures in various foliated or lineated glass templates have yielded samples with the expected ams symmetry i e k3 perpendicular to glass plates or k1 parallel to glass rods slicing of near isotropic synthetic blank and natural rock cores and interleaving with glass slides showed a progressive increase in anisotropy of the blank from 0 4 to 2 and the natural rock from 0 9 to 2 9 in both the k3 axis migrates to be perpendicular to the glass slides

an alternative emplacement model for the classic ardnamurchan cone sheet swarm nw scotland involving lateral magma supply via regional dykes abstract sub volcanic intrusive networks of which cone sheets are recognised as a major constituent partially control volcano growth and eruption style the accepted cone sheet model is that these confocally dipping intrusions originate from an unexposed central magma chamber through dip parallel magma flow however the emplacement mechanism of cone sheets remains poorly understood the 58 ma cone sheet swarm of ardnamurchan nw scotland offers an excellent opportunity to further resolve the emplacement dynamics of cone sheets through studying magma flow and their importance in volcanic edifice construction structural measurements and anisotropy of magnetic susceptibility ams analyses have constrained a lateral magma flow pattern consistently oriented nw se in the majority of the ardnamurchan cone sheets field observations also highlight the importance of host rock structure and interference between competing local and regional stress fields in controlling intrusion geometry our observations suggest cone sheet formation may be linked to laterally propagating nw se trending regional dykes sourced from laterally adjacent magmatic systems likely the palaeogene mull central complex which are deflected around the central complex by stress field interference implicitly edifice construction and potential eruption precursors observed at a volcano may instigate or result from magmatic activity in laterally adjacent volcanic systems

modeling dynamic rupture in a 3d earthquake fault model we propose a fourth order staggered grid finite difference method to study dynamic faulting in three dimensions the method uses an implementation of the boundary conditions on the fault that allows the use of general friction models including slip weakening and rate dependence because the staggered grid method defines stresses and particle velocities at different grid points we preserve symmetry by implementing a two grid row thick fault zone slip is computed between points located at the borders of the fault zone while the two components of shear traction on the fault are forced to be symmetric inside the fault zone we study the properties of the numerical method comparing our simulations with well known properties of seismic ruptures in 3d among the properties that are well modeled by our method are full elastic wave interactions frictional instability rupture initiation from a finite initial patch spontaneous rupture growth at subsonic and supersonic speeds as well as healing by either stopping phases or rate dependent friction we use this method for simulating spontaneous rupture propagation along an arbitrarily loaded planar fault starting from a localized asperity on circular and rectangular faults the shape of the rupture front is close to elliptical and is systematically elongated in the in plane direction of traction drop this elongation is due to the presence of a strong shear stress peak that moves ahead of the rupture in the in plane direction at high initial stresses the rupture front becomes unstable and jumps to super shear speeds in the direction of in plane shear another interesting effect is the development of relatively narrow rupture fronts due to the presence of rate weakening friction the solutions for the thick fault boundary conditions scale with the slip weakening distance do and are stable and reproducible for do greater than about 4 in terms of 2t t ax finally a comparison of scalar and vector boundary conditions for the friction shows that slip is dominant along the direction of the prestress with the largest deviations in slip rate direction occurring near the rupture front and the edges of the fault

effects of supershear rupture speed on the high frequency content of s waves investigated using spontaneous dynamic rupture models and isochrone theory 1 in this paper we achieve three goals 1 we demonstrate that crack tips governed by friction laws including slip weakening rate and state dependent laws and thermal pressurization of pore fluids propagating at supershear speed have slip velocity functions with reduced high frequency content compared to crack tips traveling at subshear speeds this is demonstrated using a fully dynamic spontaneous three dimensional earthquake model in which we calculate fault slip velocity at nine points locations distributed along a quarter circle on the fault where the rupture is traveling at supershear speed in the in plane direction and subshear speed in the antiplane direction this holds for a fault governed by the linear slip weakening constitutive equation by slip weakening with thermal pressurization of pore fluid and by rate and state dependent laws with thermal pressurization the same is also true even assuming a highly heterogeneous initial shear stress field on the fault 2 using isochrone theory we derive a general expression for the spectral characteristics and geometric spreading of two pulses arising from supershear rupture the well known mach wave and a second lesser known pulse caused by rupture acceleration 3 we demonstrate that the mach cone amplification of high frequencies overwhelms the de amplification of high frequency content in the slip velocity functions in supershear ruptures consequently when earthquake ruptures travel at supershear speed a net enhancement of high frequency radiation is expected and the alleged low peak accelerations observed for the 2002 denali and other large earthquakes are probably not caused by diminished high frequency content in the slip velocity function as has been speculated

rock magnetic and magnetic fabric variations across three thin alkaline dykes from sunnhordland western norway influence of initial mineralogy and secondary chemical alterations summary profiles of natural remanent magnetization nrm rock magnetic properties and anisotropy of magnetic susceptibility ams across three thin 26 to 109 cm alkaline dykes from sunnhordland western norway are presented the magnetic remanence is carried by pseudo single domain psd lowtemperature oxidized titanomagnetite with circa 6 per cent of cations substituted by aluminium systematic variations in magnetic properties across the dykes are caused by differences in titanomagnetite concentration as well as grain size lowtemperature alteration and laboratory oxidation by heating in air to a maximum temperature of 600 c do not homogenize the magnetic properties across the dykes fine grained regions near dyke margins carry the most stable remanent magnetization but are also most amenable to secondary chemical remanent magnetization crm resetting between dyke variations in magnetic properties are larger than within dyke variations a stable primary remanence is found in two of the dykes while the third is totally reset by a soft viscous remanent magnetization vrm despite similar optical grain sizes magnetic mineralogy and oxidation state the difference is attributed to a lower titanium content in the latter dyke preventing effective partitioning of magnetic grains during low temperature alteration laboratory heating affects anisotropy variations across the dykes resulting in more consistent axis orientations within each dyke and symmetric edge to edge variations in the degree of anisotropy ams directions are different from those in the country rock and vary systematically across the thinnest dyke before heating

dynamic modeling of the 1992 landers earthquake we have used observed band pass filtered accelerograms and a previously determined slip distribution to invert for the dynamic rupture propagation of the 1992 landers earthquake in our simulations dynamic rupture grows under the simultaneous control of initial stress and rupture resistance by friction which we modeled using a simple slip weakening law we used a simplified landers fault model where the fault segments were combined into a single vertical planar fault by trial and error we modified an initial stress field inferred from the kinematic slip distribution proposed by wald and heaton 1994 until dynamic rupture generated a rupture history and final slip distribution that approximately matched those determined by the kinematic inversion we found that rupture propagation was extremely sensitive to small changes in the distribution of prestress and that a delicate balance with energy release rate controls the average rupture speed for the inversion we generated synthetic 0 5 hz ground displacements using an efficient green s function propagator method axitra this method enables us to propagate the radiation generated by the dynamic rupture to distances greater than those feasible using the finite difference method the dynamic model built by trial and error inversion provides a very satisfactory fit between synthetics and strong motion data we validated this model using records from stations used in the slip inversion as well as some which were not included we also inverted for a complementary model that fits the data just as well but in which the initial stress was perfectly uniform while rupture resistance was heterogeneous this demonstrates that inversion of ground motion is nonunique

seismic imaging of the 1999 izmit turkey rupture inferred from the near fault recordings we use near fault accelerograms to infer the space time history of rupture on the fault during the izmit earthquake the records show that the ground displacement and velocity near the fault were surprisingly simple rupture propagated toward the west at a velocity of about 3 km s and toward the east at a remarkably high average velocity of 4 7 km s over a distance of about 45 km before decelerating to about 3 1 km s on the eastern segment slip on the fault is particularly large down to a depth of 20 km on the central portion of the fault where it reaches about 7 m slip is large also below 10 km on the eastern fault segment and this may have contributed to the loading of shear stress on the duzce fault on the western fault segment large slip seems confined to shallow depths

magnetic and shape fabrics of magnetite in simple shear flows abstract the magnetite fabrics measured by anisotropy of magnetic susceptibility ams and by shape preferred orientation spo optical methods are classically used as flow kinematics indicators in lava flows the development of magnetite fabrics during simple shear strains 20 was performed using a suspension of 1 volume fraction of multidomain magnetite randomly contained in a mixture of silicone and wax we measured ams fabric and spo ellipsoids by calculating a quadratic shape tensor from oriented thin sections for jeffery 1922 fabrics are usable to determine the flow kinematics including the amount of applied finite strain for 8 fabric elements foliation and lineation are stabilised closely parallel to the flow plane and the shear direction respectively two and three dimensional numerical simulations using measured aspect ratios of magnetite point out that the large scattering of aspect ratios and the initial orientation distribution of particles are together responsible for a wide ranging loss of periodicity the stable ams and spo fabrics observed at large strains in experiments are the result of these primary fabric properties combined to collisions between particles and possibly their complex three dimensional shapes in addition the constant angular relationship observed at large strains between fabrics and flow components is related to the transient collisions consequently the determination of the lava flow kinematics by using fabric properties measured either by ams or by spo analyses should be indubitably associated to a detailed study of the three dimensional shape of the solid carriers regularly shaped populations of low elongated particles will be capable to produce cyclic to oscillating fabrics while the fabric of elongated particles will be more sensitive to the shape parameters and collisions ultimately favouring stable fabrics at large strains

gravity and magnetic modelling of layered mafic ultramafic intrusions in large igneous province plume centre regions case studies from the 1 27 ga mackenzie 1 38 ga kunene kibaran 0 06 ga deccan and 0 13 0 08 ga high arctic events gravity and magnetic data from the global egm2008 and emag2 datasets are used to identify geophysical anomalies in large igneous province lip plume centre regions with the goal of characterizing mafic ultramafic intrusions linked to those lips geophysical anomalies within 18 lips distributed globally are investigated four of these lips are selected for detailed modelling the 1 27 ga mackenzie 1 38 ga kunene kibaran 0 06 ga deccan and 0 13 0 08 ga high arctic lips we recognize three spatial distribution types for intrusions in plume centre regions these are 1 intrusions emplaced along a circular fault system that circumscribes the plume centre 2 intrusions emplaced along linear rifts that in some cases converge towards the plume centre and 3 single or unclassified intrusions modelling supports that the geophysical anomalies associated with these lips tend to be produced by large radius 30 km and deep seated crustal intrusions with densities consistent with mafic ultramafic rock and

group therapy in crete a comparison between analyses by naa and thin section petrography of early minoan pottery abstract chemical analyses of minoan ceramics have divided the island of crete into three broad composition zones for the purposes of provenance work recent work by petrographic analysis has identified individual production centres within these compositional zone areas and a hitherto unsuspected broad distribution of their products also clear is the variation in choice and manipulation of raw material with traditions of paste technology over the island with a broad intra island exchange of pottery calling into question control groups sampled as being local to major bronze age sites a reassessment is needed this paper compares the chemical and mineralogical composition of stylistic and petrographic groups from the early bronze age of crete emphasizing the variety of factors which affect the grouping of pottery the primacy of the geographical scale of analysis and the interplay between provenance and technology

the aftershock signature of supershear earthquakes recent studies show that earthquake faults may rupture at speeds exceeding the shear wave velocity of rocks this supershear rupture produces in the ground a seismic shock wave similar to the sonic boom produced by a supersonic airplane this shock wave may increase the destruction caused by the earthquake we report that supershear earthquakes are characterized by a specific pattern of aftershocks the fault plane itself is remarkably quiet whereas aftershocks cluster off the fault on secondary structures that are activated by the supershear rupture the post earthquake quiescence of the fault shows that friction is relatively uniform over supershear segments whereas the activation of off fault structures is explained by the shock wave radiation which produces high stresses over a wide zone surrounding the fault

emplacement of dikes sills and crustal magma chambers at divergent plate boundaries abstract a model of dike emplacement at divergent plate boundaries is developed which predicts that below the uppermost 1 3 km of the oceanic crust the number and lengths of dikes in any particular swarm should increase but that the thickness should decrease with depth in the crust these predictions are supported by field observations it is proposed that rapid injections of dikes at divergent plate boundaries may temporarily alter the normal stress field in such a way that at a certain level in the crust subsequent dikes change into sills such sills may develop into crustal magma chambers either when many nearby sills combine into a large one or when an initial sill absorbs the magma of many dikes that enter it in a rapid succession while the sill is liquid for the fastest spreading ridges the initial thickness of a sill that is to develop into a crustal chamber need only attain 20 m whereas at the slowest spreading ridges the initial thickness must exceed 170 m because thin sills are much more likely to form than thick sills the probability of sills developing into crustal chambers increases with dike intrusion frequency hence with spreading rate furthermore the model predicts first that with increasing spreading rate the length of the chamber along the ridge axis should increase and second that the most likely site for magma chambers is at depths of 2 4 km both these predictions are supported by available data

near source ground motion from steady state dynamic rupture pulses 1 ground motion from two dimensional steady state dynamic ruptures is examined for both subshear and supershear rupture velocities synthetic seismograms demonstrate that coherent high frequency information about the source process rapidly attenuates with distance from the fault for subshear ruptures such records provide almost no resolution of the spatial extent of the stress breakdown zone at supershear speeds s waves radiate away from the fault preserving the full source spectrum and carrying an exact history of the slip velocity on both the fault parallel and fault normal components of motion whose amplitudes are given by a function of rupture speed that vanishes at times the s wave speed the energy liberated from the strain field by the passage of a supershear rupture is partitioned into fracture energy dissipated within the fault zone and far field s wave radiation the partition depends on both the rupture velocity and the size of the breakdown zone

stress breakdown time and slip weakening distance inferred from slip velocity functions on earthquake faults we estimate the critical slip weakening distance on earthquake faults by using a new approach which is independent of the estimate of fracture energy or radiated seismic energy the approach is to find a physically based relation between the breakdown time of shear stress tb the time of peak slip velocity tpv and the slip weakening distance dc from the time histories of shear stress slip and slip velocity at each point on the fault which can be obtained from dynamic rupture calculations using a simple slip weakening friction law numerical calculations are carried out for a dynamic shear crack propagating either spontaneously or at a fixed rupture velocity on a vertical fault located in a 3d half space and a more realistic horizontally layered structure with finite difference schemes the results show that tpv is well correlated with tb for faults even with a heterogeneous stress drop dis tribution except at locations near strong barriers and the fault edges we also inves tigate this relation for different types of slip weakening behavior we have applied the method to two recent strike slip earthquakes in western japan the 2000 tottori and the 1995 kobe events we integrated the slip velocity functions on the vertical fault obtained from kinematic waveform inversion of strong motion and teleseismic records from the arrival time of rupture tr to the time of the peak slip velocity tpv and we then corrected the slip obtained at tpv for the errors expected from the dynamic calculations it was found that the slip weakening dis tance d c estimated in the frequency window between 0 05 and 0 5 hz ranges between 40 and 90 cm on the two earthquake faults however if we consider the limited frequency resolution of the observed waveforms probable time errors in the slip velocity functions obtained from kinematic inversion and the uncertainty of the slip weakening behavior the above estimates may be those located between the minimum resolvable limit and the upper bound of their real values the estimated dc values do not necessarily seem to indicate larger values in the shallower part and smaller values in the deeper part of the fault but rather a spatially heterogeneous distribution that appears to be dependent on the local maximum slip this possible dependence might be interpreted by the frictional properties of the fault such as the degree of roughness or the thickness of gouge layers in addition to stress heterogeneities

magma emplacement into the lemptegy scoria cone chaine des puys france explored with structural anisotropy of magnetic susceptibility and paleomagnetic data the lemptegy volcano is a small monogenetic scoria cone located in the chaine des puys auvergne france which erupted about 32 000 years ago a first edifice lemptegy 1 formed during a trachybasalt eruption as a group of satellite vents of the puy de gouttes scoria cone a second trachyandesitic edifice lemptegy 2 formed soon after and completely covered lemptegy 1 with an 80 m high breached cone since 1946 the lemptegy volcano has been quarried for scoria and today offers unprecedented three dimensional exposure of the subvolcanic plumbing system to map the internal flow architecture of the plumbing system and to study the subvolcanic deformation of lemptegy 2 structural mapping petrographic observations anisotropy of magnetic susceptibility ams rock magnetic and paleomagnetic data were collected field structural mapping and thin section study of tension gashes riedel shears striations as well as ductile shear zones and bubbles allow the direction and sense of the magma flow to be determined twenty ams sites were established in ten dikes one to four sites in each dike with 504 specimens analyzed and 479 specimens used to infer magma flow patterns structural data the maximum susceptibility axis k 1 and the imbrication of the magnetic foliation k 1 k 2 planes indicate both upward and downward sense of flow as well as flow toward and away from the central vent rock magnetic experiments reveal that a cubic fe ti oxide phase likely low ti titanomagnetite is the principal magnetic phase carrying both the remanence and anisotropy paleomagnetic data from some sites yield statistically distinct at the 95 confidence level remanence directions while at other sites the data are indistinguishable at the 95 confidence level the paleomagnetic results observed steeply tilted scoria layers internal unconformities and faults show that as each dike was emplaced it displaced earlier dikes evidencing subvolcanic deformation the lemptegy 2 volcano shares similarities in terms of inferred eruption style and structures with other scoria cones such as cerro negro nicaragua and thus provides an excellent field laboratory to investigate active scoria cones world wide

the 1999 izmit earthquake sequence in turkey seismological and tectonic aspects the izmit duzce earthquake sequence is the largest sequence recorded by modern digital networks in the history of turkish earthquakes the m w 7 4 izmit mainshock occurred at 29 967 e and 40 729 n and 13 km depth at 00 hr 01 min 38 7 sec utc on 17 august 1999 and was preceded by at least eight immediate foreshocks of magnitude m l 3 0 aftershock locations show that a large region has become seismically very energetic following the mainshock and that most of the aftershock zone developed quickly in the first 3 days the distribution of the aftershocks illuminates numerous seismotectonic structures associated with the various branches of the north anatolian fault naf in the aftershock zone the seismicity pattern defines a rupture plane that extends about 80 km east and 75 km west of the mainshock epicenter confirming bilateral rupture during the mainshock the western aftershock activity is almost entirely located offshore extending along the axis of the izmit bay we show that the westernmost extent of the mainshock rupture reached the offshore region just 10 km south of the princes islands the dominance of strike slip motion is one of the main characteristics of the izmit sequence a number of significant aftershock clusters occur off the inferred mainshock rupture plane the spatial and temporal occurrence of these clusters show that they are triggered by the complex redistribution of stresses and dynamic strains imposed by the dislocation of the mainshock the focal mechanisms of the events located in the western end of the rupture zone show both strike slip and normal faulting mechanisms we infer that the yalova cluster was triggered by the mainshock the dominance of the normal faulting mechanisms supports the notion that dynamic extensional strains played a major role in triggering this cluster four events located off the princes islands indicate similar source mechanisms characterized by right lateral strike slip motion parallel to the strike of the northern boundary fault nbf and local trend of aftershock distribution this observation leads to a first order seismotectonic evidence that the nature of the nbf which is known to be a transtensional feature at the surface is dominated by strike slip motion at depth the coexistence of strike slip and extensional seismotectonic features in the sea of marmara highlights the importance of strain partitioning amongst shallow transtentional structures accommodating north south extension and deeper structures associated with the naf accommodating east west right lateral shear manuscript received 22 september 2000

40ar39ar geochronology of tertiary mafic intrusions along the east greenland rifted margin relation to flood basalts and the iceland hotspot track abstract the east greeland tertiary igneous province includes the largest exposed continental flood basalt sequence within the north atlantic borderlands more than ten layered gabbro complexes including the 55 ma skaergaard intrusion and a large dolerite sill complex are the plutonic equivalents of flood basalts both lavas and intrusions have been regarded as synchronous with continental breakup at 57 54 ma we report ten new ages of the mafic intrusions determined by40ar 39ar incremental heating experiments demonstrating that the mafic intrusions formed in two distinct time windows only intrusion ii of the imilik gabbro complex the skaergaard intrusion and the sorgenfri glestcher sill complex formed at 57 55 ma coeval with the eruption of regional flood basalts and continental breakup other layered gabbro intrusions at imilik intrusion iii kruuse fjord igtutarajik nordre aputite q kap edvard holm and lilloise are distinctly younger and formed between 50 and 47 ma plate kinematic models indicate the axis of the ancestral iceland mantle plume was located under central greenland at 60 ma and subsequently crossed the east greenland rifted continental margin we propose that tholeiitic magmatism along the east greenland rifted margin largely occurred in three distinct pulses at 62 59 ma lavas and dykes 57 54 ma lavas dykes sills and some gabbros and 50 47 ma gabbros dykes and rare lavas related to discrete mantle melting episodes triggered by plume impact continental breakup and passage of the plume axis respectively this model implies northwestward continental drift of greenland relative to the plume axis by 3 9 5 0 cm yr between 60 and 49 ma consistent with estimates from seismic studies of submerged flood basalts

evidence from gravity data for focused magmatic accretion along the mid atlantic ridge gravity data from the mid atlantic ridge between latitudes 27 50 n and 30 40 n show that the accretion of magma at the ridge is focused at discrete centres along the spreading axis large positive gravity anomalies equivalent to reductions of almost 50 in crustal thickness are observed over non transform discontinuities bounding spreading segments

evidence from magnetic fabric for the flow pattern of magma in the mackenzie giant radiating dyke swarm continental flood basalts and their associated swarms of mafic dykes are thought to result from melting in mantle plumes1 3 but the processes by which the resulting magma makes its way from mantle to crust remain obscure here we use magnetic anisotropy measurements to reconstruct the pattern of magma flow in the proterozoic mackenzie dyke swarm our results indicate that magma was injected vertically within 500 km of the focal point of the dyke swarm and then travelled horizontally at all further distances out to at least 2 100 km the sharp transition from vertical to horizontal flow at between 500 and 600 km may correspond to the outer boundary of melt generation in the mantle plume believed to be responsible for the mackenzie igneous events

near source ground motions from simulations of sustained intersonic and supersonic fault ruptures we examine the long period near source ground motions from simulations of m 7 4 events on a strike slip fault using kinematic ruptures with rupture speeds that range from subshear speeds through intersonic speeds to supersonic speeds the strong along strike shear wave directivity present in scenarios with subshear rupture speeds disappears in the scenarios with ruptures propagating faster than the shear wave speed furthermore the maximum horizontal displacements and velocities rotate from generally fault perpendicular orientations at subshear rupture speeds to generally fault parallel orientations at supersonic rupture speeds for rupture speeds just above the shear wave speed the orientations are spatially heterogeneous as a result of the random nature of our assumed slip model at locations within a few kilometers of the rupture the time histories of the polarization of the horizontal motion provide a better diagnostic with which to gauge the rupture speed than the orientation of the peak motion subshear ruptures are associated with significant fault perpendicular motion before fault parallel motion close to the fault supershear ruptures are associated with fault perpendicular motion after significant fault parallel motion consistent with previous studies we do not find evidence for prolonged supershear rupture in the long period 2 sec ground motions from the 1979 imperial valley earthquake however we are unable to resolve the issue of whether a limited portion of the rupture approximately 10 km in length propagated faster than the shear wave speed additionally a recording from the 2002 denali fault earthquake does appear to be qualitatively consistent with locally supershear rupture stronger evidence for supershear rupture in earthquakes may require very dense station coverage in order to capture these potentially distinguishing traits

the e greenland continental margin a transition between oceanic and continental magmatism the e greenland continental margin differs from typical margins in having a superimposed regional uplift which has exposed the early rift products for direct examination most of the igneous rocks in this region find counterparts in the oceanic setting with the proviso that the area is part of an aseismic ridge thus picrites basalts sheeted dyke swarms and layered gabbros are present and these can be used to shed light on oceanic magmatic events nephelinitic rocks are products of off axis volcanism and find no counterparts in the present oceanic rift system although they occur in marginal regions of the ocean basins and were generated in distal regions at the initiation of rifting in icelandic anomaly voluminous salic rocks in the area owe their generation largely to continental crust which has undergone extensive melting they find no direct counterparts in oceanic areas the outstanding degree of exposure in this area combined with its transitional nature and wide spectrum of rock types make it a unique laboratory for the elucidation of a variety of problems regarding both oceanic and continental volcanism together with the nature of rifted continental margins in general

magmatism at rifted continental margins the intense volcanism and uplift observed on many rifted continental margins forming basaltic seaward dipping reflector sequences is accompanied by the emplacement of a thick igneous section at depth partial melting by decompression of passively upwelling asthenosphere that is hotter than normal because it is near a hotspot explains both the thickness of the new igneous crust and the initial elevation of the rifted margins

space and time evolution of rupture and faulting during the 1999 izmit turkey earthquake we use the records of the ground motion obtained at near fault accelerometers to study the space and time evolution of rupture and faulting during the izmit earthquake we find that the rupture propagated at the sub rayleigh speed of about 3 km sec on the western and eastern segments of the fault but that the central segment izmit sapanca lake sakarya nearly 50 km long broke at the supershear speed of about 4 8 km sec this value within the range of uncertainties is the one theoretically predicted sqrt 2 v mathrm s in fracture dynamics for stable shear crack growth at intersonic speed we infer an average fault slip of about 2 9 m over a total rupture length of about 150 km with the largest values of up to 6 m occurring in the golcuk area to the west and in the sakarya region to the east the strong motion data also indicate that the slip diminished gradually to the west beyond the hersek peninsula over about 30 km whereas it stopped abruptly at depth at the termination of the eastern karadere segment the slip duration is between 2 and 4 sec except in the hypocentral area which slipped in about 1 sec the seismic moment inferred is about 2 5 10 20 n m

slip weakening distance in dynamic rupture of in slab normal faulting earthquakes we apply the above method to three large in slab normal faulting earthquakes in the mexican subduction zone to do this we calculate the spatial distribution of slip velocity functions and final slip from kinematic waveform inversion of strong motion and teleseismic records and the stress history and final stress change from dynamic rupture calculations by integrating the slip velocity functions obtained from the inversion from the rupture arrival time to the time of peak slip velocity we obtain slip dat t pv and then correct it for dc at t b through dynamic calculations we also estimate the lowest resolvable limit and probable errors of dc from the slip velocity functions and its upper bound from a theoretical constraint between the dynamic stress drop and dc we found that the slip weakening distance dc estimated in the frequency band between 0 05 and 0 5 hz ranges between 40 and 120 cm on the in slab fault of the 1999 oaxaca earthquake mw 7 5 the largest dc is detected in the central fault and in part of the deeper sections and dc in the zone around the hypocentre ranges between 50 and 70 cm the estimated dc values appear to be less depth dependent but are rather more dependent on the local maximum slip this possible slip dependence might be interpreted by the degree of fault roughness in addition to stress heterogeneities the fracture energy g in the central section and in the hypocentral zone are roughly estimated to be of the order of 10 15 and 5 8 mj m 2 respectively both of the estimated dc and g values are somewhat larger than those on the vertical fault of two recent shallow strike slip earthquakes in western japan

magmatic flow direction indicators in a giant diabase feeder dike connecticut recent studies have suggested that many large diabase dikes have been emplaced laterally rather than vertically these studies report only unidirectional flow which is usually based on measurements of anisotropy of magnetic susceptibility however in the 250 km long higganum dike of southern new england united states seven independent flow direction indicators record a complex history of forward and backward moving magma we would expect other similar size dikes to have had equally complex histories measurements within one en echelon segment of the higganum dike show that as the magma rose it spread laterally toward the ends of the segment and that backflow was in the reverse direction of the forward flow the higganum dike is interpreted to have formed from a linear array of intrusive fingers each finger spreading laterally to form an en echelon dike segment as it rose through the lithosphere

dike propagation and magma flow in a glassy rhyolite dike a structural and kinematic analysis exhumed magma conduits provide important evidence of the development and evolution of subvolcanic plumbing systems we use a 5 14 m thick flow banded rhyolite dike in arran scotland to present the first reconstruction of the directions and styles of initial propagation and subsequent magma flow based on mesoscale kinematic indicators the dike has concave inward dike margin segments with plumose like structures that record vertical and horizontal propagation of lobes which inflated and linked to form a through going sheet devitrified rhyolite zones at the dike margins show gentle to open folds in contrast glassy central parts of the dike are flow laminated and preserve folded and refolded isoclinal curvilinear folds and sheath folds that record sustained progressive deformation the inner interface between the glassy and devitrified facies is abrupt and marked by elongation lineations and mullions in the dike center fold axes plunge 27 ne along the dike and parallel to elongation lineations combined with shear sense indicators and objects sheared vesicles and asymmetric folds these features indicate that magma flow was obliquely upward to the southwest and locally 60 to the propagation direction of the dike the distribution of structures within the rhyolite indicates local accretion of the now devitrified material to the margins with localization of flow into the center of the dike we find that the initial magma flow direction was controlled by fracture propagation and interaction with the subsequent flow record controlled by accretion and flow localization in the conduit this study demonstrates that analysis of mesoscopic structural and kinematic features several of which have not previously been reported from dikes is a powerful tool that can be used to reconstruct the complex evolution of conduit initiation and magma flow processes

magnetic fabric and its significance in the 1400 ma mealy diabase dykes of labrador canada an anistropy of magnetic susceptibility ams study of the relatively fresh mealy diabase dyke swarm indicates a correlation between the ams and the deformational effects in the rocks several ams components were revealed during experimental work on 211 samples that had been heated to 640 c during a paleomagnetic study and 22 fresh samples progressively treated in alternating fields seventeen sites with slight or no deformation have a magnetic fabric composed of two major components the stable component attributed to a primary shape anisotropy has minimum and maximum principal susceptibility directions preferentially oriented vertically and along dyke strike n 60 e respectively this configuration is consistent with an origin due to a post emplacement stress field whereby horizontal extension acted on a vertical sheet of partially crystallized magma the orientation of this component with respect to the present horizontal suggests that the dykes have not been tilted since their emplacement the unstable component removed by fields up to 250 mt is attributed to induced domain anisotropy produced by mild tectonic stress during the grenville orogeny the approximate orientation of its principal components as found by the method of substitutional ellipsoids suggests an east west directed maximum stress four other sites with noticeable deformation have magnetic fabric that could not be reliably analyzed but which may be largely composed of domain and stress induced anisotropies developed during the grenville orogeny

magma flow in the east greenland dyke swarm inferred from study of anisotropy of magnetic susceptibility magmatic growth of a volcanic margin summary volcanic passive margins vpms are characterized by large volumes of melt emplaced within the lithosphere during break up processes several data and a recently proposed conceptual model of volcanic margin development suggest that vpms are fed from localized crustal zones of magma storage underlying large polygenetic volcanoes localized above diapir like instabilities of the asthenosphere we investigated the magma flow pattern within the coastparallel dyke swarm of the east greenland vpm which is the only outcropping vpm over a distance of 125 km the 44 sampled dykes are representative of the successive families of intrusions igneous petrofabrics are constrained by the measurements of the anisotropy of magnetic susceptibility the magnetic fabrics are of medium to low anisotropy p 0 flow related fabrics are recorded in 75 per cent of the sampled dykes we infer the flow directions from the imbrication geometry of the magnetic foliation planes at the dyke margins and check the results by measuring the preferred orientation of plagioclase in thin sections cut in the magnetic principal planes due to probable fabric superposition the magnetic lineation represents the zone axis for the distribution of magnetic foliation plane we obtained 23 reliable flow directions that are predominantly horizontal and directed away from identified crustal reservoirs this flow pattern supports the proposed model of vpm growth and emphasizes the localized nature of the magma sources in the mantle the entire flood basalt sequence appears to have been fed by a restricted number of crustal reservoirs and associated dyke swarms

3d finite difference dynamic rupture modeling along nonplanar faults properunderstandingofseismicemissionsassociatedwiththe growth of complexly shaped microearthquake networks and larger scale nonplanar fault ruptures both in arbitrarily heterogeneousmedia requiresaccuratemodelingoftheunderlyingdynamic processes we present a new 3d dynamic rupture finitedifference model called the finite difference fault element fdfe method it simulates the dynamic rupture of nonplanar faults subjected to regional loads in complex media fdfe is based on a 3d methodology for applying dynamic rupture boundary conditions along the fault surface the fault is discretizedbyasetofparallelepipedfaultelementsinwhichspecific boundary conditions are applied these conditions are applied to the stress tensor once transformed into a local fault reference frame numerically determined weight functions multiplying particle velocities around each element allow accurate estimates offaultkinematicparametersi e slipandsliprateindependent of faulting mechanism assuming a coulomb like slip weakeningfrictionlaw aparametricstudysuggeststhatthefdfemethod converges toward a unique solution provided that the cohesive zone behind the rupture front is well resolved i e four or more elements inside this zone solutions are free of relevant numerical artifacts for grid sizes smaller than approximately 70 m results yielded by the fdfe approach are in good quantitative agreement with those obtained by a semianalytical boundaryintegralmethodalongplanarandnonplanarparabola shaped faults thefdfemethodthusprovidesquantitative accurateresults for spontaneous rupture simulations on intricate fault geometries