backazimuthal variations of splitting parameters of teleseismic sks phases observed at the broadband stations in germany sks phases observed at broadband stations in germany show significant shear wave splitting we have analyzed sks and skks phases for shear wave splitting from 13 stations of the german regional seismic network grsn from 3 three component stations of the grafenberg array grf and from one austrian station sqta the data reveal strong differences in the splitting parameters fast direction and delay time t from a single event at various stations as well as variations at the individual stations for events with different backazimuths the backazimuthal variations of the splitting parameters at some stations can be explained by two layer anisotropy models with horizontal symmetry axes the best resolved two layer model is the gra1 model upper layer 40 t 1 15s lower layer 115 t 1 95s the upper layer can be attributed to the lithosphere because of the magnitude of the delay time of the upper layer the lower layer must lie within the asthenosphere at other stations splitting parameters are consistent with an anisotropic one layer model for the upper mantle stations near the bohemian massif show fast directions near ew throughout ne germany the directions are oriented nw se the reason for this direction is probably the nearby tornquist teisseyre line the observed fast axes are subparallel to this prominent transeuropean suture zone at stations in southern germany near the alps we observed ene wsw directions below some stations we also found indications of inclined anisotropic layers

isotropy or weak vertical transverse isotropy in d beneath the atlantic ocean 1 shear velocity properties of d beneath the central atlantic ocean are explored using predominantly european seismic recordings of intermediate and deep focus 100 km south american earthquakes broadband data are analyzed and when possible corrected for upper mantle models of receiver side anisotropic structure regional shear velocity heterogeneity in d is mapped by analysis of 306 s sks differential times that have been corrected for three dimensional seismic velocity structure above d using a whole mantle tomographic model this correction yields modest less than 2 estimates of seismic velocity heterogeneity in d with a transition from high to low seismic velocities traversing from west to east beneath the central atlantic in agreement with global tomographic models additionally shear wave splitting of s and sdiff for the same recordings was analyzed to assess seismic anisotropy in d the highest quality data provide 105 splitting times between sh and sv onsets that are mostly within the 1 s uncertainty level the few larger values generally exhibit sv delayed relative to sh assuming an anisotropy geometry involving vertical transverse isotropy vti as preferred in most regions of d that have been studied to date 0 5 anisotropy strength within a 100 km thick layer or 0 25 anisotropy within a 300 km thick layer are compatible with the data these values are low in comparison to those found in high velocity regions beneath the circum pacific ocean or in the low velocity region beneath the central pacific and many observations are in fact consistent with isotropic structure the lack of strong vti relative to other regions may be due to 1 the absence of stress from overlying midmantle downwelling 2 relatively weaker shear flow in the d boundary layer and or 3 lack of chemical heterogeneity that could develop either lattice preferred orientation or shape preferred orientation the azimuthal sampling of this region of d is quite limited thus the precise geometry and mechanism of any anisotropy are difficult to constrain it remains possible that this region may contain subtle azimuthal anisotropy that could couple the sv and sh signals however amplitude observations suggest that any such coupling is minor

d shear velocity heterogeneity anisotropy and discontinuity structure beneath the caribbean and central america the d region in the lowermost mantle beneath the caribbean and central america is investigated using shear waves from south american earthquakes recorded by seismic stations in north america we present a large scale composite study of volumetric shear velocity heterogeneity anisotropy and the possible presence of a d discontinuity in the region our data set includes 328 s sdiff sks differential travel times 300 scs s differential travel times 125 s sdiff and 120 scs shear wave splitting measurements and 297 seismograms inspected for scd the seismic phase refracted from a high velocity d layer broadband digital data are augmented by high quality digitized analog wwssn data providing extensive path coverage in our study area in all data from 61 events are utilized in some cases a given seismogram can be used for velocity heterogeneity anisotropy and discontinuity analyses significant mid mantle structure possibly associated with the ancient subducted farallon slab affects shear wave travel times and must be corrected for to prevent erroneous mapping of

double beam imaging mapping lower mantle heterogeneities using combinations of source and receiver arrays we present a new technique for imaging spatially distributed heterogeneities using information from combinations of source and receiver arrays the method is based on the single scattering assumption and is closely related to the double beam method of kruger et al 1993 1995 1996 in that it exploits amplitude delay time slowness and azimuth information in two arrays simultaneously a crucial step in the method is the application of static time corrections for a chosen reference phase here pcp one type of image is obtained as a spatial likelihood distribution for prescribed source and receiver array slowness and azimuth values and delay times with respect to a reference phase this allows the determination of the spatial origin of coherent phases in the p coda in a complementary approach we perform a double beam stack migration according to the theoretical slowness and azimuth values for candidate scatterers distributed over a three dimensional 3 d grid the scattering strength is expressed as the resulting beam power beam amplitude or semblance in a time window determined from the theoretical delay time with respect to the reference phase while likelihood mapping provides an image reflecting the probability of a region to explain observed kinematic phase properties by single scattering double beam migration provides information about the scattering strength the method is applied to nuclear explosion sources in eastern kazakhstan recorded at the yellowknife array in northern canada we identify several scattering volumes within the lower mantle below the arctic producing individual anomalous lower mantle phases in the p coda these anomalies range in depth from the core mantle boundary cmb up to 500 km into the lowermost mantle and are mainly located under the eurasian side of the arctic the joint interpretation of these results with the results of previous studies suggests a connection of these anomalies to the cenozoic and mesozoic subduction of the pacific plate and the kula plate the determination of the nature of such 3 d anomalies necessitates further array studies with better azimuthal coverage

shear wave splitting and subcontinental mantle deformation splitting observations are interpreted in terms of the strain induced lattice preferred orientation of mantle minerals especially olivine we consider three hypotheses concerning the origin of the continental anlsotropy 1 strain associated with absolute plate motion as in the oceanic upper mantle 2 crustal stress and 3 the past and present internal deformation of the subcontinental upper mantle by tectonic episodes it is found that the last hypothesis is the most successful namely that the most recent significant episode of internal deformation appears to be the best predictor of q for stable continental regions this is interpreted as fossil anisotropy whereas for presently active regions such as alaska the anisotropy reflects present day tectonic activity in the stable portion of north america there is a good correlation between delay time and lithospheric thickness this is consistent with the anisotropy being localized in the subcontinental lithosphere and suggests that intrinsic anisotropy is approximately constant the acceptance of this hypothesis has several implications for subcontinental mantle deformation first it argues for coherent deformation of the continental lithosphere crust and mantle during orogenies this implies that the anisotropic portion of the lithosphere was present since the deformational episode and rules out the addition of undeformed material to this layer by subsequent underplating or conductive growth of the thermal boundary layer one of the most important issues in the study of orogenies is the need to reconcile the formation of thickened lithosphere with the paradoxically high mantle temperatures often associated with orogenic episodes most efforts to date have focussed on modes of deformation whereby the cold lithospheric mantle is removed by convective instability or delamination and replaced by warm asthenosphere these models however are incompatible with the evidence for preserved coherent lithospheric deformation rather the deformed mantle appears to have been heated in place we suggest that the elevated mantle temperatures may be due to the strain heating accompanying the deformation

probing two low velocity regions with pkp b caustic amplitudes and scattering summary we probe low velocity regions of the lower mantle with pkp waves using south pacific earthquakes recorded at stations from the temporary ethiopia afar geoscientific lithospheric experiment array in ethiopia and the german regional seismic network grsn germany the ethiopian stations show strong scattering of pkpab waves which arrive as precursors to pkpdf they also show anomalously large amplitude diffractions from the pkp b caustic in contrast the events recorded at the grsn show scattering of pkpab as precursors to pkpdf but more normal amplitudes in the pkp b caustic diffractions to understand the cause for the observed amplitudes we test a large set of 1 d models with anomalies in p and s wave velocities and density in the lowermost mantle and outer core structures that can cause the observed large pkp b caustic amplitudes are low velocity regions just above the core mantle boundary of the order of 50 km thick for the south pacific to ethiopia path with velocity reductions of around 5 per cent for p and 25 per cent for the s wave velocity alternatively a 30 km thick lamella starting at 50 km above the core mantle boundary with 10 per cent p wave velocity reduction and 25 per cent s wave velocity reduction also produces a good fit to the data we find that the low velocity regions act as waveguides enhancing the amplitude of the pkp b caustic diffracted waves and causing superposition of pkp b diffractions with pkpdf the low velocity regions cannot be smooth layers but must show a certain degree of small scale heterogeneity to produce the high frequency scattering observed in the data the path to germany needs small scatterers in the lowermost mantle but the low velocity region is less strong compared with the south pacific ethiopia path a plausible explanation for these low velocity regions is the presence of melt inclusions perhaps entrained in small scale convective systems that act as high frequency wavefield scatterers in regions where these inclusions or convective systems are abundant they appear as a large scale low velocity zone to lower frequency waves these observations support ideas that ultra low velocity zones at the base of the mantle may be heterogeneous regions of melt rich material

p velocity heterogeneities in the lower mantle determined with the german regional seismic network improvement of previous models and results of 2d modelling abstract p waves recorded at the german regional seismic network grsn the grafenberg array the german experiment seismic system geress mini array and two czech stations are used to analyse the inhomogeneous lowermost mantle d under northern siberia and the nansen basin events with a good signal to noise ratio are used to determine the extension of the anomaly relative travel time residuals compared with a 1d model with a relfector in d give information on the shape of the anomaly improved travel time residuals are derived by applying a source or and station correction this correction which was not previously applied reduces the size of the inferred peak to peak topography at the top of d by 20 km if the sources are near lateral inhomogeneities e g near the bend of a subduction zone such corrections are necessary otherwise models of d structure can be severely biased modelling of waveforms for such laterally inhomogeneous structures in d show strong effects on the wave field at dominant periods of 1 to 5 s as multipathing of pdp the p wave reflected at the top of d travel time differences and strong variations in amplitude in agreement with the observations typical dimensions of such hills and valleys are 400 to 1200 km laterally and 10 to 100 km in topograph

variable azimuthal anisotropy in earth s lowermost mantle a persistent reversal in the expected polarity of the initiation of vertically polarized shear waves that graze the d layer the layer at the boundary between the outer core and the lower mantle of earth in some regions starts at the arrival time of horizontally polarized shear waves full waveform modeling of the split shear waves for paths beneath the caribbean requires azimuthal anisotropy at the base of the mantle models with laterally coherent patterns of transverse isotropy with the hexagonal symmetry axis of the mineral phases tilted from the vertical by as much as 20 are consistent with the data small scale convection cells within the mantle above the d layer may cause the observed variations by inducing laterally variable crystallographic or shape preferred orientation in minerals in the d layer

the effect of temperature on the seismic anisotropy of the perovskite and post perovskite polymorphs of mgsio3 using elastic constants determined by ab initio molecular dynamics over a range of temperatures we show that the seismic anisotropy of mgsio3 perovskite is significantly temperature dependent at 90 gpa the direction of greatest shear wave splitting changes from 001 at 0 k to 010 at 3500 k in addition there is no shear wave splitting in the 100 direction at 0 k but a significant amount of splitting is exhibited at 3500 k we have also calculated the first set of high temperature elastic constants of the new post perovskite phase by ab initio molecular dynamics methods we find that in contrast to mgsio3 perovskite temperature has little effect on the acoustic anisotropy of this phase we determine ln vs ln vp for the post perovskite phase and find it to be very low about 0 9 this may be an important parameter for mapping out the post perovskite phase in the lowermost mantle

evidence of a laterally variable lower mantle structure from p and s waves summary examples of grafenberg array data showing anomalous p waves which typically arrive 3 5s after the direct p wave and which have a slowness 0 7 0 8 sdeg smaller than direct p are presented this additional phase is most frequently observed for events located in the ne portion of the southern kurile island subduction zone 73 80 from grafenberg but systematically disappears for events in the sw portion of this zone because of the magnitude of the slowness difference these observations cannot be attributed to a complex source rupture process nor to multipathing through the descending slab likewise they may not be accounted for by near receiver structure because these phases are not seen for all kurile events if present they appear at all stations of the array but they follow direct p too closely to be a multiple from the moho therefore we conclude they are very likely caused by lower mantle velocity structure the most likely explanation is the presence of a p velocity jump of about 3 per cent approximately 290 km above the core mantle boundary since such a reflector in the lowermost mantle not only gives a good fit of traveltimes and slowness but is also able to model the waveform and the amplitudes of this additional p phase the distribution of bounce points on this reflector for the kurile events indicates a lateral extension of this velocity anomaly under northern siberia of about 150 km by at least 200 km the best fitting s wave model has a reflector in the same depth but the velocity contrast seems to be only about 2 per cent suggesting a different behaviour of the p and s velocity in d few events from other regions in this distance range are suitable for a definitive analysis of this kind from among this group some observations indicate a lower mantle anomaly under the lomonosow ridge and under northern greenland but since the lower mantle under western siberia northern novaya zemlya the azores islands region and the ussr afghanistan border region does not produce an additional phase in the grafenberg recordings it is very unlikely that such a velocity anomaly in the lowermost mantle is a global feature

the classification and complex thermal history of the enstatite chondrites we have carried out instrumental neutron activation analyses of 11 enstatite chondrites and electron microprobe analyses of 17 enstatite chondrites most of which were previously little described we report here the third known eh5 chondrite lew 88180 and an unusual el6 chondrite lew 87119 new data on four el3 chondrites alh 85119 eet 90299 pca 91020 and mac 88136 which is paired with mac 88180 and mac 88184 the second el5 chondrite til 91714 and an unusual metal rich and sulfide poor el3 chondrite lew 87223 the often discussed differences in mineral composition displayed by the eh and el chondrites are not as marked after the inclusion of the new samples in the database and the two classes apparently experienced a similar range of equilibration temperatures however texturally the el chondrites appear to have experienced much higher levels of metamorphic alteration than eh chondrites of similar equilibration temperatures most of the petrologic type criteria are not applicable to enstatite chondrites and unlike the ordinary chondrites texture and mineralogy reflect different aspects of the meteorite history we therefore propose that the existing petrologic type scheme not be used for enstatite chondrites we suggest that while textural type reflects peak metamorphic temperatures the mineralogical type reflects equilibration during postmetamorphic probably regolith processes unlike the ordinary chondrites and eh chondrites el chondrites experienced an extensive low temperature metamorphic episode there are now a large number of enstatite meteorite breccias and impact melts and apparently surface processes were important in determining the present nature of the enstatite chondrites

abee and related eh chondrite impact melt breccias abstract the abee eh chondrite is an impact melt breccia about 80 90 of the chondrules were melted many of those that remain have been partly resorbed euhedral enstatite phenocrysts are abundant many have nucleated on the surfaces of relict chondrules and now protrude into kamacite globules euhedral laths of graphite essentially identical to those produced magmatically on earth occur throughout the rock abee experienced two episodes of impact melting with an intervening period of brecciation the first impact melting event produced a relatively homogeneous rock with abundant melt upon cooling the rock was shattered by bombarding meteoroids igneous oldhamite rich dark inclusions were introduced at this time a second impact event produced melt which engulfed the abee fragments and dark inclusions primarily digesting the small ones heat lost to the surrounding clastic debris caused the groundmass melt to quench previously described diamonds in abee containing typical solar system isotopic compositions of n and xe were probably produced from graphite by shock several other eh chondrites may also be impact melt breccias the most extensively melted of these is adhi kot which contains shocked but largely unmelted chondrule rich clasts as well as chondrule free sulfide and silica rich clasts that also appear to be impact products

the development of shape and crystallographic preferred orientation in capto3 post perovskite deformed in pure shear knowledge of the deformation mechanism of post perovskite is important for interpreting observed seismic anisotropy in terms of mantle flow experiments on post perovskite mgsio3 and the low pressure analog material cairo3 yield different textures leaving the interpretation of the observed seismic signatures unclear here we present results of deformation experiments on capto3 post perovskite that may be a better analog to mgsio3 post perovskite capto3 deforms by glide of 100 dislocations on the 010 plane consistent with previous experimental results on cairo3 in addition samples containing a weak minority phase also display shape preferred orientation with grains elongated in the crystallographic a direction forming a planar fabric perpendicular to the compression direction this shape preferred orientation strengthens the observed crystallographic preferred orientation and results in a rapid development of texture during deformation this observation supports the recent suggestion that the d reflector might be due to a rapid generation of texture in post perovskite furthermore the role of shape preferred orientation in generating seismic anisotropy in multi phase assemblages should be considered for the d assemblage

investigating d structure beneath the north atlantic abstract data from the grsn and spiced array are used to investigate lower mantle structure beneath the north atlantic ocean just south of nova scotia we study s waves in an attempt to determine whether a d discontinuity can be observed in this region the reflected phase sds cannot be observed unambiguously so we cross correlate the traces with synthetic seismograms generated using a number of different 1d earth models that incorporate positive and negative velocity discontinuities and gradients as well as thin high velocity lamellae the cross correlation shows that a discontinuity with an s wave velocity increase of magnitude 1 3 exists which exhibits significant topography between 86 and 286 km above the cmb the study area is located on the edge of a region of anomalously high seismic velocity identified by tomographic modelling which fits with our observation of a high velocity discontinuity

metal sulfide silicate intergrowth textures in el3 meteorites origin by impact melting on the el parent body abstract we document the petrographic setting and textures of fe ni metal the mineralogy of metallic assemblages and the modal mineral abundances in the el3 meteorites asuka a 881314 a 882067 allan hills 85119 elephant moraine eet 90299 eet 90992 lapaz icefield 03930 macalpine hills mac 02635 mac 02837 mac 02839 mac 88136 northwest africa nwa 3132 pecora escarpment 91020 queen alexandra range que 93351 que 94321 que 94594 and higher petrologic type els dar al gani 1031 el4 sayh al uhaymir 188 el4 mac 02747 el4 que 94368 el4 and nwa 1222 el5 large metal assemblages often containing schreibersite and graphite only occur outside chondrules and are usually intergrown with silicate minerals euhedral to subhedral enstatite silica and feldspar sulfides troilite daubreelite and keilite are also sometimes intergrown with silicates numerous authors have shown that metal in enstatite chondrites that are interpreted to have been impact melted contains euhedral crystals of enstatite we argue that the metal sulfide silicate intergrowths in the els we studied were also formed during impact melting and that metal in el3s thus does not retain primitive i e nebular textures likewise the el4s are also impact melt breccias modal abundances of metal in the el3s and el4s range from approximately 7 to 30 wt these abundances overlap or exceed those of el6s and this is consistent either with pre existing heterogeneity in the parent body or with redistribution of metal during impact processes

seismic evidence links the subducted mongol okhotsk slab to deformation in d near the northeastern margin of the perm anomaly abstract numerical simulations have indicated that the large low shear velocity provinces llsvps dynamically interact with subducted slabs in d and slabs in turn play an important role in influencing the shape and location of the llsvps in this study we find evidence for the existence of the d discontinuity near the northeastern margin of the perm anomaly a smaller llsvp beneath eurasia by analyzing the scd the lower mantle triplication phases generated from two deep events occurring in the northwestern pacific subduction zone and recorded at stations in europe and turkey the d discontinuity in the north is characterized with a 2 velocity increase occurring at a height of 300 km above the core mantle boundary cmb whereas the discontinuity in the south is characterized with a 3 velocity increase occurring at a height of 260 km above the cmb the d discontinuity is most likely attributed to a phase transition from perovskite to post perovskite in conjunction with the presence of relic slab materials likely related to the subducted mongol okhotsk slab our observations along with previous seismic anisotropic detections provide a clue to dynamically link the subducted slab to deformation taking place near the edges of the llsvp

sr z sbnd nd isotopic record of multi stage interactions between mantle derived magmas and crustal components in a collision context the ultramafic granitoid association from vivero hercynian belt nw spain late stage epizonal calc alkaline granitoids of the hercynian belt are often associated with coeval more basic rocks such granitoids have been the focus of special attention when studying the role of mantle derived magmas in granitoid genesis these studies conclude that basic magmas are involved in granitoid production but their contribution to the chemical composition of granitoids is often controversial the vivero massif nw spain includes calc alkaline granitoids and coeval ultramafic and mafic rocks but these granitoids are syntectonic and intruded into middle crustal levels during an early stage scattered major and trace element chemistry together with heterogeneous sr nd isotopic results preclude an origin of this ultramafic to granitoid association by any closed system fractionation process ultramafic rocks include peridotites pyroxenites and homblendites with endi from 3 5 to 1 3 and 87sr86sr i ratio from 0 7059 to 0 7080 they are interpreted as cumulates formed by in situ fractional crystallization from contaminated basaltic magmas particularly by interaction with 87sr enriched fluids mafic rocks comprise biotite amphibole diorites quartz diorites and tonalites with ndi from 3 8 to 0 0 and 87sr86sr i ratio from 0 7058 to 0 7073 they represent either cumulates or more evolved magmas of calc alkaline affinity these rocks display an inverse correlation between the degree of differentiation and the degree of contamination inferred from isotope data this points to wall rock assimilation by thermal erosion during ascent of basic magmas up to the middle crust and supports the participation of mantle sources with positive end and 87sr86sr i 0 7058 granitoids range from amphibole biotite tonalites through biotite granodiorites to monzonitic granites they show an isotopic trend which diverges from that of ultramafic and mafic rocks endi 6 0 to 1 9 decreases while 87sr86sr i ratio 0 7061 to 0 7082 increases from tonalites to monzonitic granites reflecting an increasing contribution of crustal components in more differentiated rocks assimilation and fractional crystallization mixing included account for such an evolution contamination of basic magmas by this process occurring at deeper crustal level produced hybrid tonalites and related mafic magmatic enclaves only monzonitic granites are interpreted as pure crustal melts their hree fractionated patterns gd yb 4 5 suggest equilibrium with a garnet rich residuum partial melting of moderate time integrated rb sr protolith such as greywackes and or meta igneous rocks of intermediate acid composition could account for the isotopic composition of these monozogranites this example shows that at least some calc alkaline granitoids of collisional setting are not produced by pure intra crustal melting they are formed through a complex multi stage hybridization process involving mantle derived magmas crustal components and several concomitant magmatic processes

anisotropy in the deep earth abstract seismic anisotropy has been found in many regions of the earth s interior its presence in the earth s crust has been known since the 19th century and is due in part to the alignment of anisotropic crystals in rocks and in part to patterns in the distribution of fractures and pores in the upper mantle seismic anisotropy was discovered 50 years ago and can be attributed for the most part to the alignment of intrinsically anisotropic olivine crystals during large scale deformation associated with convection there is some indication for anisotropy in the transition zone particularly in the vicinity of subducted slabs here we focus on the deep earth the lower mantle and core where anisotropy is not yet mapped in detail nor is there consensus on its origin most of the lower mantle appears largely isotropic except in the last 200 300 km in the d region where evidence for seismic anisotropy has been accumulating since the late 1980s mostly from shear wave splitting measurements recently a picture has been emerging where strong anisotropy is associated with high shear velocities at the edges of the large low shear velocity provinces llsvps in the central pacific and under africa these observations are consistent with being due to the presence of highly anisotropic mgsio 3 post perovskite crystals aligned during the deformation of slabs impinging on the core mantle boundary and upwelling flow within the llsvps we also discuss mineral physics aspects such as ultrahigh pressure deformation experiments first principles calculations to obtain information about elastic properties and derivation of dislocation activity based on bonding characteristics polycrystal plasticity simulations can predict anisotropy but models are still highly idealized and neglect the complex microstructure of polyphase aggregates with strong and weak components a promising direction for future progress in understanding the origin of seismic anisotropy in the deep mantle and its relation to global mantle circulation is to link macroscopic information from seismology and microscopic information mineral physics through geodynamics modeling anisotropy in the inner core was proposed 30 years ago to explain faster p wave propagation along the direction of the earth s axis of rotation as well as anomalous splitting of core sensitive free oscillations there is still uncertainty about the origin of this anisotropy in particular it is difficult to explain its strength based on known elastic properties of iron as it would require almost perfect alignment of iron crystals indeed the strongly anomalous p travel times observed on paths from the south sandwich islands to alaska may or may not be due to inner core anisotropy and will need to be explained before consensus can be reached on the strength of anisotropy in the inner core and its origin

ar ar ages and thermal histories of enstatite meteorites abstract compared with ordinary chondrites there is a relative paucity of chronological and other data to define the early thermal histories of enstatite parent bodies in this study we report 39ar 40ar dating results for five el chondrites khairpur pillistfer hvittis blithfield and forrest five eh chondrites parsa saint marks indarch bethune and reckling peak 80259 three igneous textured enstatite meteorites that represent impact melts on enstatite chondrite parent bodies zaklodzie queen alexandra range 97348 and queen alexandra range 97289 and three aubrites norton county bishopville and cumberland falls several ar ar age spectra show unusual 39ar recoil effects possibly the result of some of the k residing in unusual sulfide minerals such as djerfisherite and rodderite and other age spectra show 40ar diffusion loss few additional ar ar ages for enstatite meteorites are available in the literature when all available ar ar data on enstatite meteorites are considered preferred ages of nine chondrites and one aubrite show a range of 4 50 4 54 ga whereas five other meteorites show only lower age limits over 4 35 4 46 ga ar ar ages of several enstatite chondrites are as old or older as the oldest ar ar ages of ordinary chondrites which suggests that enstatite chondrites may have derived from somewhat smaller parent bodies or were metamorphosed to lower temperatures compared to other chondrite types many enstatite meteorites are brecciated and or shocked and some of the younger ar ar ages may record these impact events although impact heating of ordinary chondrites within the last 1 ga is relatively common for ordinary chondrites only bethune gives any significant evidence for such a young event

thermal history of the enstatite chondrites from silica polymorphs here we report the results of our petrologic and mineralogical study of enstatite e chondrites in order to explore their thermal histories we studied silica phases in 20 e chondrites by laser micro raman spectroscopy to determine the silica polymorphs they contain silica phases are commonly present in e chondrites and their polymorphs reflect the physical conditions of formation the samples studied here include eh3 5 el3 6 e chondrite melt rocks and an anomalous e chondrite we identified quartz tridymite cristobalite and silica glass in the samples studied eh4 5 and eh melt rocks are divided into high and low temperature classes based on niningeritealabandite solid solutions eh3 el3 and some eh melt rocks of the high temperature class contain tridymite and cristobalite we suggest that tridymite and cristobalite crystallized in chondrules and e chondrite melts followed by rapid cooling leading to the survival of these silica polymorphs eh4 and el4 chondrites also contain tridymite and cristobalite in their chondrules indicating that these silica polymorphs survived low temperature metamorphism as estimated from opaque mineral geothermometers because of the sluggishness of the transition to a more stable polymorph tridymite and cristobalite in el6 chondrites reflect the high temperature processes experienced by these meteorites on the other hand some eh5 chondrites and eh melt rocks of the low temperature class contain quartz which may be a product of the transition from tridymite or cristobalite during a long period of low temperature metamorphism although the thermal history of e chondrites have been previously estimated from opaque minerals such compositions mainly reflect low temperature processes however we can reconstruct the primordial thermal processes and subsequent cooling histories of e chondrites from their silica polymorphs the e chondrites have complicated thermal histories which produced the observed variations among them

complex seismic anisotropy at the border of a very low velocity province at the base of the earth s mantle 1 we constrain the anisotropy associated with a very low velocity province vlvp at the base of the earth s mantle using the sks and skks waves sampling the region our selected high quality data sets consist of 415 sks and 111 skks waveforms for 127 deep earthquakes recorded at distances between 90 and 150 by the seismic stations in three temporary broadband passcal seismic arrays the kaapvaal seismic array 1997 1999 the tanzania seismic array 1994 1995 and the ethiopia kenya seismic array 1999 2001 as well as the permanent stations in the global seismographic network these seismic data provide good sampling coverage for some portion of the vlvp and its surrounding areas our results show when the sks or skks phases sample the regions away from the border of the vlvp inside or outside the vlvp the apparent splitting parameters inferred from the sks phases are consistent with those inferred from the skks phases and their variations strongly correlate with seismic stations but not with the exit points at the core mantle boundary of these seismic phases however when the sks or skks phases sample near the border of the vlvp the apparent splitting parameters inferred from the sks phases and skks phases are different and their variations no longer correlate with seismic stations these features indicate that part of the shear wave splitting for the seismic data sampling the border of the vlvp has to originate from deep mantle most likely near the border of the vlvp we assume that the anisotropy in the shallow mantle beneath seismic stations has a horizontal hexagonal symmetry axis and infer the splitting parameters associated with the shallow anisotropy beneath the seismic stations using the sks and skks waveforms for the seismic data sampling the regions away from the border of the vlvp we then obtain the splitting parameters associated with the lowermost mantle anisotropy using the sks and skks waveforms corrected for the inferred shallow anisotropy beneath seismic stations assuming that the medium in the lowermost mantle also has a horizontal hexagonal symmetry axis our results reveal a complex anisotropy pattern near the border of the vlvp such a complex anisotropy pattern may be explained by lattice preferred orientation of anisotropic mantle aggregates aligned by a complex mantle flow near the vlvp margins the complex flow pattern near the vlvp margins may indicate strong interactions between the vlvp a compositional anomaly and the surrounding normal mantle and may provide an explanation to the concentration of some hot spots geographically near the borders of the vlvps in the lowermost mantle

lateral variations in d thickness from long period shear wave data we explore global variations in d shear wave structure by examining 12 years of long period global digital seismograph network gdsn data to identify reflected phases from the d layer in the lowermost mantle we restrict our search to epicentral distances between 63 and 74 where a precritical d reflection will lie between the s and scs arrivals with gdsn long period records a d signal at these ranges is generally obscured by the stronger neighboring phases so to isolate the d reflections we develop a technique for stripping away the interfering s and scs waveforms we interpret travel time variations in the observed sds phases in terms of variations in d thickness and generate maps of inferred d thickness for areas of good seismic coverage the regions of observed d reflections are beneath australasia north central asia the arctic alaska and central america we find considerable variation in d structure in these regions a result consistent with previous studies which have found evidence for lateral heterogeneity in d on average the estimated d thickness is 260 km but there is a fairly uniform distribution of thicknesses between 150 km and 350 km a possible correlation is observed between d thickness and regions of predicted strong horizontal mantle flow

laterally varying reflector at the top of d beneath northern siberia summary we analyse p wave data from the french laboratoire de dctection et gcophysique l d g network of digital short period seismic stations whereas our previous paper houard nataf 1992 showed a few promising kurd events and depicted the methods we proposed to apply an extensive analysis of 32 earthquakes from the kurivkamchatka region has now been carried out by using the deconvolution technique previously described houard nataf 1992 clear and quantitative evidence of intermediate arrivals between the p and pcp waves is shown in the 75 82 distance range our approach is complementary to previous studies in the same region which used broad band data recorded at a dense narrow aperture array weber davis 1990 weber 1993 or long period data recorded at the wwssn stations gaherty lay 1992 for each station record sections of deconvolved signals from different earthquakes are built and lead us to eliminate crustal reverberations beneath l d g stations as a major explanation for these secondary arrivals a compilation of the most conclusive signals has been made it is perhaps one of the most complete pdp wave record sections since the move out with distance can be followed clearly with a good sampling over more than 5 the move out with distance of this extra pdp arrival its residual slowness value eliminates source complexities diffraction or multipathing through the kuril subduction slab as possible explanations the extra arrivals indicate the existence of a lower mantle reflector the structure of the lay discontinuity is investigated using tpdp tp residual traveltimes two different residual time versus distance tendencies are found which correlate to different geographical bounce point regions indicating the existence of lateral variations of the structure the d region is also sampled further south at distances well beyond the theoretical triplication cross over by analysing events from the honshu ryukyu islands region but only one p wave branch is observed the hodochrones and p a curves of the first p arrival are used as a complementary analysis the global prem model dziewonski anderson 1981 predicts correct p wave slownesses though its smooth lower mantle structure does not account for pdp observations this may be an indication that the lay discontinuity is a global feature however the hodochrones of the first p arrival do not allow us to test the existence of the discontinuity in the 82 85 cross over distance range

shock metamorphism of ordinary chondrites a revised petrographic classification of progressive stages of shock metamorphism of ordinary chondrites is proposed six stages of shock s1 to s6 are defined based on shock effects in olivine and plagioclase as recognized by thin section microscopy the characteristic shock effects of each shock stage are s1 unshocked sharp optical extinction of olivine s2 very weakly shocked undulatory extinction of olivine s3 weakly shocked planar fractures in olivine s4 moderately shocked mosaicism in olivine s5 strongly shocked isotropization of plagioclase maskelynite and planar deformation features in olivine and s6 very strongly shocked recrystallization of olivine sometimes combined with phase transformations ringwoodite and or phases produced by dissociation reactions s6 effects are always restricted to regions adjacent to melted portions of a sample which is otherwise only strongly shocked in stages s3 to s6 localized melting results from stress and temperature peaks which locally deviate from the equilibration shock pressure due to differences in shock impedance these melting effects are a opaque melt veins shock veins b melt pockets with interconnecting melt veins c melt dikes and d troilite metal deposits in fractures based on a critical evaluation of data from shock recovery experiments a shock pressure calibration for the six shock stages is proposed which defines the s1 s2 s2 s3 s3 s4 s4 s5 and s5 s6 transitions at 5 5 10 15 20 30 35 and 45 55 gpa respectively whole rock melting and formation of impact melt rocks or melt breccias occurs at about 75 90 gpa the symbol for the shock stage may be used in combination with the symbol for the petrologic type to abbreviate the complete classification of a chondrite e g h5 s3 we propose this new shock classification and pressure calibration system to replace previous systems which are out of date with respect to the pressure calibration nominally restricted to l chondrites and based on incomplete and in part illdefined sets of shock effects we have classified seventy six ordinary chondrites using the new classification system and conclude the following 1 1 shock effects and the sequence of progressively increasing degrees of shock metamorphism are very similar in h l and ll groups differences in the frequency distribution of shock stages are relatively minor e g l chondrites appear to have the largest fraction with stages s5 and s6 this suggests that the collisional histories of the h l and ll parent bodies were similar 2 2 petrologic type 3 chondrites are deficient in stages s4 to s6 and with increasing petrologic type the frequency of stages s4 to s6 increases we suggest that the more porous and volatile rich type 3 chondrites are subject to melting at a lower shock pressure than the nonporous chondrites of higher petrologic type volatiles trapped in pores cause shock induced dispersal of the shocked and melted material into small particles which are not expected to survive as meteorites 3 3 stage s3 is the most abundant in nearly all petrologic types 4 4 at shock pressures in excess of about 35 gpa s5 and s6 4he and 40ar are almost completely lost pressures below 10 gpa s1 and s2 do not cause noble gas losses

investigation of laterally heterogeneous shear velocity structure in d beneath eurasia seismic waves are used to explore large scale radial structure and lateral heterogeneity of shear velocity in the d region the chemical and thermal boundary layer at the base of the mantle we analyze long period horizontally polarized sh shear waveforms from earthquakes in the northwest pacific recorded in europe and the middle east for evidence of proposed laterally varying shear velocity stratification in d beneath eurasia specifically our data sample the lower mantle extending from beneath the arctic ocean near novaya zemlya ussr to under northern india using a reflectivity technique to model the sh waveforms we determine a locally stratified shear velocity model sgle that provides a good average travel time and waveform fit to over 80 of the data sensitive to d structure sgle is characterized by a rapid 2 75 increase in shear velocity approximately 290 km above the core mantle boundary similar to previously proposed models for this region unlike for other regions analyzed using similar waveform modeling techniques the percentage of data inconsistent with this average radial model is fairly substantial indicating that laterally heterogeneous d velocity structure is important in this region some waveforms lack the reflected arrival expected for a stratified model while others show evidence of two arrivals both possibly originating within d while our data are insufficient for direct inversion of the three dimensional structure we utilize several simple parameterizations of the waveforms to gain some understanding of the region both the relative amplitude and differential travel time variations of the waveforms indicate that d beneath eurasia is characterized by 1 broad regions in which the structure has large scale stratification giving coherent long period reflections with superimposed velocity variations that scatter shorter wavelength waves and 2 smaller regions which are sufficiently heterogeneous to scatter even the long period energy the scale length of the lateral heterogeneity in d that produces these variations is inferred to be somewhat less than 500 km which contrasts with weaker lateral heterogeneity that is observed in the central region of the lower mantle where heterogeneity scale lengths of approximately 500 1000 km appear to dominate these observations are consistent with a model of the lower mantle in which distributed heterogeneity is concentrated in d generating a locally stratified layer in some regions the effects of a diffuse or topographically variable d discontinuity must be further explored using two and three dimensional synthetic calculations before inversions for the detailed structure can be undertaken

lowermost mantle anisotropy beneath the north pacific from differential s scs splitting summary seismic anisotropy is an important tool for studying the nature origin and dynamics of the lowermost mantle d we introduce differential s scs splitting as a tool for removing the effect of near source and near receiver anisotropy to estimate splitting accrued in the d region this is applicable to events recorded at epicentral distances between 60 and 85 near source anisotropy has often been ignored in previous studies of lowermost mantle anisotropy we apply differential s scs splitting to records from canadian national seismic network stations of western pacific earthquakes these sample the lowermost mantle beneath the north pacific the residual splitting in scs which we attribute to d shows lag times between 1 0 and 3 9 s given the near horizontal ray path of scs in d we interpret the recovered fast directions as the orientation of the fast shear wave in the plane defined by the vertical and transverse directions and observe a clearly non vti transverse isotropy with a vertical axis of symmetry style of anisotropy the largest population of results shows an approximately southeasterly dipping symmetry axis which we speculate might be explained by descending palaeoslab material being swept horizontally across the core mantle boundary towards an upwelling region beneath the central pacific non vti symmetry and the many possible contributions to d anisotropy from lower mantle minerals melt and subducted materials suggest that our understanding of the lowermost mantle could be greatly improved by trying to resolve a more general style of anisotropy

origin and history of impact melt rocks of enstatite chondrite parentage we have conducted petrologic chemical and isotopic studies of two impact produced rocks of enstatite chondrite parentage ilafegh 009 is a total impact melt rock with no residual lithic clasts formation on the el chondrite parent body is suggested by its mineralogy and mineral compositions cooling of the impact melt was rapid at melt temperatures and decreased at subsolidus temperatures in contrast to previous studies we show that happy canyon is not a new enstatite achondrite but an impactmelt breccia of enstatite chondrite and not aubrite parentage this rock formed by impact melting and incorporation into the melt of clastic material which resulted in relatively rapid cooling at all temperatures mineralogical and bulk compositional data probably biased by the heterogeneous nature of this rock do not allow unequivocal determination of its parent body i e el vs eh although some data such as bulk total fe content seem to favor el parentage both rocks were subjected to post solidification shock which was more severe for ilafegh 009 than for happy canyon it appears that both impact melt rocks could have formed by impact melting 4 57 ga ago as is indicated by the nearly identical ixe closure ages of 1 6 and 1 4 ma before bjurbole for ilafegh 009 and happy canyon respectively an apparently younger 39ar40ar age of 4 53 ga for happy canyon may be due to small biases in the intercalibration of the ixe and 39ar40ar chronometers whereas the much younger 4 34 4 44 ga age for ilafegh 009 reflects thermal resetting during shock metamorphism shallowater which was impact derived from a different enstatite achondrite parent body has an ixe closure age 0 4 ma younger than that for ilafegh 009 and an 39ar40ar age of 4 53 ga the ancient ages of these three rocks attest to the intense early bombardment in this region of the solar system

possible impact induced refractory lithophile fractionations in el chondrites abstract literature data show that refractory lithophile elements in most chondrite groups are unfractionated relative to ci chondrites the principal exception is the el chondrite group whose observed falls all of which are type 6 are depleted in ca and light ree in contrast literature data and our new inaa data on el3 pca 91020 el3 mac 88136 and el4 grein 002 show that some replicates of these samples have nearly flat ree patterns unlike those of el6 chondrites other replicates exhibit fractionated ree patterns similar to those of el6 chondrites petrographic examination shows that many el6 and some el3 and el4 chondrites are impact melt breccias or contain impact melted portions we suggest that the same impact processes that formed these breccias and produced melt are responsible for the observed bulk compositional fractionations in refractory lithophile elements i e el6 chondrites were produced from initially unequilibrated el3 material when large amounts of impact heat were deposited plagioclase and or oldhamite cas the major ree carriers in enstatite chondrites may have been melted and then transported appreciable 10 cm distances el6 chondrites represent the residuum that is depleted in ree particularly in lree and ca unlike the case for el chondrites our new inaa data on alh 84170 eet 87746 and sah 97096 all eh3 show some scatter but are consistent with the eh group having uniform refractory lithophile abundances

a lower mantle s wave triplication and the shear velocity structure of d summary a lower mantle s wave triplication detected with short and longperiod wwssn and csn recordings indicates a substantial shear velocity discontinuity near 280 km above the core mantle boundary the triplication can be observed in rotated sh seismograms from intermediate and deep focus events throughout the distance range from 70 to 95 three distinct source region receiver array combinations that have been investigated in detail demonstrate consistent travel time and relative amplitude behaviour of the triplication with slight systematic shifts in the triplication indicating up to 40km variations in the depth of the discontinuity modelling of the observations with synthetic seismograms produced with the cagniard de hoop and reflectivity methods constrains the shear velocity increase to be 2 75 k 0 25 per cent comparable to upper mantle discontinuities short period observations indicate that the velocity increase may be a sharp first order discontinuity or may extend over a transition zone no more than 50 km thick the shear velocity gradient below the discontinuity within the d layer is not well constrained by the sh data but slightly positive or near zero velocity gradients are consistent with the long period amplitude ratios of scsh sh

partial melting of the indarch eh4 meteorite a textural chemical and phase relations view of melting and melt migration to test whether aubrites can be formed by melting of enstatite chondrites and to understand igneous processes at very low o fugacities we have conducted partial melting experiments on the indarch eh4 chondrite at 1000 1500 c silicate melting begins at 1000 c and indarch is completely melted by 1500 c the metal sulfide component melts completely at 1000 c substantial melt migration occurs at 1300 1400 c and metal migrates out of the silicate charge at 1450 c and 50 silicate partial melting as a group our experiments contain three immiscible metallic melts si p and c rich two immiscible sulfide melts fe and femgmnca rich and silicate melt our partial melting experiments on the indarch eh4 enstatite chondrite suggest that igneous processes at low fo2 exhibit several unique features the complete melting of sulfides at 1000 c suggests that aubritic sulfides are not relics aubritic oldhamite may have crystallized from ca and s complexed in the silicate melt significant metal sulfide melt migration might occur at relatively low degrees of silicate partial melting substantial elemental exchange occurred between different melts e g s between sulfide and silicate si between silicate and metal a feature not observed during experiments at higher fo2 this exchange may help explain the formation of aubrites from known enstatite chondrites

shear wave splitting and waveform complexity for lowermost mantle structures with low velocity lamellae and transverse isotropy 1 shear waves that traverse the lowermost mantle exhibit polarization anomalies and waveform complexities that indicate the presence of complex velocity structure above the core mantle boundary synthetic seismograms for horizontally and vertically polarized shear waves sh and sv respectively are computed using the reflectivity method for structures with low velocity sheets lamellae and for comb like models approximating long wavelength vertical transverse isotropy vti motivated by evidence for partial melt in the deep mantle lamella parameter ranges include 1 vp from 5 to 10 vs 3 vp 2 100 to 300 km thickness of vertical stacks of lamella 3 lamella spacing and thickness varying from 0 5 to 20 km and 4 lamellae concentrated near the top bottom or throughout the d region at the base of the mantle such lamellae represent in effect horizontally emplaced dikes within d excessively complex waveforms are produced when more than 20 of d volume is comprised of low velocity lamellae many lamellae models can match observed sdiff splitting 1 10 s delays of svdiff but typically underpredict scs splitting 1 4 s delays of scsv vti model parameters are selected to address d observations and include 1 0 5 to 3 anisotropy 2 discontinuous d shear velocity increases up to 3 3 d thicknesses from 100 to 300 km and 4 vti concentrated at the top bottom or throughout d vti models readily match observed splits of scs and sdiff we discuss lamellae and vti model attributes in relationship to waveform complexities splitting magnitude triplications from a high velocity d discontinuity and apparently reversed polarity svdiff onsets the possible presence of melt filled lamellae indicates that local chemical or thermal perturbations can produce regions that exceed the solidus within d such melt could occur in the bulk of d because the melt is either close to neutral buoyancy advective velocities exceed percolative velocities or both

shock metamorphism of enstatite chondrites abstract we have extended the shock metamorphism classification scheme of stoffler et al 1991 to account for shock effects in orthopyroxene and applied it to sixty enstatite chondrites orthopyroxene exhibits the following sequence of progressive shock effects shock stage s1 unshocked sharp optical extinction s2 very weakly shocked undulose extinction s3 weakly shocked development of clinoenstatite lamellae parallel to 100 s4 moderately shocked weak mosaicism s5 strongly shocked strong mosaicism as in the stoffler et al 1991 scheme stage s5 is characterized principally by the solid state transformation of crystalline plagioclase into maskelynite most el3 chondrites exhibit foliations caused by impact deformation of chondrules and metal particles el5 and el6 chondrites are all shock stage s2 because opaque veins are rare in s2 ordinary chondrites it is plausible that the centimeter size kamacite veins in atlanta and blithfield and the 1 6 cm long oldhamite rich vein in jajh deh kot lalu all el6 formed when their hosts were shocked to s3 s5 levels because removal of shock stage s3 s5 features including repair of shock damaged orthopyroxene requires levels of metamorphism comparable to those experienced by petrologic type 5 to 6 chondrites we infer that proto el6 material was shocked to s3 s5 levels prior to peak metamorphism and shocked again to stage s2 after metamorphism overall enstatite chondrites appear to have suffered greater shock damage than ordinary and carbonaceous chondrites

petrographic and mineralogical study of new eh melt rocks and a new enstatite chondrite grouplet enstatite chondrites are classified into eh and el groups and some of them were melted once in this paper we report petrography and mineral chemistry of five new antarctic enstatite chondrites yamato 793225 is characterized by intermediate properties between the eh and el groups and probably represents a new grouplet of enstatite chondrites the three paired meteorites y 8404 y 8414 and y 86004 may have cooled rapidly near the surface of the parent body as impact eh melt rocks yamato 82189 is also classified as eh melt rock but it experienced slow cooling as did y 793225 yamato 82189 contains the first occurrence of phlogopite in enstatite meteorites

detecting seismic anisotropy above the 410 km discontinuity using reflection coefficients of underside reflections abstract we investigate the effect of various types of deformation mechanisms on the reflection coefficients of p and s waves underside reflections off the 410 km discontinuity to develop a diagnostic tool to detect the style of deformation at boundary layers we calculate the reflection coefficient for p and sh underside reflections using velocity perturbations resulting from aligned minerals above the 410 km discontinuity for different deformation scenarios the results show that in the case of an anisotropic olivine layer above an isotropic wadsleyite layer the p wave reflection coefficient amplitudes are only slightly influenced by the joint effect of angle of incidence and the strength of imposed deformation without any polarity reversal and for all deformation styles for the sh wave underside reflections the incidence angle for which a polarity reversal occurs changes with distance for all scenarios and in addition changes with azimuth for shear deformation scenarios these differences in amplitude and polarity patterns of reflection coefficients of different deformation geometries especially for s wave at shorter distances potentially provide a possibility to distinguish between different styles of deformation at a boundary layer we also show a first test using currently available elastic constants of anisotropic wadsleyite beneath anisotropic olivine

observational evidence for diffracted sv in the shadow of the earth s core on the basis of observational evidence it has been assumed for a long time in the seismological literature that the diffracted s wave in the deep shadow of the earth s core is polarized as sh in good agreement with models that do not contain a well pronounced low velocity zone at the base of the mantle while this seems often to be the case broad band horizontal records from the geoscope network show that for some wavepaths the diffracted sv is strong enough to be clearly observed at distances between 107 and 117 moreover in some cases the diffracted sv is shifted with respect to sh by approximately a quarter period such observations should put additional constraints on the structure near the core mantle boundary among other possibilities lateral variations of the gradient at the base of the mantle allowing for regions with negative gradient and possibly azimuthal anisotropy in d should be considered in order to explain these data

azimuthal anisotropy in the d layer beneath the caribbean 1 the lowermost mantle beneath central america has anisotropic seismic velocity structure manifested in shear wave splitting of signals from south american earthquakes recorded at north american broadband recording stations prior studies of deep mantle anisotropy in this region have characterized the structure as having vertical transverse isotropy vti which is sufficient to explain a general trend of early tangential sh component arrivals however vti models cannot quantitatively match systematic waveform complexities in the onset of many of the shear waves that graze this region after accounting for splitting effects of upper mantle anisotropy beneath the recording stations we model the corrected waveform data using full wave theory for mantle velocity models with an anisotropic d layer this is the first attempt to quantitatively model a large data set including azimuthal anisotropy in d the models include transverse isotropy with either a vertical or tilted symmetry axis the latter resulting in azimuthal anisotropy for some initial shear wave polarizations tilted transverse isotropy tti produces small reversed polarity arrivals on the sv components at the arrival time of sh consistent with the data geographical variations in the azimuth of the tti symmetry axis are indicated by the data the lack of azimuthal coverage prevents unique resolution of the tti orientation and also precludes distinguishing between tti and other azimuthal anisotropy structures such as that predicted for lattice preferred orientation of minerals nonetheless our modeling demonstrates the need for laterally varying anisotropic structure of more complex form than vti for this region

evidence of midmantle anisotropy from shear wave splitting and the influence of shear coupled p waves 1 until recently the midmantle region of the earth the transition zone and the uppermost lower mantle has been considered isotropic recent work by wookey et al 2002 presented evidence of anisotropy in the midmantle region near the tonga kermadec subduction zone inferred from shear wave splitting in teleseismic s phases at australian seismic stations this data set is revisited to examine the possibility of contamination by shear coupled p waves using reflectivity modeling we explore the effects of these phases on the measurement of shear wave splitting contamination of such measurements is demonstrated for short epicentral distances 30 wave field decomposition is a simple technique which can be used to separate these phases from the main s wave arrival the tonga australia data set is reprocessed after wave field decomposition shear wave splitting ranging between 0 7 and 6 2 s is observed a range comparable to that observed by wookey et al 2002 as before the polarization of the fast shear wave is observed to be predominantly horizontal however the magnitudes of splitting for several events at short epicentral distances are significantly reduced suggesting some influence of shear coupled p waves in the original analysis the results are compared with possible models of mantle anisotropy and the results can be best explained by significant anisotropy in the midmantle region the best candidate location for this anisotropy is in the uppermost lower mantle and scenarios involving the alignment of lower mantle or subducted materials by regional dynamic processes are suggested to account for this

seismic anisotropy in the core mantle transition zone split s waves observed at hockley texas from events in the tonga fiji region of the southwest pacific show predominantly vertically polarized shear wave sv energy arriving earlier than horizontally polarized sh energy for rays propagating horizontally through d after corrections are made for the effects of upper mantle anisotropy beneath hockley a time lag of 1 5 to 2 0 s remains for the furthest events 93 9 100 6 while the time lags of the nearer observations 90 5 92 9 nearly disappear at closer distances the s waves from these same events do not penetrate as deeply into the lower mantle and are not split these observations suggest that a patch of d beneath the central pacific is anisotropic while the mantle immediately above the patch is isotropic the thickness of the anisotropic zone appears to be of the order of 100 200 km observations of shear wave splitting have previously been made for paths that traverse d under the caribbean and under alaska sh leads sv the reverse of the hockley observations but in these areas the fact that sv leads sh in the hkt data shown here suggests a different sort of anisotropy under the central pacific from that under alaska and the caribbean the case of sh travelling faster than sv is consistent with transverse isotropy with a vertical axis of symmetry vti and does not require variations with azimuth the case of sv leading sh is consistent with transverse isotropy with a horizontal axis of symmetry hti an azimuthally anisotropic medium and with a vti medium formed by a hexagonal crystal given that mg fe sio3 perovskite appears unlikely to form anisotropic fabrics on a large scale the presence of anisotropy may point to chemical heterogeneity in the lowermost mantle possibly due to mantle core interactions

upper and lower mantle anisotropy inferred from comprehensive sks and skks splitting measurements from india abstract in this study we investigate the upper mantle anisotropy beneath india using high quality sks and skks waveforms from 382 teleseismic earthquakes recorded at 119 broadband seismic stations in addition we present evidence for anisotropy in the d layer beneath southeast asia using sks and skks splitting discrepancies on the same seismogram during this exercise we obtain 200 new splitting measurements from 35 stations recently deployed in the indo gangetic plains igp central india and northeast india while the delay times between the fast and slow axes of anisotropy t range from 0 3 to 1 7 s the fast polarization azimuths at a majority of stations in the igp and central india coincide with the absolute plate motion of india implying shear at the base of the lithosphere as the dominant mechanism for forging anisotropy however stations in ne india reveal fast polarization azimuths mainly in the ene wsw direction suggestive of lithospheric strain induced by the ongoing indo eurasian collision our analysis for d anisotropy yielded a total of 100 sks skks pairs which can be categorized into those exhibiting i null measurements for one phase and significant splitting for the other phase ii null measurement for both the phases iii significant splitting for both the phases a pair is considered to be anomalous if the splitting difference between sks and skks is 0 5 s and the individual split time is 0 5 s using this criterion we obtain 12 measurements under category iii and 9 under category i that show a null measurement for sks and large splitting for the skks phase further we quantify the strength of the lower mantle anisotropy by correcting the skks measurement for the upper mantle anisotropy obtained by the sks phase on the same seismogram the sks delay times are found to be consistently less than skks times suggesting that the sks phases do not capture the lower mantle anisotropy in comparison to their skks counterparts seven coherent measurements thus obtained reveal measurable d anisotropy with fast polarization azimuths oriented mainly in the ene wsw direction these results suggest presence of a large region of deformation in the lowermost mantle beneath southeast asia a possible model for anisotropy in these regions could be the presence of slab material that pounded upon the core mantle boundary cmb and is experiencing large shear deformation resulting in lattice preferred orientation lpo of the lower mantle van der hilst and karason 1999 long 2009 the other possibility is the phase transformation from mgsio3 perovskite to a more stable post perovskite phase under favorable conditions which results in lpo of the lower mantle

small scale lateral shear velocity and anisotropy heterogeneity near the core mantle boundary beneath the central pacific imaged using broadband scs waves core reflected scs waves from 49 large m 5 1 deep tonga fiji events recorded in western north america are used to study a localized region of the core mantle boundary cmb under the central pacific ocean a total of 248 observations from the berkeley digital seismic network bdsn caltech united states geological survey usgs terrascope and the incorporated research institutions for seismology iris broadband arrays span epicentral distances of 73 85 scs reflection points sample a cmb patch southeast of the hawaiian islands at latitude 4 to 16 n and longitude 156 to 144 w near the proposed source location of the hawaiian plume highly variable scs travel times amplitudes waveforms and shearwave splitting indicate that the lowermost mantle in this region is heterogeneous both laterally and radially scsh sh differential travel times are on average 4 s larger than predicted by the preliminary reference earth model prem this is due to delayed scs arrivals and is largely accounted for by a model with a strong velocity decrease in d the scsh sh residuals also show a spatial trend not accounted for by a radial model indicating a lateral decrease in lower mantle shear velocity to the northeast this lateral velocity gradient appears to cause focusing of scs energy scsh sh amplitude ratios are larger than predicted by prem or by a model with a strong negative gradient with depth in d which enhances scs scs splitting corrected for lithospheric anisotropy beneath the receivers indicates spatial variations in d anisotropy with a systematic change in the orientation of the fast polarization direction from transverse to the ray path fast scsh to parallel to the ray path fast scsv along a northeast traverse these spatial trends suggest lateral gradients in the boundary layer shear flow on scale lengths of a few hundred kilometers which may be related to dynamical flow near the near the root of the large hawaiian plume

seismic anisotropy in the d layer we present observations of diffracted sv svd for a path between the fiji tonga islands and the eastern coast of north america at distances greater than 110 observed features of s diffracted suggest that coupling between svd and shd can be ruled out as a first order effect for this path arrivals of shd are late relative to iasp91 travel times by about 10 s and those of svd are late relative to shd by 3 s for most records the slope of the log svd shd spectral ratio is around 3hz 1 in the range 0 06 0 15 hz a transversely isotropic low velocity layer in the lowermost mantle with a thickness of 200 300 km may account for most of the observed properties of svd

experimental study and tem characterization of dusty olivines in chondrites evidence for formation by in situ reduction an analytical transmission electron microscopy atem study was undertaken in order to better understand the formation conditions of dusty olivines i e olivines containing abundant tiny inclusions of fe ni metal in primitive meteorites dusty olivines from type i chondrules in the bishunpur chondrite ll3 1 and from synthetic samples obtained by reduction of san carlos olivines were examined in both natural and experimental samples micron size metal blebs observed in the dusty region often show preferential alignments along crystallographic directions of the olivine grains have low ni contents typically 2 wt and are frequently surrounded by a silica rich glass layer these features suggest that dusty olivines are formed by a sub solidus reduction of initially fayalitic olivines according to the following reaction fe2sio4 in olivine 2fe metal xsio2 in melt 1 x sioin gas 3 x 2o2 in gas some volatilization of siogas may account for the apparent excess of metal relative to silica rich glass observed in both experimental and natural samples comparison with experimentally produced dusty olivines suggests that time scales of the order of minutes usually inferred for chondrule formation are also adequate for the formation of dusty olivines these observations are in agreement with the hypothesis that at least part of the metal phase in chondrites originated from reduction during chondrule formation

small scatterers in the lower mantle observed at german broadband arrays seismograms of earthquakes from the south pacific recorded at a german broadband array and network show precursors to pkpdf these precursors mainly originate from off path scattering of pkpab or a nearby pkpbc to p for receiver side scattering or from scattering of p to pkpab or pkpbc on the pkpdf path for source side scattering standard array processing techniques based on plane wave approximations such as vespagram or frequency wavenumber analysis are inadequate for investigating these precursors since scattered waves cannot be approximated as plane waves for arrays and networks larger than 300 300 km for short period waves we therefore develop a migration method to estimate the location of scatterers in the mantle at the core mantle boundary and at the top of the outer core with our method we are able to find isolated scatterers at the source side and the receiver side although the depth of the scatterer is not well constrained however from looking at the first possible arrival time of precursors at different depth and the region where scattering can take place scattering volume we believe that the location of the scatterers is in the lowermost mantle since we have detected scatterers in regions where ultralow velocity zones have been discovered recently we think that the precursor energy possibly originates from scattering at partial melt at the base of the mantle comparing results from broadband and band pass filtered data the detection of small scale structure of the ultralow velocity zones becomes possible

the lowermost mantle beneath northern asia ii evidence for lower mantle anisotropy summary we have analysed prediffracted s waves with turning points beneath northern siberia in a study of anisotropy in the lowermost mantle deep focus earthquakes beneath the marianas izu bonin and the sea of japan recorded at stations in western europe are used a correction for upper mantle anisotropy is applied to the data comparisons of the data with synthetic data for models with and without a high velocity d layer suggest that there is a velocity discontinuity at the top of the d region and that the style of anisotropy is transversely isotropic in this region time separations between s waves on the radial and transverse component show a weak trend where the separation increases with epicentral distance a normalization of this separation with the travel distance within d 300 km thick in this region suggests that the anisotropy is uniformly distributed within this layer and has an average value of 0 5 per cent a combination of different studies which investigate the structure of the lowermost mantle beneath europe and northern siberia reveals a complicated picture tomographic models from this area and evidence of d anisotropy lower mantle scatterers reflections from a d discontinuity and ultralow velocity zones suggest two distinct regions one exhibits high velocities d anisotropy a d discontinuity and no evidence of scatterers or ultra low velocity zones these features are likely associated with the palaeosubduction of the izanagi plate well into the lowermost mantle the other region has a lower overall velocity and shows evidence of scatterers and ultra lowvelocity zones perhaps suggesting the presence of partial melt these results suggest dramatic lateral variations in the nature of the lowermost mantle beneath northern asia over a length scale of roughly 30 degrees

shdiff svdiff splitting in an isotropic earth 1 the vertically polarized svdiff seismic phase sometimes arrives later than the horizontally polarized shdiff both phases propagate in the d layer at the base of the mantle and the difference in their arrival times is usually interpreted in terms of seismic anisotropy in d using numerical simulations we demonstrate that a significant 1 s delay of svdiff relative to shdiff is present even in the iasp91 isotropic reference earth model and in a model with a high s velocity in d relative to iasp91 the delay is accompanied by strong amplitude decay of svdiff with distance this relationship resembles the effect of anelasticity in a model with a low s velocity in d the waveform of svdiff can be very different from shdiff and this difference can be mistaken for a delay of shdiff relative to svdiff in a laterally heterogeneous isotropic d svdiff can also be delayed because unlike shdiff svdiff has strong amplitude decay along fast paths its traveltime is thus dominated by slow propagation paths and tends to be positively biased for d models with an s velocity difference of 3 between two hemispheres there are regions where shdiff arrives up to 15 s earlier than the first clearly visible svdiff this is observed in a corridor several hundred kilometers wide if these effects are misinterpreted in terms of d anisotropy in actual seismic data estimates of anisotropy are biased and may contribute to the idea of seismic anisotropy being different in high and low s velocity regions of d

lowermost mantle anisotropy beneath the northwestern pacific evidence from pcs scs sks and skks phases 1 in this study we utilize data from broadband seismic stations of the japanese f net network to investigate anisotropic structure in the lowermost mantle beneath the northwestern pacific the comparison of shear wave splitting from phases that have similar paths in the upper mantle but different paths in the lowermost mantle such as pcs scs or sks skks can yield constraints on anisotropy in the d region we measured splitting for sks skks pcs and scs phases at f net stations and compared the measured splitting parameters to previously published estimates of upper mantle anisotropy beneath the region we observed many examples of discrepant sks skks splitting and splitting of scs pcs sks and skks phases that does not agree well with the known upper mantle anisotropy beneath individual stations the most likely explanation for these discrepancies is a contribution to splitting from anisotropy in the lowermost mantle in particular for sks skks phases that sample the lowermost mantle just to the south and east of the kamchatka peninsula we observed generally n s fast directions with delay times between 0 5 and 1 5 s these data suggest the presence of a fairly large coherent region of deformation in the lowermost mantle beneath the northwestern pacific our preferred model for these observations is that solid state flow at the base of the mantle induces a lattice preferred orientation of lowermost mantle minerals generating a seismically anisotropic fabric

deformation crystal preferred orientations and seismic anisotropy in the earth s d layer abstract we use a forward multiscale model that couples atomistic modeling of intracrystalline plasticity mechanisms dislocation glide twinning in mgsio3 post perovskite ppv and periclase mgo at lower mantle pressures and temperatures to polycrystal plasticity simulations to predict crystal preferred orientations cpo development and seismic anisotropy in d we model the cpo evolution in aggregates of 70 ppv and 30 mgo submitted to simple shear axial shortening and along corner flow streamlines which simulate changes in flow orientation similar to those expected at the transition between a downwelling and flow parallel to the core mantle boundary cmb within d or between cmb parallel flow and upwelling at the borders of the large low shear wave velocity provinces llsvp in the lowermost mantle axial shortening results in alignment of ppv 010 axes with the shortening direction simple shear produces ppv cpo with a monoclinic symmetry that rapidly rotates towards parallelism between the dominant 100 010 slip system and the macroscopic shear these predictions differ from mgsio3 post perovskite textures formed in diamond anvil cell experiments but agree with those obtained in simple shear and compression experiments using cairo3 post perovskite development of cpo in ppv and mgo results in seismic anisotropy in d for shear parallel to the cmb at low strain the inclination of scs sdiff and skks fast polarizations and delay times vary depending on the propagation direction at moderate and high shear strains all s waves are polarized nearly horizontally downwelling flow produces sdiff scs and skks fast polarization directions and birefringence that vary gradually as a function of the back azimuth from nearly parallel to inclined by up to 70 to cmb and from null to 5 change in the flow to shear parallel to the cmb results in dispersion of the cpo weakening of the anisotropy and strong azimuthal variation of the s wave splitting up to 250 km from the corner transition from horizontal shear to upwelling also produces weakening of the cpo and complex seismic anisotropy patterns with dominantly inclined fast scs and skks polarizations over most of the upwelling path models that take into account twinning in ppv explain most observations of seismic anisotropy in d but heterogeneity of the flow at scales

formation of feldspathic and metallic melts by shock in enstatite chondrite reckling peak a80259 the enstatite chondrite reckling peak rkp a80259 contains feldspathic glass kamacite troilite and unusual sets of parallel fine grained enstatite prisms that formed by rapid cooling of shock melts metallic fe ni and troilite occur as spherical inclusions in feldspathic glass reflecting the immiscible fe ni s and feldspathic melts generated during the impact the fe ni s and feldspathic liquids were injected into fractures in coarse grained enstatite and cooled rapidly resulting in thin 10 m semicontinuous to discontinuous veins and inclusion trails in host enstatite whole rock melt veins characteristic of heavily shocked ordinary chondrites are conspicuously absent raman spectroscopy shows that the feldspathic material is a glass elevated mgo and sio2 contents of the glass indicate that some enstatite and silica were incorporated in the feldspathic melt metallic fe ni globules are enclosed by sulfide and exhibit nienrichment along their margins characteristic of rapid crystallization from a fe ni s liquid metal enclosed by sulfide is higher in si and p than metal in feldspathic glass and enstatite possibly indicating lower o fugacities in metal sulfide than in silicate domains fine grained elongate enstatite prisms in troilite or feldspathic glass crystallized from local pyroxene melts that formed along precursor grain boundaries but most of the enstatite in the target rock remained solid during the impact and occurs as deformed coarsegrained crystals with lower cao al2o3 and feo than the fine grained enstatite reckling peak a80259 represents an intermediate stage of shock melting between unmelted e chondrites and whole rock shock melts and melt breccias documented by previous workers the shock petrogenesis of rkpa80259 reflects the extensive impact processing of the enstatite chondrite parent bodies relative to those of other chondrite types

properties of chondrules in el3 chondrites comparison with eh3 chondrites and the implications for the formation of enstatite chondrites the study of chondrules provides information about processes occurring in the early solar system in order to ascertain to what extent these processes played a role in determining the properties of the enstatite chondrites the physical and chemical properties of chondrules from three el3 chondrites and three eh3 chondrites have been examined by optical cathodoluminescence cl and electron microprobe techniques properties examined include size texture cl and composition of both individual phases and bulk chondrules the textures distribution of textures and composition of silicates of the el3 chondrules resemble those of eh3 chondrules however the chondrules from the two classes differ in that 1 the size distribution of the el chondrules is skewed to larger values than eh chondrules 2 the enstatite in el chondrules displays varying shades of red cl due to the presence of fine grained sulfides and metal in the silicates and 3 the mesostasis of eh chondrules is enriched in na relative to that of el chondrules the similarities between the chondrules of the two classes suggest similar precursor materials while the differences suggest that there was not a single reservoir of meteoritic chondrules but that their origin was fairly local the differences in the size distribution of chondrules in eh and el chondrites may be explained by aerodynamic and gravitational sorting during accumulation of the meteoric material while differences in cl and mesostasis properties may reflect differences in formation conditions and cooling rate following chondrule formation we argue that our observations are consistent with the formation of enstatite chondrites in a thick dynamic regolith on their parent body

apparent splitting of s waves propagating through an isotropic lowermost mantle observations of shear wave anisotropy are key for understanding the mineralogical structure and flow in the mantle several researchers have reported the presence of seismic anisotropy in the lowermost 150 250 km of the mantle i e d urn x wiley jgrb media jgrb52636 jgrb52636 math 0002 layer based on differences in the arrival times of vertically sv and horizontally sh polarized shear waves by computing waveforms at a period 6 s for a wide range of 1 d and 3 d earth structures we illustrate that a time shift i e apparent splitting between sv and sh may appear in purely isotropic simulations this may be misinterpreted as shear wave anisotropy for near surface earthquakes apparent shear wave splitting can result from the interference of s with the surface reflection ss for deep earthquakes apparent splitting can be due to the s wave triplication in d urn x wiley jgrb media jgrb52636 jgrb52636 math 0003 reflections off discontinuities in the upper mantle and 3 d heterogeneity the wave effects due to anomalous isotropic structure may not be easily distinguished from purely anisotropic effects if the analysis does not involve full waveform simulations

a very slow basal layer underlying large scale low velocity anomalies in the lower mantle beneath the pacific evidence from core phases a multi phase analysis using long period world wide standardized seismograph network and canadian network data has been conducted using core phases for deep focus events from the southwest pacific these include sks s2ks svd iff and sp dks the last phase emerges from sks near 106 and is associated with a p wave diffracting along the bottom of the mantle patterns in s2ks sks differential travel times t s2ks sks correlate with those in sp dks sks t spdks sks t s2ks sks values strongly depend on variations in v s structure in the lower third of the mantle whereas t spdks sks values mainly depend on v p structure and variations in a thin zone 100 km or less at the very base of the mantle anomalously large t s2ks sks and t spdks sks values relative to the preliminary reference earth model prem are present for fiji tonga and kermadec events recorded in north and south america along with anomalously large sv diff amplitudes well into the core s shadow more northerly paths beneath the pacific to north america for indonesian and solomon events display both prem like and anomalous times a model compatible with the observations is presented and contains a thin very low velocity layer at the base of the mantle that underlies the large volumetric lower mantle low velocity regions in the southwest pacific a low velocity layer of 20 100 km thickness with reductions of up to 5 10 relative to prem can reproduce t spdks sks as well as sv diff amplitudes large scale more than 1000 km lower mantle v s heterogeneity 2 4 can explain long wavelength trends in t s2ks sks the exact thickness and velocity reduction in the basal layer is uncertain owing to difficulties in resolving whether anomalous structure occurs on the source and or receiver side of wavepaths at the cmb

p and s wave reflections from anomalies in the lowermost mantle summary seismic p and s waves recorded at the grf array in germany are used to study the inhomogeneous structure of the boundary layer d at the base of the mantle the use of the seismic array allows the detection of small scale anomalies in the lowermost mantle the lateral resolution attainable is about 100 to 400 km i e more than 10 times better than with tomographic methods the analysis of 13 years of grf broad band array data yields 255 events with high signal to noise ratio which are used to map the lowermost mantle 74 of these events show anomalous p waves pdp which arrive 3 6 s after the direct p wave and have a slowness 0 7 0 8 s smaller than the slowness of direct p other events close to the events with pdp do not have such an anomalous phase using slowness traveltime amplitude and waveform information it is demonstrated that pdp is caused by an anomalous lower mantle velocity structure below the turning point of the p wave if the pp phase i e a p wave first reflected at the free surface near the source is used together with the p phase distinct and well sepaarted p velocity anomalies can be determined under the nansen basin the kara sea and northern siberia the areas of the bounce points of pdp in the lower mantle have a lateral extension of about 100 by 200 km but this is not the size of the anomaly since the resolution of p waves at 1 hz in this depth is 130 km by 260 km fresnel zone the accuracy to which the depth of the reflector 2612 km under the nansen basin and 2605 km under northern siberia respectively can be determined for 1 d models is 10 to 20 km the velocity contrast at the lower mantle discontinuity is about 3 per cent 1 to 1 5 per cent areas which have velocity fluctuations smaller than about 1 per cent can not be detected as anomalous areas 2 d models of the anomalies reveal the range of adequate models and possible trade offs if the lateral extension of the anomaly is about 7 the reflector has to be 40 km deeper than in the 1 d model if the dip of the reflector in the lowermost mantle is only about 1 5 it is difficult to resolve if the reflector is tilted towards the source or towards the receiver for the anomaly under the nansen basin deviations from the great circle path are observed for pdp indicating 3 d effects below the three anomalies the core mantle boundary cmb can be located by pcp with a depth in agreement with standard earth models under northern siberia has a region where p and s velocity anomalies coincide but also a region where only an s velocity anomaly but no p velocity anomaly is observed this results in changes of the poisson ratio from 5 9 per cent to 4 8 per cent 1 5 per cent across the discontinuity in the lowermost mantle for the regions studied the analysis of the s waves s sds and scs reveals two s velocity anomalies in the lowermost mantle under northern siberia the first anomaly coincides with the p velocity anomaly under northern siberia and can be explained by a 1 d model with a reflector depth of 2610 km 15 km and a velocity contrast of 2 3 per cent the second s velocity anomaly is in an area where no p velocity anomaly can be detected the corresponding 1 d s velocity model has a reflector depth of 2575 km 15 km and an s velocity contrast of 2 6 per cent the smallest structures that can be resolved with the s waves in this depth are about twice as large as for the p waves i e 230 km by 460 km fresnel zone the joint analysis of p and s waves therefore shows a region with a p velocity anomaly together with weak indications for an s velocity anomaly nansen basin and a second region with a p but no s anomaly kara sea the lowermost mantle

tomographic imaging of the lowermost mantle with differential times of refracted and diffracted core phases pkp pdiff the mapping of variations in p wave speed in the deep mantle is restricted by the uneven sampling of p waves in particular beneath the southern hemisphere to enhance data coverage we augmented the 1 6 million summary rays of p pp and pwp that we used in previous studies with differential travel times of diffracted and refracted core phases for the core refracted differential travel time residuals pkpab pkpdf and pkpab pkpbc we used 1383 cross correlated digital waveforms as well as 27 000 routinely processed bulletin data we used the waveform data to define quality criteria for the selection and reduction of the bulletin pkp data for pkpdf pdiff we only considered 543 records derived from waveform cross correlation no pcp data were used in this study we used optical ray theory to calculate the ray paths associated with the p pp pwp and pkp data which are measured at 1 hz however to account for the large fresnel zones of the low frequency 50 mhz pkpdf pdiff data we estimated the three dimensional shape of the frechet sensitivity kernels from kernels calculated by normal mode summation the use of these kernels allows us to properly distribute the sensitivity for a given seismic phase over a large mantle volume while allowing the high frequency data to constrain small scale structure the differential times are relatively insensitive to source mislocation and to structure in the shallow mantle beneath source and receivers and they have previously been interpreted exclusively in terms of lateral structure directly above the core mantle boundary cmb however images thus obtained can be contaminated by effects of small scale structure elsewhere in the mantle here we do not make a priori assumptions about the mantle source of anomalous time differentials from test inversions we conclude that both upper and lower mantle structures that are poorly resolved by p data can be mapped into the core along pkp paths but that the effect of outer core structures if any on the mantle model is small compared to the inversion of the p pp and pwp alone the inclusion of the pkpab pkpdf and pkpab pkpbc and pkpdf pdiff data improves the resolution of structure beneath 2200 km depth in particular the joint inversion puts better constraints on the long wavelength variations in the very deep mantle and yields an increase in the amplitude of velocity perturbations near the cmb that is in agreement with but still smaller than inferences from shear wave studies resolution tests indicate that in some regions the enhanced definition of structure is significant but in most regions the improvements are subtle and structure remains poorly resolved in large regions of the mantle

seismicity of the bora tullu moye volcanic field 2016 2017 the bora tullu moye tm volcanic field is a geothermal energy prospect in the central main ethiopian rift but little is known about the seismicity of the region here we document seismic activity between february 2016 and october 2017 locating more than a 1 000 seismic events of local magnitude 0 to 2 7 this provides new insights into fluid movement and deformation beneath what we only now realize is a complicated volcanic system a discrete cluster of events lies beneath tm but surprisingly most of the seismicity lies in two clusters that are beneath neither the bora or tm edifices in these regions we use earthquake cluster orientations fault plane solutions and fast seismic shear wave orientations to show that seismicity is triggered by hydrothermal circulation of fluids along preexisting fractures the fractures trend in multiple directions and are in general not parallel to rifting related structures instead the fractures are parallel to structures created during previous caldera forming eruptions at both bora and tm highly fractured regions such as this could be attractive targets for geothermal power generation we estimate a minimum depth for a magmatic body beneath tm to be 6 5 km using the mapped brittle ductile transition frequency analysis of the earthquake waveforms reveal the bulk of the events to be volcano tectonic but some low frequency seismicity is present at a depth of 5 km beneath the tm edifice triggered by high pore fluid pressures

a strong lateral shear velocity gradient and anisotropy heterogeneity in the lowermost mantle beneath the southern pacific 1 velocity heterogeneity and seismic anisotropy in the d region beneath the southern pacific ocean is investigated using broadband shear waves from southwest pacific and south american earthquakes recorded by permanent and temporary stations in antarctica australia south america and on islands in the south pacific shear velocity perturbations vs are inferred from the differential travel times between hand picked horizontally polarized sh direct or diffracted s arrivals and vertically polarized sv sks arrivals derived patterns in vs roughly agree with global shear velocity models but reveal a stronger shear velocity lateral gradient north of approximately 52 s where vs transitions from approximately 0 5 to 1 0 relative to radially averaged reference models over less than 600 km along the base of d waveform analyses provide an even stronger constraint on the transitional region where there is a change in waveform behavior occurring over length scales less than 300 km differential travel times between skks and sks calculated by cross correlating sks with hilbert transformed skks support large vs amplitudes and help to constrain the location of the transitional region anisotropy in d ks is inferred from handpicked differential travel times between the sv and sh components of direct or diffracted s and display a slight spatial trend where the sh component precedes the sv in the north and east of the study area but arrives after sv in the center and southwest of it however there is high ks variability in the center of the study region there is no apparent correlation between vs and ks abrupt changes in the character of velocity and anisotropy in this region of d may be related to chemical heterogeneity at the boundary of the pacific superswell as well as small scale convection in the deep mantle

the lowermost mantle beneath northern asia i multi azimuth studies of a d heterogeneity summary in most areas of the globe it is only possible to sample the lowermost mantle or d region with a limited range of azimuths this is due to limitations in source receiver geometries an exception is the region beneath northern siberia where we are able to sample d along orthogonal great circle paths using northern pacific earthquakes recorded in europe and hindu kush earthquakes recorded in canada receiver vespagrams slant stacks of one source recorded at many receivers for a kurile event recorded at the german regional seismic network grsn and the grafenberg array and source vespagrams many sources one receiver of hindu kush events recorded at stations of the canadian national seismograph network cnsn show reflections from a discontinuity at the top of d vespagrams for real data are compared to vespagrams for synthetic data in order to guide their interpretation additional support for our interpretation comes from frequency wavenumber f k analysis which provides a measure of the slowness and azimuth of a given phase the traveltimes of the d reflected phases are then used to provide estimates of d thickness in regions where the roughly orthogonal paths overlap we estimate the thickness to be 310 15 km which is in agreement with existing models for this region we also estimate a thickness of 210 15 km in a previously unsampled region east of this area using the slowness backazimuth and traveltime values for a third source receiver combination we find a reflection point at a depth of 2800 km thickness of d of 91 km in a region several degrees to the west these results and those of previous studies suggest considerable topographic relief 100 km over distances of less than 5 approximately 300 km at the cmb comparing the results from this study with results from other seismic studies of d structure in this area we find this to be a complicated region which exhibits heterogeneity on many scales this region may be the lower mantle graveyard of the subducted izanagi plate and one interpretation of this complexity is that it is due to subduction related deformation at the base of the mantle

inferences on flow at the base of earth s mantle based on seismic anisotropy we applied global waveform tomography to model radial anisotropy in the whole mantle we found that in the last few hundred kilometers near the core mantle boundary horizontally polarized s wave velocities vsh are on average faster by 1 than vertically polarized s wave velocities vsv suggesting a large scale predominance of horizontal shear this confirms that the d region at the base of the mantle is also a mechanical boundary layer for mantle convection a notable exception to this average signature can be found at the base of the two broad low velocity regions under the pacific ocean and under africa often referred to as superplumes where the anisotropic pattern indicates the onset of vertical flow

characterising hydrothermal fluid pathways beneath aluto volcano main ethiopian rift using shear wave splitting abstract geothermal resources are frequently associated with silicic calderas which show evidence of geologically recent activity hence development of geothermal sites requires both an understanding of the hydrothermal system of these volcanoes as well as the deeper magmatic processes which drive them here we use shear wave splitting to investigate the hydrothermal system at the silicic peralkaline volcano aluto in the main ethiopian rift which has experienced repeated uplift and subsidence since at least 2004 we make over 370 robust observations of splitting showing that anisotropy is confined mainly to the top 3 km of the volcanic edifice we find up to 10 shear wave anisotropy swa is present with a maximum centred at the geothermal reservoir fast shear wave orientations away from the reservoir align nne ssw parallel to the present day minimum compressive stress orientations on the edifice however are rotated ne sw in a manner we predict from field observations of faults at the surface providing fluid pressures are sufficient to hold two fracture sets open these fracture sets may be due to the repeated deformation experienced at aluto and initiated in caldera formation we therefore attribute the observed anisotropy to aligned cracks held open by over pressurised gas rich fluids within and above the reservoir this study demonstrates that shear wave splitting can be used to map the extent and style of fracturing in volcanic hydrothermal systems it also lends support to the hypothesis that deformation at aluto arises from variations of fluid pressures in the hydrothermal system these constraints will be crucial for future characterisation of other volcanic and geothermal systems in rift systems and elsewhere

evaluation of image classification routines for determining modal mineralogy of rocks from x ray maps the mineral proportions of rocks or rock fragments in thin section can be retrieved from electron microprobe element x ray maps using multispectral image classification software however different image classification algorithms can yield different inferred mineral proportions and thus deviate from the true modal mineralogy several image classification algorithms implemented in the programs erdas imagine and multispec were evaluated on a very simple rock a fragment of noritic anorthositic granulite in a thin section of lunar meteorite alha81005 this fragment contains only plagioclase orthopyroxene olivine clinopyroxene and chromite each has a constant composition across the fragment the true reference mineral proportions of the fragment were measured manually from the x ray maps results of classification algorithms varied widely in accuracy from 1 to more than 15 misclassified pixels compared to the reference the best match for this fragment comes from an unsupervised classification using the isodata algorithm most minerals were placed into multiple classes which had to be manually recombined and the second closest match was a supervised classification using a euclidean distance classifier the results of supervised classifications depend strongly on the user s selection of training areas and small errors in defined training areas e g a training area for one phase contains an inclusion of another can produce large errors the most common source of misclassification is mixed pixels those that contain signatures from multiple phases in this example pixels including both plagioclase and olivine tend to be classified as clinopyroxene thus mineral proportions calculated from classification of x ray maps must be considered critically even for simple rocks in more complex samples like those with more phases or with chemically zoned phases mineral classification of x ray maps are likely to require care and may be difficult to validate

in situ observation of texture development in olivine ringwoodite magnesiowustite and silicate perovskite at high pressure magnesium silicates are the dominant minerals in the earth s mantle their preferred orientation is important for understanding the rheology and seismic anisotropy in the deep earth here we report results of radial synchrotron diffraction diamond anvil cell dac experiments on san carlos olivine axially compressed to 50 gpa experiments were performed at room temperature except for brief laser heating to induce phase transformations high stresses and development of preferred orientation were observed in diffraction images quantitative texture information was obtained by analyzing the images with the rietveld method with increasing pressure between 9 and 43 gpa olivine develops a texture with 001 axes perpendicular to the compression direction that is compatible with hk0 001 pencil glide ringwoodite between 20 and 50 gpa develops weak preferred orientation with 011 lattice planes perpendicular to the compression direction after the phase transformation perovskite and magnesiowqstite display transformation textures that are then modified by continuing deformation magnesiowqstite has a weak h111i maximum parallel to the compression direction that changes towards h001i with increasing deformation perovskite transforming from olivine has a pronounced 100 transformation texture and with increasing deformation a 012 maximum develops the transformation texture is probably produced by mechanical 110 twinning and nucleation in orientations that minimize elastic strain energy the deformation texture of perovskite is due to slip d 2004 elsevier b v all rights reserved

lateral variations in lowermost mantle shear wave anisotropy beneath the north pacific and alaska s waves recorded by long period world wide standardized seismograph network and broadband stations in north america for deep northwest pacific subduction zone earthquakes provide evidence for anisotropy in the lowermost mantle shear velocity structure beneath the north pacific and alaska systematic delays of up to 4 s are observed between longitudinal components sv and transverse components sh of motion for core reflected scs waves as well as for core grazing and diffracted s waves the absence of significant splitting for s waves that have turning points more than a few hundred kilometers above the core mantle boundary indicates that anisotropy is localized within the d region the lowermost portion of the mantle sv sh differential arrival times for both scs and sdiff along with path length estimates assuming a 250 km thick d region indicate spatial variations in the strength of shear wave anisotropy the strongest anisotropy 1 1 5 is found in the eastern part of the study area with systematic reduction in magnitude toward the west a transverse isotropy model can explain the data with the velocity structure for horizontally polarized waves vsh having a 2 3 discontinuous shear velocity increase at the top of d as proposed in earlier studies of the region and a similar structure for s wave particle motion in the direction normal to the core mantle boundary vsv but with the velocity jump at the top of d and the velocity within d being reduced from that for vsh by 0 5 1 5 km s large uncertainties exist for velocity gradients above and below the velocity jump but the requirement of a reduced vsv relative to vsh in d is clear synthetic waveforms calculated by using separate isotropic structures for sh and sv match the observations well and constrain the basic anisotropic structure because the shear waves all traverse the region with near grazing geometries the study area exhibits strong lateral variations in lower mantle shear velocity structure and variable thickness of the d layer topography of the d layer is not well resolved because of trade off with volumetric heterogeneity

constraints from seismic anisotropy on the nature of the lowermost mantle the d layer lies at the bottom of the earth s rocky mantle and separates it from the liquid metal alloy core this region extending from the core mantle boundary to a few hundred kilometres above fig 1 is geodynamically analogous to the more easily studied lithosphere at the top of the mantle the structure of d may reflect the style of lower mantle convection the nature of core mantle interaction and perhaps even the fate of subducting lithosphere1 observations of lithospheric seismic anisotropy have provided valuable insight into the nature of the upper mantle boundary layer but discussion of lower mantle seismic anisotropy has been somewhat contentious2 5 here we present evidence from seismic waves that have traversed the lowermost mantle beneath the caribbean region for a zone of seismic anisotropy below the d discontinuity which in this region lies 250 km above the core mantle boundary the anisotropy is most probably due to horizontal layering or aligned inclusions of a material with differing shear wave velocity if d is a graveyard for subducted lithosphere a plausible explanation of the anisotropy may be the contrast between cold lithospheric mantle and material that formerly constituted the oceanic crust which may have lower shear wave velocity owing to the presence of melt

evidence for anisotropy in the deep mantle beneath alaska we consider the possibility that the velocity structure of d is anisotropic the data we examined consist of seismograms from 9 deep japanese earthquakes recorded at wwssn receiver stations in north america the source receiver combinations span distances of 70 106 with associated s waves passing through d beneath alaska differential travel times of the s scd scs and sks phases are used to constrain the velocity structure in d shear waves refracted by d are observed beyond 72 2 and provide a sensitive measurement of the velocity structure in d beyond 93 the vertically polarized sv and horizontally polarized sh shear waves often appear distinctly split although at distances less than 89 the components are more nearly synchronous near 94 sh occurs as a double arrival sv in this range however remains a single arrival roughly synchronous with the second sh arrival we have been unable to reproduce these effects in isotropic model synthetics synthetics for transversely isotropic models have been computed that do match these waveforms the anisotropy was constrained to be only within d with a vertical symmetry axis we conclude that these observations may be explained by an anisotropic d layer the d discontinuity may be due to a transition to anisotropic mantle a few hundred kilometers above the core mantle boundary