Isostatic response of the Laccadive Ridge from admittance analysis of gravity and bathymetry data

The Laccadive Ridge (L-R), trending roughly parallel to the west coast of India, is an intriguing segment of the northernmost Chagos-Laccadive Ridge (C-L-R) system. Although crustal nature and isostatic response of the southern C-L-R is well known, there are no similar studies on the L-R. In the present study, the isostatic response of the lithosphere beneath the L-R is estimated so as to characterize its crustal nature, total crustal as well as effective elastic plate thickness and mode of compensation. Twelve gravity and bathymetry profiles across the ridge were analyzed using linear transfer function and forward model techniques. The observed admittance function within the diagnostic waveband of 250 < λ > 80 km (0.025 < k > 0.080 km⁻¹) fits well with (i) the Airy model whose average crustal thickness (T_c) and density are 17 ± 2 km and 2.7 × 10³ kg m⁻³, respectively, and (ii) the thin plate flexure model of isostasy with an effective elastic plate thickness (T_e) of 2–3 km. The estimated average crustal thickness and density are in good agreement with published seismic refraction results over the ridge. The results of the present study support an Airy model of isostasy for the L-R. The low T_e value, in view of other published results in the study area, suggests stretched and loaded continental lithosphere of the L-R during the evolution of the western continental margin of India.     

The migration history of the Nazca Ridge along the Peruvian active margin: a re-evaluation

The collision zone of the 200 km wide and 1.5 km high Nazca Ridge and the Peruvian segment of the convergent South American margin between 14°S and 17°S is characterized by deformation of the upper plate and several hundred meters of uplift of the forearc. This is evident by a narrowing of the shelf, a westward shift of the coastline and the presence of marine terraces. As the Nazca Ridge is oblique with respect to both trench and convergence direction of the Nazca Plate, it migrates southward along the active plate boundary. For reconstructing the migration history of the Nazca Ridge, this study uses updated plate motion data, resulting from a revision of the geomagnetic time scale. The new model suggests that the ridge crest moved laterally parallel to the margin at a decreasing velocity of ∼75 mm/a (before 10.8 Ma), ∼61 mm/a (10.8-4.9 Ma), and ∼43 mm/a (4.9 Ma to present). Intra-plate deformation associated with mountain building in the Peruvian Andes since the Miocene reduces the relative convergence rate between Nazca Plate and Peruvian forearc. Taking an intra-plate deformation at a rate of ∼10 mm/a, estimated from space-geodetic and geological data, into account, does not significantly reduce these lateral migration velocities. Constraining the length of the original Nazca Ridge by its conjugate feature on the Pacific Plate yields a length of 900 km for the subducted portion of the ridge. Using this constraint, ridge subduction began ∼11.2 Ma ago at 11°S. Therefore, the Nazca Ridge did not affect the northern sites of Ocean Drilling Program (ODP) Leg 112 located at 9°S. This is supported by benthic foraminiferal assemblages in ODP Leg 112 cores, indicating more than 1000 m of subsidence since at least Middle Miocene time, and by continuous shale deposition on the shelf from 18 to 7 Ma, recorded in the Ballena industrial well. At 11.5°S, the model predicts the passage of the ridge crest ∼9.5 Ma ago. This agrees with the sedimentary facies and benthic foraminiferal stratigraphy of ODP Leg 112 cores, which argue for deposition on the shelf in the Middle and Late Miocene with subsequent subsidence of a minimum of several hundred meters. Onshore at 12°S, the sedimentary record shows at least 500 m uplift prior to the end of the Miocene, also in agreement with the model.     

Oblique subduction, collision of microcontinents and subduction of oceanic ridge: Their implications on the Cretaceous tectonics of Japan

Abstract The Izanagi plate subducted rapidly and obliquely under the accretionary terrane of Japan in the Cretaceous before 85 Ma. A chain of microcontinents collided with it at about 140 Ma. In southwest Japan the major part of it subducted thereafter, but in northeast Japan it accreted and the trench jumped oceanward, resulting in a curved volcanic front. The oblique subduction and the underplated microcon-tinent caused uplifting of high-pressure (high-P) metamorphic rocks and large scale crustal shortening in southwest Japan. The oblique subduction caused left-lateral faulting and ductile shearing in northeast Japan. The arc sliver crossed over the high-temperature (high-T) zone of arc magmatism, resulting in a wide high-T metamorphosed belt. At about 85 Ma, the subduction mode changed from oblique to normal and the tectonic mode changed drastically. Just after this the Kula/Pacific ridge subducted and the subduction rate of the Pacific plate decreased gradually, causing the intrusion of huge amounts of granite magma and the eruption of acidic volcanics from large cauldrons. The oblique subduction of the Pacific plate resumed at 53 Ma and the left-lateral faults were reactivated.     

The interaction between regional and local tectonics during resurgent doming: the case of the island of Ischia, Italy

The volcanic island of Ischia is located on the Tyrrhenian margin of Central Italy, characterized by Plio-Quaternary NW–SE- and NE–SW-trending extensional fractures. Ischia displays a resurgent dome uplifted by at least 800 m in the last 33 ka. Remote sensing and field data have been collected to study the structural setting of the island, the deformation pattern associated with resurgence and the superimposition of the regional and the resurgence-induced stress fields. NW–SE and NE–SW extensional fracture systems predominate throughout the island and around the resurgent block, suggesting a relationship with the regional extensional structures. These systems were formed before resurgence and were partly reactivated during resurgence. The reactivation of pre-existing regional systems during resurgence confined the extent of the uplifted area. N–S- and E–W-trending systems have been found exclusively at the borders of the dome and are interpreted as being induced by resurgence. The topmost resurgent block shows an octagonal shape in map view and is tilted at an angle of 15° around a NE–SW-trending horizontal axis; the block is partly bordered by high-angle, inward-dipping regional faults. More than 90% in volume of the volcanic products coeval with resurgence on Ischia have been erupted outside the resurgent block area, suggesting that the resurgence process locally replaced volcanic activity in the last 33 ka.     

南极洲东部普里兹湾海域重磁场特征及地壳结构

普里兹湾位于南极洲东部大陆边缘,其深部地壳结构特征对认识白垩纪冈瓦纳古陆裂解和新生代大陆边缘形成具有重要意义.本文利用重磁、多道反射地震、声纳浮标折射地震和ODP钻井数据对普里兹湾海域的深部地壳结构进行了研究.研究结果显示,普里兹凹陷表现为典型的盆地负重力异常特征,其沉积基底较深,而在四夫人浅滩为高幅重力正异常,其沉积基底普遍抬升.在大陆架中部存在SW-NE向条带状基底的抬升,且呈朝NE向逐渐变深的趋势.在中大陆架外侧,均衡残余重力异常呈V字形负异常条带状分布,其两翼分别与四夫人浅滩和弗拉姆浅滩外的大陆坡相连.该异常带在大陆架中部向陆的偏移可能是由于古大陆架边缘的地形影响,推测其与普里兹冲积扇同属于洋陆过渡带向陆的部分,在重力模拟剖面表现为地壳向海逐渐减薄.普里兹冲积扇的地壳厚度较薄,平均为6 km,最薄处可达4.6 km,并且根据洋陆过渡带向海端的位置,推测可能属于接近洋壳厚度的过渡壳.重力异常分区的走向与兰伯特地堑在普里兹湾的构造走向基本一致,可能主要反映了二叠纪—三叠纪超级地幔柱对普里兹湾的裂谷作用的影响.该区域的自由空间重力异常和均衡残余异常均表现为超过100×10-5m/s2的高幅正异常特征,可能由位于大陆架边缘的巨厚沉积体负载在高强度岩石圈之上的区域挠曲均衡作用所导致,可能与该区域第二期裂谷期之后的沉积间断以及快速进积加厚的演化过程有关.普里兹湾磁力异常的走向与重力异常明显不同,大致可分为东北高幅正异常区和西南低幅异常区.重磁异常在走向上的差异反映高磁异常主要来源于岩浆作用形成的铁镁质火成岩的影响,并且岩浆作用的时代不同于基底隆升的时代,而可能形成于前寒武纪或者南极洲和印度板块裂谷期间(白垩纪).     

Morphology and magnetic survey of the Rivera-Cocos plate boundary of Colima,Mexico

The propagation of the Pacific-Cocos Segment of the East Pacific Rise (EPR-PCS) has significantly altered the plate configuration at the north end of the Middle America Trench. This ridge propagation, the collision of the EPR-PCS with the Middle America Trench, the separation of the Rivera and Cocos plates and the formation of the Rivera Transform have produced a complex arrangement of morphotectonic elements in the area of Rivera-Cocos plate boundary, atypical of an oceanic transform boundary. Existing marine magnetic and bathymetric data has proved inadequate to unravel this complexity, thus, a dense grid of total field magnetic data were collected during campaigns MARTIC-04 and MARTIC-05 of the B/O EL PUMA in 2004 and 2006. These data have greatly clarified the magnetic lineation pattern adjacent to the Middle America trench, and have revealed an interesting en echelon, NE-SW oriented magnetic high offshore of the Manzanillo Graben. We interpret these new data to indicate that the EPR-PCS ridge segment reached the latitude (~18.3°N) of the present day Rivera Transform at about Chron 2A₃ (~3.5Ma) and propagated further northward, intersecting the Middle America Trench at about 1.7 Ma (Chron 2). At 1.5 Ma spreading ceased along the EPR north of 18.3°N and the EPR-PCS has since retreated southward in association with a southward propagation of the Moctezuma Spreading Segment. North of 18.3°N the seafloor near the trench has been broken into small, uplifted blocks, perhaps due to the subduction of the young lithosphere generated by the EPR-PCS.     

南海深部构造特征及其地质意义:来自重磁位场反演的认识

莫霍面和居里面是认识深部过程重要的地质与地球物理界面.为了进一步理解南海深部构造活动与洋盆扩张的关系,本文以OBS剖面和深反射地震剖面作为约束,对卫星测高重力异常进行海水、沉积层影响校正,采用最小曲率位场分离方法消除局部密度体的重力影响,获取了反映莫霍面起伏的重力异常,并利用双界面模型重力场快速反演方法计算得到了南海地区莫霍面深度值.通过与居里面起伏的对比研究,发现南海莫霍面和居里面整体均表现为"洋盆浅、周缘深"的菱形特征,两者在洋陆转换区呈现明显的窄梯级带特征,反映了南海扩张期岩石圈的强烈伸展减薄、南北向构造拉张作用等深部构造过程.洋盆莫霍面和居里面的西南向楔形形态是对南海由东向西渐进式扩张的深部构造响应.洋盆南部莫霍面浅于北部,这与扩张中心逐渐向南迁移的特征一致,而洋盆居里面南深北浅的特征则可能与洋盆的简单剪切扩张方式以及洋盆北部的岩浆活动更活跃有关.南海地区莫霍面和居里面呈现交错叠置关系,南、北陆缘表现为明显的深部构造差异,说明南海为非对称式扩张.北部陆缘区居里面深度浅于莫霍面,而洋盆区和南部陆缘区居里面深于莫霍面,这与南、北陆缘性质的差异和南部陆缘复杂的中-新生代俯冲碰撞等构造演化相关,而洋盆区居里面深于莫霍面的现象推测与大洋上地幔橄榄岩蛇纹石化导致的岩石磁性增强有关.     

Fast mapping of magnetic basement depth,structure and nature using aeromagnetic and gravity data: combined methods and their application in the Paris Basin

Assessment of deep buried basin/basement relationships using geophysical data is a challenge for the energy and mining industries as well as for geothermal or CO₂ storage purposes. In deep environments, few methods can provide geological information; magnetic and gravity data remain among the most informative and cost‐effective methods. Here, in order to derive fast first‐order information on the basement/basin interface, we propose a combination of existing and original approaches devoted to potential field data analysis. Namely, we investigate the geometry (i.e., depth and structure) and the nature of a deep buried basement through a case study SW of the Paris Basin. Joint processing of new high‐resolution magnetic data and up‐to‐date gravity data provides an updated overview of the deep basin. First, the main structures of the magnetic basement are highlighted using Euler deconvolution and are interpreted in a structural sketch map. The new high‐resolution aeromagnetic map actually offers a continuous view of regional basement structures and reveals poorly known and complex deformation at the junction between major domains of the Variscan collision belt. Second, Werner deconvolution and an ad hoc post‐processing analysis allow the extraction of a set of magnetic sources at (or close to) the basin/basement interface. Interpolation of these sources together with the magnetic structural sketch provides a Werner magnetic basement map displaying realistic 3D patterns and basement depths consistent with data available in deep petroleum boreholes. The last step of processing was designed as a way to quickly combine gravity and magnetic information and to simply visualize first‐order petrophysical patterns of the basement lithology. This is achieved through unsupervised classification of suitably selected gravity and magnetic maps and, as compared to previous work, provides a realistic and updated overview of the cartographic distribution of density/magnetization of basement rocks. Altogether, the three steps of processing proposed in this paper quickly provide relevant information on a deep buried basement in terms of structure, geometry and nature (through petrophysics). Notwithstanding, limitations of the proposed procedure are raised: in the case of the Paris Basin for instance, this study does not provide proper information on Pre‐Mesozoic basins, some of which have been sampled in deep boreholes.     

The behavior of magneto-acoustic-gravity waves near the cusp resonance in a lossless,compressible, isothermal,stratified, electrically conducting and uniformly magnetized atmosphere. II: The valve effect

Abstract

In the present paper the behavior of (internal) magneto-acoustic-gravity waves near the cusp resonance in a lossless, compressible, isothermal, stratified, electrically conducting atmosphere that is permeated by a uniform, nearly horizontal magnetic field is re-addressed (Kamp, 1989). The previously analyzed linear conversion of long acoustic-gravity waves into short magneto-acoustic waves that carry off the energy from the resonance region along the magnetic field, is re-analyzed with boundary layer techniques that are based on the smallness of the vertical component of the magnetic field. More specifically the existence of the so-called valve effect for the generated magneto-acoustic mode near the critical level is explicitly demonstrated and shown to be governed by two rivalling effects.     

The Whitsunday Volcanic Province, Central Queensland, Australia: lithological and stratigraphic investigations of a silicic-dominated large igneous province

Contrary to general belief, not all large igneous provinces (LIPs) are characterised by rocks of basaltic composition. Silicic-dominated LIPs, such as the Whitsunday Volcanic Province of NE Australia, are being increasingly recognised in the rock record. These silicic LIPs are consistent in being: (1) volumetrically dominated by ignimbrite; (2) active over prolonged periods (40–50 m.y.), based on available age data; and (3) spatially and temporally associated with plate break-up. This silicic-dominated LIP, related to the break-up of eastern continental Gondwana, is also significant for being the source of >1.4×10⁶ km³ of coeval volcanogenic sediment preserved in adjacent sedimentary basins of eastern Australia.The Whitsunday Volcanic Province is volumetrically dominated by medium- to high-grade, dacitic to rhyolitic lithic ignimbrites. Individual ignimbrite units are commonly between 10 and 100 m thick, and the ignimbrite-dominated sequences exceed 1 km in thickness. Coarse lithic lag breccias containing clasts up to 6 m diameter are associated with the ignimbrites in proximal sections. Pyroclastic surge and fallout deposits, subordinate basaltic to rhyolitic lavas, phreatomagmatic deposits, and locally significant thicknesses of coarse-grained volcanogenic conglomerate and sandstone are interbedded with the ignimbrites. The volcanic sequences are intruded by gabbro/dolerite to rhyolite dykes (up to 50 m in width), sills and comagmatic granite. Dyke orientations are primarily from NW to NNE.The volcanic sequences are characterised by the interstratification of proximal/near-vent lithofacies such as rhyolite domes and lavas, and basaltic agglomerate, with medial to distal facies of ignimbrite. The burial of these near-vent lithofacies by ignimbrites, coupled with the paucity of mass wastage products such as debris-flow deposits indicates a low-relief depositional environment. Furthermore, the volcanic succession records a temporal change in: (1) eruptive styles; (2) the nature of source vents; and (3) erupted compositions. An early explosive dacitic pyroclastic phase was succeeded by a later mixed pyroclastic-effusive phase producing an essentially bimodal suite of lavas and rhyolitic ignimbrite. From the nature and distribution of volcanic lithofacies, the volcanic sequences are interpreted to record the evolution of a multiple vent, low-relief volcanic region, dominated by several large caldera centres.     

A new ellipsoidal gravimetric, satellite altimetry and astronomic boundary value problem, a case study: The geoid of Iran

A new gravimetric, satellite altimetry, astronomical ellipsoidal boundary value problem for geoid computations has been developed and successfully tested. This boundary value problem has been constructed for gravity observables of the type (i) gravity potential, (ii) gravity intensity (i.e. modulus of gravity acceleration), (iii) astronomical longitude, (iv) astronomical latitude and (v) satellite altimetry observations. The ellipsoidal coordinates of the observation points have been considered as known quantities in the set-up of the problem in the light of availability of GPS coordinates. The developed boundary value problem is ellipsoidal by nature and as such takes advantage of high precision GPS observations in the set-up. The algorithmic steps of the solution of the boundary value problem are as follows:

- Application of the ellipsoidal harmonic expansion complete up to degree and order 360 and of the ellipsoidal centrifugal field for the removal of the effect of global gravity and the isostasy field from the gravity intensity and the astronomical observations at the surface of the Earth.

- Application of the ellipsoidal Newton integral on the multi-cylindrical equal-area map projection surface for the removal from the gravity intensity and the astronomical observations at the surface of the Earth the effect of the residual masses at the radius of up to 55 km from the computational point.

- Application of the ellipsoidal harmonic expansion complete up to degree and order 360 and ellipsoidal centrifugal field for the removal from the geoidal undulations derived from satellite altimetry the effect of the global gravity and isostasy on the geoidal undulations.

- Application of the ellipsoidal Newton integral on the multi-cylindrical equal-area map projection surface for the removal from the geoidal undulations derived from satellite altimetry the effect of the water masses outside the reference ellipsoid within a radius of 55 km around the computational point.

- Least squares solution of the observation equations of the incremental quantities derived from aforementioned steps in order to obtain the incremental gravity potential at the surface of the reference ellipsoid.

- The removed effects at the application points are restored on the surface of reference ellipsoid.

- Application of the ellipsoidal Bruns’ formula for converting the potential values on the surface of the reference ellipsoid into the geoidal heights with respect to the reference ellipsoid.

- Computation of the geoid of Iran has successfully tested this new methodology.

The internal structure of lava flows—insights from AMS measurements II: Hawaiian pahoehoe, toothpaste lava and 'a'

We studied the anisotropy of magnetic susceptibility (AMS) of 22 basaltic flow units, including S-type pahoehoe, P-type pahoehoe, toothpaste lava and 'a' emplaced over different slopes in two Hawaiian islands. Systematic differences occur in several aspects of AMS (mean susceptibility, degree of anisotropy, magnetic fabric and orientation of the principal susceptibilities) among the morphological types that can be related to different modes of lava emplacement. AMS also detects systematic changes in the rate of shear with position in a unit, allowing us to infer local flow direction and some other aspects of the velocity field of each unit. 'A' flows are subject to stronger deformation than pahoehoe, and also their internal parts behave more like a unit. According to AMS, the central part of pahoehoe commonly reveals a different deformation history than the upper and lower extremes, probably resulting from endogenous growth.     

The transition from Alboran to Algerian basins (Western Mediterranean Sea): Chronostratigraphy,deep crustal structure and tectonic evolution at the rear of a narrow slab rollback system

The eastern Alboran basin and its transition to the Algerian basin is a key area in the Mediterranean realm where controversial kinematic and geodynamical models are proposed. Models imply striking differences regarding the nature of the crust, the prevalence of brittle faulting and ductile shear, the origin of magmatism, the style of Miocene deformation and the driving mechanisms of the Alboran plate kinematics. Combining a new chronostratigraphic chart of the Alboran and Algerian basins based on the Habibas (HBB-1) core drill, deep seismic sections striking WSW-ENE and SSE-NNW, and potential field data, we re-assess the tectonic evolution that controlled the sedimentation and basement deformation of the westernmost limit of the Algerian basin and its transition with the Alboran domain. A WSW-directed extensional tectonic phase has shaped a stretched continental crust with typical tilted blocks along ∼100 km from Burdigalian to Tortonian times, which is assumed to result from the WSW-directed migration of the Alboran block driven by a narrow slab rollback. In the Algerian basin, this event was followed by the emplacement of an oceanic-type crust. Potential field signatures of the deep basin as well as geometrical correlations with onland outcrops of inner zones suggest a minimum WSW-directed displacement of the Alboran terrane of ∼200 km. At the southern foot of the Algerian basin, the continent-ocean transition is sharp and may result from the westward propagation of a slab tear at depth, forming two segments of STEP (Subduction-Transform Edge Propagator) margins. Our results support models of intense shear tractions at the base of an overriding plate governed by slab rollback-induced mantle flow. Finally, Messinian salt tectonics affected overlying deposits until today. A late Tortonian to Quaternary dominantly transpressive tectonic episode linked to the Africa-Iberia convergence post-dates previous events, deforming the whole margin.     

The last 40 ka tephrostratigraphic record of Lake Ohrid,Albania and Macedonia: a very distal archive for ash dispersal from Italian volcanoes

A 1075 cm long core (Lz1120) was recovered in the south-eastern part of the Lake Ohrid (Republics of Macedonia and Albania) and sampled for identification of tephra layers. Magnetic susceptibility investigations show rather high magnetic values throughout the core, with peaks unrelated to the occurrence of tephra layers but instead to the relative abundance of detrital magnetic minerals in the sediment. Naked-eye inspection of the core allowed us to identify of two tephra layers, at 896–897 cm and 1070–1075 cm. Laboratory inspection of the grain-size fraction > 125 μm allowed for the identification of a third cryptotephra at 310–315 cm. Major element analyses on glass shards of the tephra layers at 896–897 cm and 1070–1075 cm show a trachytic composition, and indicate a correlation with the regionally dispersed Y-3 and Y-5 tephra layers, dated at ca 30 and 39 cal ka BP. The cryptotephra at 310–315 cm has a mugearitic–benmoreitic composition, and was correlated with the FL eruption of Mt. Etna, dated at 3370 ± 70 cal yr BP. These ages are in agreement with five ¹⁴C AMS measurements carried out on plant remains and macrofossils from the lake sediments at different depths along the core.     

The phosphoric acid leak from the wreck of the MV Ece in the English Channel in 2006: Assessment with a ship of opportunity, an operational ecosystem model and historical data

This study evaluates the ship of opportunity (Ferrybox) concept for both sustained monitoring of UK shelf sea waters and numerical model validation. Release of phosphate from the wreck of a chemical tanker (MV Ece) in the western English Channel (49.73 degrees N, 3.25 degrees W) in March 2006 is used to demonstrate the importance of sustained observations in decision support systems and policy development. The Ferrybox system continuously collects sea surface (5m) data from a suite of autonomous electronic sensors installed on a passenger ferry operating year-round between Portsmouth (UK) and Bilbao (Spain). The detection of anomalously high concentrations of phosphate (1.54mmolm(-3), four times the usual level) and onset of phytoplankton growth close to the wreck site in March 2006 was placed in the context of multiple years of measurements (phosphate, nitrate, silicate and chlorophyll) collected from the Ferrybox system (2003-2006) and the long-term time series station E1 (50.03 degrees N, 4.65 degrees W, 1930-1987) in the English Channel. With regard to decision support, release of phosphate from the tanker is unlikely to pose a threat as phytoplankton growth at the end of winter is not unusual in this region and dissolved inorganic nitrogen rather than phosphate (DIN:DIP=10-18) is likely to ultimately limit algal growth in spring 2006. With regard to policy development, the Oslo and Paris (OSPAR) commissions recommendation of sampling every three years in "non-problem areas" is likely to provide statistically inadequate data, given the interannual and decadal variability identified in the Ferrybox and E1 data: the Ferrybox data show that oceanic winter nutrient concentrations varied by 35-50% between 2003/2004 and 2005/2006 due to deeper mixing of water off-shelf in early 2005/2006 and comparisons between the Ferrybox and E1 years show that the western English Channel is currently experiencing a low in phosphate concentrations similar to those in the 1960s. The importance of Ferrybox data in evaluating the reliability of predictive operational models needed in decision support is also demonstrated, by highlighting both strengths and weaknesses in a state-of-the-art ecosystem model designed for UK shelf waters.