Crustal motion and deformation in Greece from a decade of GPS measurements, 1993–2003

The Hellenic plate boundary region, located in the collision zone between the Nubian/Arabian and Eurasian lithospheric plates, is one of the seismo-tectonically most active areas of Europe. During the last 15 years, GPS measurements have been used to determine the crustal motion in the area of Greece with the aim to better understand the geodynamical processes of this region. An extended reoccupation network covering whole Greece has been measured periodically in numerous GPS campaigns since the late eighties, and a continuous GPS network has been operated in the region of the Ionian Sea since 1995. In this paper, we present a new detailed high-quality solution of continuous and campaign-type measurements acquired between 1993 and 2003. During the GPS processing, a special effort was made to obtain consistent results with highest possible accuracies and reliabilities. Data of 54 mainly European IGS and EUREF sites were included in the GPS processing in order to obtain results which are internally consistent with the European kinematic field and order to allow for a regional interpretation. After an overview of the results of the IGS/EUREF sites, the results from more than 80 stations in Greece are presented in terms of velocities, time series, trajectories and strain rates. Previous geodetic, geological and seismological findings are generally confirmed and substantially refined. New important results include the observation of deformation zones to the north and to the south of the North Aegean Trough and in the West Hellenic arc region, arc-parallel extension of about 19 mm/yr along the Hellenic arc, and compression between the Ionian islands and the Greek mainland. Due to continuous long-term observations of 4–8 years, it was possible to extract height changes from the GPS time series. In Greece, we observe a differential subsidence of the order of 2 mm/yr between the northern and central Ionian islands across the Kefalonia fault zone. The differential subsidence of the central Ionian islands with respect to the northwestern Greek mainland amounts to 4 mm/yr.     

Active crustal deformation in two seismogenic zones of the Pannonian region — GPS versus seismological observations

The seismicity and the associated seismic hazard in the central part of the Pannonian region is moderate, however the vulnerability is high, as three capital cities are located near the most active seismic zones. In our analysis two seismically active areas, the Central Pannonian and Mur-Mürz zones, have been considered in order to assess the style and rate of crustal deformation using Global Positioning System (GPS) and earthquake data.We processed data of continuous and campaign GPS measurements obtained during the years 1991–2007. Velocities relative to the stable Eurasia have been computed at HGRN, CEGRN and EPN GPS sites in and around the Pannonian basin. Uniform strain rates and relative displacements were calculated for the investigated regions. GPS data confirm the mostly left lateral strike slip character of the Mur-Mürz–Vienna basin fault system and suggest a contraction between the eastward moving Alpine-North Pannonian unit and the Carpathians.The computation of the seismic strain rate was based on the Kostrov summation. The averaged unit norm seismic moment tensor, which describes the characteristic style of deformation, has been obtained from the available focal mechanism solutions, whereas the annual seismic moment release showing the rate of the deformation was estimated using the catalogues of historical and recent earthquakes.Our analysis reveals that in the Central Pannonian zone the geodetic strain rate is significantly larger than the seismic strain rate. Based on the weakness of the lithosphere, the stress magnitudes and the regional features of seismicity, we suggest that the low value of the seismic/geodetic strain rate ratio can be attributed to the aseismic release of the prevailing compressive stress and not to an overdue major earthquake. In the Mur-Mürz zone, although the uncertainty of the seismic/geodetic strain rate ratio is high, the seismic part of the deformation seems to be notably larger than in the case of the Central Pannonian zone. These results reflect the different deformation mechanism, rheology and tectonic style of the investigated zones.     

Change in seismic activity in the Tokai region related to weakening and strengthening of the interplate coupling

We investigate the relationship between changes in seismicity and crustal deformations in the Tokai region. We describe how seismicity in the subducted slab increased remarkably in the fall of 2000 and decreased in the fall of 2001, while in contrast, the crust seismicity decreased in the fall of 2000 and increased in the fall of 2001. We note that the trend of horizontal displacement at GPS stations changed coincidentally and we propose interpreting the increase and decrease in seismic activities and the changes in crustal deformations in a unified way based on changes in the state of the interplate coupling, i.e., the back-slip rate was reduced in the fall of 2000 and was partially restored in the fall of 2001. We explain why reduction of the back-slip rate increases seismogenic stress in the slab and decreases stress in the crust. We also describe the substantial positive dilatation observed in the region around Mt. Fuji in the fall of 2000 and suggest that the remarkable increase of low-frequency earthquakes beneath Mt. Fuji in October 2000 may have been caused by deceleration of the converging motion of the Izu micro-plate with the Eurasian plate. The decrease of the subduction velocity of the Izu micro-plate on the Suruga Trough in late 2000 would also have contributed to weakening of the interplate coupling beneath the Tokai region, since reduction of the relative velocity between overriding and subducting plates produces the same effect on the plate interface as a diminishing back-slip rate. However, subduction of the Izu micro-plate on the Suruga Trough was accelerated in early 2003, which may have caused increases in both slab and crust seismicities in that period.     

Crustal deformation in the Central Ionian Islands (Greece): Results from DGPS and DInSAR analyses (1995–2006)

Ground deformation studies based on Differential GPS (DGPS) and Differential Interferometric SAR (DInSAR) analyses have been conducted in the seismically active area of the Central Ionian Islands. Local GPS networks were installed in Cephallonia (2001) and Zakynthos (2005). The Cephallonian network has been remeasured five times and Zakynthos' once as of July 2006. The studies have yielded detailed information regarding both local and regional deformations that are occurring in the area.For Lefkas Island, DInSAR analysis (March to September 2003) revealed 56 mm of uplift in the central and western parts and is attributed to the August 2003 earthquake (M_w = 6.3) that occurred offshore to the west. Synthetic DInSAR modelling of the magnitude and extent of deformation is consistent with the seismologically deduced parameters for the ruptured segment along the Lefkas Transform Fault. Subsidence (< 28 mm) along the northern part of the island is attributed to local conditions unrelated to the earthquake. For Zakynthos Island, large-magnitude earthquakes that occurred offshore to the south in October 2005 and April 2006 most likely contributed to the observed deformation as deduced from DGPS measurements for an encompassing period (August 2005 to July 2006). The largest amount of horizontal deformation occurred in the south, where its western part moved in a W–NW direction, while the eastern part moved towards the NE, with magnitudes ranging from 15 to 26 mm. The southern part of the island uplifted a maximum of 65 mm whereas the north subsided from 12 to 28 mm.For Cephallonia Island, DInSAR analysis (1995 to 1998) indicated ground deformation up to 28 mm located in small sections of the island. Further interferometric analysis for the period 2003 to 2004, encompassing the occurrence of the Lefkas earthquake in August 2003, indicated 28 mm of uplift in the northern part, while during the next two years (2004 to 2005), further uplift of at least 56 mm had taken place at the western and northern part of the island.DGPS measurements for the period 2001 to 2006 revealed a clockwise rotation of the island with respect to a centrally located station on Aenos Mt. The horizontal component of deformation generally ranged from 6 to 34 mm, with the largest values at the western and northern parts of the island. Considering the vertical deformation, two periods are distinguished. The first one (2001 to 2003) is consistent with anticipated motions associated with the main geological and tectonic features of the island. The second one (2003 to 2006) has been tentatively attributed to dilatancy in which relatively small uplift (20–40 mm) occurred along the southern and southeastern parts of the island, while larger magnitudes (> 50 mm) happened at the western part (Paliki Peninsula). These large magnitudes of uplift over an extended area (> 50 km), in conjunction with an accelerated Benioff strain determined from the analysis of the seismicity in the broader region, are consistent with dilatancy. This effect commenced some time after 2003 and is probably centered in the area between Zakynthos and Cephallonia. If this interpretation is correct, it may foreshadow the occurrence of a very strong earthquake(s) sometime during 2007 to 2008 in the above designated region.     

Pressure–temperature conditions of ongoing regional metamorphism beneath the Japanese Islands

We evaluate the pressure–temperature (P–T) conditions of ongoing regional metamorphism at the top of the oceanic crust of the subducted Pacific and Philippine Sea plates through a combination of phase diagrams and hypocenter distribution and based on the dehydration-induced earthquake hypothesis. The brute-force method was employed to find the best match thermal structure to link the hypocenter distribution and dehydration. The estimated thermal structure varies far from the values obtained from numerical simulation. Our estimates are consistent with the qualitative physical prediction for the variation of temperature in different subduction zones and provide quantitative constraints for the models.In northeastern Japan, the P–T path for the top of the oceanic crust turns to the high-T side at a depth of around 90 km. The depth corresponds to the location of the volcanic front and an active convection of the wedge mantle below this depth is suggested. Our computations also reveal the effect of an exceptional scenario beneath the Kanto region. The temperature in the Kanto region, where the cold lid of the Philippine Sea plate prevents heating by the return-flow of mantle wedge above, is much lower than that of northeastern Japan. The subduction of younger Philippine Sea plate leads to a higher-temperature in the oceanic crust. In the central Shikoku region, the thermal structure exhibits high-T/P nature. Heating by shear deformation can explain the high-T/P path in the depth range from 20 to 35 km. The Kyushu area shows moderate type T/P path reaching up to eclogite facies conditions. In the Kii and central Shikoku region, the thermal structure exhibits high-T/P nature. However, the absolute values for the areas seem to have problem in physical context. Our results has risen the significance of sediment subduction in the southwest Japan and requirement for further improvements in this technique including the aspect of variation of the bulk composition of the subducted material.     

Space–time distribution of interplate moment release including slow earthquakes and the seismo-geodetic coupling in the Sanriku-oki region along the Japan trench

Along the Japan trench where some Mw8 class interplate earthquakes occurred in the past century such as the 1896 Sanriku tsunami earthquake (M6.8, Mt8.6, 12×10²⁰ N m) and the 1968 Tokachi-oki earthquake (Mw8.2, 28×10²⁰ N m), the Pacific plate is subducting under northeast Japan at a rate of around 8 cm/year. The seismic coupling coefficient in this region has been estimated to be 20–40%. In the past decade, three ultra-slow earthquakes have occurred in the Sanriku-oki region (39°N–42°N): the 1989 Sanriku-oki (Mw7.4), the 1992 Sanriku-oki (Mw6.9), and the 1994 Sanriku-oki (Mw7.7) earthquakes. Integrating their interplate moments released both seismically and aseismically, we have the following conclusions. (1) The sum of the seismic moments of the three ultra-slow earthquakes was (4.8–6.6)×10²⁰ N m, which was 20–35% of the accumulated moment (18.6–23.0)×10²⁰ N m, in the region (39°N–40.6°N, 142°E–144°E) for the 21–26 years since the 1968 Mw8.2 Tokachi-oki earthquake. This is consistent with the previous estimates of the seismic coupling coefficient of 20–40%. On the other hand, the sum of the interplate moments including aseismic faulting is (11–16)×10²⁰ N m, leading to a “seismo-geodetic coupling coefficient” of 50–85%, which is an extension of the seismic coupling coefficient to include slow events. (2) The time constants showed a large range from 1 min (10² s) for the 1968 Tokachi-oki earthquake to 10–20 min (10³ s) for the 1896 Sanriku tsunami earthquake, to one day (10⁵ s) for the 1992 Sanriku-oki ultra-slow earthquake, to on the order of one year (10⁷ s) for the 1994 Sanriku-oki ultra-slow earthquakes. (3) Based on the space–time distribution, three “gaps of moment release,” (40.6°N–42°N, 142°E–144°E) 39°N–40°N, 142°E–143°E) and (39°N–40°N, 142°E–144°E), are identified, instead of the gaps of seismicity.     

Local seismic array observations at north Evoikos, central Greece, delineate crustal deformation between the North Aegean Trough and Corinthiakos Rift

In order to better constrain and define the microseismic activity at the north Evoikos Gulf and its surrounding area we deployed an onshore/offshore seismic array consisting of 31 three-component seismic digital stations. The array was active from 30 June to 24 October 2003, and covered an area of 2500 km². We located more than 2000 seismic events ranging from 0.7 to 4.5 M_L by using six stations as a minimum in order to define the foci parameters. Recorded seismicity delineated three major zones of deformation: from south to north, the Eretria–Parnis–eastern Corinthiakos zone, the Psachna–Viotia zone, and the Northern Sporades–North Evia–Bralos zone. Alignments of the recorded seismicity follow the tectonic trends and their orientation in the above zones. The whole area accommodates the stress field between the North Aegean Trough and the Corinthiakos Gulf. Rate of deformation intensifies from north to south, as revealed also by historical and instrumental seismicity. The successive change of orientation between the two stress fields fragments the crust in relatively small units and the fault systems developed do not permit the generation of major earthquakes in the north Evoikos area and its immediate vicinity. This is also supported by the instrumental seismicity of the last century. Larger events reported in historical times are probably overestimated.Most seismic activity is crustal. Subcrustal events were recorded mainly below the Lichades area and are interpreted as the consequence of the subduction of the Ionian oceanic lithosphere below the Hellenides. The Lichades volcano is the most northern end of the Hellenic volcanic arc.At present the highest seismic activity is associated with the Psachna region of north Evia that has been continuously active since 2001. Considering, however, the development of the seismic activity during the last decade, there has been a sequence of large events, i.e., Parnis in 1999, Skyros in 2001 and Psachna in 2001–2003. This demonstrates the fact that the tectonic deformation in all this area is intense and important for the accommodation of the stress field of the North Aegean Trough to that of the Corinthiakos Rift.     

Robust and exploratory analysis of active mesoscale tectonic zones in Japan utilizing the nationwide GPS array

A monitoring GPS array recently developed in Japan can yield nationwide maps of active inland tectonic zones (ATZs) on a mesoscale, approximately 70 to several hundred kilometers in lateral extent. But it has been difficult to characterize ATZs in Japan, as they are in fact operational on multiple scales and our efforts are often hindered by various irregularities in the data. The key to overcoming these problems would be to gain an insight into the available data before any precise kinematic modeling is performed with indefinite assumptions. In this study, horizontal velocity fields, deduced from the nationwide GPS array, were treated with a set of techniques in robust smoothing and exploratory data analysis that brought out exceptionally powerful mesoscale ATZs, and made them easier to characterize. The resolved ATZs were then retrospectively monitored to study their regional and temporal variations, using a set of approx. 840 observation stations, about 30 km apart, for a 4-year series of fixed observation time-intervals, 810 days each. The smoothing operation involved three steps: (1) imputation of the velocity fields for the purpose of anti-aliasing, (2) robust smoothing of the velocity fields with the median operative, and (3) visualization of deformation-rate distributions in several coordinate independent parameters, and post-filtering. The geometrical resolvability of mesoscale ATZs was confirmed by calibrating the smoothing scheme against synthetic tectonic boundary models before it was applied to the case study in Japan. ATZs in Japan, which are essentially visible as systematic deviations in the velocity fields on the International Terrestrial Reference Frame (ITRF) and as strain rate anomalies, were highlighted sharply along some known tectonic zones, chains of active volcanoes, and areas above low seismic velocity anomalies in the crust and upper mantle, all of which generally paralleled the offshore trench axes. The geometrical agreements among the mapped ATZs and the physical anomalies in the crust are presumably due to their common structural weakness on the mesoscale. In the four main islands of Japan, all but 30–40% of the strain rate anomalies persisted during the entire 6 years of the case study period, while the rest sporadically appeared or disappeared in a period from several months to a few years. The transient shifts in the deformation rates were remarkably synchronous with some nearby major tectonic episodes: large earthquakes and slow events. Differential plate coupling strengths along the subduction zones can also be inferred from the persistent pattern of rotational strain rate anomalies forming clockwise and counterclockwise pairs along the Pacific. Our empirical observations suggest that the first-order features of interseismic crustal deformations in Japan can be characterized as collateral processes behaving in response to fluctuations of the tectonic stresses on multiple scales, likely influenced by changes of plate coupling strengths on the contiguous subduction faults.     

Upper and middle crustal deformation of an arc–arc collision across Hokkaido, Japan, inferred from seismic refraction/wide-angle reflection experiments

The Hidaka Collision Zone (HCZ), central Hokkaido, Japan, is a good target for studies of crustal evolution and deformation processes associated with an arc–arc collision. The collision of the Kuril Arc (KA) with the Northeast Japan Arc (NJA), which started in the middle Miocene, is considered to be a controlling factor for the formation of the Hidaka Mountains, the westward obduction of middle/lower crustal rocks of the KA (the Hidaka Metamorphic Belt (HMB)) and the development of the foreland fold-and-thrust belt on the NJA side. The “Hokkaido Transect” project undertaken from 1998 to 2000 was a multidisciplinary effort intended to reveal structural heterogeneity across this collision zone by integrated geophysical/geological research including seismic refraction/reflection surveys and earthquake observations. An E–W trending 227 km-long refraction/wide-angle reflection profile found a complicated structural variation from the KA to the NJA across the HCZ. In the east of the HCZ, the hinterland region is covered with 4–4.5 km thick highly undulated Neogene sedimentary layers, beneath which two eastward dipping reflectors were imaged in a depth range of 10–25 km, probably representing the layer boundaries of the obducting middle/lower crust of the KA. The HMB crops out on the westward extension of these reflectors with relatively high Vp (>6.0 km/s) and Vp/Vs (>1.80) consistent with middle/lower crustal rocks. Beneath these reflectors, more flat and westward dipping reflector sequences are situated at the 25–27 km depth, forming a wedge-like geometry. This distribution pattern indicates that the KA crust has been delaminated into more than two segments under our profile. In the western part of the transect, the structure of the fold-and-thrust belt is characterized by a very thick (5–8 km) sedimentary package with a velocity of 2.5–4.8 km/s. This package exhibits one or two velocity reversals in Paleogene sedimentary layers, probably formed by imbrication associated with the collision process. From the horizontal distribution of these velocity reversals and other geophysical/geological data, the rate of crustal shortening in this area is estimated to be greater than 3–4 mm/year, which corresponds to 40–50% of the total convergence rate between the NJA and the Eurasian Plate. This means that the fold-and-thrust belt west of the HCZ is absorbing a large amount of crustal deformation associated with plate interaction across Hokkaido Island.     

Tectonic features of the Lipari–Vulcano complex (Aeolian archipelago,Italy) from 10 years (1996–2006) of GPS data

The tectonic deformation of the Lipari–Vulcano complex, one of the most important active volcanic areas of the Mediterranean region, is studied here through the analysis of 10 years (1996–2006) of GPS data from both three permanent and 13 non‐permanent stations. This area can be considered crucial for the understanding of the interaction between the Eurasian and African plates in the Mediterranean area, and, in general, this work emphasizes a methodological approach, already applied in other areas worldwide ( J. Geophys. Res., 1996, 101 , 27 957 ; J. Geodyn., 1999, 27 , 213 ) where geodetic data and strain parameters maps of critical areas can help to improve our understanding of their geodynamical aspects. In this framework, this study is aimed at providing a kinematic deformation model on the basis of the dense geodetically estimated velocities of the Lipari–Vulcano complex. In particular, the observed deformation pattern can be described by a combination of (1) the main N–S regional compression and (2) a NNE–SSW compression with a small right‐lateral strike slip component acting along a tectonic structure trending N°40W between the two islands. This pattern was inspected through a simplified synthetic model.     

Seismotectonics of strike–slip earthquakes within the deep crust of southern Italy: Geometry, kinematics, stress field and crustal rheology of the Potenza 1990–1991 seismic sequences (Mmax 5.7)

We present a revision and a seismotectonic interpretation of deep crust strike–slip earthquake sequences that occurred in 1990–1991 in the Southern Apennines (Potenza area). The revision is motivated by: i) the striking similarity to a seismic sequence that occurred in 2002  140 km NNW, in an analogous tectonic context (Molise area), suggesting a common seismotectonic environment of regional importance; ii) the close proximity of such deep strike–slip seismicity with shallow extensional seismicity (Apennine area); and iii) the lack of knowledge about the mechanical properties of the crust that might justify the observed crustal seismicity. A comparison between the revised 1990–1991 earthquakes and the 2002 earthquakes, as well as the integration of seismological data with a rheological analysis offer new constraints on the regional seismotectonic context of crustal seismicity in the Southern Apennines. The seismological revision consists of a relocation of the aftershock sequences based on newly constrained velocity models. New focal mechanisms of the aftershocks are computed and the active state of stress is constrained via the use of a stress inversion technique. The relationships among the observed seismicity, the crustal structure of the Southern Apennines, and the rheological layering are analysed along a crustal section crossing southern Italy, by computing geotherms and two-mechanism (brittle frictional vs. ductile plastic strength) rheological profiles. The 1990–1991 seismicity is concentrated in a well-defined depth range (mostly between 15 and 23 km depths). This depth range corresponds to the upper pat of the middle crust underlying the Apulian sedimentary cover, in the footwall of the easternmost Apennine thrust system. The 3D distribution of the aftershocks, the fault kinematics, and the stress inversion indicate the activation of a right-lateral strike–slip fault striking N100°E under a stress field characterized by a sub-horizontal N142°-trending σ1 and a sub-horizontal N232°-trending σ3, very similar to the known stress field of the Gargano seismic zone in the Apulian foreland. The apparent anomalous depths of the earthquakes (> 15 km) and the confinement within a relatively narrow depth range are explained by the crustal rheology, which consists of a strong brittle layer at mid crustal depths sandwiched between two plastic horizons. This articulated rheological stratification is typical of the central part of the Southern Apennine crust, where the Apulian crust is overthrusted by Apennine units. Both the Potenza 1990–1991 and the Molise 2002 seismic sequences can be interpreted to be due to crustal E–W fault zones within the Apulian crust inherited from previous tectonic phases and overthrusted by Apennine units during the Late Pliocene–Middle Pleistocene. The present strike–slip tectonic regime reactivated these fault zones and caused them to move with an uneven mechanical behaviour; brittle seismogenic faulting is confined to the strong brittle part of the middle crust. This strong brittle layer might also act as a stress guide able to laterally transmit the deviatoric stresses responsible for the strike–slip regime in the Apulian crust and may explain the close proximity (nearly overlapping) of the strike–slip and normal faulting regimes in the Southern Apennines. From a methodological point of view, it seems that rather simple two-mechanism rheological profiles, though affected by uncertainties, are still a useful tool for estimating the rheological properties and likely seismogenic behaviour of the crust.     

Pb and Nd isotopic composition of the Jigongshan granite: constraints on crustal structure of Tongbaishan in the middle part of the Qinling–Tongbai–Dabie orogenic belt, Central China

The Qinling–Dabie–Sulu belt is the world's largest ultrahigh pressure (UHP) metamorphic belt. The UHP metamorphism is well dated at 220–245 Ma in the Dabie–Sulu belt but at 507 Ma in the Qinling belt. The Tongbaishan is located between the Qinling orogenic belt to the west and the Dabie–Sulu UHP metamorphic belt to the east. It is the key area for studying the tectonic relation between the Qinling and Dabie–Sulu belts and the diachronous UHP metamorphism. The Jigongshan granitic pluton (t=128 Ma) with a total area of 1200 km², composed of monzogranite, was mostly emplaced into the Tongbai complex, an exposed basement in the Tongbaishan. The Jigongshan granites have SiO₂=69.85–72.35%, K₂O/Na₂O=0.87–1.13, A/CNK=0.91–1.03, Rb/Sr=0.14–0.25 and Th/U=3.3–12. Their REE compositions show strongly fractionated patterns with (La/Yb)_N=14–58 and Eu*/Eu=0.79–1.05. The granites are characterized by low radiogenic Pb isotopic composition. The present-day whole-rock Pb isotopic ratios are ²⁰⁶Pb/²⁰⁴Pb=16.707–17.055, ²⁰⁷Pb/²⁰⁴Pb=15.239–15.326 and ²⁰⁸Pb/²⁰⁴Pb=37.587–37.853, which are similar to that of the continental lower crust. Their _Nd(t) values range from −16 to −20, and depleted-mantle Nd model ages (T_DM) from 1.8 to 2.2 Ga. The above evidence indicates that the magma of the Jigongshan granites was derived from the partial melting of the continental crust. The Pb and Nd isotopic compositions of the Jigongshan granites resemble those of the Dabie core complex in the Dabieshan but are distinct from those of the Tongbai complex in the Tongbaishan. Thus, the Dabie core complex would be the magma source of the Jigongshan granites. The result implies that the Dabie core complex is extended to the west and constitutes the unexposed basement underlaying the Tongbai complex in the Tongbaishan.     

U–Pb ages of Neoarchean granitoids from the Black Hills, South Dakota, USA: implications for crustal evolution in the Archean Wyoming province

Zircons in basement rocks from the eastern Wyoming province (Black Hills, South Dakota, USA) have been analyzed by ion microprobe (SHRIMP) in order to determine precise ages of Archean tectonomagmatic events. In the northern Black Hills (NBH) near Nemo, Phanerozoic and Proterozoic (meta)sedimentary rocks are nonconformably underlain by Archean biotite–feldspar gneiss (BFG) and Little Elk gneissic granite (LEG), both of which intrude older schists. The Archean granitoid gneisses exhibit a pervasive NW–SE-trending fabric, whereas an earlier NE–SW-trending fabric occurs sporadically only in the BFG, which is intruded by the somewhat younger LEG. Zircon crystals obtained from the LEG and BFG exhibit double terminations, oscillatory zoning, and Th/U ratios of 0.6±0.3—thereby confirming a magmatic origin for both lithologies. In situ analysis of the most U–Pb concordant domains yields equivalent ²⁰⁷Pb/²⁰⁶Pb ages (upper intercept, U–Pb concordia) of 2559±6 and 2563±6 Ma (both ±2σ) for the LEG and BFG, respectively, which constrains a late Neoarchean age for sequential pulses of magmatism in the NBH. Unzoned (in BSE) patches of 2560 Ma zircon commonly truncate coeval zonation in the same crystals with no change in Th/U ratio, suggesting that deuteric, fluid-assisted recrystallization accompanied post-magmatic cooling. A xenocrystic core of magmatic zircon observed in one LEG zircon yields a concordant age of 2894±6 Ma (±2σ). This xenocryst represents the oldest crustal material reported thus far in the Black Hills. Whether this older zircon originated as unmelted residue of 2900 Ma crust that potentially underlies the Black Hills or as detritus derived from 2900 Ma crustal sources in the Wyoming province cannot be discerned. In the southern Black Hills (SBH), the peraluminous granite at Bear Mountain (BMG) of previously unknown age intrudes biotite–plagioclase schist. Zircon crystals from the BMG are highly metamict and altered, but locally preserve small domains suitable for in situ analysis. A U–Pb concordia upper intercept age of 2596±11 Ma (±2σ) obtained for zircon confirms both the late Neoarchean magmatic age of the BMG and a minimum age for the schist it intrudes. Taken together, these data indicate that the Neoarchean basement granitoids were emplaced at 2590–2600 Ma (SBH) and 2560 Ma (NBH), most likely in response to subduction associated with plate convergence (final assembly of supercontinent Kenorland?). In contrast, thin rims present on some LEG–BFG zircons exhibit strong U–Pb discordance, high common Pb, and low Th/U ratios—suggesting growth or modification under hydrothermal conditions, as previously suggested for similar zircons from SE Wyoming. The LEG–BFG zircon rims yield a nominal upper intercept date of 1940–2180 Ma, which may represent a composite of multiple rifting events known to have affected the Nemo area between 2480 and 1960 Ma. Together, these observations confirm the existence of a Paleoproterozoic rift margin along the easternmost Wyoming craton. Moreover, the 2480–1960 Ma time frame inferred for rifting in the Black Hills (Nemo area) corresponds closely to a 2450–2100 Ma time frame previously inferred for the fragmentation of supercontinent Kenorland.     

Kinematics of the Iberia–Maghreb plate contact from seismic moment tensors and GPS observations

The Iberian Peninsula and the Maghreb experience moderate earthquake activity and oblique,  NW–SE convergence between Africa and Eurasia at a rate of  5 mm/yr. Coeval extension in the Alboran Basin and a N35°E trending band of active, left-lateral shear deformation in the Alboran–Betic region are not straightforward to understand in the context of regional shortening, and evidence complexity of deformation at the plate contact. We estimate 86 seismic moment tensors (M_W 3.3 to 6.9) from time domain inversion of near-regional waveforms in an intermediate period band. Those and previous moment tensors are used to describe regional faulting style and calculate average stress tensors. The solutions associated to the Trans-Alboran shear zone show predominantly strike-slip faulting, and indicate a clockwise rotation of the largest principal stress orientation compared to the regional convergence direction (σ₁ at N350°E). At the N-Algerian and SW-Iberian margins, reverse faulting solutions dominate, corresponding to N350°E and N310°E compression, respectively. Over most of the Betic range and intraplate Iberia, we observe predominately normal faulting, and WSW–ENE extension (σ₃ at N240°E). From GPS observations we estimate that more than 3 mm/yr of African (Nubian)–Eurasian plate convergence are currently accommodated at the N-Algerian margin,  2 mm/yr in the Moroccan Atlas, and  2 mm/yr at the SW-Iberian margin. 2 mm/yr is a reasonable estimate for convergence within the Alboran region, while Alboran extension can be quantified as  2.5 mm/yr along the stretching direction (N240°E). Superposition of both motions explains the observed left-lateral transtensional regime in the Trans-Alboran shear zone. Two potential driving mechanisms of differential motion of the Alboran–Betic–Gibraltar domain may coexist in the region: a secondary stress source other than plate convergence, related to regional-scale dynamic processes in the upper mantle of the Alboran region, as well as drag from the continental-scale motion of the Nubian plate along the southern limit of the region. In the Atlantic Ocean, the  3.5 mm/yr, westward motion of the Gibraltar Arc relative to intraplate Iberia can be accommodated at the transpressive SW-Iberian margin, while available GPS observations do not support an active subduction process in this area.     

Imaging the ramp–décollement geometry of the Chelungpu fault using coseismic GPS displacements from the 1999 Chi-Chi, Taiwan earthquake

We use coseismic GPS data from the 1999 Chi-Chi, Taiwan earthquake to estimate the subsurface shape of the Chelungpu fault that ruptured during the earthquake. Studies prior to the earthquake suggest a ramp–décollement geometry for the Chelungpu fault, yet many finite source inversions using GPS and seismic data assume slip occurred on the down-dip extension of the Chelungpu ramp, rather than on a sub-horizontal décollement. We test whether slip occurred on the décollement or the down-dip extension of the ramp using well-established methods of inverting GPS data for geometry and slip on faults represented as elastic dislocations. We find that a significant portion of the coseismic slip did indeed occur on a sub-horizontal décollement located at 8 km depth. The slip on the décollement contributes 21% of the total modeled moment release. We estimate the fault geometry assuming several different models for the distribution of elastic properties in the earth: homogeneous, layered, and layered with lateral material contrast across the fault. It is shown, however, that heterogeneity has little influence on our estimated fault geometry. We also investigate several competing interpretations of deformation within the E/W trending rupture zone at the northern end of the 1999 ground ruptures. We demonstrate that the GPS data require a 22- to 35-km-long lateral ramp at the northern end, contradicting other investigations that propose deformation is concentrated within 10 km of the Chelungpu fault. Lastly, we propose a simple tectonic model for the development of the lateral ramp.     

Sources of Phanerozoic granitoids in the transect Bayanhongor–Ulaan Baatar, Mongolia: geochemical and Nd isotopic evidence, and implications for Phanerozoic crustal growth

The Central Asian Orogenic Belt (CAOB) is renowned for massive generation of juvenile crust in the Phanerozoic. Mongolia is the heartland of the CAOB and it has been subject to numerous investigations, particularly in metallogenesis and tectonic evolution. We present new petrographic, geochemical and Sr–Nd isotopic analyses on Phanerozoic granitoids emplaced in west-central Mongolia. The data are used to delineate their source characteristics and to discuss implications for the Phanerozoic crustal growth in Central Asia. Our samples come from a transect from Bayanhongor to Ulaan Baatar, including three tectonic units: the Baydrag cratonic block (late Archean to middle Proterozoic), the Eo-Cambrian Bayanhongor ophiolite complex and the Hangay–Hentey Basin of controversial origin. The intrusive granitoids have ages ranging from ca. 540 to 120 Ma. The majority of the samples are slightly peraluminous and can be classified as granite (s.s.), including monzogranite, syenogranite and alkali feldspar granite. Most of the rocks have initial ⁸⁷Sr/⁸⁶Sr ratios between 0.705 and 0.707. Late Paleozoic to Mesozoic granitoids (≤250 Ma) are characterized by near-zero _Nd(T) values (0 to −2), whereas older granitoids show lower _Nd(T) values (−1.5 to −7). The data confirm the earlier observation of Kovalenko et al. [Geochemistry International 34 (1996) 628] who showed that granitoids emplaced outside of the Pre-Riphean basement rocks are characterized by juvenile positive _Nd(T) values, whereas those within the Pre-Riphean domain and the Baydrag cratonic block, as for the present case, show a significant effect of ‘contamination’ by Precambrian basement rocks. Nevertheless, mass balance calculation suggests that the granitoids were derived from sources composed of at least 80% juvenile mantle-derived component. Despite our small set of new data, the present study reinforces the general scenario of massive juvenile crust production in the CAOB with limited influence of old microcontinents in the genesis of Phanerozoic granitoids.     

Zircon U–Pb ages, Hf and O isotopes constrain the crustal architecture of the ultrahigh-pressure Dabie orogen in China

The crustal structure of the Dabie orogen was reconstructed by a combined study of U–Pb ages, Hf and O isotope compositions of zircons from granitic gneiss from North Dabie, the largest lithotectonic unit in the orogen. The results were deciphered from metamorphic history to protolith origin with respect to continental subduction and exhumation. Zircon U–Pb dating provides consistent ages of 751 ± 7 Ma for protolith crystallization, and two group ages of 213 ± 4 to 245 ± 17 Ma and 126 ± 4 to 131 ± 36 Ma for regional metamorphism. Majority of zircon Hf isotope analyses displays negative ε_Hf(t) values of − 5.1 to − 2.9 with crust Hf model ages of 1.84 to 1.99 Ga, indicating protolith origin from reworking of middle Paleoproterozoic crust. The remaining analyses exhibit positive ε_Hf(t) values of 5.3 to 14.5 with mantle Hf model ages of 0.74 to 1.11 Ga, suggesting prompt reworking of Late Mesoproterozoic to Early Neoproterozoic juvenile crust. Zircon O isotope analyses yield δ¹⁸O values of − 3.26 to 2.79‰, indicating differential involvement of meteoric water in protolith magma by remelting of hydrothermally altered low δ¹⁸O rocks. North Dabie shares the same age of Neoproterozoic low δ¹⁸O protolith with Central Dabie experiencing the Triassic UHP metamorphism, but it was significantly reworked at Early Cretaceous in association with contemporaneous magma emplacement. The Rodinia breakup at about 750 Ma would lead to not only the reworking of juvenile crust in an active rift zone for bimodal protolith of Central Dabie, but also reworking of ancient crust in an arc-continent collision zone for the North Dabie protolith. The spatial difference in the metamorphic age (Triassic vs. Cretaceous) between the northern and southern parts of North Dabie suggests intra-crustal detachment during the continental subduction. Furthermore, the Dabie orogen would have a three-layer structure prior to the Early Cretaceous magmatism: Central Dabie in the upper, North Dabie in the middle, and the source region of Cretaceous magmas in the lower.     

Crustal lineaments, distribution of lower crustal intrusives and structural evolution of the Vøring Margin, NE Atlantic; new insight from wide-angle seismic models

Five lineaments on the volcanic Vøring Margin, NE Atlantic, have been identified in crustal scale models derived from Ocean Bottom Seismograph (OBS) data. It is suggested that the Vøring Basin can be divided in four compartments bounded by the Jan Mayen Fracture Zone/Lineament, a new lineament defined from this study, the Gleipne Lineament, the Surt Lineament and the Bivrost Lineament. The NW–SE trending Jan Mayen-, Gleipne- and Bivrost lineaments probably represent old zones of weakness controlling the onset of the early Eocene seafloor spreading, whereas the Surt- and New lineaments, rotated ca. 30° symmetrically from the azimuth of the Gleipne Lineament, may represent adjustment features related to the early Cretaceous/early Tertiary rifting. The longest landward extent of a lower crustal high-velocity body, assumed to represent intrusions related to the last phase of rifting, is found between the New Lineament and the Gleipne Lineament, where the body extends across the Helland Hansen Arch. Northeastwards in the Vøring Basin, the landward limit of the body steps gradually seawards, closely related to the interpreted lineaments. Northeast of the Gleipne Lineament, the body terminates close to the Fles Fault Complex, north of the Surt Lineament, it extends across the Nyk High, and northeast of the Bivrost Lineament the intrusions terminate around the Vøring Escarpment. Evidence for an interplay between active and passive rifting components is found on regional and local scales on the margin. The active component is evident through the decrease in magmatism with increased distance from the Icelandic plume, and the passive component is documented through the fact that all found crustal lineaments to a certain degree acted as barriers to magma emplacement. The increased thickness of the continental crust on the seaward side of the Vøring Escarpment, the upwarping of Moho and thinning of the lower crustal high-velocity layer in the western part of the Vøring Basin, as well as a strong shallowing of the Moho observed in parts of the area between the Jan Mayen Fracture Zone/Lineament and the New Lineament, can be explained by lithospheric delamination models.     

Age and kinematics of ductile deformation in the Cerro Durazno area, NW Argentina: Significance for orogenic processes operating at the western margin of Gondwana during Ordovician – Silurian times

The Cerro Durazno Pluton belongs to a suite of Paleozoic granitoid intrusions in NW-Argentina, that are central for understanding the tectonic setting of the western margin of Gondwana in Ordovician and Silurian times. The pluton and its host rocks were tectonically overprinted by metamorphic mineral shape fabrics formed under middle greenschist-facies metamorphic conditions and associated with the nearby Agua Rosada Shear Zone. Kinematic analysis of the shear zone based on the geometric relationship between individual segments of the shear plane and principal axes of mineral fabric ellipsoids indicates reverse-sense of shear with a minor component of left-lateral displacement. This is compatible with the kinematics of other ductile deformation zones in this area, collectively forming a network, which accomplished orogen-parallel extension in addition to vertical thickening. Using the Rb–Sr isotopic system, an undeformed pegmatite dike of the Cerro Durazno Pluton was dated at 455.8 ± 3.6 Ma and mineral fabrics of the Agua Rosada Shear Zone formed at middle greenschist-facies metamorphism gave deformation ages of 437.0 ± 3.8 Ma and 428.4 ± 4.5 Ma. Thus, tectonic overprint at low metamorphic grade occurred about 20–30 Ma after terminal magmatism in the Cerro Durazno area. Our data from the Cerro Durazno area and regional considerations suggest that the western margin of Gondwana was characterized by orogen-parallel extension in addition to crustal thickening as well as episodes of magmatism and ductile deformation that varied greatly in time and space.