导体边缘附近虚感应矢量性质的数值模型研究

本文用E极化的二维有限元法和有限差分法研究导体层边缘附近的虚感应矢量.首先,研究有一定埋深的导体--围岩界面处的情况,其中包括虚感应矢量的频率响应.其次,研究了海岸效应的虚感应矢量的频率响应.在海岸效应中不仅考虑了不同海水层厚度的影响,还考虑了海水与陆地的电导率比值大小对虚感应矢量的影响.最后,研究了被动陆缘（passive continental margin）处的海岸效应,检查磁感应矢量的测量是否能揭示海岸下面深部电导率的横向变化（即与岩石层厚度变化有关联的电导率变化）.     

海洋板块俯冲作用下上覆大陆岩石层减薄机制的动力学模拟

几乎所有大陆岩石层的减薄现象,可能都与海洋板块的俯冲作用相关,但是两者之间的内在联系迄今仍不十分明确,为此,我们设计了一系列包含洋-陆俯冲系统的二维数值模型,来探讨海洋板块的俯冲作用对上覆大陆岩石层变形行为的影响,尤其对大陆岩石层减薄效应的制约.模型结果表明,海洋板块俯冲过程中的地幔楔熔体对大陆岩石层地幔的热侵蚀以及由熔体上升所诱发的地幔局部对流的强烈扰动会导致上覆大陆岩石层的减薄效应.这种效应不仅表现在横向上的向陆内蔓延,还表现在垂向上的向浅部发展.且多类动力学参数都能制约大陆岩石层的减薄效应.具体地,随着汇聚速率和洋壳厚度的增加,上覆大陆岩石层在横向上的减薄范围越大,在垂向上的减薄程度也越深;而随着俯冲海洋板块年龄的增加,上覆大陆岩石层在横向上的减薄范围增大,但在垂向上的减薄程度会减小;随着上覆大陆岩石层厚度的增加,其横向减薄范围会减小,但在垂向上的减薄程度会加深.本文研究成果能为揭示华北克拉通减薄/破坏的动力学过程提供一定的理论参考依据.     

Electromagnetic evidence for lateral inhomogeneities within the Earth's upper mantle

The composition of the upper mantle is of great significance to our understanding of plate tectonics and global evolution. Information about the physical properties of the Earth at upper mantle depths, including lateral variations in electrical conductivity, can be deduced from measurements of the electric and magnetic fields at the Earth's surface. Electromagnetic methods appear to give poorer resolution than do some other methods, for example seismics, but as they are sensitive to quite different properties of a medium they provide a different and complementary class of information.The basic theory of electromagnetic sounding methods is briefly reviewed below, and evidence regarding lateral conductivity inhomogeneities in the Earth's upper mantle is examined. While lateral electrical conductivity inhomogeneities appear to be the rule rather than the exception, the interpretation of electromagnetic data still presents difficulties and the results from many regions are not as yet unambiguous. Where the data are of sufficient resolution, a rapid increase in electrical conductivity can usually be identified within the upper mantle. The depth to this highly conductive zone is different in different tectonic environments, but is broadly consistent between analogous but widely separated tectonic environments. A comparatively shallow conducting region is found beneath the ocean lithosphere. The depth of this region is dependent on lithospheric age. Many of the more shallow conducting regions in both continental and oceanic environments are associated with high heat flow values and seismic low velocity zones. These highly conducting regions may be zones of partial melt.     

南海岩石层及边界构造的地球物理特征

南海经历了中生代主动大陆边缘到新生代被动大陆边缘的转换,其岩石层地球物理场具有明显的块、带特征.本文通过综合分析南海地区深地震探测、面波层析成像、重磁异常以及地热与岩石层流变学等各种地质地球物理资料,对南海地壳及岩石层的综合地球物理特征进行了深入总结,发现深地震探测剖面所确定的洋、陆壳转换位置与空间重力异常梯级带位置较为一致,据此拟定了南海洋、陆壳的转换边界;依据多条地壳结构剖面中拉张减薄的程度确定了正常减薄陆壳、洋陆壳过渡带及洋壳等属性特征,并初步圈定了南海下地壳高速层的分布范围.对比分析了南、北陆缘地壳结构及其拉张减薄的变化特征,从综合地球物理特征的相似性上推测了北部陆缘的中西沙陆块与南部陆缘的南沙礼乐滩陆块具有共轭对称性.依据S波速度梯度变化确定了南海岩石层厚度分布情况,揭示出南海北部陆缘存在一条岩石层厚度的减薄带,且该减薄带与高热流带具有较好的一致性.在综合分析的基础上,以深地震探测剖面与重、磁异常变化的对应性为基础,划定了南海边界构造的位置.     

Lithospheric model with thick oceanic crust at the continental boundary: A mechanism for shallow spreading ridges in young oceans

In a general lithospheric model of a simple divergent ocean and continental margin that satisfies the constraints of isostasy and gravity anomalies, the free-air gravity anomaly at the margin is modelled by an oceanic crust that thickens exponentially toward the margin from its common value of 6.4 km about 600 km from the margin to 17.7 km at the margin; this postulated thickening is supported empirically by seismic refraction measurements made near continental margins. The thickness of the oceanic crust matches that of the continental lithosphere at breakup, as observed today in Afar and East Africa, and is interpreted as the initial oceanic surface layer chilled against the continental lithosphere. With continued plate accretion, the chilled oceanic crust thins exponentially to a steadystate thickness, which is achieved about 40 m.y. after breakup. These findings contrast with the generally held view that the oceanic crust has a uniform thickness.During the first 40 m.y. of spreading, the thicker oceanic crust, of density 2.86 g/cm³, displaces the denser (3.32 g/cm³) subjacent material; by isostasy, the spreading ridge and the rest of the seafloor thus stand higher in younger( <40m.y.) oceans than they do in older(>40m.y.) oceans. This is postulated to be the cause of the empirical relationship between the crestal depth of spreading ridges and the age (or half-width) of ocean basins.     

地幔热导率的选取对动力学数值模拟的影响——以岩石圈张裂过程为例

目前存在有多种地幔热导率模型,不同模型在数值和随温压变化的特征上有明显的差异.为探究不同热导率模型对动力学数值模拟结果的影响,本文对不同模型下的岩石圈张裂过程进行模拟研究,探讨地幔热导率对岩石圈热传输、变形和熔融过程的影响及其作用机理.结果显示,不同热导率模型下,岩石圈的变形和熔融特征表现出明显差异.高热导率模型下,岩石圈破裂较晚,形成陆缘较为宽阔,地壳熔融强烈而地幔熔融较弱;低热导率模型下,岩石圈破裂较早,形成陆缘较为狭窄,地幔熔融强烈而地壳熔融较弱.这种差异源于不同地幔热导率下岩石圈和地幔热状态的变化及相应力学性质的改变.高热导率下,热传导的增温效应显著,岩石圈呈现较热的状态,其强度整体较低,壳幔耦合减弱;而低热导率下,热对流的增温效应显著,岩石圈呈较冷的状态,其强度整体较高,壳幔耦合增强.基于模拟结果,本文认为地幔热导率的选取对动力学模拟的结果有着较为显著的影响,相对于随温压的变化,热导率数值的差异对动力学数值模拟的结果影响更大,尤其是对于地幔熔融过程的影响.     

Experimental models of extension of continental lithosphere weakened by percolation of asthenospheric melts

Analogue models are used to investigate extension of a continental lithosphere weakened by asthenospheric melts percolating through the upper mantle, a process that has been hypothesised to control the opening of the Ligurian Tethys. Models were performed in a centrifuge apparatus and reproduced, by using materials such as sand and viscous mixtures, extension of 60-km thick, three-layer continental lithosphere floating above the asthenosphere. The percolated lithospheric mantle was assumed to be characterised by a rheological behaviour similar to that of the asthenosphere. Two sets of experiments investigated the influence on deformation of (1) the thickness of the percolated mantle and the associated strength contrast between the normal and weakened lithosphere, and (2) the lateral width of the weakened zone. Model results suggest that mantle percolation by asthenospheric melts is able to promote strong localised thinning of the continental lithosphere, provided that a significant thickness of the lithospheric mantle is weakened by migrating melts within a narrow region. Strain localisation is maximised for percolation of the whole lithospheric mantle and strong strength contrast between the normal and weakened lithosphere. Under these conditions, the thickness of the lithosphere may be reduced to less than 12 km in 3 Ma of extension. Conversely, localised thinning is strongly reduced if the thickness of the percolated zone is ≤1/3 of the thickness of the whole lithospheric mantle and/or the lithosphere is weakened over wide regions. Overall, model results support the working hypothesis that mantle percolation by asthenospheric melts is a controlling factor in the transition from distributed continental deformation to localised oceanic spreading.     

Thermal anomalies and magmatism due to lithospheric doubling and shifting

We present some thermal and magmatic consequences of the processes of lithospheric doubling and lithospheric shifting.Lithospheric doubling concerns the obduction of a cold continental or old oceanic lithospheric plate over a young and hot oceanic lithosphere/upper mantle system, including an oceanic ridge.Lithospheric shifting concerns the translation and rotation of a lithospheric plate relative to the upper mantle.In both cases the resulting thermal state of the upper mantle below the obducting or shifting lithosphere may be perturbed relative to a “normal” continental or oceanic geothermal situation.The perturbed geothermal state gives rise to a density inversion at depth and thus induces a vertical gravitational instability which favours magmatism.We speculate about the magmatic consequences of this situation and infer that in the case of lithospheric doubling our model may account for the petrology and geochemistry of the resulting magma.The original layering and composition of the overridden young oceanic lithosphere may strongly influence magmatic processes.We dwell shortly on the genesis of kimberlites within the framework of our lithospheric doubling model and on magmatism in general. Lithospheric recycling is inherent to the mechanism of lithospheric doubling.     

Marine electromagnetic research

All aspects of Marine Electromagnetic Research have made important advances over the last few years: theoretical studies, instrument design and data from equipment on the bottom of the ocean. The seafloor results show that the depth to the conducting asthenosphere is greater under older lithosphere and thus the thickness of the lithosphere increases with age. To obtain greater resolution of the electrical conductivity structure of the upper layers, several controlled source systems were developed. The first observations indicate a decrease in conductivity a few kilometers below the seafloor. Improved theoretical response curves for the electromagnetic fields at an ocean-continent boundary are now available. The theoretical curves, combined with land and seafloor data from coastal regions, allow the effects of the electric currents induced in the seawater to be separated from those caused by currents in the tectonic structure at the continental margin. Of growing interest is the application of electromagnetic methods to determine oceanographic wave parameters. Recent studies have investigated meanders in ocean current patterns, the response of internal edge waves and the tidal effect. With the forseeable improvements in seafloor instrumentation it will be possible to investigate the conductivity structure of offshore basins and the hydrothermal deposits associated with spreading ridges.     

Electromagnetic study of lithospheric structure in the marginal zone of East European Craton in NW Poland

The marginal zone of the East European Platform, an area of key importance for our understanding of the geotectonic history of Europe, has been a challenge for geophysicists for many years. The basic research method is seismic survey, but many important data on physical properties and structure of the lithosphere may also be provided by the electromagnetic methods. In this paper, results of deep basement study by electromagnetic methods performed in Poland since the mid-1960s are presented. Over this time, several hundred long-period soundings have been executed providing an assessment of the electric conductivity distribution in the crust and upper mantle. Numerous 1D, 2D, and pseudo-3D electric conductivity models were constructed, and a new interpretation method based on Horizontal Magnetic Tensor analysis has been applied recently. The results show that the contact zone is of lithospheric discontinuity character and there are distinct differences in geoelectric structures between the Precambrian Platform, transitional zone (TESZ), and the Paleozoic Platform. The wide-spread conducting complexes in the crust with integral conductivity values reaching 10 000 S at 20–30 km depths are most spectacular. They are most likely consequences of geological processes related to Caledonian and Variscan orogenesis. The upper mantle conductivity is also variable, the thickness of high-resistive lithospheric plates ranging from 120–140 km under the Paleozoic Platform to 220–240 km under the East European Platform.     

Softening of the subcontinental lithospheric mantle by asthenosphere melts and the continental extension/oceanic spreading transition

The majority of ophiolitic peridotites in the Alpine–Apennine system show evidence of extensive interaction between subcontinental lithospheric mantle and fractional melts of asthenospheric origin. This interaction led to petrological, structural, and geochemical changes in the lithospheric mantle, and was accompanied by a temperature increase to near-asthenospheric values, resulting in the thermomechanical erosion of the lithosphere. We term the parts of mantle lithosphere thus affected the asthenospherized lithospheric mantle or ALM.The thermal and rheological consequences of thermomechanical erosion are explored by modelling the temperature and rheological properties of the thinned lithosphere as a function of thickness of ALM and time since asthenospherization (i.e., since the beginning of thermal relaxation). Results are given both in terms of rheological profiles (strength envelopes) and total lithospheric strength (TLS) for different lower crustal rheologies. The TLS decreases as a consequence of thermomechanical erosion. This decrease is a non-linear function of the thickness of the ALM. While practically negligible if less than 50% of lithospheric mantle is affected, it becomes significant (up to almost one order of magnitude) if thermomechanical erosion approaches the Moho. The maximum decrease in TLS is achieved within a short time span (1–2 Ma) after the end of the heating episode.As a working hypothesis, we propose that thermomechanical erosion of the lithospheric mantle, related to lithosphere/asthenospheric melts interaction, can be an important factor in a geologically rapid decrease in TLS. This softening could lead to whole lithospheric failure and consequently to a transition from continental extension to oceanic spreading.     

南海东北部及其邻近地区的Pn波速度结构与各向异性

利用中国地震台网和ISC台站1980～2004年的地震数据,反演了南海东北部及其邻近地区的Pn波速度结构和各向异性.上地幔顶部的速度变化揭示出区域地质构造的深部特征：华南地区速度较高并且变化平缓,具有构造稳定地区的岩石层地幔特征;华南沿海尤其是滨海断裂带附近出现低速异常,表明该断裂可能穿过壳幔边界深达上地幔顶部.南海北部至台湾海峡较高的速度与华南地区类似,反映出大陆边缘和陆架地区的岩石层地幔性质;西沙海槽附近较高的速度不仅反映了华南大陆向南的延伸,而且与海槽裂谷拉张引起的地幔上拱有关,整个南海北部没有发现大规模地幔热流的活动痕迹.相比之下,南海东部次海盆的上地幔顶部存在明显的低速异常,对应于海底扩张中心的地幔上涌区,表明岩石层地幔强烈减薄甚至缺失;台湾东部－吕宋－菲律宾北部的低速异常与地震、火山活动以及岩浆作用紧密相关,揭示了西太平洋岛弧俯冲带的活动特征;南海东北部的洋－陆边界清晰,南海东部和菲律宾海西部较高的速度代表了海洋岩石层地幔的性质.Pn波各向异性反映出区域性构造应力状态及岩石层地幔的变形痕迹：华南地区的各向异性较小,说明这一构造稳定地区的岩石层地幔变形程度较弱;南海北部的快波方向与地壳浅表层构造的伸展方向一致,主要反映了中、新生代以来的大陆边缘张裂和剪切作用对岩石层地幔结构的影响;琉球－台湾－吕宋岛弧两侧各向异性十分强烈,平行于海沟的快波方向表明菲律宾海板块和欧亚大陆的相互作用导致俯冲板块前缘的岩石层地幔强烈变形;台湾东南海域快波方向的变化可能与欧亚大陆和菲律宾海板块俯冲机制的转换以及岩石层被撕裂有关.     

Dynamics of thinning and destruction of the continental cratonic lithosphere: Numerical modeling

Thinning of the cratonic lithosphere is common in nature, but its destruction is not. In either case, the mechanisms for both thinning and destruction are still widely under debate. In this study, we have made a review on the processes and mechanisms of thinning and destruction of cratonic lithosphere according to previous studies of geological/geophysical observations and numerical simulations, with specific application to the North China Craton (NCC). Two main models are suggested for the thinning and destruction of the NCC, both of which are related to subduction of the oceanic lithosphere. One is the “bottom-up” model, in which the deeply subducting slab perturbs and induces upwelling from the hydrous mantle transition zone (MTZ). The upwelling produces mantle convection and erodes the bottom of the overriding lithosphere by the fluid-melt-peridotite reaction. Mineral compositions and rheological properties of the overriding lithospheric mantle are changed, allowing downward dripping of lithospheric components into the asthenosphere. Consequently, lithospheric thinning or even destruction occurs. The other is the “top-down” model, characterized by the flat subduction of oceanic slab beneath the overriding cratonic lithosphere. Dehydration reactions from the subducting slab would significantly hydrate the lithospheric mantle and decrease its rheological strength. Then the subduction angle may be changed from shallow to steep, inducing lateral upwelling of the asthenosphere. This upwelling would heat and weaken the overriding lithospheric mantle, which led to the weakened lithospheric mantle dripping into the asthenosphere. These two models have some similarities, in that both take the subducting oceanic slab and relevant fluid migration as the major driving mechanism for thinning or destruction of the overriding cratonic lithosphere. The key difference between the two models is the effective depth of the subducting oceanic slab. One is stagnation and flattening in the MTZ, whereas the other is flat subduction at the bottom of the cratonic lithosphere. In the NCC, the eastern lithosphere was likely affected by subduction of the Izanagi slab during the Mesozoic, which would have perturbed the asthenosphere and the MTZ, and induced fluid migration beneath the NCC lithosphere. The upwelling fluid may largely have controlled the reworking of the NCC lithosphere. In order to discuss and analyze these two models further, it is crucial to understand the role of fluids in the subduction zone and the MTZ. Here, we systematically discuss phase transformations of hydrous minerals and the transport processes of water in the subduction system. Furthermore, we analyze possible modes of fluid activity and the problems to explore the applied feasibility of each model. In order to achieve a comprehensive understanding of the mechanisms for thinning and destruction of cratonic lithosphere, we also consider four additional possible dynamic models: extension-induced lithospheric thinning, compression-induced lithospheric thickening and delamination, large-scale mantle convection and thermal erosion, and mantle plume erosion. Compared to the subduction-related models presented here, these four models are primarily controlled by the relatively simple and single process and mechanism (extension, compression, convection, and mantle plume, respectively), which could be the secondary driving mechanisms for the thinning and destruction of lithosphere.     

基于变参数板块模型的大洋岩石圈参数反演与岩石圈伸展成盆模拟

本文利用热导率、热容、热膨胀系数等参数随温度变化的经验表达式,在板块模型的基础上用隐式有限差分方法解非线性热传导方程,并利用北太平洋和北大西洋海底年龄与水深数据反演了大洋岩石圈厚度与底界温度等参数,结果表明大洋岩石圈的厚度在105 km左右,岩石圈底界温度在1450℃左右,这与Stein等用全球大量数据反演的结果一致.将变参数模型用到岩石圈拉张成盆的模拟中,结果表明当考虑岩石圈热参数随温度变化之后,预测的地表热沉降要大于常参数均匀伸展模型的预测量.由此我们指出：McKenzie的均匀伸展模型预测的初始沉降偏大而热沉降偏小,可能与该模型没有考虑热参数随深度(即温度)变化有关.     

Tectonic interpretation of regional conductivity anomalies

This review paper selects key results from electromagnetic induction studies of a variety of distinctive tectonic phenomena in the top 200 km of the Earth. Its main theme is that electromagnetic data are essential for an understanding of tectonism involving partial melting, recycling of large volumes of fluids (CO₂ and H₂O) and underthrusting of metasedimentary rocks. The wide variety of tectonic regimes in which these processes are known to be important is reflected in the choice of case studies. A discussion of conductivity models for young oceanic lithosphere and asthenosphere is followed by results from induction studies across the S.E. Australian passive margin, the North American active margin, the Ryukyu Island-Arc and the Oregon Cascades continental arc. The importance of partial melting and free fluid movement i apparent in these regions. Terrain accretion and/or continent-continent collisions recorded at palaeosuture zones in Ireland, Germany and Scandinavia have left distinctive conductivity structures. These are often associated with grain-boundary graphite either in weakly-metamorphosed black shales in underthrust sedimentary basins or precipitated from CO₂-rich fluids. They are discussed in the context of the evolution of mature continental crust. All of the case studies are based on experiments published since 1989 in which the electromagnetic results have been central to an integrated geophysical and geological interpretation.     

Age-dependent subduction of oceanic lithosphere beneath western South America

A recently established relation between the penetration depth of oceanic lithosphere and the lithospheric age appears to be of special interest to the understanding of the South American subduction zone. The main characteristics of this complicated zone, such as the absence of deep-focus earthquakes south of 30°S, the variations in the dip angle of the descending Nazca plate and the gap in seismic activity between depths of approximately 300 and 525 km, can be understood if the spatial and temporal variations in the age of the descending oceanic lithosphere are taken into account. In view of the significance of local aspects of the subduction process the South American-Nazca plate interaction cannot simply be considered as a type-example of the interaction between a continental and an oceanic plate.     

上通流对大陆岩石圈地幔-地壳热结构模式的潜在影响

本文以多孔介质中大尺度传热问题为基础,结合热平衡理论分析与数值计算,探讨了上通流对大陆岩石圈地幔-地壳热结构模式的潜在影响.根据大陆岩石圈中孔隙波传热概念模型的初步理论分析结果,指出了采用理论分析和数值模拟相结合的方法在研究大陆岩石圈地幔-地壳热结构模式时的重要性.理论分析方法可用来确定岩石圈尺度范围内大陆岩石圈的厚度和大陆地壳相关的边界条件,从而为地壳范围内数值模型的建立提供一些重要信息.数值模拟方法可以用来模拟地壳尺度范围内地壳的详细结构和复杂几何形状.如果地壳内的热分布是所考虑的主要因素,采用具有地壳尺度的合理数值模型可以有效减少计算机工作量.利用理论分析方法求出的岩石圈尺度范围内大陆岩石圈厚度与地幔传导热流之间关系的理论解,不仅可以用来验证模拟大陆岩石圈内传热问题所采用的数值方法, 而且可以用来初步研究大陆岩石圈内热分布的基本规律,为研究岩石圈地幔热事件中大陆岩石圈热减薄过程提供相应的边界条件.本文从理论分析的观点初步探讨了中国大陆不同构造背景下大陆岩石圈的热结构模式,其结果与从地球物理和地质资料中获得的大陆岩石圈热结构模式十分吻合.研究结果表明由大陆岩石圈中孔隙波传播所导致的上通流是影响大陆岩石圈地幔-地壳热结构模式及大陆岩石圈地幔与地壳之间物质和能量交换的可能机制之一.     

三河-平谷8级大震区地壳上地幔电性结构特征研究

用电磁阵列剖面法 (EMAP)、大地电磁测深方法 (MT) ,沿三河 -平谷 8级大震震源区 ,作了 31 8km长的EMAP探测和两条总长 150 0 5km共 36个点的MT探测。获得了研究范围内的地壳上地幔电性结构、高导层特征和陡变带、高导异常体、断裂展布、岩石圈结构等结果 ,为搞清地震危险区的深浅构造关系、从电性结构特征推测发震模式和预测未来强震的可能地点提供了介质电性的多种参数     

How is the European Lithosphere Imaged by Magnetotellurics?

I review recent investigations on the electrical conductivity of the lithosphere and asthenosphere in Europe. The principal method in the reviewed studies is the magnetotelluric method, but in many cases other electromagnetic methods (e.g., magnetovariational profilings and geomagnetic depth soundings) have provided additional information on subsurface conductivity or have been the primary method. The review shows that the magnetotelluric method has been used, and is being used, in all kinds of environments and for many different processes shaping the crust and lithosphere. The crust is very heterogeneous, both with respect to the scale of conductive/resistive features and interpretations: research targets vary from Archaean palaeostructures to ongoing processes. The European database of the depth to the lithosphere-asthenosphere boundary (LAB) in Europe is updated, and a new map showing lateral variations of the depth of LAB is provided. The compilation shows that (1) the Phanerozoic European lithosphere, with considerable variations (45–150 km), is much thinner than the Precambrian European lithosphere, (2) the Trans-European Suture Zone is a major electrical border in Europe separating electrically (as well as geophysically and geologically in general) two quite different settings, (3) the thinnest lithosphere is found under the extensional Pannonian Basin (45–90 km), (4) in most of the East European Craton there are no indications of a high conductivity zone in upper mantle. In many regions there is no information at all on upper mantle conductivity, which calls for pan-European projects to operate arrays of simultaneously recording instruments with long recording periods (2–8 months) and dense spatial sampling (20–50 km).     

Coast effect distortion of marine magnetotelluric data: Insights from a pilot study offshore northeastern Japan

We report on strong coast effect distortions observed for broadband marine magnetotelluric (MT) data collected on the forearc offshore northeastern Japan. Eight days of horizontal electric and magnetic fields recorded at eight seafloor stations and the horizontal magnetic fields from a land remote station were processed with a robust multiple-station algorithm, yielding good MT responses and inter-station transfer functions at periods of 7–10,000 s. Transverse electric (TE) mode responses have cusps in apparent resistivity and negative phases at periods around 1000 s, while the transverse magnetic (TM) mode responses are galvanically depressed below the TE responses. An analysis of inter-station transfer functions confirms that the apparent resistivity cusps are a magnetic field, rather than electric field, phenomenon, consisting of an amplitude minimum and rapid phase change around a characteristic frequency. Poynting vectors for a TE coast effect model study illustrate that the anomalous phases are associated with energy diffusing back up to the seafloor from below, after being turned around from its usual downward propagating trajectory by inductive coupling between the conductive ocean and the resistive seafloor along the continental margin. We show that the characteristic frequency and position of the TE mode apparent resistivity cusps are determined by a relatively simple combination of the electrical resistivity of the seafloor, the depth of the ocean, and the distance from the coastline. By including coastlines and bathymetry in 2D inversion, we recover the seafloor conductivity structure along the forearc, demonstrating that broadband data can constrain the thickness of conductive forearc sediments and the underlying high resistivity associated with the mantle wedge and subducting oceanic lithosphere.