Andreas Wüstefeld ⁎

Götz Bokelmann

Guilhem Barruol

Evidence for ancient lithospheric deformation in the East European Craton based on mantle seismic anisotropy and crustal magnetics

Keywords:

We present new shear wave splitting measurements performed at 16 stations on the East European Craton, and discuss their implications in terms of upper mantle anisotropy for this geophysically poorly known region. Previous investigations of mantle anisotropy in Central Europe have shown fast directions aligning smoothly with the craton's margin and various suggestions have been proposed to explain their origin such as asthenospheric flow or lithospheric frozen in deformation. Here, we aim at investigating the continuation of this shear wave splitting pattern further to the East, into the East European Craton. For the craton, the interpretation appears to be less ambiguous than for central Europe since several arguments support lithospheric anisotropy in this region: 1) The large scale coherence within either of the four constituting blocks and the significant variations between the blocks on a small scale, 2) the weak correlation with absolute plate motion vectors, and 3) the good correlation between anisotropy and crustal features, for which we use magnetic field alignments as a proxy. Rather good correlation of these magnetic features with seismic fast orientations strongly supports the idea of vertically coherent deformation throughout upper mantle and crust. The observed splitting orientations thus reflect the last tectonic events of each block, frozen in into the lithosphere for hundreds of millions of years.

Introduction

Deformation in the Earth's mantle can be constrained by seismic anisotropy, which is caused by deformation induced alignment of anisotropic minerals, especially olivine, into "crystal preferred orienta tions" (CPO) [START_REF] Nicolas | Formation of anisotropy in upper mantle peridotites -A review[END_REF][START_REF] Tommasi | Forward modeling of the development of seismic anisotropy in the upper mantle[END_REF][START_REF] Ben Ismaïl | An olivine fabric database: an overview of upper mantle fabrics and seismic anisotropy[END_REF][START_REF] Mainprice | The seismic anisotropy of the earths mantle: from single crystal to polycrystal[END_REF]. A change in stress regime (and thus deformation) may erase any existing CPO relatively quickly by developing new fabric or by recrystallisation processes [START_REF] Nicolas | Mechanisms of flow in naturally and experimentally deformed Peridotites[END_REF][START_REF] Mainprice | Interpretation of SKS-waves using samples from the subcontinental lithosphere[END_REF]. Therefore, any anisotropy in the lithosphere can be assumed to reflect the last period of significant deformation, or the last thermal event. In stable tectonic regimes, the anisotropy caused by the last period of deformation may remain "frozen" in the lithospheric rocks during post tectonic thermal relaxa tion [START_REF] Vauchez | Mountain building: strike parallel motion and mantle anisotropy[END_REF][START_REF] James | Tectonic implications of S-wave anisotropy beneath SE Brazil[END_REF][START_REF] Barruol | Seismic anisotropy in the eastern US: deep structure of a complex continental plate[END_REF][START_REF] Barruol | Lithospheric anisotropy beneath the Pyrenees from shear wave splitting[END_REF][START_REF] Heintz | The apparently isotropic Australian upper mantle[END_REF].

Seismic anisotropy can be detected with a broad range of seismic waves: Surface waves (e.g., [START_REF] Montagner | Dimensional structure of the Indian Ocean inferred from long period surface waves[END_REF], P waves (e.g., [START_REF] Babuška | Seismic Anisotropy in the Earth[END_REF][START_REF] Bokelmann | P-wave array polarization analysis and effective anisotropy of the brittle crust[END_REF][START_REF] Bear | Multi-wavelet analysis of three-component seismic arrays: application to measure effective anisotropy at Pinon Flats, California[END_REF][START_REF] Schulte-Pelkum | Upper mantle anisotropy from long-period P polarization[END_REF]Bokelmann, 2002a), by means of Receiver Functions [START_REF] Levin | Crustal anisotropy in the Ural Mtns foredeep from teleseismic receiver functions[END_REF] and particularly by using the splitting of shear waves [START_REF] Fukao | Evidence from core-reflected shear waves for anisotropy in the earths mantle[END_REF] for reviews see [START_REF] Silver | Seismic anisotropy beneath the continents: probing the depths of geology[END_REF][START_REF] Savage | Seismic anisotropy and mantle deformation: what have we learned from shear wave splitting[END_REF], which uses the birefringence of seismic shear waves: When passing through an anisotropic medium, a shear wave is split into two separate waves that are polarized along the anisotropic symmetry axes and travel at different speeds. That velocity difference leads to a delay between the arrival times of the two waves, which is proportional to the length of the travel path in the anisotropic layer and the strength of the anisotropy. Shear wave splitting techniques generally try to reverse the splitting effect by grid searching over all possible anisotropy orientations and delay times. The best removal of the splitting is based on different criteria, such as cross correlation of fast and slow wave [START_REF] Bowman | Shear-wave splitting in the upper-mantle wedge above the Tonga subduction zone[END_REF], eigenvalue criteria [START_REF] Silver | Implications for continental structure and evolution from seismic anisotropy[END_REF], or minimal energy on the transverse component [START_REF] Silver | Shear wave splitting and subcontinental mantle deformation[END_REF]. The method proposed by [START_REF] Chevrot | Multichannel analysis of shear wave splitting[END_REF], on the other hand, uses the backazimuthal variation of splitting strength as indicator for anisotropic orientation and strength. These methods are compared by [START_REF] Long | Estimating shear-wave splitting parameters from broadband recordings in Japan: a comparison of three methods[END_REF], [START_REF] Evans | Automated SKS splitting and upper-mantle anisotropy beneath Canadian seismic stations[END_REF], and [START_REF] Wüstefeld | Null detection in shear-wave splitting measurements[END_REF]. For backazimuthal variations of splitting and 3D modeling of fabrics see [START_REF] Sileny | Inversion of shear-wave splitting parameters to retrieve 3-D orientation of anisotropy in continental lithosphere[END_REF].

Over the past two decades the method of shear wave splitting has been widely applied in several active tectonic settings: Subduction zones (e.g., [START_REF] Margheriti | SKS splitting measurements in the Apenninic-Tyrrhenian domain (Italy) and their relation with lithospheric subduction and mantle convection[END_REF][START_REF] Levin | Detailed mapping of seismic anisotropy with local shear waves in southeastern Kamchatka[END_REF][START_REF] Nakajima | Shear-wave polarization anisotropy and subductioninduced flow in the mantle wedge of northeastern Japan[END_REF][START_REF] Plomerová | Seismic anisotropy beneath the Northern Apennines (Italy): mantle flow or lithosphere fabric? Earth Planet[END_REF], rifts [START_REF] Kendall | Teleseismic arrivals at a mid-ocean ridge: effect of mantle melt and anisotropy[END_REF][START_REF] Gao | SKS splitting beneath continental rift zones[END_REF][START_REF] Vauchez | Comment on "SKS splitting beneath rift zones[END_REF][START_REF] Walker | On the relationship between extension and anisotropy: constraints from shear wave splitting across the East African Plateau[END_REF][START_REF] Kendall | Magma-assisted rifting in Ethiopia[END_REF], hotspots [START_REF] Walker | Shear-wave splitting to test mantle deformation models around Hawaii[END_REF][START_REF] Barruol | A Tertiary asthenospheric flow beneath the southern French Massif Central indicated by upper mantle seismic anisotropy and related to the west Mediterranean extension[END_REF][START_REF] Walker | Shear-wave splitting around the Eifel hotspot: evidence for a mantle upwelling[END_REF], oceanic islands [START_REF] Behn | Detection of upper mantle flow associated with the African Superplume[END_REF][START_REF] Fontaine | Shear-wave splitting beneath the Galápagos archipelago[END_REF][START_REF] Fontaine | Upper-mantle flow beneath French Polynesia from shear wave splitting[END_REF], orogens (e.g., [START_REF] Barruol | Lithospheric anisotropy beneath the Pyrenees from shear wave splitting[END_REF][START_REF] Flesch | Constraining the extent of crust-mantle coupling in central Asia using GPS, geologic, and shear wave splitting data[END_REF].

In contrast to those active zones, cratons constitute the old, stable parts of continents. Studies of cratonic areas include North America (e.g., [START_REF] Silver | Constraints on mantle anisotropy beneath Precambrian North America from a transportable teleseismic experiment[END_REF][START_REF] Barruol | Seismic anisotropy in the eastern US: deep structure of a complex continental plate[END_REF][START_REF] Bokelmann | Mantle variation within the Canadian Shield: travel times from the APT89 portable broadband transect[END_REF][START_REF] Fouch | Shear wave splitting, continental keels, and patterns of mantle flow[END_REF][START_REF] Snyder | Two anisotropic layers in the Slave craton[END_REF], Kaapvaal craton [START_REF] Vinnik | Seismic constraints on dynamics of the mantle of the Kaapvaal craton[END_REF][START_REF] Silver | Mantle deformation beneath southern Africa[END_REF][START_REF] Gao | Mantle discontinuities beneath Southern Africa[END_REF] South American Craton [START_REF] Heintz | Shear wave splitting in SE Brazil: an effect of active or fossil upper mantle flow, or both?[END_REF][START_REF] Assumpção | Upper mantle anisotropy in SE and Central Brazil from SKS splitting: evidence of asthenospheric flow around a cratonic keel[END_REF], Australian shield (Heintz andKennett, 2005, 2006). See [START_REF] Fouch | Seismic anisotropy beneath stable continental interiors[END_REF] for a review. The information on upper mantle deformation that they supply is particularly valuable for cratonic regions, since other observational constraints on depth, such as lithospheric flexure, are lacking there. The particularity of cratonic regions is that they are generally associated with thick lithospheric roots that move coherently with the plates. The surrounding of lithospheric roots can deform more easily than the stable, thicker and colder inner parts. This raises the question of how mantle deformation due to the relative motion between surface plate and deeper mantle is accommodated. The thick lithospheric roots are also thought to act as keels that could deviate mantle flow around them [START_REF] Mckenzie | Finite deformation during fluid flow[END_REF][START_REF] Bormann | Teleseismic shear-wave splitting and deformations in Central Europe[END_REF][START_REF] Fouch | Shear wave splitting, continental keels, and patterns of mantle flow[END_REF]. Seismic anisotropy in such environments may help to distinguish present day deformation at depth, associated with plate motion [START_REF] Mckenzie | Finite deformation during fluid flow[END_REF], from fossil deformation (e.g., [START_REF] Bokelmann | Shear stress at the base of shield lithosphere[END_REF][START_REF] Vecsey | Shear wave splitting as a diagnostic of variable anisotropic structure of the upper mantle beneath central Fennoscandia[END_REF]. In this study, we focus on anisotropy beneath the East European Craton (EEC).

Fig. 1 shows shear wave splitting measurements from 28 studies available in Central Europe. Relevant references are given below in the corresponding sections. The data are taken from the interactive shear wave splitting database (Wüstefeld et al., submitted for publication; http://www.gm.univ montp2.fr/splitting/DB/), which is a freely accessible repository of teleseismic shear wave splitting publications. Researchers are invited to add their results to provide a comprehen sive, up to date database. Although shear wave splitting has been investigated around the EEC in detail, the anisotropy of the craton itself is poorly known. In Central Europe, shear wave splitting observations show orientations of the fast shear wave polarization plane Ф sub parallel to the edge of the craton. They may thus indicate an asthenospheric mantle flow deflected by the craton's edge at depth [START_REF] Bormann | Teleseismic shear-wave splitting and deformations in Central Europe[END_REF]. However, these directions are also close to the trend of the Hercynian fold belt in Central Europe, which might on the other hand suggest, that the anisotropy is frozen into the subcrustal lithosphere (e.g., Babuška et al., 1993;[START_REF] Plomerová | Seismic anisotropy and velocity variations in the mantle beneath the Saxothuringicum-Moldanubicum contact in Central Europe[END_REF][START_REF] Brechner | Backazimuthal variations of splitting parameters of teleseismic SKS phases observed at the broadband stations in Germany[END_REF][START_REF] Plenefisch | Upper mantle anisotropy at the transition zone of the Saxothuringicum and Moldanubicum in southeast Germany revealed by shear wave splitting[END_REF][START_REF] Babuška | European mantle lithosphere assembled from rigid microplates with inherited seismic anisotropy[END_REF]. The debate on the nature of the anisotropy in Central Europe is not completely resolved. Perhaps additional data will help (e.g., [START_REF] Walther | Automated analysis of SKS shear-wave splitting for regional seismic networks[END_REF]. If one could determine the sense of shear at depth from anisotropy observations, this question might be resolved (e.g. Bokelmann, 2002a). On the other hand, the depth extent of anisotropy would be indicative as well [START_REF] Montagner | Where can seismic anisotropy be detected in the earths mantle? In boundary layers[END_REF][START_REF] Mainprice | The seismic anisotropy of the earths mantle: from single crystal to polycrystal[END_REF][START_REF] Kendall | Seismic anisotropy in the boundary layers of the mantle[END_REF], but shear wave splitting gives only relatively weak constraints on depth. [START_REF] Sieminski | Finite-frequency sensitivity of body waves to anisotropy based upon adjoint methods[END_REF] showed that SKS shear wave splitting is mostly sensitive to anisotropy in the upper mantle. The depth at which splitting occurs can be further constrained by comparing observations from different azimuths, and at different stations [START_REF] Silver | The interpretation of shear-wave splitting parameters in the presence of two anisotropic layers[END_REF][START_REF] Rümpker | Apparent shear-wave splitting parameters in the presence of vertically varying anisotropy[END_REF][START_REF] Saltzer | How are vertical shear wave splitting measurements affected by variations in the orientation of azimuthal anisotropy with depth?[END_REF], although this approach assumes horizontal symmetry axes of anisotropy. We may also resort to surface waves which have better vertical resolution, but weaker lateral resolution.

Studying seismic anisotropy in cratonic regions gives an interest ing opportunity to analyze deformation processes in the mantle and to address the following geodynamics questions, amongst others:

• How did the lithosphere develop and deform?

• Is the whole lithosphere affected by the deformation, or only part of the keel (in the case of lithospheric anisotropy)? • Is there evidence for multi layer (e.g., lithospheric plus astheno spheric) anisotropy? • How is mantle flow affected by a deep and more resistant lithospheric root? • How does the topography at the base of the craton affect mantle flow? • How and to what extent do abandoned rifts affect the splitting measurements in Eastern Europe?

In recent years, there has been growing interest in studying the geodynamics of the stable portion of Europe (e.g., [START_REF] Bogdanova | Europrobe 1996 -Lithosphere dynamics[END_REF]Fig. 1. Lithospheric blocks in Eastern Europe (after Zonenshain et al., 1990). Also shown (black bars) are previously-determined shear-wave splitting observations of numerous authors, taken from the splitting database (http://www.gm.univ-montp2.fr/splitting). White arrows indicate plate motion vectors relative to a hotspot reference system HS3-Nuvel1A [START_REF] Gripp | Young tracks of hotspots and current plate velocities[END_REF]. The average velocity is 24 mm/yr. PLBT indicates the Polish-Lithuanian-Belarus terrane, whose south-eastern extent remains unclear. DDA: Dnieper-Donets Aulacogen, VOT: Volhyn-Orsha Trough. The gray triangles represent the stations included in this study. [START_REF] Grad | Crustal structure of the Trans-European suture zone region along POLONAISE97 seismic profile P4[END_REF][START_REF] Bruneton | Complex lithospheric structure under the central Baltic Shield from surface wave tomography[END_REF][START_REF] Bogdanova | EUROBRIDGE: new insight into the geodynamic evolution of the East European Craton[END_REF][START_REF] Pedersen | Lithospheric and sublithospheric anisotropy beneath the Baltic shield from surface-wave analysis[END_REF][START_REF] Vecsey | Shear wave splitting as a diagnostic of variable anisotropic structure of the upper mantle beneath central Fennoscandia[END_REF][START_REF] Artemieva | Dynamic topography of the East European craton: shedding light upon lithospheric structure, composition and mantle dynamics[END_REF]. However, only a few of the available seismological stations on the EEC have been analysed so far in terms of mantle anisotropy (e.g., [START_REF] Silver | Shear wave splitting and subcontinental mantle deformation[END_REF][START_REF] Makeyeva | Shear wave splitting and small-scale convection in the continental upper mantle[END_REF][START_REF] Helffrich | Shear-wave splitting variation over short spatial scales on continents[END_REF][START_REF] Dricker | Upper mantle flow in eastern Europe[END_REF], and the craton has not yet been discussed as a whole, motivating this study of the anisotropy and deformation of the EEC.

Geology of the East European Craton

The ancient core of Europe is composed of three major crustal segments: Fennoscandia, Sarmatia and Volgo Uralia (Fig. 1). The EEC is today mostly covered by Phanerozoic sediments. Only the north west (Fennoscandian Shield) and parts in the south (Ukrainian Shield and Voronezh Massif) display Proterozoic and Archean rocks [START_REF] Gorbatschev | Frontiers in the Baltic shield[END_REF][START_REF] Bogdanova | Riphean rifting and major Palaeoproterozoic crustal boundaries in the basement of the East European Craton: geology and geophysics[END_REF].

Between 2.1 and 2.0 Ga the oceanic domain, which separated Sarmatia from Volgo Uralia, closed. Simultaneously, subduction started in the (present day) northern edge of Sarmatia and lead to the collision with Fennoscandia at around 1.8 Ga. The accretion of the two blocks was completed at around 1.75 Ga. The central part of Fennoscandia remained stable for at least 1.2 Ga. The Sarmatian block is a stable Archean Craton, created by late Archean / early Palaeopro terozoic welding of several terranes with different ages ranging from 3.8 to 2.8 Ga [START_REF] Shcherback | The Precambrian Subdivision of the Ukrainian Shield and their Comparison with USSR and International Scales Geologicheskiy Zhurnal[END_REF]. The Archean crust of Volgo Uralia underwent apparently reorganization into a system of dome like structures [START_REF] Bogdanova | Europrobe 1996 -Lithosphere dynamics[END_REF].

The Central Russian Rift Systems is a consequence of successive rifting episodes during (1.25 0.8 Ga) the assemblage of Fennoscandia, Sarmatia and Volgo Uralia. The riftings occurred mainly along the old suture zones [START_REF] Bogdanova | Europrobe 1996 -Lithosphere dynamics[END_REF]. A suite of small failed arm rifts accompany the major rifts, namely: The late Riphean Volhyn Orsha trough divides Fennoscandia from Sarmatia and shows purely extensional features. In contrast, the Central Russian rift system is ca. 250 km wide and comprises a series of en echelon troughs. Finally, the Pachelma rift system is offset by ca 150 km to the northwest of the surface trace of the suture. [START_REF] Bogdanova | EUROBRIDGE: new insight into the geodynamic evolution of the East European Craton[END_REF] suggest that during that period Volgo Uralia acted as a wedge between Fennos candia and Sarmatia. As such Volgo Uralia moved to the present day west, while Fennoscandia escaped to the present day north.

The Trans European suture zone (TESZ) marks the contact between the Precambrian East European Craton and the Phanerozoic collage of central Europe terranes. It is a broad, complex zone resulting from the collision of continental blocks which lasted until the Variscan orogeny [START_REF] Ziegler | Geological Atlas of Western and Central Europe[END_REF][START_REF] Berthelsen | Mobile Europe[END_REF][START_REF] Pharaoh | Geological and geophysical studies in the Trans-European Suture zone[END_REF][START_REF] Grad | Crustal structure of the Trans-European suture zone region along POLONAISE97 seismic profile P4[END_REF]. In its southern portion, the TESZ roughly corresponds to an upper mantle transition from the high (S wave) velocity EEC to lower velocities beneath the younger terranes to the west [START_REF] Zielhuis | Shear-wave velocity in the upper mantle beneath central Europe[END_REF]. Furthermore, [START_REF] Babuška | Age-dependent large-scale fabric of the mantle lithosphere as derived from surface-wave velocity anisotropy[END_REF] found different anisotropic structure on either side of the TESZ.

The Ural Mountains are the topographic and tectonic feature separating Europe from Asia and therefore represent the eastern boundary of the EEC. The carboniferous (ca. 350 Ma) closure of the Ural paleo ocean accreted the East European Craton to the Siberian Craton, while the convergence continued into early Triassic (~250 Ma). The Urals extend in North South direction more than 2500 km, from the Aral Sea to the Arctic Ocean. The final subduction was directed towards the present day east [START_REF] Perez-Estaun | Uralidesa key to understanding collisional orogeny[END_REF]. The eastern flank is a mosaic of accreted terranes and obducted slabs [START_REF] Brown | The structural architecture of the footwall to the main Uralian Fault[END_REF].

As a fourth constituent of the EEC, Bogdanova (2005, and papers therein) points out the evidence for a terrane in the area NE of the TESZ and Fennoscandia, comprising NE Poland, parts of Belarus and the Baltic states. According to [START_REF] Bogdanova | The East-European craton: some aspects of the Proterozoic evolution in its south-west[END_REF] geophysical data and drillcore materials in this area identify several independent tectonic blocks between Fennoscandia and Sarmatia (Fig. 1). Some features continue through NE Poland all the way to the TESZ. The evolutionary histories of these belts allow for grouping them into three different subterranes though tectonic reconstruction and age determination remains imprecise.

Geophysical characteristics of the EEC

Tomography

Surface wave tomography provides the unique possibility to coherently study the lithospheric structure of the vast area covered by the EEC and to provide images of upper mantle structures, such as the presence of high velocity areas that may constrain the thickness of the lithosphere. Fig. 2 shows the EEC as seen in the global model of [START_REF] Debayle | Global azimuthal seismic anisotropy and the unique plate motion deformation of Australia[END_REF]. The EEC can be well identified as a region of faster than average material down to depths well below 200 km. This high velocity body extends from the Urals to the TESZ and from the Arctic Sea and Finland towards the Black Sea. Note the thinner than average lithosphere beneath Sarmatia. The other constituting blocks cannot be identified in surface wave topographic models. This is probably due to the resolution of the method and the low ray coverage in this region.

Present day motion of the EEC

Fig. 1 shows the present day plate motion according to the hot spot reference frame model HS3 [START_REF] Gripp | Young tracks of hotspots and current plate velocities[END_REF]. The previously available results from Central Europe are in good alignment with absolute motion and shear wave splitting only for a few places. Close to the boundary of the craton, fast orientations clearly deviate from plate motion direction and become tangential to the craton edge. This might be an effect of mantle flow deviated by the thick lithosphere of the EEC. On the other hand, this pattern also follows the general Hercynian structures, which alternatively might indicate lithospheric origins of anisotropy. Quantifying the motion of the plates is a key point in discussing and interpreting upper mantle anisotropy measurements since part of the signal may be acquired in the asthenosphere, due to plate drag. In that case, minerals should align parallel to the direction of the relative motion between moving plate and underlying mantle [START_REF] Nicolas | Formation of anisotropy in upper mantle peridotites -A review[END_REF]. This effect is expected to be independent of plate velocity. Slow moving plates would simply require a longer time of uniform plate motion to align the minerals. This 'simple asthenospheric flow' model [START_REF] Vinnik | Anisotropiya litosfery po nablyudeniyam voln SKS and SKKS[END_REF][START_REF] Bormann | Teleseismic shear-wave splitting and deformations in Central Europe[END_REF][START_REF] Silver | Seismic anisotropy beneath the continents: probing the depths of geology[END_REF] predicts that anisotropic orientations should be parallel to the absolute plate motion direction, assuming that mantle convection is much slower than plate motions.

Plate kinematic models describe the velocities of points at the Earth's surface due to plate tectonic motions. These absolute plate motion vectors (APM) are either derived from (a) geodetic space techniques such as VLBI, SLR and GPS (e.g. [START_REF] Drewes | Combination of VLBI, SLR and GPS determined station velocities for actual plate kinematic and crustal deformation models[END_REF][START_REF] Sella | REVEL: a model for recent plate velocities from space geodesy[END_REF], or from (b) geophysical observations such as sea floor spreading rates, transform fault and earthquake slip azimuths (e.g., [START_REF] Gripp | Current plate velocities relative to the hotspots incorporating the NUVEL-1 global plate motion model[END_REF][START_REF] Argus | No-net-rotation model of current plate velocities incorporating plate motion model NUVEL-1[END_REF][START_REF] Demets | Effect of recent revisions to the geomagnetic reversal time scale on estimates of current plate motions[END_REF][START_REF] Gripp | Young tracks of hotspots and current plate velocities[END_REF][START_REF] Kreemer | An integrated global model of present-day plate motions and plate boundary deformation[END_REF]. Hotspots have been thought to be the surface manifestation of deep mantle plumes [START_REF] Morgan | Convection plumes in the lower mantle[END_REF] but it remains a matter of debate if the hotspots can be assumed as stationary [START_REF] Steinberger | Advection of plumes in mantle flow: implications for hotspot motion, mantle viscosity and plume distribution[END_REF][START_REF] Tarduno | The Emperor Seamounts: southward motion of the Hawaiian hotspot plume in earths mantle[END_REF][START_REF] Andrews | Uncertainties in plate reconstructions relative to the hotspots, Pacific-hotspot rotations and uncertainties for the past 68 million years[END_REF]. Nevertheless, the hotspot frame remains the appropriate choice for geodynamic interpretations. The direction of HS3 plate motion vectors in Europe bend from W in the Urals to SW in the Iberian Peninsula.

All of the above mentioned models agree that the APM velocity of Europe is relatively slow. No Net Rotation (NNR) solutions give an average velocity of 25 mm/yr and 24.6 mm/yr for the REVEL2000 [START_REF] Sella | REVEL: a model for recent plate velocities from space geodesy[END_REF] and APKIM2000 [START_REF] Drewes | Combination of VLBI, SLR and GPS determined station velocities for actual plate kinematic and crustal deformation models[END_REF]) models, respectively. In the older hotspot reference frame model HS2 Nuvel1A [START_REF] Gripp | Current plate velocities relative to the hotspots incorporating the NUVEL-1 global plate motion model[END_REF], average APM is 10.6 mm/yr. This is half the velocity of the latest model, HS3 Nuvel1A, which results in 20 mm/yr. However, the directions of the NNR and hotspot reference frames result in opposite directions of plate motion for Europe. In the NNR, Europe moves generally towards NE and E, while in a hotspot frame this motion is SW to west. We will later in this study compare the splitting observations with the different plate motion models. ). Some geological features, such as the Pachelma Rift, are well represented as lineaments with short scale variations. Linear features also mark the Polish Lithuanian Belarus terrane, curving from NNE SSW in the south to E W directions in the northern part of the terrane. The Ural Mountains show patches of positive anomalies largely aligning parallel to the trend of the belt. The boundary between the EEC and Urals is sharply defined as a negative magnetic anomaly. At 58°N latitude, just north of station ARU, a lineament of alternating positive and negative anomalies is E W oriented.

Magnetic anomalies

For the shield region, it shows frequently a 'checkerboard pattern' of magnetic intensity. A zone of intense non segmented positive field can be distinguished in the central part of the Kola Peninsula, close to station LVZ. The central part of the Russian platform can be characterized by a complex structure of the magnetic field, by a combination of systems of checkerboard field with various orienta tions, and linear anomalies. Sections of the magnetic field correspond ing to the Ukrainian Shield, and Voronezh Massif, amongst others, can be distinguished within its limits. However, these limits do not always coincide with their surface expression [START_REF] Zonenshain | Geology of the USSR: A Plate-Tectonic Synthesis[END_REF]. The magnetic field of the Ukrainian shield can be characterized by great differentiation, by the presence of an anomalous zone with steep horizontal gradients, sometimes as much as 10 μT/km.

Magnetic anomalies represent structural and/or compositional differences at depths where temperatures are below the Curie Temperature. At higher temperatures, the magnetization of the material is lost. For rock magnetism, the most important mineral is magnetite (Fe 3 O 4 ) which has a Curie temperature of 578 °C. Depending on the geothermal gradient, this temperature is reached in the Earth at depth as shallow as 15 km at mid ocean ridges and 30 km for continents (e.g., [START_REF] Turcotte | Geodynamics[END_REF]. In cratonic regions, the Curie isotherm can be significantly deeper and may reach depths of 70 km (e.g., [START_REF] Artemieva | Global 1°× 1°thermal model TC1 for the continental lithosphere: implications for lithosphere secular evolution[END_REF]. Furthermore, the magnetic intensity is strongly distance dependant. Therefore, the observed magnetic anomalies represent essentially only crustal features. This renders magnetic data particularly interesting to us, since they might give constraints on crustal structure and fabric. We are particularly interested in magnetic lineaments, for developing a proxy for crustal fabric, which is otherwise difficult to obtain for the EEC, since much of surface extent is covered by sediments. Fabric of the subsedimentary basement is therefore hardly known so far.

The comparison of magnetic lineaments with anisotropic Φ can thus yield valuable information of crust mantle interaction [START_REF] Bokelmann | Comparing crustal and mantle fabric from the North American craton using magnetics and seismic anisotropy[END_REF]. In particular, parallelism of these two indepen dent datasets allows addressing the question whether deformation in crust and lithosphere is vertically coherent (VCD, cf. [START_REF] Silver | Seismic anisotropy beneath the continents: probing the depths of geology[END_REF].

Seismic anisotropy: data origin and processing

The EEC is covered by a number of permanent broad band seismographs, maintained by several networks. Figs. 1 and3 give an overview of the 16 stations analysed in this study. The stations are distributed irregularly over the area. However, they cover all major geologic regions of the EEC. The data quality is also variable, where the best signal to noise ratios (SNR) were found at stations KEV and ARU.

Unfortunately, the NARS network stations [START_REF] Muyzert | A seismic cross-section through the east European continent[END_REF] provided only few events for shear wave splitting measurements.

Most of the events used for shear wave splitting analysis have backazimuth between 45°and 100°, since they occur in the western Pacific subduction zones. Another set of events originate from the Andean subduction zone, having backazimuths of 240°to 310°(Fig. 4).

Handling of data from several data centres can easily been achieved by using the SplitLab environment [START_REF] Wüstefeld | SplitLab: a shear-wave splitting environment in Matlab[END_REF]. In total, we were able to identify 663 SKS shear wave splitting events. With SplitLab, we were able to quickly manage the entire splitting process, and it provides a comprehensive interface to simultaneously compare different splitting techniques (Fig. 5) and to test various filter time window combinations. This proved helpful since small delay times make careful data and filter selection necessary. Where possible, we used either raw data or a broad band pass filter (0.01 1 Hz). The narrowest filter we applied on an individual seismogram was a 0.02 0.2 Hz band pass. The dominant frequency of an SKS phase is usually around 0.125 Hz. Narrow filtering usually results in a Null event (i.e., no energy on the initial transverse component) and thus has a strong influence on the resulting splitting parameter estimates. Note, that care has to be taken to find a good trade off between signal fidelity and noise suppression. [START_REF] Levin | Shear wave splitting in the Appalachians and the Urals: a case for multilayered anisotropy[END_REF] tested the influence of filtering on synthetic seismograms modelled for a multilayer case. They found that the apparent splitting parameters are most variable for the lowest band pass filter applied, in their case 0.05 0.15 Hz. We remark that low frequent narrow filters (e.g., 0.01 0.15 Hz) may in general be applied for shear wave splitting measurements. These require however a very good signal to noise ratio (N10), defined as the ratio of the maximum amplitude on the radial component to the standard deviation of the transverse components, both after the removal of splitting (Restivo and Helffrich, Manual and automatic quality and Null assignment are given: g=good, f=fair, p=poor, fN=fair Null, gN = good Null.

a More complex geometry possible/beyond resolvability. b For 2 s dominant frequency.

1999). For a detailed discussion on shear wave splitting processing we refer to [START_REF] Vecsey | Shear-wave splitting measurementsproblems and solutions[END_REF].

The quality of our measurements has been determined using the automatic method defined by [START_REF] Wüstefeld | Null detection in shear-wave splitting measurements[END_REF]. This method reliably detects Null measurements and assigns a quality indicator, as reported in Table 1. An automatic detection is important to produce a homogeneous dataset for all stations by eliminating possible bias for specific stations. Such bias can be caused by a low signal to noise ratio, in which case the seismologist tends to assign a too high quality. Furthermore, subjective biases can be caused by a sequence of Null events followed by a rather poor event (or a sequence of good events followed by a fair one). A dataset with automatic Null and quality detection is therefore reproducible with fixed criteria for each assignment, and thus objective.

Results

We analysed the shear wave splitting of 16 broad band seismo graph stations on the East European Craton. Fig. 6 and Table 1 give an overview of the results (see also supplemental material). A detailed discussion) of each station can be found in Wüstefeld (2007, http:// www.gm.univ montp2.fr/splitting/). Station LVZ (Fig. 7) has also been discussed in [START_REF] Wüstefeld | Null detection in shear-wave splitting measurements[END_REF]. The observed splitting on the EEC is generally low, with delay times varying between 0.4 and 1.1 s.

Wherever backazimuthal variation could not simply be explained by a single layer of anisotropy, we have tested for a two layer case [START_REF] Silver | The interpretation of shear-wave splitting parameters in the presence of two anisotropic layers[END_REF]. Only at OBN and ARU, a two layer model explained the observations better than a single layer of horizontal anisotropy with sufficient confidence. At OBN, located close to the central cratonic triple junction of rifts (Fig. 1), the observed overall anisotropy is relatively small. The mean of the good non Null measurements is 0.6 s. Therefore, any complex anisotropic structure is hard to distinguish. Our two layer model for OBN represents the best fit with the obtained splitting parameters. More years of measurements might eventually lead to a closer insight into the anisotropic structure of these two stations in particular, and the EEC as a whole. At 70 km east of OBN, the station MHV only observations from events with easterly backazimuths are available. These show results parallel to events from the same backazimuthal range at OBN. However, the lack of backazimuthal coverage renders the interpreta tion of this more complex anisotropy difficult.

At ARU, in the central Uralian foredeep, Levin et al. (1999) identified a three layer case with the lower two layers dipping in opposite directions. A simpler two layer model explains our SKS splitting observations within similar error margins. Station AKTK, in the southern Ural Mountains (Fig. 1), has only a quarter of the number of events available for ARU (Table 1), with most results only available for events from easterly backazimuths. The proposed (one layer model) anisotropy orientation is parallel to the upper layer of ARU.

Anisotropy at station KIEV proved to be particularly difficult to detect. Despite 17 years of 3 component digital data available, only few events showed reasonably well the effect of shear wave splitting. We Fig. 6. Lower hemisphere plots of good and fair quality measurements for the stations of the EEC. Note, that we also used Null measurements (marked as circles) to constrain the seismic fast orientation Ф (cf. Fig. 6). Gridlines are in 5°steps from the centre. observed many Nulls, often from similar backazimuths as non Nulls. The automated Null detection [START_REF] Wüstefeld | Null detection in shear-wave splitting measurements[END_REF], based on characteristic differences between the cross correlation and minimum energy splitting methods, proved to be rather helpful for KIEV (Table 1). The manual selection of events showed splitting directions close to the actual backazimuth of the corresponding event.

The automated Null detection method identified many of those events as Nulls. Apparently, subjective Null detection tends to prefer non Nulls. The seismologist interprets small energy on the transverse component as evidence for a non Null event, especially in cases when the separation of fast and slow component is below resolution of splitting methods.

Discussion

On the EEC, teleseismic shear waves from a broad range of backazimuths can be recorded, notably from the Andean and South West Pacific subduction zones. Interestingly, this good coverage lead only to a relatively small selection of splitting measurements per station. This must be attributed to a low signal to noise ratio, probably due to large sediment thickness, accompanied by a generally low level of anisotropy. Furthermore, complex anisotropy patterns, as may be suggested for some stations are at the limit of resolution of the method applied here.

How can we interpret the results?

Before we can address the major questions posed in the introduction, we need to know the nature of the seismic anisotropy that is observed on the EEC. Does it arise from ongoing deformation in the asthenosphere, or does it rather represent the fossil deformation of the lithosphere? Responding to this question would be easier, if we had stronger constraints on the depth range of the anisotropy. We can nevertheless address the nature of the anisotropy by correlating with expected anisotropy orientations that would result either from asthenospheric flow, or from frozen lithospheric deformation. In the first case, anisotropy can be expected to be more or less parallel to the absolute plate motion direction. In the second case, it would be parallel to crustal fabric, at least in regions where crust and upper mantle have been deformed coherently [START_REF] Silver | Seismic anisotropy beneath the continents: probing the depths of geology[END_REF]. In the following we thus test these two hypotheses.

Correlation with plate motion

The plate velocity of the EEC, relative to the most recently proposed hotspot reference frame HS3, is about 20 mm/yr [START_REF] Gripp | Young tracks of hotspots and current plate velocities[END_REF], and thus relatively small, and it was even smaller relative to the earlier hotspot reference frame HS2, about 10 mm/yr [START_REF] Gripp | Current plate velocities relative to the hotspots incorporating the NUVEL-1 global plate motion model[END_REF]. It is however significantly different from zero and might thus in principle be able to produce a CPO. NoNetRotation frames would suggest similar plate velocities, but pointing in the opposite direction [START_REF] Drewes | Combination of VLBI, SLR and GPS determined station velocities for actual plate kinematic and crustal deformation models[END_REF][START_REF] Sella | REVEL: a model for recent plate velocities from space geodesy[END_REF]. This would produce a similar CPO.

Fig. 8 shows the angular difference at each station between the plate motion direction and the observed Φ, within 0°and 90°, which reflect parallel and perpendicular orientations, respectively. Indepen dent of the considered model, the correlation between plate motion direction and Φ is poor.

Good coherence between the observed fast axis and the plate motion vector is only observed at ARU, NE51/PUL and NE54. Since plate motion is uniform and splitting orientations change over short spatial distances, particularly in the west, this coherence seems to be only coincidental. At the other stations, plate motion vectors and Φ are almost perpendicular. The absence of coherence between splitting and plate motion indicates only a small, if any, asthenospheric contribution to shear wave splitting. Topography at the base of the lithosphere might also cause small scale variations [START_REF] Bormann | Teleseismic shear-wave splitting and deformations in Central Europe[END_REF].

Correlation with crustal fabric via magnetic lineations

Lithospheric anisotropy, on the other hand, corresponds to the frozen record of past tectonic processes. If the accompanying deformation is vertically coherent [START_REF] Silver | Seismic anisotropy beneath the continents: probing the depths of geology[END_REF], one should expect a parallelism between the fast S wave orientation and geological fabric at the surface. Studying North America, [START_REF] Bokelmann | Comparing crustal and mantle fabric from the North American craton using magnetics and seismic anisotropy[END_REF] show indeed an excellent agreement of the two quantities for most of the Canadian Shield, strongly suggesting vertically coherent deformation for that craton. For the EEC, thick sediments cover much of the surface extent [START_REF] Zonenshain | Geology of the USSR: A Plate-Tectonic Synthesis[END_REF], impeding a geological determination of crustal fabric, except for the few areas of exposed precambrian crust. For an ancient region as the EEC, topography is mostly eroded. This requires indirect observations to determine crustal structure and fabric. [START_REF] Bogdanova | Riphean rifting and major Palaeoproterozoic crustal boundaries in the basement of the East European Craton: geology and geophysics[END_REF] pointed out that the long wavelength magnetic anomalies of T MAGSAT data correspond to first order to the three major units of the EEC. Two large positive anomalies in the south coincide with Sarmatia and Volgo Uralia, respectively. In contrast, Fennoscandia is characterized by more magnetic inhomogeneity, possibly related to the complex mosaic of crustal provinces. High gradient zones coincide remarkably well with the ancient suture zones combining the three segments of the EEC.

The fact that splitting orientations are rather consistent within each tectonic unit while being different from one tectonic unit to another (Fig. 3) already represents a basic argument for a lithospheric origin of the anisotropy. We also note that a higher resolution experiment in Finland [START_REF] Vecsey | Shear wave splitting as a diagnostic of variable anisotropic structure of the upper mantle beneath central Fennoscandia[END_REF] has shown a variation across that portion of the Fennoscandia block, that agrees well with the nearly North South fast directions in the North (at stations KEV and LVZ) and more NE SW fast directions to the South (at stations PUL and NE51).

The subsurface imprints of tectonic activities in the crust are clearly visible in magnetic maps (e.g. Fig. 3) and can thus be used to constrain crustal fabric. Alternatively, gravity anomalies might be used [START_REF] Simons | Seismic and mechanical anisotropy and the past and present deformation of the Australian lithosphere[END_REF], but they offer lower resolution than magnetic data. The longer wavelength portion of the gravity field may also still contain information from the mantle, which is not desirable for this study, since we wish to compare mantle fabric (via seismic anisotropy) with crustal fabric.

An important issue is how magnetic structures are related with the tectonic regime, and more specifically, with seismic anisotropy. In transpressional regimes, folding and faulting yield magnetic (and other types of) structures parallel to the collision front. As has been shown, fast seismic S wave orientations in orogens are also oriented along the main trend of the mountain chain. We can thus expect a parallelism between magnetic lineaments and seismic fast axes.

Visual inspection of map data for (linear) features does not necessarily reveal all structures, since potential variations at smaller scales are not visible or hidden by the colormap. We prefer using an objective and reproducible method for extracting lineaments from the magnetic maps. Applying a Radon Transform to test regions around the seismic stations provides such objective measure (see Appendix A). [START_REF] Bokelmann | Comparing crustal and mantle fabric from the North American craton using magnetics and seismic anisotropy[END_REF] used that technique to identify very good correlation between magnetic and shear wave splitting data in parts of North America. The technique uses a Radon transformation, performed on circular portions of the residual magnetic map of Korhonen et al. (2007, see Fig. 3) around the stations. The test radii are 50, 80 and 100 km, respectively (cf. Fig. A1c andd). The smallest circle thus represents approximately the Fresnel zone in the uppermost mantle for shear waves with 10 s dominant period. This ensures that similar regions of crust and mantle are examined. Similar azimuths recovered from several test regions indicate the robustness of the result. Fig. A1d shows the result of a Radon Transform around station TRTE for an 80 km radius. The linear features of the magnetic anomalies are well resolved and align well with the observed fast S wave direction.

Magnetic anomalies reflect compositional and/or structural con trasts in the crust. A parallelism between splitting orientations (associated with mantle processes) and crustal magnetic lineaments thus suggests vertically coherent deformation (VCD). The results of that comparison for stations of the EEC are given in Table 2.

Fig. 9 shows the misfit between the fast splitting orientation and the magnetic lineament for all the 16 stations at different radii, grouped for each tectonic unit. In the PLTB and Fenno Scandia the fit is generally fairly good. This is striking, since fast direction vary strongly between the different stations. In fact, there is a much better correlation between seismic and magnetic directions than between the seismic directions at the various stations. The good correlation for all those stations suggests vertically coherent deformation of crust and lithosphere for Fenno Scandia and the PLTB. Interestingly, stations at the western edge of the EEC show only weak correlation (KEV, SUW). Perhaps the location of these stations over the transition from thick lithosphere to thinner and more mobile lithosphere to the west (see Fig. 2) helps to explain this anomalous behaviour [START_REF] Wylegalla | Anisotropy across the Sorgenfrei-Tornquist Zone from shear wave splitting[END_REF][START_REF] Wiejacz | Shear wave splitting across Tornquist-Teisseyre zone in Poland[END_REF].

For Sarmatia, however, there is no correlation between seismic fast axes and magnetic lineations. The correlation is nearly antiparallel for the station MHV, and oblique at KIEV and OBN. MHV and OBN are located close to the mutually perpendicular Central Russian Rift and the Pachelma Rift (Fig. 2), both of which showing clear magnetic signatures. In fact, the presence of a rift can strongly perturb the magnetic signal, since they are associated with mafic material [START_REF] Bogdanova | EUROBRIDGE: new insight into the geodynamic evolution of the East European Craton[END_REF]. This can also be seen for stations on the McKenzie dyke swarm in Northern Canada: at those stations there were two sets of magnetic lineations present, one associated with crustal fabric, and one with the dykes in the area [START_REF] Bokelmann | Comparing crustal and mantle fabric from the North American craton using magnetics and seismic anisotropy[END_REF]. The anisotropy was also difficult to determine for those stations, including KIEV. The observed splitting orientations are less well constrained, showing a complex backazimuthal pattern and many Nulls, which is especially true for KIEV. It is therefore plausible to assume that low anisotropy in Central and North Eastern Sarmatia exist and poorly constrained fast directions thus causing the poor correlation with magnetic trends.

We note also that seismic anisotropy is not yet well understood in a rift setting. While the extension of the VCD model would suggest rift normal fast directions [START_REF] Silver | Seismic anisotropy beneath the continents: probing the depths of geology[END_REF], both rift normal and rift parallel orientations have been found in rift zones (e.g. [START_REF] Gao | SKS splitting beneath continental rift zones[END_REF][START_REF] Walker | On the relationship between extension and anisotropy: constraints from shear wave splitting across the East African Plateau[END_REF][START_REF] Kendall | Mantle upwellings, melt migration and magma assisted rifting in Africa: insights from seismic anisotropy[END_REF]. The effect on the anisotropy in a stable shield caused by two superimposing rifting episodes is expected to be complex and can probably not be determined from a single station.

The still acceptable misfits of 30°and 15°for the Ural stations AKTK and ARU, respectively, might be explained by (a) the distance from the deformation front and (b) the strong magnetic anomaly north of ARU, which trends ENE WSW.

Nature of the anisotropy under the Eastern European Craton

The weak correlation of the shear wave splitting fast orientations Φ with absolute plate motion vectors of various reference systems does not favour a model of asthenospheric mantle flow beneath the East European Craton. Such a model might, however, be valid for the adjacent areas [START_REF] Bormann | Teleseismic shear-wave splitting and deformations in Central Europe[END_REF], though rejected by a growing number of studies (e.g., Babuška et al., 1993;[START_REF] Plomerová | Seismic anisotropy and velocity variations in the mantle beneath the Saxothuringicum-Moldanubicum contact in Central Europe[END_REF][START_REF] Brechner | Backazimuthal variations of splitting parameters of teleseismic SKS phases observed at the broadband stations in Germany[END_REF][START_REF] Babuška | European mantle lithosphere assembled from rigid microplates with inherited seismic anisotropy[END_REF]. Directly beneath the EEC however, our data clearly favour anisotropy frozen into the lithosphere. The relative weakness of anisotropy, (delay times between 0.4 and 1.1 s), does not allow to constrain a two layer model, which might account for asthenospheric contributions, if present. For the EEC, the coherence of anisotropy within tectonic units, and variation between units suggested that the anisotropy is primarily of lithospheric ('fossil') origin, and this was confirmed by the correlation with magnetic lineaments for Fenno Scandia and the PLTB.

An asthenospheric contribution is apparently present for large parts of the Earth. [START_REF] Gung | Global anisotropy and the thickness of continents[END_REF] interpret the presence of V SH N V SV anisotropy beneath continents at depths between 200 km and 400 km as an indicator of present day simple shear. This would be consistent with some shear wave splitting studies, which propose a two layer aniso tropy, one frozen in into the upper layer, and a deeper layer, related to ongoing processes (e.g. [START_REF] Levin | Shear wave splitting in the Appalachians and the Urals: a case for multilayered anisotropy[END_REF][START_REF] Fouch | Shear wave splitting, continental keels, and patterns of mantle flow[END_REF][START_REF] Heintz | The apparently isotropic Australian upper mantle[END_REF]. For Fennoscandia, [START_REF] Pedersen | Lithospheric and sublithospheric anisotropy beneath the Baltic shield from surface-wave analysis[END_REF] used Rayleigh and Love phase velocities to conclude that the lithosphere beneath the Baltic Shield contains little azimuthal anisotropy, but that coherent azimuthal anisotropy does exist below the lithosphere. In contrast, a joint inversion of body wave anisotropy in the same area by [START_REF] Vecsey | Shear wave splitting as a diagnostic of variable anisotropic structure of the upper mantle beneath central Fennoscandia[END_REF] interprets the existence of several anisotropic domains as fossil olivine fabric frozen in the lithosphere before assemblage. Our findings support this latter interpretation, which is also in line with [START_REF] Fouch | Seismic anisotropy beneath stable continental interiors[END_REF]. They suggest that the oceanic asthenospheric flow model might not be appropriate beneath continents. Furthermore, [START_REF] Becker | Stochastic analysis of shear-wave splitting length scales[END_REF] proposed a stochastic model of anisotropy, which shows a correlation length of splitting parameters of L ~1600 km for stable continental regions, that the authors relate to large scale tectonic processes (lithospheric anisotropy). In oceanic regions the correlation length determined from surface waves is much longer (L ~4500 km), probably reflecting asthenospheric anisotropy under oceans.

The question of the presence of an asthenosphere is especially interesting for cratons, e.g., to know whether they are firmly coupled with mantle convection or not. For the Canadian shield, a second layer which has the characteristics of 'asthenospheric flow' anisotropy was indeed required to explain the 'relative variation paradox' between P and S wave delays [START_REF] Bokelmann | Mantle variation within the Canadian Shield: travel times from the APT89 portable broadband transect[END_REF]. This was later used to study plate mantle interaction under the craton in more detail (Bokelmann, 2002a,b). Surface wave studies also suggest presence of asthenospheric anisotropy under certain cratons (e.g., [START_REF] Debayle | Global azimuthal seismic anisotropy and the unique plate motion deformation of Australia[END_REF].

The availability of information about the geology of the EEC in international literature is only growing by now and the station coverage of the vast area of the East European Craton is extremely sparse at this moment. The findings of this study should however provide a framework for planning further studies, focused on specific parts of the EEC. Suitable regions for such seismic experiment would be for instance an array over the Polish Belarus Lithuanian Terrane, the Pachelma Rift/Central Russian Rift system junction close to OBN and MHV and across the Ural Mountains. The EEC represents one of the thickest blocks of lithosphere on Earth, and it is probably a place of strong plate mantle interaction. Better understanding of the aniso tropy under Eastern Europe is therefore essential for understanding convection processes in the deep Earth (e.g., Bokelmann, 2002a), and well as their effect on the motion of plates at the Earth's surface, and their tectonic evolution.

Conclusions

We analyzed the shear wave splitting of 16 broad band seismo graph stations on the East European Craton. The observed splitting is generally low, with delay times varying between 0.4 and 1.1 s. This is in agreement with general findings in several other shield regions [START_REF] Fouch | Seismic anisotropy beneath stable continental interiors[END_REF], although some of the highest splitting delay times (N2 s) on the globe are observed on the Canadian Shield.

In order to constrain the origin of anisotropy we compared the observed fast split azimuths Φ with several datasets: the Φ are quite homogeneously oriented inside each tectonic unit but display rather strong variations between the tectonic units, suggesting "frozen in", lithospheric anisotropy inside each block. Such interpretation is supported by a lack of correlation of Φ with plate motion direction. Finally, the trend of (crustal) magnetic structures aligns relatively well with the observed splitting orientations. Of 16 analysed stations, 7 agree within 15°, 9 within 30°. This is remarkable given the sparse station distribution and varying tectonic setting. We interpret this result as vertically coherent deformation throughout the crust and the upper most mantle, supporting the idea that shear wave splitting under cratons may be dominated by the signature of the last tectonic event. It appears clear that magnetic data are useful as a proxy for crustal fabric orientation, which is especially valuable in inaccessible regions or regions with thick sediments covering the underlying basement.

There appears to be no strong indicator of present day mantle flow causing the anisotropy beneath the EEC. The spatial coherence of Φ across the TESZ is however an intriguing observation. If caused by present day mantle flow, this would require that the same mechanism applies across the TTZ. Such mechanism could involve lateral mechanical erosion at the edge of the EEC lithosphere that would be consistent with the thinner lithosphere observed beneath Sarmatia. Such interpretation is beyond the limit of resolution of surface wave tomography and seismic anisotropy studies available at this moment. For the remaining stations of the EEC, no indicator of mantle flow deviated by the keel of the EEC can be observed in splitting measure ments and the observed delay times are too small to resolve any dipping anisotropies.

Fig. 2 .

 2 Fig. 2. Tomographic surface-wave velocity model of Eastern Europe (after Debayle et al., 2005), in two depth sections, at 150 (left) and 200 km (right) depth. The velocity perturbation is given in percent relative to PREM at the corresponding depth level. Previous studies of shear wave splitting are given as white bars, triangles represent the stations included in this study.

Fig. 3

 3 Fig.3shows the residual magnetic field at an altitude of 5 km based on a recent compilation of global magnetic data(Korhonen et al.,

  2007

Fig. 4 .

 4 Fig. 4. Earthquake data distribution for splitting analysis, shown for station LVZ on the Kola Peninsula (epicentral distance between 90 and 130 degrees). The event locations are given by circles whose sizes indicate event magnitude (M W ). Shading represents the epicentral depth.

Fig. 5 .

 5 Fig.5. Example splitting analysis for an event recorded at station LVZ (output of the SplitLab code of[START_REF] Wüstefeld | SplitLab: a shear-wave splitting environment in Matlab[END_REF]. The upper panel displays a) the observed seismograms (radial/transverse components as dashed/solid lines) b) shows the splitting result on a lower hemisphere. The shaded area represents the selected time window. The center panel displays the results for the Rotation-Correlation (RC) technique: c) normalized components after rotation in RC-anisotropy system, d) Radial (Q) and transverse (T) seismogram components after RC-correction, e) surface particle motion in geographical coordinates before (dashed) and after (solid) RC correction and f) map of correlation. Lower panel displays the results for the minimum energy (SC) technique: g) normalized components after rotation in SC-anisotropy system, h) SC-corrected radial and transverse seismogram component, i) SC particle motion before and after correction and j) map of minimum energy on transverse component.

Fig. 8 .

 8 Fig. 8. a) Angular difference between the seismic fast orientation Ф and the absolute plate motion direction, as calculated from various models. The splitting directions of the 16 analyzed stations on the EEC show no correlation with any of the 4 plate motion models. b) Histogram of the angular difference between the seismic fast orientation Ф and magnetic preferred orientations for three search radii.

Fig. 7 .

 7 Fig. 7. Shear-wave splitting estimates from 41 measurements at station LVZ for the Minimum Energy method. The left panel displays the fast orientation, the right panel the delay time estimates. The solid horizontal lines indicates our interpretation of the LVZ with fast axis at 15°and 1.1 sec delay time, based on the mean of the good splitting measurements, shown by + signs. Note that the placement of the Nulls agrees with that interpretation.

Fig. 9 .

 9 Fig. 9. Misfit between magnetic lineaments and seismic fast orientation for different test radii (50, 80, and 100 km) around each station, grouped by tectonic unit. In the PLTB and Fenno-Scandia there is generally good agreement, except for stations on the western edge of the craton (KEV, SUW). There is no correlation for Sarmatia, the reasons for which are discussed in the text.

Table 1

 1 Splitting results of the stations on the East European Craton.

	Station	Lat	Long	Network	Total	Automatic / manual	Φ		dt
							Upper	Lower	Upper	Lower
	AKTK	50.435	58.02	Kazachstan	35	6g + 1f + 10p + 2fN + 16gN	56°0.5
						3g + 6f + 6p + 11fN + 9gN
	ARU	56.43	58.56	GSN	139	34g + 19f + 35p + 20fN + 31gN	55°78°0.5	0.4
						49g + 32f + 11p + 26fN + 21gN
	KEV	69.76	27.01	GSN	44	4g + 3f + 12p + 9fN + 16gN	0°0.4
						3g + 4f + 0p + 19fN + 18gN
	KIEV a	50.69	29.21	GSN	61	3g + 4f + 24p + 3fN + 27gN	75°1.1
						19g + 14f + 5p + 13fN + 10gN
	LVZ	67.9	34.65	GSN	41	10g + 8f + 6p + 7fN + 10gN	10°1.1
						16g + 7f + 1p + 9fN + 8gN
	MHV	54.96	37.77	GEOFON	61	4g + 2f + 10p + 2fN + 13gN	35°1.1
						5g + 6f + 3p + 8fN + 9gN
	NE51	59.881	29.826	NARS	12	4g + 0f + 3p + 2fN + 3gN	47°0.7
						5g + 4f + 0p + 0fN + 3gN
	NE52 a	57.819	28.39	NARS	8	1g + 0f + 4p + 0fN + 3gN	60°a		0.8
						1g + 0f + 0p + 5fN + 2gN
	NE53	54.904	26.793	NARS	18	0g + 1f + 5p + 3fN + 9gN	NS		b 0.5
						2g + 1f + 3p + 6fN + 6gN
	NE54	52.568	23.861	NARS	3	1g + 0f + 1p + 0fN + 1gN	60°0.6
						1g + 1f + 0p + 1fN + 0gN
	NE55	49.716	29.656	NARS	6	0g + 1f + 2p + 0fN + 3gN	37°0.8
						1g + 0f + 0p + 1fN + 4gN
	NE56	46.676	30.899	NARS	6	4g + 0f + 1p + 0fN + 1gN	45°0.6
						3g + 1f + 0p + 0 fN + 2gN
	OBN b	55.11	36.57	GSN	115	12g + 16f + 42p + 8fN + 37gN	30°70°0.8	0.5
						26g + 22f + 13p + 24fN + 30gN
	PUL	59.77	30.32	GEOFON	20	2g + 1f + 5p + 0fN + 12gN	41°0.7
						4g + 4f + 0p + 1fN + 11gN
	SUW	54.012	23.18	GEOFON	72	13g + 11f + 20p + 8fN + 20gN	70°0.8
						19g + 21f + 2p + 15fN + 15gN
	TRTE	58.38	26.72	GEOFON	22	1g + 5f + 7p + 3fN + 6gN	35°0.7
						6g + 7f + 0p + 5fN + 4gN

Table 2

 2 Difference (in degrees) between upper-layer seismic fast-axis orientation Φ and dominant magnetic trend Ψ for each stations of the EEC for different radii.

	Station Fast orientation Magnetic Orientation Ψ	Misfit
	Ф		50 km 80 km 100 km 50 km 80 km 100 km
	AKTK	56	24	22	24	32	34	32
	ARU	55	69	71	71	14	16	16
	K E V	0	4 7	6 8	7 7	4 7	6 8	7 7
	KIEV	75	35	28	29	40	47	46
	LVZ	10	59	30	25	69	20	15
	MHV	35	51	53	74	86	88	71
	NE51	47	67	57	49	20	10	2
	NE52	60	66	59	58	6	1	2
	NE53	0	4	1	9	4	1	9
	NE54	60	10	48	46	50	12	14
	NE55	37	40	30	18	3	7	19
	NE56	45	21	11	11	24	34	34
	OBN	30	66	89	70	84	59	40
	PUL	41	60	49	48	19	8	7
	SUW	70	29	22	26	81	88	84
	TRTE	35	44	44	46	9	9	11

Acknowledgements

We would like to thank the IRIS, NARS, and Geofon networks for providing data, E. Debayle for the tomography model and A. Vauchez for fruitful discussions. This work was performed using funding of the European Community (contract MIRG CT 2004 003953). The sugges tions of two anonymous reviewers helped to improve the manuscript. The splitting results of the 16 stations on the EEC analysed in this study can be found in the supplemental material. A detailed discussion of each station is provided in Wüstefeld (2007), available at http:// www.gm.univ montp2.fr/splitting.

Appendix A

An objective method to identify linear features of a complex map can be based on the Radon Transform, as has been done in other recent approaches [START_REF] Zhang | Automatic lineament extraction from potentialfield images using the Radon transform and gradient calculation[END_REF][START_REF] Hansen | Linear feature analysis for aeromagnetic data[END_REF][START_REF] Bokelmann | Comparing crustal and mantle fabric from the North American craton using magnetics and seismic anisotropy[END_REF]:

The transform corresponds to a slant stack along a line in direction θ of the magnetic field f(x,y). Linear features in f(x,y) are projected into a point in R[f] at angle θ and offset s, whereas points in f(x,y) are represented as sinusoids in R[f]. Fig. A1a presents a synthetic dataset with two linear features with 30°and 60°azimuth, with positive and negative amplitude, respectively. Also shown are two point features of positive and negative amplitude. Fig. A1b shows the Radon Transform of this test. Note that the azimuths are well recovered at 30°and 60°as well as 210°and 240°. The two point anomalies in a) are shown as sinoids in b), whose maximum amplitude represents their distances from the centre of the map.

The Radon Transform therefore gives a basis for either detecting individual lineaments in a map, or for characterizing preferred orientations by summing over all offsets s in an appropriate way. Here, we are primarily interested in the overall preferred orientations rather individual lineaments, both of positive and negative amplitude. This "polarization" of the magnetic map can be identified by the measure

The exponent λ may be chosen larger than 1 to enhance the effect. In this study, we have chosen λ = 2.

We perform this transformation on residual magnetic data, based on aeromagnetic surveys and ship cruises relative to satellite magne tometer measurements and geomagnetic observatories (Fig. 3, [START_REF] Korhonen | Magnetic Anomaly Map of the World[END_REF]. This compilation has a resolution of 5 km and we downsampled the data to 1 km. The orientations are calculated for several circular test regions around each station with radii of 50, 80 and 100 km, respectively (cf. Fig. A1c andd). Similar azimuths recovered from several test regions indicate the robustness of the result. Fig. A1d shows the results of a Radon Transform around station TRTE for an 80 km radius. The linear features of the magnetic anomalies are well resolved and align well with the fast S wave direction observed at that station (see main text), thus confirming that the method is applicable to the real magnetic dataset.

Appendix A. Supplementary data

Supplementary data associated with this article can be found, in the online version, at doi:10.1016/j.tecto.2009.01.010.