Anelasticity of the crust and upper mantle beneath north and central India from the inversion of observed Love and Rayleigh wave attenuation data

The fundamental mode Love and Rayleigh waves generated by earthquakes occurring in Kashmir, Nepal Himalaya, northeast India and Burma and recorded at Hyderabad, New Delhi and Kodaikanal seismic stations are analysed. Love and Rayleigh wave attenuation coefficients are obtained at time periods of 15–100 seconds, using the spectral amplitude of these waves for 23 different paths along northern (across Burma to New Delhi) and central (across Kashmir, Nepal Himalaya and northeast India to Hyderabad and Kodaikanal) India. Love wave attenuation coefficients are found to vary from 0.0003 to 0.0022 km^–1 for northern India and 0.00003 km^–1 to 0.00016 km^–1 for central India. Similarly, Rayleigh wave attenuation coefficients vary from 0.0002 km^–1 to 0.0016 km^–1 for northern India and 0.00001 km^–1 to 0.0009 km^–1 for central India. Backus and Gilbert inversion theory is applied to these surface wave attenuation data to obtainQ ^–1 models for the crust and uppermost mantle beneath northern and central India. Inversion of Love and Rayleigh wave attenuation data shows a highly attenuating zone centred at a depth of 20–80 km with lowQ for northern India. Similarly, inversion of Love and Rayleigh wave attenuation data shows a high attenuation zone below a depth of 100 km. The inferred lowQ value at mid-crustal depth (high attenuating zone) in the model for northern India can be by underthrusting of the Indian plate beneath the Eurasian plate which has caused a low velocity zone at this shallow depth. The gradual increase ofQ ^–1 from shallow to deeper depth shows that the lithosphere-asthenosphere boundary is not sharply defined beneath central India, but rather it represents a gradual transformation, which starts beneath the uppermost mantle. The lithospheric thickness is 100 km beneath central India and below that the asthenosphere shows higher attenuation, a factor of about two greater than that in the lithosphere. The very lowQ can be explained by changes in the chemical constitution taking place in the uppermost mantle.     

印度板块与欧亚板块在兴都库什-帕米尔地区相互俯冲的动力作用分析

兴都库什-帕米尔地区是印度板块与欧亚板块相互碰撞的强烈变形区域,也是中国大陆与周边板块动力传递的关键部位,明确该地区两大板块俯冲接触的几何形态和动力作用对研究区域动力环境具有实际意义.本文首先基于Hayes等在2009和2010年提出的Slab1.0的研究思路,利用地震参数准定量地给出两大板块在兴都库什-帕米尔地区碰撞接触的几何形态.结果表明,印度板块在兴都库什地区呈现自南往北的俯冲;欧亚板块在帕米尔地区呈现由北往南的俯冲;同时在兴都库什和帕米尔之间存在俯冲交汇区,在该区印度板块以北西方向、欧亚板块以南东方向相互俯冲.其次基于哈佛大学提供的震源机制解,对不同接触部位进行了应力张量反演,结果显示在兴都库什俯冲区域主要表现为逆冲性质,帕米尔弧西段主要表现为走滑性质,且均具有较好的一致性;而在俯冲交汇区域,走滑、逆冲性质并存,表现为震源机制一致性紊乱.结合两大板块接触的几何形态和区域应力场反演结果,认为印度板块在兴都库什地区主动往北俯冲,而欧亚板块在帕米尔地区被动往南东-南向俯冲,形成两大板块的相互俯冲.本文从几何形态和应力场反演分析两大板块在兴都库什-帕米尔地区碰撞的动力作用方式,可为该区域地球动力学相关研究提供基础资料.     

Eastward subduction of the Indian plate beneath the Indo-Myanmese arc revealed by teleseismic P-wave tomography

The deep structure of the eastward-subducting Indian plate can provide new information on the dynamics of the India-Eurasia collision. We collected and processed waveform data from temporary seismic arrays (networks) on the eastern Tibetan Plateau, seismic arrays in Northeast India and Myanmar, and permanent stations of the China Digital Seismic Network in Tibet, Gansu, Qinghai, Yunnan, and Sichuan. We combined these data with phase reports from observation stations of the International Seismological Center on the Indian plate and selected 124,808 high-quality P-wave relative travel-time residuals. Next, we used these data to invert the 3-D P-wave velocity structure of the upper mantle to a depth of 800 km beneath the eastern segment of the arcuate Himalayan orogen, at the southeastern margin of the Tibetan Plateau. The results reveal a high-angle, easterly dipping subducting plate extending more than 200 km beneath the Indo-Myanmese arc. The plate breaks off at roughly 96°E; its fragments have passed through the 410-km discontinuity (D410) into the mantle transition zone (MTZ). The MTZ beneath the Tengchong volcanic area contains a high-velocity anomaly, which does not exceed the Red River fault to the east. No other large-scale continuous subducted plates were observed in the MTZ. However, a horizontally spreading high-velocity anomaly was identified on the D410 in some regions. The anomaly may represent the negatively buoyant 90°E Ridge plate or a thickened and delaminated lithospheric block experiencing collision and compression at the southeastern margin of the Tibetan Plateau. The Tengchong volcano may originate from the mantle upwelling through the slab window formed by the break-off of the subducting Indian continental plate and oceanic plate in the upper mantle. Low-velocity upper mantle materials on the west side of the Indo-Myanmese arc may have supplemented materials to the Tengchong volcano.     

Regionality of deep low-frequency earthquakes associated with subduction of the Philippine Sea plate along the Nankai Trough, southwest Japan

The Fukuoka District Meteorological Observatory recently logged three possible deep low-frequency earthquakes (LFEs) beneath eastern Kyushu, Japan, a region in which LFEs and low-frequency tremors have never before been identified. To assess these data, we analyzed band-pass filtered velocity seismograms and relocated LFEs and regular earthquakes using the double-difference method. The results strongly suggest that the three events were authentic LFEs, each at a depth of about 50 km. We also performed relocation analysis on LFEs recorded beneath the Kii Peninsula and found that these LFEs occurred near the northwest-dipping plate interface at depths of approximately 29–38 km. These results indicate that LFEs in southwest Japan occur near the upper surface of the subducting Philippine Sea (PHS) plate. To investigate the origin of regional differences in the occurrence frequency of LFEs in western Shikoku, the Kii Peninsula, and eastern Kyushu, we calculated temperature distributions associated with PHS plate subduction. Then, using the calculated thermal structures and a phase diagram of water dehydration for oceanic basalt, the water dehydration rate (wt.%/km), which was newly defined in this study, was determined to be 0.19, 0.12, and 0.08 in western Shikoku, the Kii Peninsula, and eastern Kyushu, respectively; that is, the region beneath eastern Kyushu has the lowest water dehydration rate value. Considering that the Kyushu–Palau Ridge that is subducting beneath eastern Kyushu is composed of tonalite, which is low in hydrous minerals, this finding suggests that the regionality may be related to the amount of water dehydration associated with subduction of the PHS plate and/or differences in LFE depths. Notable dehydration reactions take place beneath western Shikoku and the Kii Peninsula, where the depth ranges for dehydration estimated by thermal modeling agree well with those for the relocated LFEs. The temperature range in which LFEs occur in these regions is estimated to be 400–500 °C.     

Rotational overthrusting of the northwestern Himalaya: further palaeomagnetic evidence from the Riasi thrust sheet, Jammu foothills, India

Thermal demagnetization results (316 samples) are presented for the Tertiary succession of the Riasi thrust sheet (Jammu foothills, northwestern Himalaya). Primary and secondary magnetization directions of Murree Group red beds (Miocene to Upper Eocene) sampled northeast of Jammu indicate, for this part of the Riasi thrust sheet, a clockwise rotation over about 45° with respect to the Indian shield since Late Eocene/Early Miocene time. This accords with clockwise rotations of similar magnitude observed in the Panjal Nappe and the Krol Belt, and is interpreted as representative for the northwestern Himalaya. Results from the western part of the Kalakot inlier, sampled northwest of Jammu, i.e. basal Murree claystone (Middle Eocene) and carbonate from the Subathu Group (lower Middle to Lower Eocene), indicate an aberrant 20–25° counterclockwise rotation which is of local importance only. Available observations on rotation of Himalayan thrust sheets with respect to the Indian shield, indicate that the Himalayan Arc has formed through oroclinal bending. This supports Powell and Conaghan's and Veevers et al.'s model of Greater India with large-scale intracontinental underthrusting along the Main Central Thrust beneath the Tibetan Plateau. Minimal magnitudes of underthrusting of 550 km in the Krol Belt and 650 km in the Thakkhola region are concluded. Palaeolatitude observations (herein and in [1[) agree with absolute positioning of the Indian plate based on India-Africa relative movement data fixed to a hotspot frame in the Atlantic Ocean, and with palaeolatitude observations from DSDP cores on the Indian plate. Collision-related secondary magnetic components observed both to the north and to the south of the Indus-Tsangpo Suture zone show palaeolatitudes between the equator and 7°N. Comparison of both datasets indicates that initial contact between Greater India and south-central Asia had been established in the Hindu Kush—Karakorum region by about 60 Ma ago whereas eastwards progressive suturing had advanced to the Lhasa Block segment of the Indus-Tsangpo Suture zone before 50 Ma ago.     

帕米尔—兴都库什地区构造应力场反演及拆离板片应力形因子特征研究

从Global CMT目录搜集了1976年1月至2016年6月之间的震源深度大于70km的255个震源机制解,用阻尼应力反演方法,分70～160km和170～310km两个深度,计算了帕米尔—兴都库什地区的构造应力场;同时以10km为间隔计算了兴都库什地区深度介于70～310km之间的应力形因子.得到以下初步结论:兴都库什板片向下俯冲和帕米尔地区断裂带的横向拉张,可能是导致应力场不同的原因.兴都库什俯冲带与帕米尔俯冲带碰撞,导致交汇地区(37°N—37.5°N)的应力场参数突变.兴都库什俯冲板片受到深部温度、压力等因素,出现薄弱面进而形成拆离板片.其脱离了主俯冲板片的束缚后,重力的上下拉张作用导致空区附近张轴倾伏角接近90°,拆离板片俯冲至上地幔不连续面,导致板片部分熔融进而应力形因子随着深度变小.而拆离板片受到地幔挤压其内部发生破碎,其压应力轴由西部的NS到东部NW-SE方向偏转及纵向张应力轴倾伏角变小.     

中国西部及邻区岩石圈S波速度结构面波层析成像

本文利用瑞利波群速度频散资料和层析成像方法,研究了中国西部及邻近区域(20°N—55°N,65°E—110°E)的岩石圈S波速度结构.结果表明这一地区存在三个以低速地壳/上地幔为特征的构造活动区域:西蒙古高原—贝加尔地区,青藏高原,印支地区.西蒙古高原岩石圈厚度约为80 km,上地幔低速层向下延伸至300 km深度,说明存在源自地幔深部的热流活动.缅甸弧后的上地幔低速层下至200 km深度,显然与印度板块向东俯冲引起俯冲板片上方的热/化学活动有关.青藏高原地壳厚达70 km,边缘地区厚度也在50 km以上并且具有很大的水平变化梯度,与高原平顶陡边的地形特征一致.中下地壳的平均S波速度明显低于正常大陆地壳,在中地壳20~40 km深度范围广泛存在速度逆转的低速层,这一低速层的展布范围与高原的范围相符.这些特征说明青藏高原中下地壳的变形是在印度板块的北向挤压下发生塑性增厚和侧向流动.地幔的速度结构呈现与地壳显著不同的特点.在高原主体和川滇西部地区上地幔顶部存在较大范围的低速,低速区范围随深度迅速减小;100 km以下滇西低速消失,150 km以下基本完全消失.青藏高原上地幔速度结构沿东西方向表现出显著的分段变化.在大约84°E以西的喀喇昆仑—帕米尔—兴都库什地区,印度板块的北向和亚洲板块的南向俯冲造成上地幔显著高速;84°E—94°E之间上地幔顶部速度较低,在大约150~220 km深度范围存在高速板片,有可能是俯冲的印度岩石圈,其前缘到达昆仑—巴颜喀拉之下;在喜马拉雅东构造结以北区域,存在显著的上地幔高速区,可能阻碍上地幔物质的东向运动.川滇西部岩石圈底界深度与扬子克拉通相似,约为180 km,但上地幔顶部速度较低.这些现象表明青藏高原岩石圈地幔的变形/运动方式可能与地壳有本质的区别.     

The geometry of the Burmese-Andaman subducting lithosphere

The gross seismotectonic features for the Burmese-Andaman arc system which defines the northeast margin of the Indian plate are rather well known but variations in the subduction zone geometry along and across the arc and fault pattern within the subducting Indian plate have not been studied. Present workaims to study these by using seismicity data whose results are presented in the form of: (a) Lithospheric across-the-arc sections at about every 100–120 km (approximately one degree latitude apart) covering the 3500 km longBurmese-Andaman arc system, (b) a structure contour map showing the depth tothe top surface of the seismically active lithosphere and (c) interpretationof focal mechanism solutions for 148 Benioff zone earthquakes. Both penetrationdepth and the dip of the Benioff zone vary considerably along the arc in correspondence to the curvature of the fold-thrust belt which varies from concave to convex in different sectors of the arc. Several extensive `Hinge faults' that abut at high angles to the arc orientation, are inferred from aninterpretation of the structure contour map. Active nature of the hinge faultsis established in several areas by their association with earthquakes andcorroborated through fault plane solutions. At shallow level of the Benioffzone along these faults, focal mechanism solutions display left lateral strikeslip movement while at deeper levels reverse fault solutions are common.     

Shear wave velocity structure of the upper mantle under the NW Indian Ocean

Shear wave velocity structure of the NW Indian ocean is analysed by using fundamental mode Rayleigh wave dispersion data of 67 events occurred during 1990–98 at the central Indian Ridge and Carlsberg Ridge and recorded at Hyderabad Geoscope station (HYB). These events provide a dense coverage of the NW Indian ocean and Chagos-Laccadive Ridge (CLR) in the back-azimuthal range of 192–253° with respect to HYB. The dispersion curves, corrected for continental and young ocean paths, indicate large variations in the shear wave velocity structure of the region. The group velocities along the CLR path support a typical aseismic ridge-type structure. However, the central region bounded between the Central Indian Ridge and India in the back-azimuth of 206–234° indicates a decrease in the group velocity by 0.1 km/s. Inversion of these data sets indicates presence of aseismic-ridge type lithospheric structure for CLR, a thin lithosphere and high velocity block in the depth range of 125–200 km for the central region, and a continental-type lithospheric structure for the northern-most part of the Indian ocean. It is inferred that the dynamic state of the upper mantle in this region has been significantly perturbed during the recent geological past.     

Historical 1942 Ecuador and 1942 Peru subduction earthquakes and earthquake cycles along Colombia-Ecuador and Peru subduction segments

Two large shallow earthquakes occurred in 1942 along the South American subduction zone inclose proximity to subducting oceanic ridges: The 14 May event occurred near the subducting Carnegie ridge off the coast of Ecuador, and the 24 August event occurred off the coast of southwestern Peru near the southern flank of the subducting Nazca ridge. Source parameters for these for these two historic events have been determined using long-periodP waveforms,P-wave first motions, intensities and local tsunami data.We have analyzed theP waves for these two earthquakes to constrain the focal mechanism, depth, source complexity and seismic moment. Modeling of theP waveform for both events yields a range of acceptable focal mechanisms and depths, all of which are consistent with underthrusting of the Nazca plate beneath the South American plate. The source time function for the 1942 Ecuador event has one simple pulse of moment release with a duration of 22 suconds, suggesting that most of the moment release occurred near the epicenter. The seismic moment determined from theP waves is 6–8×10²⁰N·m, corresponding ot a moment magnitude of 7.8–7.9. The reported location of the maximum intensities (IX) for this event is south of the main shock epicenter. The relocated aftershcks are in an area that is approximately 200 km by 90 km (elongated parallel to the trench) with the majority of aftershocks north of the epicenter. In contrast, the 1942 Peru event has a much longer duration and higher degree of complexity than the Ecuador earthquake, suggesting a heterogeneous rupture. Seismic moment is released in three distinct pulses over approximately 74 seconds; the largest moment release occurs 32 seconds after rupture initiation. the seismic moment as determined from theP waves for the 1942 Peru event is 10–25×10²⁰N·m, corresponding to a moment magnitude of 7.9–8.2. Aftershock locations reported by the ISS occur over a broad area surrounding the main shock. The reported locations of the maximum intensities (IX) are concentrated south of the epicenter, suggesting that at least part of the rupture was to the south.We have also examined great historic earthquakes along the Colombia-Ecuador and Peru segments of the South American subduction zone. We find that the size and rupture length of the underthrusting earthquakes vary between successive earthquake cycles. This suggests that the segmentation of the plate boundary as defined by earthquakes this century is not constant.     

Moment Tensor Solutions of Earthquakes in the Indian Ocean from Surface Waves and their Tectonic Implications

—Rayleigh and Love waves generated by sixteen earthquakes which occurred in the Indian Ocean and were recorded at 13 WWSSN stations of Asia, Africa and Australia are used to determine the moment tensor solution of these earthquakes. A combination of thrust and strike-slip faulting is obtained for earthquakes occurring in the Bay of Bengal. Thrust, strike slip or normal faulting (or either of the combination) is obtained for earthquakes occurring in the Arabian Sea and the Indian Ocean. The resultant compressive and tensional stress directions are estimated from more than 300 centroid moment tensor (CMT) solution of earthquakes occurring in different parts of the Indian Ocean. The resultant compressive stress directions are changing from north-south to east-west and the resultant tensional stress directions from east-west to north-south in different parts of the Indian Ocean. The results infer the counterclockwise movement of the region (0°–33°S and 64°E–94°E), stretching from the Rodriguez triple junction to the intense deformation zone of the central Indian Ocean and the formation of a new subduction zone (island arc) beneath the intense deformation zone of the central Indian Ocean and another at the southern part of the central Indian basin. The compressive stress direction is along the ridge axis and the extensional stress manifests across the ridge axis. The north-south to northeast-south west compression and east-west to northwest-southeast extension in the Indian Ocean suggest the northward underthrusting of the Indian plate beneath the Eurasian plate and the subduction beneath the Sunda arc region in the eastern part. The focal depth of earthquakes is estimated to be shallow, varying from 4 to 20 km and increasing gradually in the age of the oceanic lithosphere with the focal depth of earthquakes in the Indian Ocean.     

Gravity anomaly, lithospheric structure and seismicity of Western Himalayan Syntaxis

A compiled gravity anomaly map of the Western Himalayan Syntaxis is analysed to understand the tectonics of the region around the epicentre of Kashmir earthquake of October 8, 2005 (Mw = 7.6). Isostatic gravity anomalies and effective elastic thickness (EET) of lithosphere are assessed from coherence analysis between Bouguer anomaly and topography. The isostatic residual gravity high and gravity low correspond to the two main seismic zones in this region, viz. Indus–Kohistan Seismic Zone (IKSZ) and Hindu Kush Seismic Zones (HKSZ), respectively, suggesting a connection between siesmicity and gravity anomalies. The gravity high originates from the high-density thrusted rocks along the syntaxial bend of the Main Boundary Thrust and coincides with the region of the crustal thrust earthquakes, including the Kashmir earthquake of 2005. The gravity low of HKSZ coincides with the region of intermediate–deep-focus earthquakes, where crustal rocks are underthrusting with a higher speed to create low density cold mantle. Comparable EET (∼55 km) to the focal depth of crustal earthquakes suggests that whole crust is seismogenic and brittle. An integrated lithospheric model along a profile provides the crustal structure of the boundary zones with crustal thickness of about 60 km under the Karakoram–Pamir regions and suggests continental subduction from either sides (Indian and Eurasian) leading to a complex compressional environment for large earthquakes.     

帕米尔—兴都库什深俯冲残留体对410km间断面起伏形态的影响

受俯冲残留体影响的410km间断面起伏形态的研究对于确定地球内部物质构成及地球动力学过程具有重要作用.帕米尔—兴都库什俯冲区域拥有全球少有的中、深源地震,为研究410km间断面起伏提供了良好的资源.利用日本Hi-net地震台网和美国TA台阵记录的帕米尔—兴都库什俯冲区域的6个震源深度为154.0~220.9km、震级为Mb5.6~6.4的中、深源地震的短周期/宽频带波形资料,经过4次根倾斜叠加处理,获得了36组Hi-net子台网和TA记录资料的倾斜叠加灰度图,从中提取了与410km间断面相关的次生转换震相SdP,发现受俯冲残留体影响下的410km间断面的深度位于372~398km.较之持续俯冲的西太平洋地区海洋岩石圈,研究区域俯冲滞留体对于410km间断面的相变线的影响要小得多.     

Tectonics and paleomagnetism of structural arcs of the Pamir-Punjab syntaxis

The Pamir-Punjab syntaxis consists of two structural arcs with the crests facing northward—the Pamir and Hindu-Kush-Karakorum arcs. These arcs are mutually disharmonious, and the exterior (Pamir) arc is more tight, as compared to the inner one. Paleomagnetic study of the Pamir arc has shown that the structures of the future Pamirs had a northeastern strike in the Paleogene and Cretaceous, and they occurred on the eastern limb of the Darvaz-Kopetdag structural arc, whose crest faces south. The Pamir arc originated after the Paleogene during the process of the formation of the Pamir-Punjab syntaxis. Knowledge of the kinematics of the Pamir arc, combined with data on the geometry of the syntaxis and the character of its boundaries, enable us to choose a model of the development of the syntaxis. The process of “plastic flow” of crustal masses around the underthrusting segment of the Indian plate was likely of paramount importance.     

Eruption characteristics and cycles at Pavlof Volcano, Alaska, and their relation to regional earthquake activity

Pavlof Volcano (55° 25′N, 161° 54′W) exhibits two eruption styles: magmatic eruptions of one-to-two-days duration, and phreatic-phreatomagmatic activity lasting several days to two months. Thirty-four eruptions have occurred in historic times; of these the largest are Volcano Explosivity Index=3. Nine magmatic and 13 phreatomagmatic eruptions occurred between 1973–1983. All the magmatic eruptions occurred in the fall, between Sept. 9–Nov. 20. Four magmatic eruptions occurred during November 11–15, but in four different years. A 3-year-long period of eruptive activity between 1973–1976 bears striking resemblance to a period of activity between 1980–1983. No locatable shallow earthquakes (<50 km) have occurred within 30 km of Pavlof since 1973, which is quite unusual for an active island-arc volcano. Shallow events in the adjacent are segments have focal mechanisms with P-axes perpendicular to the arc (and parallel to plate convergence). Deep earthquakes (> 100 km) are clustered beneath Pavlof and several other volcanoes. Their T-axes show downdip tension within the slab. Deep teleseisms (> 160 km) mostly occurred between 1977–1979 when the volcano was not erupting. Catalogued volcanic activity throughout the Alaska/Aleutian arc shows a weak tendency to increase around the time of great (M > 7.8) earthquakes.     

Initiation of subduction by post-collision foreland thrusting and back-thrusting

While postulated causes of initial subduction and trench formation include underthrusting, controls on its location and age have not been determined. Consideration of the age of subduction zones bordering five collisional orogens suggests that subduction may have been initiated by foreland thrusts and back-thrusts. Foreland thrusts develop within a continental foreland on the subducting plate mostly within 50 my of collision with an arc system; where the foreland is narrow the thrusts may intersect the continent-ocean crust boundary. Back-thrusts develop in the fore-arc or back-arc area on the overriding plate within 10 to 20 my of collision, and can result in tectonic burial of the magmatic arc; where the arc system is oceanic the back-thrusts may intersect the arc-ocean crust boundary. Possible examples of subduction initiated by foreland thrusts are the start of subduction in the late Jurassic beneath the northern Sunda Arc, and at the end-Miocene in the Negros Trench. Examples of back-thrusts which have initiated or may initiate subduction are the late Cenozoic eastward translation of Taiwan over the Philippine Sea plate, the incipient southward subduction of the Banda Sea beneath Timor, and the W-dipping back-thrust comprising the Highland Boundary Fault zone and postulated early Ordovician thrusts to the SE in Scotland. The suggested relationship of subduction to collision helps to explain the persistence of Wilson cycles in the still-active late Mesozoic to Cenozoic orogenic belts and implies that orogeny will cease only with collision between major continents.     

Geophysical modeling of the northern Appalachian Brompton-Cameron, Central Maine, and Avalon terranes under the New Jersey Coastal Plain

A regional terrane map of the New Jersey Coastal Plain basement was constructed using seismic, drilling, gravity and magnetic data. The Brompton-Cameron and Central Maine terranes were coalesced as one volcanic island arc terrane before obducting onto Laurentian, Grenville age, continental crust in the Taconian orogeny [Rankin, D.W., 1994. Continental margin of the eastern United States: past and present. In: Speed, R.C., (Ed.), Phanerozoic Evolution of North American Continent-Ocean Transitions. DNAG Continent-Ocean Transect Volume. Geological Society of America, Boulder, Colorado, pp. 129–218]. Volcanic island-arc rocks of the Avalon terrane are in contact with Central Maine terrane rocks in southern Connecticut where the latter are overthrust onto the Brompton-Cameron terrane, which is thrust over Laurentian basement. Similarities of these allochthonous island arc terranes (Brompton-Cameron, Central Maine, Avalon) in lithology, fauna and age suggest that they are faulted segments of the margin of one major late Precambrian to early Paleozoic, high latitude peri-Gondwana island arc designated as “Avalonia”, which collided with Laurentia in the early to middle Paleozoic. The Brompton Cameron, Central Maine, and Avalon terranes are projected as the basement under the eastern New Jersey Coastal Plain based on drill core samples of metamorphic rocks of active margin/magmatic arc origin. A seismic reflection profile across the New York Bight traces the gentle dipping (approximately 20 degrees) Cameron's Line Taconian suture southeast beneath allochthonous Avalon and other terranes to a 4 sec TWTT depth (approximately 9 km) where the Avalonian rocks are over Laurentian crust. Gentle up-plunge (approximately 5 degrees) projections to the southwest bring the Laurentian Grenville age basement and the drift-stage early Paleozoic cover rocks to windows in Burlington Co. at approximately 1 km depth and Cape May Co. at approximately 2 km depths. The antiformal Shellburne Falls and Chester domes and Chain Lakes-Pelham dome-Bronson Hill structural trends, and the synformal Connecticut Valley-Gaspe structural trend can be traced southwest into the New Jersey Coastal Plain basement. A Mesozoic rift basin, the “Sandy Hook basin”, and associated eastern boundary fault is identified, based upon gravity modeling, in the vicinity of Sandy Hook, New Jersey. The thickness of the rift-basin sedimentary rocks contained within the “Sandy Hook basin” is approximately 4.7 km, with the basin extending offshore to the east of the New Jersey coast. Gravity modeling indicates a deep rift basin and the magnetic data indicates a shallow magnetic basement caused by magnetic diabase sills and/or basalt flows contained within the rift-basin sedimentary rocks. The igneous sills and/or flows may be the eastward continuation of the Watchung and Palisades bodies.     

Source depletion and extent of melting in the Tongan sub-arc mantle

The fluid immobile High Field Strength Elements (HFSE) Nb and Ta can be used to distinguish between the effects of variable extents of melting and prior source depletion of the Tongan sub-arc mantle. Melting of spinel lherzolite beneath the Lau Basin back-arc spreading centres has the ability to fractionate Nb from Ta due to the greater compatibility of the latter in clinopyroxene. The identified spatial variation in plate velocities and separation of melt extraction zones, combined with extremely depleted lavas make Tonga an ideal setting in which to test models for arc melt generation and the role of back-arc magmatism.We present new data acquired by laser ablation-ICPMS of fused sample glasses produced without the use of a melt fluxing agent. The results show an arc trend towards strongly sub-chondritic Nb/Ta (< 17) with values as low as 7.2. Melting models show that large degree melts of depleted MORB mantle fail to reproduce the observed Nb/Ta. Alternatively, incorporation of residual back-arc mantle that has undergone less than 1% melting into the sub-arc melting regime reproduces arc values. However, the extent of partial melting required to produce the composition of the Lau Basin back-arc basalts averages 7%. This apparent discrepancy can be explained if only the lowermost 4 km of the residua from the mantle melt column beneath the back-arc is added to the source of arc magmas. We have identified that the degree of arc/back-arc coupling displayed in the rock record provides an index of the depth of hydrous melting beneath the arc. In this case, this would imply a depth of ~ 75 km for generation of arc magmas, indicating that hydrous melting in the mantle wedge is triggered by the breakdown of hydrous phases in the subducting slab.     

Thermal state, rheology and seismicity in the pannonian basin, Hungary

On the basis of data on crustal structure and terrestrial heat flow, a 3-D geothermal model for the lithosphere in the Pannonian basin, Hungary, has been calculated. This model, together with information on crustal composition, laboratory data on rock friction, and certain assumptions about fluid conditions and strain-rate levels within the lithosphere, has been used to construct a rheological model of the area.The results obtained show a layered rheological structure where an aseismic part of the crust is “sandwiched” between an upper and a lower seismogenic crustal layers. According to the proposed rheological model, seismic activity in the upper crust may be expected down to depths of 10–12 km, which is confirmed well by the observed depth distribution of seismicity. The model also predicts a lower crustal seismogenic layer down to 20–22 km. Because of infrequent occurrences of deep earthquakes and/or a generally small number of reliable hypocenter depth determinations in the study area, this seismogenic zone is less constrained by observations.The depth of the different rheologic horizons within the crust is governed mainly by thermal conditions. The lower boundary of both seismogenic layers appears isothermal. Brittle-ductile transition in the upper crust coincides with the ˜200 °C isotherm, while in the lower crust it coincides with the ˜ 375 °C isotherm. The lowermost crust and the upper mantle beneath Hungary show ductile behavior, thus the possibility of siesmic activity at these horizons can be excluded.     

Ukinrek Maars, Alaska, I. April 1977 eruption sequence, petrology and tectonic setting

During ten days of phreatomagmatic activity in early April 1977, two maars formed 13 km behind the Aleutian arc near Peulik volcano on the Alaska Peninsula. They have been named “Ukinrek Maars”, meaning “two holes in the ground” in Yupik Eskimo. The western maar formed at the northwestern end of a low ridge within the first three days and is up to 170 m in diameter and 35 m in depth. The eastern maar formed during the next seven days 600 m east of West Maar at a lower elevation in a shallow saddle on the same ridge and is more circular, up to 300 m in diameter and 70 m in depth. The maars formed in terrain that was heavily glaciated in Pleistocene times. The groundwater contained in the underlying till and silicic volcanics from nearby Peulik volcano controlled the dominantly phreatomagmatic course of the eruption.During the eruptions, steam and ash clouds reached maximum heights of about 6 km and a thin blanket of fine ash was deposited north and east of the vents up to a distance of at least 160 km. Magma started to pool on the floor of East Maar after four days of intense phreatomagmatic activity.The new melt is a weakly undersaturated alkali olivine basalt (Ne = 1.2%) showing some transitional character toward high-alumina basalts. The chemistry, an anomaly in the tholeitic basalt-andesite-dominated Aleutian arc, suggests that the new melt is primitive, generated at a depth of 80 km or greater by a low degree of partial melting of garnet peridotite mantle with little subsequent fractionization during transport.The Pacific plate subduction zone lies at a depth of 150 km beneath the maars. Their position appears to be tectonically controlled by a major regional fault, the Bruin Bay fault, and its intersection with cross-arc structural features. We favor a model for the emplacement of the Ukinrek Maars that does not link the Ukinrek conduit to the plumbing system of nearby Peulik volcano. The Ukinrek eruptions probably represent a genetically distinct magma pulse originating at asthenospheric depths beneath the continental lithosphere.