mobility and fractionation of rare earth elements in argillaceous sediments implications for dating diagenesis and low grade metamorphism abstract we report sm nd and rb sr data for the fine fractions of lower paleozoic argillaceous rocks from wales uk and new york usa spanning the range of low grade metamorphic conditions from the diagenetic zone zeolite facies to the epizone greenschist facies in all cases leaching of the fine fractions results in a high 147 sm 144 nd 0 09 0 29 acid soluble component and a complementary low 147 sm 144 nd 0 05 0 14 residual component the observed fractionation is an ancient feature related to diagenesis burial and metamorphism the magnitude of sm nd fractionation between leachates and residues as well as the resulting sm nd ages vary as a function of grain size and metamorphic grade uncleaved welsh mudrocks of the diagenetic zone yield sm nd leachate residue ages of 453 484 ma in agreement with their llanviian to caradocian biostratigraphic ages whereas higher grade rocks of the anchizone and epizone yield sm nd ages as young as 413 ma these ages are transitional between the time of deposition and the time of regional deformation related to the acadian orogeny at 390 ma distinct convex upward rare earth element ree patterns of the leachates suggest that the precipitation of early diagenetic apatite controls the trace element budget of the rock forcing a depletion of middle rees on the subsequently formed diagenetic phyllosilicates the amount of organic matter present and the extent of later prograde reactions are probable modifiers of this fractionation process ordovician and devonian clastic rocks associated with the trenton and onondaga limestones of new york yield single sample and multi sample sm nd isochron ages that agree well with their biostratigraphic ages of 454 ma and 390 ma respectively the ree fractionation observed in shale leachates of the ordovician utica formation is related to ca mg of the bulk rock and hence to the composition of the diagenetic carbonate cement in all cases the sm nd system remained closed subsequent to the peak of diagenesis or metamorphism including the north american rocks that show no evidence of being isotopically reset during widespread remagnetization of the subjacent limestone units in the late paleozoic

intercalibration of 40ar39ar dating standards a number of standards used for k ar and 40ar39ar dating have been intercalibrated multiple splits 15 40 mg each of mmhb 1 sb 3 biotite lp 6 biotite 40 60 ga1550 biotite fish canyon tuff 3 biotite taylor creek rhyolite sanidine bern 4 muscovite bern 4 biotite and pennsylvania state university orthoclase 1a were irradiated in capsules fused and analyzed in replicate analyses most of these standards proved to be homogeneous in 40ar 39ark ratio age k ca and 36ar contents mmhb 1 is inhomogeneous in age at the 15 mg level indicating it is unsuitable for use as an interlaboratory standard using material taken straight from the bottle the orthoclase sample is also inhomogeneous at the 30 mg level quoted relative to an age of 162 9 ma for the sb 3 biotite standard total fusion ages are lp 6 biotite 40 60 128 1 ma ga1550 biotite 97 8 ma fish canyon tuff 3 biotite 27 95 ma taylor creek rhyolite sanidine 28 0 ma bern 4 muscovite 18 51 ma bern 4 biotite 17 21 ma these ages are in good agreement with the absolute k ar ages of these minerals different preparations of biotite and sanidine standards from a welded tuff such as the fish canyon tuff may not be of identical ages tests on submilligram splits reveals that sb 3 biotite ga1550 biotite fish canyon tuff 3 biotite and taylor creek rhyolite sanidine are suitable for use as monitors in laser work the sanidine specimen appears to be the most suitable for use as the primary monitor for 40ar39ar dating

resetting of rb sr whole rock ages during acadian low grade metamorphism in north wales ordovician volcanic rocks and associated sub volcanic intrusions in north wales give devonian metamorphic rb sr whole rock ages isotopic resetting is not confined to acid volcanic rocks three intrusions preserve ordovician ages the closed system resetting is not directly related to metamorphic grade but appears to correlate with the extent to which the secondary metamorphic mineral assemblage is developed rocks giving devonian ages have a well developed secondary mineralogy whereas those retaining ordovician ages are less extensively altered the metamorphic ages combine to give a lower devonian mean age and weighted 2a error of 399 f 9 ma metamorphism is therefore synchronous with acadian deformation in the north wales area

short paper low grade metamorphism of the welsh basin lower palaeozoic succession an example of diastathermal metamorphism the lower palaeozoic succession of the welsh basin has suffered the effects of low grade metamorphism which shows features that do not accord with classical models of regional or burial metamorphism textural evidence demonstrates the development of greenschist facies metamorphism prior to deformation and cleavage development while mineralogical characters indicate a low pressure facies series in order to accommodate these features it is suggested that the area has suffered diastathermal metamorphism in which an enhanced thermal flux develops early in the tectonic cycle in response to an extensional setting temperatures remained relatively high during final closure associated with deformation of the basin sediments and consequent cleavage development such a model adequately resolves previously conflicting interpretations of burial and deformation related metamorphism

tem and xrd determination of crystallite size and lattice strain as a function of illite crystallinity in pelitic rocks average crystallite size and mean square strain of illite in rock specimens and clay separates were measured independently in tem images and by single line fourier voigt method profile analysis of the c 1 nm peak of xrd patterns for a prograde sequence of pelitic rocks illite crystallinity indices 0 17 0 58 d2h from the gaspepeninsula quebec the tem determined crystallite sizes in clay separates approxi mate those determined by fourier profile analyses and those calculated from illite crystallinity indices by the scherrer equation with the exception of the diagenetic sample the crystallite sizes and mean square strains of illite in rock samples exhibit a trend similar to that determined by profile analyses but the average crystallite sizes are up to five times larger than those measured for clay separates tem images show that all rock samples have a wide range of crystallite sizes and the proportions of larger crystallites increase with metamorphic grade the diagenetic illite is defect rich fine grained mean thickness by volume c 70 nm 1m d material anchizonal illite tends to occur as separate aggregates of small 1m d and larger 2m 1 crystals c 200 nm comprising arrays of subparallel coalescing packets the epizone sample has thick c 400 nm defect free crystals of muscovite occurring in stacks of parallel layers or subhedral crystals intergrown with large angle boundaries cleaved crystals that are free of intracrystalline layer terminations are dominant in clay separates of all samples having ranges of smaller sizes with volume average thicknesses of c 43 43 and 81 nm c 14 28 67 nm by the voigt method respectively for the three zones the results suggest that illite crystallinity indices do not provide a direct measure of a single microstruc tural state of illite in rocks although they yield consistent limits for average crystallite sizes for the anchizone 23 48 nm here therefore they serve as a general parameter of the degree of recrystallization on a relative basis in part because the contributions of all peak broadening variables mixed layering size and strain decrease regularly with prograde regional metamorphism of pelites the microstructural changes caused by rock disaggregation are probably a function of those variables as well the data collectively demonstrate a trend from metastable defect rich small crystals towards a stable assemblage of larger defect free crystals through dissolution of strained crystals and neocrystallization consistent with the ostwald step rule

shear wave splitting in the upper mantle wedge above the tonga subduction zone summary intermediate period seismograms of shear waves from deep earthquakes are examined for shear wave splitting diagnostic of seismic velocity anisotropy in the upper mantle wedge between the subducting tonga slab and stations in the fiji islands earthquakes with near vertical incidence angles are selected to minimize complications of converted phases and the free surface horizontal component seismograms for 25 deep earthquakes recorded at station lak on the lau ridge show splitting of up to 0 9 s with an average of 0 45 0 25 s between fast and slow s wave components anisotropy in the upper mantle is the most likely cause of the splitting the ese orientation of the fast velocity direction at station lak corresponds to an axis of symmetry of the anisotropic mantle and is approximately perpendicular to magnetic lineations in the lau basin and the northeastern portion of the south iiji basin the observations are consistent with a model of flow alignment of olivine 100 and orthopyroxene 001 axes parallel to the direction of basin extension seismograms for eight deep earthquakes at station sva on viti levu show up to 1 0 s splitting with an average of 0 65 0 29 s and a weak nne preferred orientation of fast velocity polarizations of shear wave particle motions corrected for anisotropy are within 36 5 26 of polarixations predicted for nine moment tensor solutions the magnitude of shearwave splitting and estimates of the thickness of the lithosphere in the fiji region imply that the anisotropy cannot be contained in the lithosphere alone either the anisotropy is distributed in both the lithosphere and asthenosphere or the asthenosphere alone is anisotropic

spatial variation in the crustal anisotropy and its temporal variation associated with a moderate sized earthquake in the tokai region central japan summary spatial and temporal variations in seismic anisotropy in the crust are investigated using earthquakes in the crust and the subducting philippine sea plate beneath the tokai region central japan we use a total of 351 high quality waveform data recorded from 1986 december until 1999 august by the microearthquake observation network of the research center for seismology and volcanology nagoya university the leading shear waves are polarized in an approximately e w direction independent of the focal depth showing consistency with the direction of regional horizontal maximum compressive stress in the tokai region the time delay increases for greater focal depth down to 30 km these results indicate that the regional compressive stress controls anisotropy not only in the upper crust but also in the lower crust assuming a uniform distribution of anisotropy the degree of anisotropy is estimated to be 0 6 per cent in the lower crust an increase in time delay is found after the aichi ken tobu earthquake m 5 7 in 1997 at station stn this variation is statistically significant and is not an apparent change due to the variation in hypocentre distribution or any kind of artificial effect on the other hand no temporal variation in time delay is found at station inu the spatial pattern of the temporal variation in time delay is consistent with the change in the coulomb failure function cff caused by both coseismic and post seismic fault slips we therefore conclude that the observed temporal change in time delays might be caused by increasing crack density as well as pore fluid pressure due to the static stress change accompanied with the aichi ken tobu earthquake

40ar 39ar illite dating of late caledonian acadian metamorphism and cooling of k bentonites and slates from the welsh basin u k the vacuum encapsulation laser 40ar39ar technique allows extremely small 10 6 g samples of fine grained materials such as diagenetic clays to be dated here we show that the method can be extended to higher grade clay minerals the integration of transmission electron microscopic tem characterization with 40ar39ar dating of vacuum encapsulated samples permits the resolution of the timing of metamorphic growth cooling from the time of diagenesis we have applied this technique to well characterized lower paleozoic slates and k bentonites from the welsh basin which span the transition from anchizonal to epizonal grade which had been previously studied using rbsr and smnd dating tem observations of epizonal k bentonites and slate showed that illite in these samples is of 2m1 polytype of muscovite like composition and oriented parallel to cleavage suggesting that they are of metamorphic origin total gas ages equivalent to conventional kar ages for encapsulated epizonal k bentonites and slate 340 408 ma are considerably variable the ar retention ages calculated from 39ar and 40ar atoms retained in the sample after irradiation are more consistent 383 411 ma the 39ar recoil losses are minor for illites from whole rock samples of epizonal k bentonites but very significant for clay separates of epizonal slate plateaus in age spectra were observed in epizonal k bentonites and slate the plateau ages 414 421 ma and retention ages 383 411 ma can be correlated with the onset of acadian metamorphism and culmination of uplift and inversion of the welsh basin respectively these ages are significantly younger than the 450 ma ages previously reported for diagenetic clays using the same method suggesting that diagenetic history has been lost in these epizonal k bentonites and slate tem observations of anchizonal slates showed that there are two modes of illite the first mode is similar to that observed in epizonal samples suggesting a metamorphic origin the second mode consists of the 1md polytype has typical diagenetic illite composition and is oriented parallel to bedding suggesting a diagenetic origin total gas ages for encapsulated anchizonal slates vary considerably 361 422 ma the retention ages are more consistent 413 432 ma than the total gas ages the 39ar recoil losses are more significant than those for epizonal k bentonites and slate plateaus in age spectra are generally not observed however the consistent retention ages for the anchizonal slates correspond to the plateau ages for the epizonal samples and are inferred to represent the onset of acadian metamorphism these data when combined with our previously published results for diagenetic shales suggest that thermal conditions near the boundary of anchizonal and epizonal grades are necessary to completely reset ar systems in shales and slates

shear wave polarization anisotropy and subduction induced flow in the mantle wedge of northeastern japan abstract shear wave polarization anisotropy in the central part of northeastern ne japan was investigated we analyzed s phases of 293 intermediate depth earthquakes recorded at 77 stations and obtained 1286 splitting parameters the obtained results show the systematic variation in splitting parameters across the arc the leading shear wave polarization directions fast direction obtained at stations in the western side of the study area are oriented nearly e w which is sub parallel to the direction of relative plate motion between the pacific plate and the north american okhotsk plate in contrast fast directions obtained at stations in the eastern side of the study area are oriented approximately n s stations in the western side have a larger average delay time compared to those in the eastern side these variations suggest that the nature of anisotropy is quite different between the eastern and the western sides both experimental studies on the deformation of olivine and numerical simulations of plate driven flow predict that the lattice preferred orientation of olivine causes the fast direction sub parallel to the flow direction for the conditions in the mantle wedge of ne japan we thus infer that the lattice preferred orientation of olivine generated by flow induced strain is the most likely candidate for anisotropy that produces splitting with e w fast direction three possible causes of anisotropy in the eastern side are considered deformation of water rich olivine in the mantle wedge trench parallel flow in the mantle wedge due to along strike dip variations in the slab and anisotropy in the crust and the slab we could not however determine which causes are dominant further systematic study is required to reveal the cause of anisotropy in the eastern side

short paper a note on rb sr whole rock ages from cleaved mudrocks in the welsh basin two cleaved epizone grade palaeozoic mudrock suites from snowdonia yield a mean rb sr whole rock age of 409 11 ma this is consistent with resetting of the rb sr whole rock isotope systems of the mudrocks during cleavage development in end silurian early devonian acadian deformation in contrast cleaved palaeozoic anchizone mudrocks from north and central wales give a silurian mean age of 430 9 ma this dates either a pre tectonic isotopic event preserved through acadian anchizone metamorphism or mid silurian deformation and metamorphism at anchizone grade unaffected by the acadian event the precise timing of the main caledonian deformation in wales has never been established firmly powell phillips 1985 suggest that it occurred in north and south wales in latest silurian time whereas ashgill and end llandovery tectonic events were dominant in central wales see also george 1963 however woodcock 1984 1987 and soper et al 1987 have argued in favour of a single tectonic event of early to middle devonian age this chronological study was undertaken to address this controversy the only previous attempt to define precisely the timing of cleavage development in wales was the k ar dating of cleavage parallel chlorite from a volcanic greywacke from north wales fitch et al 1969 when recalculated using modern decay constants these data give an age of 399 10 ma bath 1974 demonstrated the feasibility of obtaining rb sr ages from cleaved mudrocks of the welsh caledonides and this work produced caledonian ages of 409 20ma and 393 22ma from north wales and late

subduction zone anisotropy beneath corvallis oregon a serpentinite skid mark of trench parallel terrane migration 1 we studied the back azimuth dependence of ps converted phases at gsn station cor corvallis oregon using broadband p receiver functions computed from 602 teleseismic earthquakes the amplitudes and polarities of the transverse ps phases are largely two lobed which indicates anisotropy with a tilted symmetry axis a double peaked ps conversion at 4 5 6 5 s delay has the moveout of a dipping slab but is not consistent with simple deflection of the ps converted phase by the dipping interface a polarity flip on the transverse receiver function near north south back azimuth indicates an anisotropic symmetry axis aligned north south far from the convergence direction n68 e the ps phase is modeled using reflectivity synthetics with a highly anisotropic layer of depressed wave speed vp 6 0 6 5 km s near 40 km depth at the slab interface with the overriding north american plate one dimensional modeling suggests 10 anisotropy in the supraslab layer with a slow symmetry axis oriented n5 w at a 60 tilt from the vertical adjustments for the effect of slab interface dip on ps amplitude suggest a somewhat lower 7 anisotropy we infer a thin 7 km anisotropic detachment zone for the northward slippage of the siletz forearc terrane along the top of the descending slab a motion consistent with gps measurements and models of regional lithospheric dynamics serpentinite is a likely constituent for the deep anisotropic layer owing to its reduced vp and high poisson ratio the ductile rheology and hydrated composition of serpentinite make it a plausible lithology for a mechanical detachment zone

geological and industrial implications of extensive dilatancy anisotropy extensive dilatancy anisotropy is the hypothesized distribution of stress aligned fluid filled microcracks pervading most rocks in the earth s crust the geometry of the cracks and the aligning stress field can be monitored by analysing the waveforms of shear waves propagating through the rockmass the ability to estimate some of the parameters of the crack and stress geometry by analysing shear wave particle displacements has widespread implications

first successful 40ar 39ar dating of glauconies argon recoil in single grains of cryptocrystalline material potassium argon dating of glaucony one of the few widely occurring authigenic materials in the sedimentary record has provided a large percentage of the dates used for the calibration of the geologic time scale this material has been considered unsuitable for 40 ar 39 ar dating because of the large losses of neutron induced 39 ar during or subsequent to irradiation of samples we have measured recoiled 39 ar and 37 ar from single grains of glaucony 100 350 m in diameter irradiations were done both in air and in vacuo by using a microampoule technique of encapsulation the observed in vacuo behavior of recoiled atoms can be readily understood by a simple model of their interaction with the network of very fine laminae composing the grains when the microampoule technique is used the integrated 40 ar 39 ar ages and their associated precision compare very favorably with conventional k ar results because glauconies are inherently heterogeneous populations owing to their complex evolution it is anticipated that single grain 40 ar 39 ar dating may lead to significant refinements in the geologic time scale especially for the cenozoic

kinematic vorticity flow analysis and 40ar 39ar geochronology related to inclined extrusion of the hp lt metamorphic rocks along the zagros accretionary prism iran abstract the zagros accretionary prism in southwestern iran is exposed along the nw se trending of the zagros thrust system and inclined zagros transpression zone this accretionary prism consists of two units the upper sedimentary melange unit on top and the high pressure metamorphic melange unit at the bottom both units show characteristics of a tectonic wedge the upper unit consists of type i ii and iii melanges which display s c and c type shear band cleavages quartz ribbons and rectangular or fish head boudins the lower unit fabrics display and type porphyroclasts and quartz ribbon mylonites these fabrics formed from a combination of 60 5 simple shear and 39 5 pure shear both components were involved in a lateral exhumation of the high pressure low temperature metamorphic rocks in an inclined transpression wedge shaped geometry the estimated kinematic vorticity number wk was calculated from quartz c axis patterns rotation of porphyroclasts and orientation of finite strain with respect to shear zone boundaries using the mean estimated wm value of 0 84 the inclination angle for the thrust wedge on top of the ne subducting neo tethyan oceanic lithosphere is 18 the 40ar 39ar plateau ages of early generations of biotite from the lower metamorphic melange are 119 95 0 88 ma and 112 58 0 66 ma this late aptian age is related to early thrusting and formation of hp lt metamorphic rocks the dating of two amphibole samples from the amphibolite yields a weighted mean age of 91 23 0 89 ma this turonian cenomanian age suggests a later metamorphic event associated with subduction and obduction of the neyriz ophiolite and later lateral extrusion of hp lt metamorphic rocks along the inclined zagros accretionary prism

regional variations in the uppermost 100 km of the earth s inner core 1 the structure of the uppermost 100 km of the inner core was examined from pkikp and pkikp waveforms in the distance range of 118 140 we found evidence of a low velocity layer in the uppermost inner core in the equatorial region predominantly located between longitude 20 w to 140 e in the latitudinal direction the anomaly is detectable from 35 s beneath the indian ocean to 60 n underneath asia the maximum thickness of the low velocity layer inferred from waveform modeling is 40 km with velocity jump of about 3 we speculate that this layer may represent newly solidified core in the area where vigorous compositional convection in the outer core coincides with new crystal growth in the inner core

subduction of the philippine sea slab in view of p wave anisotropy abstract three dimensional 3d structure of isotropic and anisotropic velocities in the chugoku and shikoku districts southwest japan is retrieved from the first p arrival times from local earthquakes by assuming p wave azimuthal anisotropy with hexagonal symmetry axis the main result of isotropic velocity structure is that the philippine sea phs slab is represented by a high velocity body with the thickness of 20 30 km lying at the around depth of 40 km beneath the shikoku and southern chugoku districts the phs slab subduction appears to be interrupted by a small scale low velocity body under middle chugoku and is not clear under north chugoku these features of 3d isotropic velocity structure are basically similar to those obtained in previous tomographic studies using recent seismic data this result suggests that seismic anisotropy has no significant effect on estimation of isotropic velocity structure from the p arrival time data on the other hand fast propagation directions of p wave azimuthal anisotropy in the crust are basically in close agreement with fast polarization directions estimated from polarization anomaly of shear waves azimuthal anisotropy in the crust of shikoku district is characterized by the fast propagation direction with the e w orientation trend which is nearly parallel to trajectory of maximum principal stress acting on the southwest japan arc and strike directions of zonal geological structure crustal anisotropy in the chugoku district is consistent with the orientation of the geological structure and not with the maximum principal stress direction these results indicate that the crustal anisotropy is attributable not only to alignment of the stress induced cracks but also to lineament of geological structures we could also point out the difference in the upper mantle anisotropy between the shikoku and chugoku districts i e p wave velocity tends to be fast in the n s directions beneath shikoku and southern chugoku and in e w directions under northern chugoku since the n s axis trend of p wave anisotropy is highly correlated with high velocity anomaly of phs slab it is regarded as slab anisotropy the reason why the p wave anisotropy shows the high propagation velocity in the e w direction beneath northern chugoku is discussed on the example of two cases in which the phs slab is absent and present in the upper mantle beneath the region although we could not locate exactly the north leading edge of the phs slab in view of anisotropic velocity structure if the phs slab is defined as a high velocity body with fast propagation axes in the directions of n s system it turns out to exist at least beneath shikoku to south chugoku

laser probe 40ar 39ar and conventional k ar dating of illites associated with the mcclean unconformity related uranium deposits north saskatchewan canada the mcclean group of uranium deposits consists of elongate pods of high grade uranium mineralization width 15 40 m tightly confined to within 40 m of the basal unconformity uraninite coffinite sulphide arsenide chlorite siderite mineralization at mcclean is surrounded by a muscovite illite haematite hydrothermal alteration halo which can contain coffinite pararammelsbergite nias2 muscovite illite nodules ten laser probe 40ar 39ar dates two of which are step heat runs showing good plateaus and 10 conventional k ar dates for this material show a distribution with asharp beginning at 1320 ma a marked peak in the 1250 1200 ma class interval and a tail to dates as young as 1002 33 1 ma these determinations are in complete agreement with direct u pb and sm nd dates on uraninite at the midwest e g 1328 9 and 1110 28 ma key lake e g 1350 4 and 1281 6 ma and collins bay b deposits e g 1281 80 ma since estimated depositional ages for the athabasca sedimentary s

shear wave anisotropy of the upper and lower crusts estimated by stripping analysis of ps converted waves abstract the anisotropic structure of the crust is investigated by measuring the shear wave polarization anisotropy of the ps phases converted at the conrad and moho discontinuities the ps converted phases are identified on the p wave receiver functions which are constructed from teleseismic records at the f net and hi net stations in the kinki chugoku shikoku and kyushu districts in southwest japan the ps phase traveling through an anisotropic layer undergoes shear wave splitting which is described by splitting parameters fast polarization direction fpd and split time the upper and lower crustal anisotropies are estimated from the splitting parameters of the conrad and moho ps phases respectively which are determined by a waveform cross correlation method to reliably estimate the lower crustal anisotropy the moho ps phase is corrected for the shear wave splitting effect in the upper crust by a stripping method because the upper crustal anisotropy masks the anisotropic behavior of the ps phase in the lower crust the split times estimated for the upper and lower crusts are less than 0 2 s but their fpd characteristics differ from each other the fpds in the upper crust are predominantly in the e w direction the fpds in the lower crust range from nw se to ne sw in the region with the philippine sea phs slab subduction but from ne sw to se nw in the region without slab subduction thus seismic anisotropy exists in both the lower and upper crusts and slab subduction significantly influences the lower crustal anisotropy the most likely cause for the upper crustal anisotropy is the alignment of vertical cracks induced in the upper crust by the tectonic stress while the lower crustal anisotropy is attributed to the lattice preferred orientation of rock forming minerals along with ductile deformation of the lower crust

receiver functions of seismic waves in layered anisotropic media application to the estimate of seismic anisotropy abstract we investigate the effect of seismic anisotropy on p wave receiver functions calculating synthetic seismograms for p wave incidence on multilayered anisotropic structure with hexagonal symmetry the main characteristics of the receiver functions affected by the anisotropy are summarized as 1 appearance of seismic energy on radial and transverse receiver functions 2 systematic change of p to s ps converted waveforms on receiver functions as ray back azimuth increases and 3 reversal of the ps phase polarity on the radial receiver function in a range of the back azimuth another important influence is shear wave splitting of the ps converted waves and other later phases reverberated as s wave by numerical experiments using synthetic receiver functions we demonstrate that the waveform cross correlation analysis is applicable to splitting ps phases on receiver functions to estimate the seismic anisotropy of layer structure advantages to utilizing the ps phases are 1 they appear more clearly on receiver functions than on seismograms and 2 they inform us about what place along the seismic ray path is anisotropic real analysis of shear wave splitting is executed to the moho generated ps phases that are identified on receiver functions at six seismic stations in the chugoku district southwest japan the time lags between the two arrivals of the split ps phases are estimated at 0 2 0 7 sec and the polarization directions of the fast arrival components are from north south to northeast southwest this result is consistent with recent results of shear wave splitting measurements and the trend of linear epicenter distributions of crustal earthquakes and active fault strikes in the chugoku district

regional variations of the shear wave polarization anisotropy in the crust and mantle wedge beneath the tohoku district abstract we investigated the regional variations of the shear wave polarization anisotropy in the upper crust lower crust and mantle wedge beneath the tohoku district in the northeast japan arc the shear wave splitting parameters fast polarization direction fpd and split time for the upper and lower crusts are estimated by shear wave splitting analysis of the ps phases which are generated by the conversion of a p wave into an s wave at the conrad and moho discontinuity interfaces the splitting parameters for the mantle wedge are determined by splitting analysis of the slab ps phase converted at the oceanic crust on the subducting pacific slab and the direct s waves from deep focus earthquakes located just below seismic stations the ps phases are identified on the p wave receiver functions constructed from the teleseismic records to accurately determine the splitting parameters for the lower crust from the moho ps phase the ps phase should be corrected for the shear wave splitting effect in the upper crust similarly to estimate the mantle wedge anisotropy the direct s wave and slab ps phase must be corrected for the splitting effect in the crust the anisotropies of the crust and mantle wedge are estimated from the corrected ps phases and direct s waves the anisotropy of the upper crust exhibits a regional variation where fpd of the split shear wave is predominant in the n s direction in the pacific coast area and in the e w direction in the rest of the tohoku district the split time is less than 0 2 s the upper crustal anisotropy is attributed to the alignment of vertical cracks induced in the upper crust by tectonic stress in the lower crust fpd is predominant in the e w direction with a split time similar to that in the upper crust and the anisotropy of the lower crust is due to the lattice preferred orientation of rock forming minerals in the mantle wedge fpd is predominant in the n s direction trench parallel in the fore arc side of the volcanic front but in the e w direction trench perpendicular in the back arc side the split time in the mantle wedge is similar to that of the upper crust the most likely cause for the trench perpendicular anisotropy in the back arc mantle wedge is the lattice preferred orientation of a type olivine along flows of the secondary mantle convection for trench parallel anisotropy in the fore arc mantle wedge we considered two possible causes the preferred orientation of b type olivine along the secondary mantle convection and that of dry olivine along the trench parallel flow produced by the along strike dip variation of the pacific slab however we could not determine which one is the origin

p sh conversions in layered media with hexagonally symmetric anisotropy a cookbook reflectivity synthetic seismograms demonstrate that the type layering and orientation of 1 d anisotropy influences strongly the coda of teleseismic p waves at periods t 1 sec particularly p sh converted waves we assume the simplest form of anisotropy described by an elastic tensor with a symmetry axis of arbitrary orientation the resulting phase velocities vary as cos 2 with respect to that axis using three families of simple crustal models we compare the effects of an anisotropic surface layer with reverberations caused by both thick and thin layers of anisotropy at depth if anisotropy in the surface layer is significant the polarization of direct p can be distorted to generate a transverse component followed by ps and a prominent shear reverberation converted from direct p at the free surface if the anisotropic layer is buried the first and often the most prominent arrival on the transverse component is the p to sh conversion at its upper surface if the anisotropic layer is sufficiently thin p to sh conversions from its boundaries interfere to form a derivative pulse shape on the transverse component which could be mistaken as the signature of shear wave splitting if is horizontal compressional p and shear s anisotropy both produce similar waveform perturbations with four lobed azimuthal patterns suggesting that a weighted stack of p coda from different back azimuths would improve signal to noise for tilted between the horizontal and vertical however the effects of p and s anisotropy differ greatly the influence of p anisotropy on p to s conversion is greatest for a symmetry axis tilted at 45 to the vertical where its azimuthal pattern has two lobes rather than four combinations of p and s anisotropy typically lead to a composite azimuthal dependence in the p coda reverberations

mechanisms of argon retention in clays revealed by laser 40ar 39ar dating a method for dating clays is important for studies of weathering diagenesis hydrocarbon migration and the formation of major metalliferous deposits however many attempts have produced imprecise or inaccurate results data from shales show that contrary to expectations the 40ar 39ar dating technique can be successfully used to determine the diagenetic age of ancient sediments because 39ar losses during irradiation are controlled by release from low retentivity sites in illite equivalent to those that have lost radiogenic 40ar in nature rather than by direct recoil as is generally assumed

precise dating of low temperature deformation strain fringe analysis by 40ar 39ar laser microprobe pyritized graptolites from the welsh basin united kingdom slate belt acted as rigid bodies during cleavage formation and epizonal white micas formed within the resulting strain shadows orthogonal to the principal stress orientation although the quantities of mica are small they are a pure synkinematic mineral and have been dated by ar 40 ar 39 infrared laser microprobe as a means to dating cleavage four samples of strain fringe mica from different hand samples yielded ages ranging from 394 4 3 1 to 397 8 1 8 ma 2sigma with a mean age of 396 1 1 4 ma 2sigma by focusing on minerals that are unequivocally synkinematic this technique provides a novel solution to the problems of isotopically dating slaty cleavage previous studies have predominantly relied on dating whole rock slate samples or separated illite grains by ar 40 ar 39 techniques problems encountered included 1 separating the effects of isotopic contamination by detrital phases 2 ar 39 loss during the irradiation of illite mineral separates and 3 thermally induced ar 40 loss in nature from fine grained minerals by circumventing these problems this new method provides the first unequivocal and high precision age data for acadian deformation in the well characterized welsh basin slate belt with such precision the method may afford geologists the opportunity to track tectonic fronts across orogens and assess the rates of accretion processes in areas that are peripheral to sites of continent continent collision

coastline extraction from repeat high resolution satellite imagery abstract this paper presents a new coastline extraction method that improves water classification accuracy by benefitting from an ever increasing volume of repeated measurements from commercial satellite missions the widely used normalized difference water index ndwi method is tested on a sample of around 12 600 satellite images for statistical analysis the core of the new water classification method is the use of a water probability algorithm based on the stacking of repeat measurements which can mitigate the effects of translational offsets of images and the classification errors caused by clouds and cloud shadows by integrating quickbird worldview 2 and worldview 3 multispectral images the final data product provides a 2 m resolution coastline as well as a 2 m water probability map and a repeat count measurement map improvements on the existing coastline gshhs the global self consistent hierarchical high resolution shoreline database 50 m 5000 m in terms of resolution 2 m is substantial thanks to the combination of multiple data sources

regional variation in shear wave polarization anisotropy of the crust in southwest japan as estimated by splitting analysis of ps converted waves on receiver functions abstract regional variation in shear wave polarization anisotropy of the crust is investigated in the kyushu chugoku shikoku and kinki districts southwest japan using ps converted phases generated at the moho discontinuity the ps phases on radial and transverse receiver functions are analyzed to determine shear wave splitting parameters of fast polarization direction and time lag by waveform cross correction analysis as a result the fast polarization directions are shown to be obviously different between the chugoku district and the kyushu and shikoku districts which are separated from each other by the seto inland sea in the kyushu shikoku and south kinki districts the ps phases exhibit fast polarization directions ranging from e w to nw se the polarization directions are not only consistent with those of s waves from shallow crustal earthquakes but also parallel to the direction of maximum horizontal compression acting on southwest japan and the strike directions of metamorphic belts accretionary zones and rift zones thus the polarization anisotropy is understood to reflect mainly the upper crustal anisotropy that is attributable to stress induced vertical cracks and geological lineament structure on the other hand in the chugoku and north kinki districts fast polarization directions are in a range from n s to ne sw they are not consistent with the tectonic stress field and the s wave splitting observed for shallow crustal earthquakes the disagreement between the fast polarization directions estimated from the s waves and the ps converted phases is interpreted as arising because the lower crustal anisotropy attributed to mineral preferred orientation has a significant influence on the ps phase splitting the regional variation in ps phase polarization anisotropy of the crust in southwest japan is dependent on which anisotropy of the upper and lower crust is prominently reflected in the moho ps phase splitting

seismic velocities and anisotropy of the lower continental crust a review seismic anisotropy is often neglected in seismic studies of the earth s crust since anisotropy is a common property of many typically deep crustal rocks its potential contribution to solving questions of the deep crust is evaluated the anisotropic seismic velocities obtained from laboratory measurements can be verified by computations based on the elastic constants and on numerical data pertaining to the texture of rock forming minerals for typical lower crustal rocks the influence of layering is significantly less important than the influence of rock texture surprisingly most natural lower crustal rocks show a hexagonal type of anisotropy maximum anisotropy is observed for rocks with a high content of aligned mica it seems possible to distinguish between layered intrusives and metasediments on the basis of in situ measurements of anisotropy which can thus be used to validate different scenarios of crustal evolution

origin of crustal anisotropy shear wave splitting studies in japan shear wave splitting manifested as leading shear wave polarization that is parallel alignment of leading shear wave particle motions from a variety of sources has been observed at a number of seismograph stations in japan detected on shear wave seismograms from crustal earthquakes over a wide range of source zones and source receiver azimuths the shear wave splitting can be attributed to crustal anisotropy this paper discusses the relation between leading shear wave polarization directions and tectonic features of japan to explain the observed shear wave splitting the author proposes that at least three phenomena should be taken into account stress induced microcracks primarily aligned in vertical or subvertical planes cracks or fractures in the vicinity of active faults having their orientation parallel to the fault planes and intrinsic rock anisotropy resulting from preferred orientation of minerals travel time differences between leading and slower split shear waves from crustal and upper mantle earthquakes analyzed for about one third of the stations suggest that the crustal anisotropy which causes the observed shear wave splitting may be limited to the upper 15 25 km this implies that the density of nonhorizontally aligned cracks or fractures below 15 25 km and into the upper mantle is muchmore smaller than that in the crust above 15 25 km less

40ar39ar dating of very fine grained samples an encapsulated vial procedure to overcome the problem of 39ar recoil loss 40ar39ar dating of very fine grained materials is compromised by the loss of 39ar by recoil and also sometimes other ar isotopes during irradiation an encapsulated sample 40ar39ar procedure is described which overcomes this problem and which produces a k ar equivalent date measurements indicate that the procedure yields reliable and reproducible results it offers the potential for dating small amounts of material of many types for example very fine grained authigenic sheet silicates from rock cores which may suffer ar loss during irradiation

shear wave anisotropy in the crust mantle wedge and subducting pacific slab under northeast japan to study the anisotropic structure beneath northeast ne japan we made 4366 shear wave splitting measurements using high quality seismograms of many earthquakes occurring in the crust and the subducting pacific slab our results provide important new information on the s wave anisotropy in the upper crust lower crust mantle wedge and subducting pacific slab in the upper crust the anisotropy is mainly caused by the stress aligned fluid saturated microcracks the measured delay times dts increase to 0 10 s at 10 11 km depth the fast velocity directions fvds are parallel to either the tectonic stress or the strike of active faults the maximum dts for the low frequency earthquakes near the moho are 0 15 0 17 s suggesting strong anisotropy at the base of the crust or in the uppermost mantle the measurements for the intermediate depth earthquakes in the pacific slab show dominant e w trench normal fvds in the back arc area and n s trench parallel fvds in the fore arc area the trench normal fvds in the back arc area are caused by the corner flow in the mantle wedge as a result of the subduction of the pacific plate the maximum dts for the slab earthquakes reach 0 30 0 32 s at 100 km depth but only half of the total dts are produced in the mantle wedge the small dts in the mantle wedge may result from an isotropic or weak anisotropic zone in the middle of the mantle wedge in the fore arc the dominant trench parallel fvds for the slab earthquakes are consistent with those for the upper crust earthquakes and 80 of the total dts can be accounted for by the anisotropy in the crust in the subducting pacific slab the trench parallel fvds may reflect either the original fossil anisotropy in the pacific plate when the plate was produced in the mid ocean ridge or the preferred orientations of the crystals and cracks in the upper part of the subducting slab

seismic anisotropy of oceanic upper mantle shear wave splitting methodologies and observations we develop methodologies to obtain accurate measurements of shear wave splitting and apply these techniques to examine the pattern of oceanic upper mantle anisotropy to obtain high quality estimates of receiver splitting at island stations we devise a stacking method that finds the optimum splitting parameters and 95 error region for teleseismic shear wave phases from a suite of earthquake events to obtain additional measurements in oceanic regions we develop techniques to measure splitting parameters and errors for ss phases that sample upper mantle anisotropy at their bounce points however we find the data are often of low resolution and anomalous characteristics are sometimes found that make splitting difficult to interpret ten second splitting of ss is observed across the banjo seismic array but we cannot unambiguously attribute this signal to mantle anisotropy at the bounce the receiver splitting methods are used to assess the adequacy of a two anisotropic layer model for the pacific region with fast polarization azimuth in the lithosphere oriented in the fossil spreading direction and in the asthenosphere oriented in the absolute plate motion direction this model has been proposed to explain surface wave data in the pacific ocean but our splitting results demonstrate that oceanic anisotropy patterns are more heterogeneous than would be predicted while the island splitting measurements could reflect the influence of individual hotspot upwellings hotspot effects do not appear to be universally dominant we propose that splitting observations alternatively indicate broad scale differences in the underlying character of oceanic upper mantle anisotropy associated with coherant patters of lithospheric structure and asthenospheric flow in particular splitting surface wave models and regional studies all support a model where lithospheric anisotropy throughout the south pacific has been erased or reoriented toward the absolute plate motion direction whereas more limited observations in the north pacific indicate that the fossil lithospheric signature appears to be preserved

40 ar 39 ar dating of cleavage formation in tuffs during anchizonal metamorphism 40ar 39ar incremental release analyses were carried out on whole rock and constituent white mica illite rich size fractions 0 63 1 to 6 3 20 m within two very low grade penetratively cleaved metatuffs of contrasting anchizonal metamorphic grade northeastern rheinisches schiefergebirge federal republic of germany one sample from the upper anchizone displays internally concordant 40ar 39ar spectra with plateau ages ranging between ca 316 and 325 ma these are similar to conventional k ar ages determined for the whole rock and size fractions together the isotopic results suggest that cleavage formed at ca 320 ma during a concomitant very low grade metamorphism this is consistent with biostratigraphic controls which suggest that metamorphism and cleavage formation occurred during the westphalian

the influence of strain lithology and stratigraphical depth on white mica illite crystallinity in mudrocks from the vicinity of the corris slate belt wales implications for the timing of metamorphism in the welsh basin low and very low grade metamorphism in a mudrock dominated succession ranging from upper cambrian to wenlock is characterized by the white mica illite crystallinity technique kubler indices indicate that grade ranges from late diagenetic to low epizonal the prograde reaction most easily observed by xrd study is the progressive ordering of k na micas diagnostic 2m 1 polytypic reflections become recognizable at 0 4 2 and strengthen as grade advances grade generally correlates with thickness of overburden in the 9 km thick succession however strain exercises additionally a profound influence on grade and the distribution of strain was related to decollement resulting from inhomogeneities and ductility contrasts within the succession strain energy is seen as an additional and important variable in the development of white mica crystallinity in slates the distribution of metamorphic grade implies that the lower palaeozoic succession was augmented by several kilometres of lower devonian old red sandstone strata at the metamorphic culmination the culmination was syntectonic and represents an acadian rather than an early diastathermal event subsequent dextral movement on the bala fault displaced the regional isocryst pattern by c 4 km retrogressive smectite appears to be localized to the vicinity of the bala fault and is probably of variscan age

azimuthal shear wave velocity anisotropy in the basin and range province using moho ps converted phases an eight station large aperture seismometer array was deployed in the basin and range province of west central nevada for the 1988 1989 program for array seismic studies of the continental lithosphere passcal passive source seismic experiment during the 10 months of data collection a total of 100 teleseismic events were recorded 30 source equalized p wave receiver functions show waveform variations in timing amplitude and polarity in the moho ps phase that is diagnostic of shear wave velocity anisotropy information regarding the crustal anisotropic component can be recovered through the ps phase since it is generated at the moho and is sensitive to shear wave splitting in the crust only we utilize single event stacked and array beam formed receiver functions to minimize the influence of structural variations and scattering on the observed converted phases analysis of shear wave splitting of the ps conversion from the crust mantle boundary indicates a fast azimuth of anisotropy oriented approximately nw se with observed time delays of 0 20 s this strongly suggests coherent crustal anisotropy with a fast direction perpendicular to the strike of fault block mountain ranges this fabric is most likely due to the preferred crystallographic orientation of seismically anisotropic minerals in the middle to lower crust rather than the distribution of basin and range upper crustal fracture systems

on the detection and identification of converted and reflected phases from receiver functions summary we develop a simple methodology to extract information related to crustal structure such as the mean vp vs ratio the magnitude of lateral variations and the presence of discontinuities from a selection of receiver functions we show that crustal phases can be detected and exploited on individual receiver functions two filtering techniques based on singular value decomposition and spectral matrix filtering are proposed to enhance the coherent phases in receiver functions additionally the relative traveltimes and slownesses between the different phases are utilized for a more reliable phase identification we apply these methods to data of station pru czech republic and compare the results with a classical time domain waveform inversion we find a good qualitative agreement between the simple crustal models inferred from the two approaches waveform modelling suggests that the identification of converted phases can still be ambiguous and other geophysical or geological constraints are needed

40ar 39ar thermochronology of the kampa dome southern tibet implications for tectonic evolution of the north himalayan gneiss domes structural and thermochronological studies of the kampa dome provide constraints on timing and mechanisms of gneiss dome formation in southern tibet the core of kampa dome contains the kampa granite a cambrian orthogneiss that was deformed under high temperature sub solidus conditions during himalayan orogenesis the kampa granite is intruded by syn tectonic leucogranite dikes and sills of probable oligocene to miocene age overlying paleozoic to mesozoic metasedimentary rocks decrease in peak metamorphic grade from kyanite staurolite grade at the base of the sequence to unmetamorphosed at the top the kampa shear zone traverses the kampa granite metasediment contact and contains evidence for high temperature to lowtemperature ductile deformation and brittle faulting the shear zone is interpreted to represent an exhumed portion of the south tibetan detachment system biotite and muscovite 40 ar 39 ar thermochronology from the metasedimentary sequence yields disturbed spectra with 14 22 0 18 to 15 54 0 39 ma cooling ages and concordant spectra with 14 64 0 15 to 14 68 0 07 ma cooling ages petrographic investigations suggest disturbed samples are associated with excess argon intracrystalline deformation mineral and fluid inclusions and or chloritization that led to variations in argon systematics we conclude that the entire metasedimentary sequence cooled rapidly through mica closure temperatures at 14 6 ma the kampa granite yields the youngest biotite 40 ar 39 ar ages of 13 7 ma immediately below the granite metasediment contact we suggest that this age variation reflects either varying mica closure temperatures re heating of the kampa granite biotites above closure temperatures between 14 6 ma and 13 7 ma or juxtaposition of rocks with different thermal histories our data do not corroborate the inverse mica cooling gradient observed in adjacent north himalayan gneiss domes instead we infer that mica cooling occurred in response to exhumation and conduction related to top to north normal faulting in the overlying sequence top to south thrusting at depth and coeval surface denudation

crustal growth and late precambrian early palaeozoic plate tectonic evolution of england and wales comparison of the basement rocks of britain reveals major contrasts to the n and s of the iapetus suture solway shannon line north of the suture polymetamorphic high grade basement schists and gneisses and igneous rocks date from about 2900 ma to c 400 ma the latter age corresponding to closure of the iapetus ocean by contrast s of the suture in england and wales the basement is dominantly low grade metamorphic rock together with intrusive and extrusive igneous and sedimentary rocks there is little geological or geochemical evidence for the presence of basement older than c 900 ma below this area and the predominant expression of crustal growth in england and wales is indicated by radiometric ages between c 700 ma and 450 ma reflecting late precambrian palaeozoic magmatism and metamorphism palaeomagnetic data from scotland demonstrate that the crust in northern britain developed as an integral part of the laurentian shield during mid proterozoic times 1600 1000 ma the palaeomagnetic record from the pre ordovician rocks of england and wales shows that this area was unrelated to the laurentian shield prior to closure of the iapetus ocean the apparent polar wander a p w path is that of a distinct microplate involved in rapid growth of iapetus during early cambrian times and emplaced in its present configuration adjacent to northern britain as a result of subduction of the ocean basin between cambrian and silurian times we conclude that whereas the basement to the north of the iapetus suture formed an essentially continental unit from 2900 to 400 ma the basement of england and wales is a result of crustal growth from about 900 ma to 400 ma and formed by accretion of island arcs associated accretionary prisms and fore arc basin sediments within the iapetus ocean

shear wave polarization anisotropy in the upper mantle beneath honshu japan shear wave polarization anisotropy in the wedge portion of the upper mantle between a subducting plate and the earth s surface is investigated using three component seismograms of intermediate depth and deep earthquakes recorded at 14 local stations in honshu japan eighty nine high quality seismograms were selected from a period of 3 years the data used in this study are restricted such that incidence angles are smaller than the critical angle of 30 to the earth s surface in order to avoid phase shifts in the shear wave train to find directions of the maximum and minimum velocities in split shear waves where shear waves are resolved into two phases with the maximum time separation each set of the two horizontal component seismograms is rotated in the horizontal plane the split shear waves thus obtained are again recombined after the correction of anisotropy and the anisotropy corrected particle motion is compared with the focal mechanism for a cross check of the observed anisotropy directions of the maximum axes are plotted on azimuthincidence angle stereograms at each station the stereograms and the cross sections of seismic ray paths show that 1 the anisotropic material is distributed at intermediate location between earthquake sources and receiving stations and 2 the anisotropic region is separated into two parts one in the north of the present study area with the polarization of the maximum velocity shear wave trending 0 to 30 from the north north anisotropy and the other in the south with it trending 90 to 120 south anisotropy the maximum time delays between the two shear waves along a vertical seismic ray is about l s for both the anisotropic regions the horizontal extent of the anisotropic area in the north is 50 km at depths of 50 to 150 km perhaps prevalent in west honshu however the depth of anisotropic material is not well constrained because of insufficient coverage of seismic ray paths and angles if we assume the vertical extent of the anisotropic material to be 100 km the maximum velocity contrast would be 4 if we adopt a crack alignment model for the observed anisotropy cracks are inferred to be distributed vertically trending 0 to 30 in the north and 90 to 120 in the south if we assume an olivine alignment model olivine crystals are inferred to be aligned in the north with the a axis of the slip direction trending 0 to 30 and the slip plane normal to the b axis being vertical and in the south anisotropy with the axis trending 90 to 120 and the slip plane being horizontal

excess argon 40are uptake during slate formation a 40ar 39ar uv laserprobe study of muscovite strain fringes from the palaeozoic welsh basin uk syn kinematic pure muscovite strain fringes grew around pyritized graptolites in response to the early devonian acadian closure and uplift of the peri gondwanan welsh basin they record ar isotope information on the timing of low temperature deformation and the mudrock to slate transition here we have applied three different 40ar 39ar laserprobe dating methods to strain fringes from the central welsh basin to test for variations in age and to explore the technical limitations of the method infrared laserprobe step heating analyses yielded only occasional plateau ages varying between ca 403 ma and 422 ma and strong evidence for discordance in the ar system likewise ages from infrared laserprobe total fusion analysis are scattered and range from ca 390 ma to 651 ma by contrast ultraviolet laserprobe intra fringe age profiles revealed a complex pattern of ages across individual strain fringes this is the highest spatial resolution technique available and the resulting age profiles consistently showed an oscillatory pattern comprising troughs or baseline ages that are all within error of ca 397 ma and peaks or apparent ages ranging up to ca 550 ma these apparent ages we interpret as the results of repeated periodic assimilation of excess 40ar 40are during deformation and the mudrock to slate transition here basinal fluids carrying 40are are driven through bulk mudrock during contractional phases of the acadian orogen and 40are is taken up by muscovite strain fringes the concentration of 40are varies with metamorphic grade and is inversely proportional to the temperature of metamorphism and we suggest is controlled by the solubility of ar in water which has potentially quite severe implications for 40ar 39ar dating of k bearing minerals in any mudrock these results provide a first record of ar isotope information in strain fringes and a first record of orogenically driven fluid flow through bulk deforming mudrock they also demonstrate how the high spatial resolution ultraviolet laserprobe 40ar 39ar method can access ar isotope information that is obscured by lower spatial resolution laser analysis or bulk sampling methods

evidence from core reflected shear waves for anisotropy in the earth s mantle seismologists have begun to recognize that the upper mantle at depth is anisotropic and that seismological data are not fully explicable if anisotropy is ignored1 5 on entering an anisotropic zone a plane shear wave splits into two orthogonally polarized phases6 an effect observed by several researchers7 notably by ando et al 5 8 9 the near vertical scs wave an s wave reflected once at the core mantle boundary is one of the most suitable phases to study the phenomenon of shear wave birefringence because it is isolated from other phases and because its particle motion is nearly horizontal and thus free from s to p conversion the split scs shear waves travel with different velocities in the anisotropic zone and consequently there is a time delay between their arrival at stations my aim was to determine this time delay and the azimuth of polarization using the scs records of two horizontal components from the japanese islands for a deep shock beneath the sea of okhotsk from the total of 41 stations 34 show a polarization of the maximum velocity phase in the nnw sse direction with an average time advance of 0 8s from the minimum velocity phase the relevant anisotropic zone seems to be deep seated

topography of the seismic velocity discontinuities beneath the chugoku and shikoku districts southwest japan the first p arrival time data from 513 local earthquakes were analyzed to study lateral variation of the depth to the conrad and moho discontinuities beneath the chugoku and shikoku districts southwest japan as well as to determine earthquake hypocenters and p wave station corrections the depth to the discontinuity was estimated by minimizing the travel time residuals of more than 8700 first p arrivals observed at 55 seismic stations the conrad and moho discontinuities are located within depth ranges of 15 25 km and 30 40 km respectively the moho is deeper under the mountain area than under the seto inland sea area and especially deep under the pacific coast of the shikoku district and the mountain area in the chugoku district the depth variation of the moho is quite similar to the bouguer gravity anomaly distribution and the lateral variations of the p wave velocity the deep moho under the southern shikoku is located at the portion in which the continental moho under the island arc meets the oceanic moho that is the boundary interface between the oceanic crust and the philippine sea phs plate dipping toward the back arc although there are high mountains in the northern and middle shikoku the moho is not so deep because subduction of the phs plate prevents the moho from getting deep while the moho is deep due to isostatic balance under the mountain area in the chugoku district in addition we indicated the possibility that the upper boundary of the oceanic crust just above the high velocity phs plate is in contact with the deep moho under the western chugoku the contact of the moho with the oceanic crust can explain the markedly negative gravity anomaly observed in the western chugoku and the later phase that appears just after the first p arrival from local earthquakes

synchronous very low grade metamorphism contraction and inversion in the central part of the welsh lower palaeozoic basin a metamorphic map white mica crystallinity based on 663 data points has been generated for the recently published british geological survey 1 50 000 geological sheets 178 llanilar and 179 rhayader epizone upper anchizone lower anchizone and sub anchizone are delineated metamorphic cross sections are drawn and superimposed on structural cross sections assumptions made in order to draw metamorphic cross sections are 1 the anchizone spans 100 c 2 the metamorphic field gradient was 36 ckm 1 3 metamorphic surfaces were initially generated parallel to the contemporary ground surface restoration of the sections indicates that at the time of metamorphism a contractional phase had begun because reverse movement on earlier extensional faults development of back thrusts footwall short cuts and limited growth of precursor folds had taken place metamorphism continued locally during later stages of the contractional phase as a strain related phenomenon finally tensional faulting followed the cessation of folding and metamorphism

the evolution of illite to muscovite mineralogical and isotopic data from the glarus alps switzerland thirty five illite and muscovite concentrates were extracted from triassic and permian claystones shales slates and phyllites along a cross section from the diagenetic alpine foreland tabular jura and borehole samples beneath the molasse basin to the anchi and epimetamorphic helvetic zone of the central alps concentrates and thin sections were investigated by microscopic x ray infrared mossbauer thermal dta and tg wet chemical electron microprobe k ar rb sr 40ar 39ar and stable isotope methods with increasing metamorphic grade based on illite crystallinity data xrd and ir the following continuous changes are observed i the 1md 2m1 polymorph transformation is completed in the higher grade anchizone ii k2o increases from 6 8 wt diagenetic zone to 8 5 10 anchizone to 10 11 5 epizone reflecting an increase in the total negative layer charge from 1 2 to 2 0 iii a decrease of the chemical variation of the mica population with detrital muscovite surviving up to the anchizone epizone boundary iv a shift of an endothermic peak in differential thermal curves from 500 to 750 c v k ar and rb sr apparent ages of the fraction 2 m decrease from the diagenetic zone to the epizone k ar ages being generally lower than rb sr ages the critical temperature for total ar resetting is estimated to be 260 30 c k ar and rb sr ages become concordant when the anchizone epizone boundary is approached the stable isotope data on the other hand show no change with metamorphic grade but are dependent on stratigraphic age these results suggest that the prograde evolution from 1 md illite to 2m1 muscovite involves a continuous lattice restructuration without rupture of the tetrahedral and octahedral bonds and change of the hydroxyl radicals however this is not a recrystallization process this restructuration is completed approximately at the anchizone epizone boundary the isotopic data indicate significant diffusive loss of 40ar and 87sr prior to any observable lattice reorganization the restructuration progressively introduces a consistent repartition of ar and k in the mineral lattices and is outlined by the 40ar 39ar age spectra concordant k ar and rb sr ages of around 35 30 ma with concomitant concordant 40ar 39ar release spectra are representative for the main phase of alpine metamorphism calanda phase in the glarus alps a second age group between 25 and 20 ma can probably be attributed to movements along the glarus thrust ruchi phase while values down to 9 ma in regions with higher metamorphic conditions suggest thermal conditions persisting at least until the middle tortonian

azimuthal anisotropy in the earth from observations of sks at geoscope and nars broadband stations abstract we have analyzed observations of low frequency sks waves recorded on horizontal components at several broadband stations of the geoscope and nars arrays splitting of sks is consistently observed and in good agreement with a hypothesis of azimuthal anisotropy due to transverse anisotropy with a horizontal axis in the upper mantle beneath the receivers we have designed an inversion method which takes advantage of observations made simultaneously at different azimuths and which permits us to obtain estimates of the direction of fast axis of anisotropy as well as the time delay between the orthogonal quasi shear waves in general the directions of fast velocity obtained both for stations around the atlantic and around the pacific appear to be related to directions of past and present plate motions

split s waveforms observed in northern taiwan implications for crustal anisotropy short period stations ent and nsk in northern taiwan have frequently recorded split s waveforms from earthquakes in the underlying subduction zone to determine splitting parameters together with objective straightforward error estimations we employ a waveform cross correlation method the error for each measurement is estimated by translating the lower limit of the 95 confidence interval of the cross correlation coefficient into azimuth and time only measurements with 95 confidence regions excluding zero delay time are accepted as splitting data we also use an energy minimization method to constrain the solutions with minor splitting we interpret the splitting as caused by stress aligned cracks in the upper crust the average delay times are 0 03 0 04 s for ent and 0 07 s for nsk the fast directions at the two stations are predominantly n s which roughly match the stress trajectories predicted by a model simulating the arc continent collision in the taiwan area however at ent there is a subset of data with fast directions oriented in roughly e w the nearly 90 rotation of the inferred crack orientation occurs in a horizontal distance of 10 15 km in the upper crust both the predicted heterogeneous stress regime and the local structural setting in northeastern taiwan favor a lateral variation with a similar length scale

clay quantification and ar ar dating of synthetic and natural gouge application to the miocene sierra mazatan detachment fault sonora mexico direct dating of brittle fault rocks has been predicated on the assumption that illite in fault gouge is a mixture of two populations of clays one detrital derived from the wall rock and the other authigenic forming in the brittle fault zone during faulting owing to complex diagenetic histories of wallrock shales in previous fault dating studies this assumption has remained largely untested we demonstrate the validity of our clay quantification technique using calculated wildfire patterns to accurately model both artificial mixtures of 2m1 and 1md illite and natural mixtures in young fault gouges using our quantification approach we obtain well defined detrital and authigenic ages from the sierra mazatan metamorphic core complex in sonora mexico because the age population of the wall rock is well constrained by cooling ages we confirm that the detrital component in the gouge is the age of the footwall granite while the age of the authigenic component is the time at which faulting ceased which is in agreement with field constraints based on the ar ar age and the thermal conditions for the growth of the authigenic illite we also conclude that the sierra mazatan detachment continued slipping in the brittle regime at strain rates similar to those for the plastic regime on the order of 10 14 s 1

the distinction between caledonian burial and regional metamorphism in metapelites from north wales an analysis of isocryst patterns in snowdonia and ll n three prograde stages in metapelites gwna group to wenlock are recognized stage i with white mica crystallinity values 0 43 2 stage ii 0 43 0 26 2 stage iii isocrysts at 0 03 2 intervals are drawn separately for pelites above and below the sub arenig unconformity approximately concordant isocryst patterns characterize the unconformity in the northwestern and southeastern limbs of the snowdonia syncline whereas discordant patterns are developed northwest of the aber dinlle fault and at the closure of the ynyscynhaiarn anticline the contrasting patterns are the result of a pre arenig burial metamorphism in post tremadoc pelites belts of steep metamorphic gradient coincide with belts of high strain and a broad area of higher grade caradoc pelites coincides with the core of the snowdonia syncline the later metamorphism is related to caledonian deformation rather than to burial pelite alteration adjacent to small intermediate to acid plutons suggests emplacement before the caledonian regional metamorphism a close correlation exists between metabasite zones and metapelite stages

a transmission electron microscope study of white mica crystallite size distribution in a mudstone to slate transitional sequence north wales uk high resolution transmission electron microscopy hrtem measurements of the thickness of white mica crystallites were made on three pelite samples that represented a prograde transition from diagenetic mudstone though anchizonal slate to epizonal slate crystallite thickness measured normal to 001 increases as grade increases whereas the xrd measured 10 a peak profile the kubler index decreases the mode of the tem measured size population can be correlated with the effective crystallite size n 001 determined by xrd the results indicate that the kubler index of white mica crystallinity measures changes in the crystallite size population that result from prograde increases in the size of coherent x ray scattering domains these changes conform to the scherrer relationship between xrd peak broadening and small crystallite size lattice strain broadening is relatively unimportant and is confined to white mica populations in the diagenetic mudstone rapid increases in crystallite size occur in the anchizone coincident with cleavage development changes in the distribution of crystallite thickness with advancing grade and cleavage development are characteristic of grain growth by ostwald ripening the kubler index rapidly loses sensitivity as an indicator of metapelitic grade within the epizone

kar ages of clay size concentrates from the mineralisation of the pedroches batholith spain and evidence for mesozoic hydrothermal activity associated with the break up of pangaea abstract the k ar ages of 32 clay concentrates extracted from samples of ore gangue and wallrock associated with mineralisation in the pedroches batholith lie in the range 119 285 ma although some of the mineralisation is of permian age more than half of the ages lie between 210 and 230 ma and indicate a hydrothermal event at this triassic time a comparison with age data for mineralisation and certain anorogenic magmatism in other areas of the north atlantic reveals a consistent pattern of a major event at ca 210 230 ma and a minor event at ca 160 ma with little cretaceous or tertiary activity it is proposed that the ca 210 230 ma event was related to an initial rapid fracturing of the crust associated with the break up of pangaea which was related with an increase in the geothermal gradient and penetration of the deep crust by surface waters which returned via both new and rejuvenated fissure systems

receiver function study of the cascadia megathrust evidence for localized serpentinization 1 we characterize the crustal and upper mantle structure of the cascadian fore arc using receiver function analysis of data from long running observatories located along the north american margin in oregon and washington we identify coherent sv polarized phases consistent with p s conversion from the top of the subducting juan de fuca plate throughout the area west of the volcanic arc in many instances we also find evidence for deeper features likely associated with the subducted oceanic crust at many locations sh polarized energy in receiver functions is associated with the surface of the subducting plate these sh polarized phases display significant directional variation and show clear polarity reversals consistent with the presence of seismic anisotropy previously a layer of serpentinized material within the megathrust has been proposed to explain an anisotropic low velocity high poisson s ratio layer identified by receiver function analysis at corvallis oregon global seismic network station cor our new results suggest that this feature is widespread in cascadia but the specific receiver function signature varies this anisotropic supraslab layer is supported by receiver function estimates beneath seismic stations which similar to cor are 100 km from the coast and 40 km above the slab surface however the orientation of rock fabric implied by the anisotropy varies greatly from station to station modeling with synthetic seismograms supports the existence of serpentinite beneath station gnw green mountain washington closer to the coast where the slab is shallower an anisotropic layer is likely beneath some stations but the evidence is less conclusive we propose that the anisotropic signature beneath the inland stations is related to serpentinization of the mantle wedge at greater depths best explained by the presence of the polymorph antigorite hypothetically suctioned upward by slab rollback a deformed sliver of serpentinite rich rock above the slab interface would likely display the observed large spatial variation in anisotropic symmetry axis within the layer consistent with our observations

seismic lamination and anisotropy of the lower continental crust abstract seismic lamination in the lower crust associated with marked anisotropy has been observed at various locations three of these locations were investigated by specially designed experiments in the near vertical and in the wide angle range that is the urach and the black forrest area both belonging to the moldanubian a collapsed variscan terrane in southern germany and in the donbas basin a rift inside the east european ukrainian craton in these three cases a firm relationship between lower crust seismic lamination and anisotropy is found there are more cases of lower crustal lamination and anisotropy e g from the basin and range province western us and from central tibet not revealed by seismic wide angle measurements but by teleseismic receiver function studies with a p s conversion at the moho other cases of lamination and anisotropy are from exhumed lower crustal rocks in calabria southern italy and val sesia and val strona ivrea area northern italy we demonstrate that rocks in the lower continental crust apart from differing in composition differ from the upper mantle both in terms of seismic lamination observed in the near vertical range and in the type of anisotropy compared to upper mantle rocks exhibiting mainly orthorhombic symmetry the symmetry of the rocks constituting the lower crust is either axial or orthorhombic and basically a result of preferred crystallographic orientation of major minerals biotite muscovite hornblende we argue that the generation of seismic lamination and anisotropy in the lower crust is a consequence of the same tectonic process that is ductile deformation in a warm and low viscosity lower crust this process takes place preferably in areas of extension heterogeneous rock units are formed that are generally felsic in composition but that contain intercalations of mafic intrusions the latter have acted as heat sources and provide the necessary seismic impedance contrasts the observed seismic anisotropy is attributed to lattice preferred orientation lpo of major minerals in particular of mica and hornblende but also of olivine a transversely isotropic symmetry system such as expected for sub horizontal layering is found in only half of the field studies azimuthal anisotropy is encountered in the rest of the cases this indicates differences in the horizontal components of tectonic strain which finally give rise to differences in the evolution of the rock fabric

transition from diagenesis to metamorphism in extensional and collision settings areas in which there is a transition from diagenesis to metamorphism occur in a variety of tectonic settings some of which are of an extensional character rather than the collision setting normally associated with regional metamorphism in collision settings the diagenetic to metamorphic transition is typically followed by a continuum into the highest grades of metamorphism the baric type is usually of an intermediate pressure variety and metamorphism occurs at a late stage in the tectonic cycle in extensional settings such as marginal basins a restricted facies series is developed that ranges from diagenetic levels up to a low greenschist facies maximum the metamorphism is normally of a low pressure type and develops early in the tectonic cycle these characteristics suggest that it is possible to subdivide broadly the field of low grade regional metamorphism into compressional orogenic and extensional settings

tem and aem constraints on the origin and significance of chlorite mica stacks in slates an example from central wales u k chlorite mica stacks grains of intergrown chlorite and white mica in a matrix of fine grained white mica and chlorite have been studied using xrd sem empa tem and aem methods the stacks occur in a weakly cleaved llandoverian mudstone central wales which has a white mica illite crystallinity index of 0 35 a20 corresponding to the lower anchizone white mica occurs as packets 100 to 8 m thick interleaved with dominant chlorite packets in stacks with both apparent coherent interfaces or cross cutting small angle boundaries with chlorite layers it is well crystallized 2m polytype with phengitic composition and low paragonite component na na k 0 0 07 chlorite in stacks is fe rich and relatively homogeneous tio2 crystals surround stacks and occur within chlorite and white mica is ti rich compared to matrix white mica fine grained white mica and chlorite in the matrix have two different orientations one parallel to bedding and one parallel to cleavage which is approximately 30 50 to bedding matrix white mica is predominantly a 2m1 polytype but some cleavage parallel white mica is 3 t and some bedding parallel white mica is 1md it is na rich na na k 0 14 0 42 and relatively heterogeneous some discrete paragonite and phengitic muscovite are observed to coexist in the cleavage orientation matrix white mica and chlorite each contain less fe than corresponding white mica and chlorite in stacks both matrix and stack chlorite are one layer polytypes the data imply that chlorite in the stacks is largely derived from the replacement of volcanogenic biotite and other ferromagnesian minerals probably via intermediate expandable trioctahedral phyllosilicates most intergrowths of chlorite and mica in stacks formed by mica replacement of chlorite and altered biotite along cleavage fissures subsequent deformation caused further modification of pre existing chlorite mica stacks whereas partial dissolution of stacks and bedding parallel matrix phyllosilicates resulted in the formation of cleavage parallel phyllosilicates

three dimensional seismic velocity structure and configuration of the philippine sea slab in southwestern japan estimated by double difference tomography 1 three dimensional seismic velocity structure in and around the philippine sea plate subducting beneath southwestern sw japan is determined by applying double difference tomography method to arrival time data for earthquakes obtained by a dense nationwide seismic network in japan a region of low s wave velocity and high vp vs of several kilometers in thickness is recognized immediately above the region of intraslab seismicity in a wide area from tokai to kyushu this characteristic layer dips shallowly in the direction of slab subduction compared with the upper surface of the philippine sea slab based on seismic reflection and refraction surveys on seven survey lines we interpret that the low vs and high vp vs layer corresponds to the oceanic crust of the philippine sea slab on the basis of the position of the low vs and high vp vs layer and the precisely relocated hypocenter distribution of intraslab earthquakes the upper surface of the philippine sea slab is reliably determined for the entire area of sw japan nonvolcanic deep low frequency earthquakes that occurred associated with the subduction of the philippine sea slab are distributed along the isodepth contour of 30 km in sw japan except for the tokai district where the depth of deep low frequency earthquakes becomes gradually deeper toward northeast

40ar 39ar muscovite dating of thrust activity a case study from the axial zone of the pyrenees abstract the pyrenees belt features well exposed deformation structures that make it ideally suited to study thrust faults in this paper we present the first 40ar 39ar age of the pic de port vieux thrust ppvt located in the south western part of the pyrenean axial zone the ppvt is a secondary thrust related to the major gavarnie thrust gt with a footwall comprising upper cretaceous limestones and a hanging wall made of lower triassic pelites and sandstones a vertical transect in the fault hanging wall was investigated to characterize the deformation mechanism and to date the thrust activity by using 40ar 39ar step heating technique on muscovite petrographic observations confirm that pressure solution and recrystallization dissolution processes in presence of fluids are the main mechanisms that controlled the chemical textural changes in the fault core zone three groups of muscovite are identified in the fault zone i e detrital diagenetic and newly formed muscovites the detrital and diagenetic muscovites are mainly present in the damage zone while the newly formed muscovite is abundant in the core zone and is synkinematic to the fault activity the fractions 50 16 m 16 2 m 2 0 5 m

mapping seismic anisotropy using harmonic decomposition of receiver functions an application to northern apennines italy 1 isotropic and anisotropic seismic structures across the northern apennines italy subduction zone are imaged using a new method for the analysis of teleseismic receiver functions rfs more than 13 000 p wave coda of teleseismic records from the 2003 2007 retreating trench extension and accretion tectonics retreat experiment are used to provide new insights into a peculiar subduction zone between two continental plates that is considered a focal point of mediterranean evolution a new methodology for the analysis of receiver functions is developed which combines both migration and harmonic decomposition of the receiver function data set the migration technique follows a classical common conversion point scheme and helps to focus on a crucial depth range 20 70 km where the mantle wedge develops harmonic decomposition of a receiver function data set is a novel and less explored approach to the analysis of p to s converted phases the separation of the back azimuth harmonics is achieved through a numerical regression of the stacked radial and transverse receiver functions from which we obtain independent constraints on both isotropic and anisotropic seismic structures the application of our method to the retreat data set succeeds both in confirming previous knowledge about seismic structure in the area and in highlighting new structures beneath the northern apennines chain where previous studies failed to clearly retrieve the geometry of the subducted interfaces we present our results in closely spaced profiles across and along the northern apennines chain to highlight the convergence of the tyrrhenian and the adriatic microplates which differ in their crustal structure where the adriatic microplate subducts beneath tuscany and the tyrrhenian sea a signature of the dipping adriatic moho is clearly observed beneath the tyrrhenian moho in a large portion of the forearc region in the area where the two mohos overlap our new analysis reveals the presence of an anisotropic body above the subducted moho there is a strong ps converted phase with anisotropic characteristics from the top of the adriatic plate to a depth of at least 80 km because the ps conversion occurs much deeper than similar ps phases in cascadia and japan dehydration of oceanic crust seems unlikely as a causative factor rather the existence of this body trapped between the two interfaces supports the hypothesis of lower crustal delamination in a postsubduction tectonic setting

fault gouge dating in the southern appalachians usa illite age analysis iaa the method of comparing radiometric ages of successive size fractions with varying percentages of detrital illite has been successfully applied to several rock types including fault gouge shales and argillaceous limestones iaa results are presented for fi ve fault rocks including four clay gouges and one cataclasite from the exhumed southern appalachian foreland fold thrust belt eastern united states determining detrital versus authigenic illite is now an established procedure utilizing x ray analysis to quantify illite polytypes less established however is how to interpret 40 ar 39 ar ages as a function of diagenetic grade both a total gas age incorporating the recoiled argon fraction after irradiation and a retention age omitting the recoiled fraction are obtained for a sample we relate their respective use to diagenetic grade and specifi cally the crystallite thickness of illite when measured crystallite sizes are equal to or smaller than that calculated from ar recoil we use total gas ages here 5 nm while for crystallites greater than 10 nm we use retention ages the four clay gouge ages are all the same within error and much younger than that of the cataclasite we conclude that a major period of frontal appalachian fault activity occurred during the early permian cisuralian 276 280 ma and that these foreland faults were active simultaneously as part of an internally deforming regional thrust wedge

dating petroleum emplacement by illite 40ar 39ar laser stepwise heating the timing of petroleum emplacement is crucial in understanding petroleum forming processes and predicting new prospecting regions this study investigates the possibility of determining the age of petroleum emplacement by an illite 40ar 39ar laser stepwise heating technique a specially designed apparatus has been used to clean up the organic gases in samples contaminated with residues from oil or gas fields resulting in credible 40ar 39ar data the sandstone samples are from the zhuhai and zhujiang formations in the huizhou sag in the pearl river mouth basin south china sea illite 40ar 39ar laser stepwise heating consistently yields gradually rising age spectra at the low temperature steps until reaching age plateaus at high temperature steps the first step apparent ages of the six samples from the zhujiang formation give a weighted mean age of 11 4 0 5 ma 1 in agreement with an age of 12 1 1 1 ma in the first three steps of the sample from the zhuhai formation these dates are interpreted as being caused by the formation of fine authigenic illite and therefore indicative of the maximum age of petroleum emplacement in the huizhou sag because illite is commonly one of the latest forming minerals prior to hydrocarbon accumulation the plateau ages ranging from 58 to 276 ma in high laser output steps are interpreted as being caused by variable contributions from detrital k feldspar in the sandstones the older total ages based on the illite 40ar 39ar stepwise heating results are equal to the k ar ages and are not usable in delineating petroleum emplacement because they represent mixtures with variable contributions from authigenic illite and detrital k feldspar the 40ar 39ar stepwise heating technique thus can be applied to resolve important petroleum emplacement age information from illite samples that otherwise cannot be provided by traditional k ar methods

upb pbpb and smnd dating of authigenic monazite implications for the diagenetic evolution of the welsh basin abstract authigenic monazite nodules from silurian turbiditic rocks in the central welsh basin are zoned with respect to nd isotope composition the cores are relatively radiogenic and the rim to core trend projects toward the composition of associated very early authigenic apatite concretions these authigenic phases therefore equilibrated with changing fluid composition as they formed the apatites and the monazite nodule cores may have included significant nd from organic matter or seawater whereas the rees in the outer parts of the monazite nodules were derived from clastic detritus the apparent sm nd age obtained from these phosphates is a mixing line which gives an ordovician age for these silurian rocks the u pb isotope system does not show similar trends and shows promise as an alternative method for dating such phosphates a date of 417 11 ma 2 mswd 107 agrees well with recent estimates of the time of the diagenetic smectite illite transition in the rocks of this study a pb pb isochron gave an age of 422 24 ma 2 mswd 0 9 this shows that the pb is in isotopic equilibrium and is undisturbed and that the elevated mswd of the u pb regression is caused by slight subsequent mobility in the u pb system the event dated by the monazite is probably compactional dewatering during accelerated mid silurian sediment loading in the central welsh basin

incipient metamorphism in the lower palaeozoic marginal basin of wales a diagenctic through anchizone to epizone transition is demonstrated in pelitic rocks of the lower palaeozoic marginal basin of wales by examination of variations in phyllo silicate mineralogy illite crystallinity and bo parameter of white micas this transition represents a temperature range from 150 c to 400 c and the metamorphism is of a low pressure facies series type with a geothermal gradient of 40 ckm 1 variations in grade can be correlated largely with the original basin and shelf form suggesting a depth related metamorphism however in areas closer to the site of caledonian plate collision an increasingly syn tectonic metamorphic event is apparent correlation of pelite data with metabasite assemblages is variable the most consistent relationship being between epizone crystallinity values andepidote actinolite greenschist facies assemblages diagenetic clay mineral assemblages are found associated with prehnite pumpellyite assemblages in metabasites and it is suggested that the latter represent non buffered and therefore non diagnostic assemblages