A composite trans-Atlantic heat flow profile between 20°S and 35°S

Sixty new measurements together with published heat flow data in the South Atlantic between 20°S and 35°S latitude have been analyzed. Heat flux is greater through the eastern Mid-Atlantic Ridge flank and basin than their counterparts on the west but the standard deviation or spatial variation is greater on the west, contrary to expectation based on sediment thickness. The variance in the data indicates that this asymmetry in mean heat flux is statistically significant at 87% confidence level. A pair of ridge-flank minima appear in a composite trans-Atlantic profile of heat flux versus sea-floor age, suggesting hydrothermal circulation in the young oceanic crust. The Walvis Ridge has a mean excess heat flux of 28 mW m^?2 (0.7 μcal cm^?2 sec^?1) above the surrounding Cape and Angola Basins, and decreases along the ridge towards the northeast. Consistent with the apparent asymmetric distribution in the South Atlantic, it is also significantly higher than that of the Rio Grande Rise. We hypothesize that the trend and larger mean heat flux of the Walvis Ridge is best explained by a hot-spot origin, perhaps combined with higher radioactivity in the crust. However, the morphologic and heat flow differences between the Walvis Ridge and Rio Grande Rise suggest that these features have different geologic histories.     

New terrestrial heat flow measurements on the Nazca Plate

Sixty-seven new heat flow measurements on the Nazca Plate are reported, and the thermal regimes of three specific areas on the plate are examined. The Nazca Ridge is an aseismic ridge which may have been generated as an “island trail” from the Easter Island “hot spot” and/or may be a fossil transform fault. The Nazca Ridge has lower heat flow than the surrounding sea floor implying that the ridge might have low “effective” thermal conductivity causing heat to preferentially flow or refract to surrounding ocean crust which has higher conductivity, or, the low heat flow values may be caused by hydrothermal circulation on the ridge. The Carnegie Plateau is an elevated region south of the Carnegie Ridge on the northeastern Nazca Plate with high heat flow and shallow topography consistent with an age of less than 20 m.y. B.P. The central Nazca Plate is an area of highly variable heat flow which is possibly related to thin sediment and to rough regional topography.     

The early opening between Broken Ridge and Kerguelen Plateau

Reconstructions to total closure of the Australia-Antarctic continents causes an unacceptable overlap of Broken Ridge and Kerguelen Plateau. This has been partially resolved in the past by supposing that the northern part of Kerguelen, that is principally involved in the overlap, is younger than the remainder. We have revised the early reconstructions using a newly proposed breakup chronology of Australia and Antarctica which suggests that opening began at least 90 m.y. B.P. at an initially slow rate. This eliminates the overlap problem without invoking major age differences within the Kerguelen Plateau. We also suggest that the northeastern flank of Kerguelen Plateau may be underlain by the “missing” westward continuation of the Diamantina Zone. It may have been isolated on the Antarctic plate by a ridge crest jump at about anomaly 24 time that also formed the Ob Trench.     

The evolution of the Parece Vela Basin,eastern Philippine Sea

The Parece Vela Basin is a back-arc basin. It is approximately 5000 m deep and is divided into two topographic provinces by the north-trending Parece Vela Rift. The western province is thinly sedimented and topographically rough. The eastern province is blanketed by a thick apron of volcaniclastic sediments which were derived from the West Mariana Ridge. The Parece Vela Rift is composed of a series of discrete deeps and troughs with depths commonly of 6 km and locally exceeding 7 km.Petrologic and seismic refraction data indicate that the Parece Vela Basin is of oceanic character.Low-amplitude, nort-trending, lineated magnetic anomalies are present in the basin and appear symmetric about a line near the Parece Vela Rift. In the central latitudes of the basin seafloor spreading anomalies 10 (30 m.y. B.P.) to 5E or 5D (18 or 17 m.y. B.P.) can be identified. The uncertainty in identifying the youngest anomaly may be due to ridge jumps near the end of spreading. Spreading may have started slightly later in the northern end of the basin. Anomalies in the eastern province are disrupted and are difficult to correlate. DSDP results indicate post-spreading volcanism on the eastern side of the basin and this may have degraded the anomalies. The age obtained in the western province of the basin at DSDP Site 449 (～25m.y. B.P.) is in close agreement with that obtained from the magnetic data (～26m.y. B.P.).It is hypothesized that subduction was occurring at a west-dipping subduction zone east of the Palau-Kyushu Ridge in the Early Oligocene. This volcanic arc split about 31 or 32 m.y. ago and interarc spreading was initiated between the Palau-Kyushu Ridge (which then became a remnant arc) and the West Mariana Ridge. The Parece Vela Basin formed between the ridges by two-limb seafloor spreading. Spreading stopped about 17 or 18 m.y. ago.Like certain other marginal basins, the Parece Vela Basin is deeper than predicted from depth vs. age curves. The average heat flow for the Parece Vela Basin is in agreement with that predicted from heat flow vs. age curves.The origin of the Parece Vela Rift is unclear. It may represent the extinct spreading center or may be a postspreading feature.     

Massive sulfides in a sedimented rift valley, northern Juan de Fuca Ridge

A number of mounds each several hundred meters across and up to sixty meters high have been observed with SeaMARC II acoustic imagery and Seabeam bathymetry in the sediment-filled axial valley at the northern end of the Juan de Fuca Ridge. The mounds are located a few kilometers west of the eastern valley-bounding normal fault scarp where the local sediment fill is approximately 300 m thick. All of the mounds are believed to be of hydrothermal origin, and one is associated with anomalously high heat flow in excess of 1 W m⁻². A piston core collected from that mound comprises coarse clastic sulfide units interbedded with sulfidic muds. Hydrothermal minerals present in the 2.3 m section include pyrrhotite, pyrite, marcasite, sphalerite, chalcopyrite, iss (intermediate solid solution in the CuFeZnS system), chalcopyrrhotite, galena, talc, barite, and amorphous silica. Mineral fabrics of the clasts indicate that the material was precipitated at or near the sea floor by mixing of hot hydrothermal fluids with cold seawater. Low concentrations of Zn, Cu, Cd, and Ag relative to those found in unsedimented ridge hydrothermal deposits, and the presence of pyrrhotite as an early phase mineral indicates that the vent fluids have been modified by reaction with sediments beneath the mound. Rapid sedimentation in a rift valley is clearly conducive to the formation of large hydrothermal mineral deposits like those believed to be present within and beneath these mounds. The relatively impermeable sediment cover insulates the crust, inhibits groundwater recharge, promotes long-lived discharge at a restricted number of sites, provides a substrate for the efficient subsurface precipitation of minerals, and through continued sedimentation, protects surficial deposits from the corrosive effects of seawater. No reliable estimate of the bulk composition of the mounds can be made with existing data, but their size is comparable to major hydrothermal mineral deposits found on land; ancient settings in which many land deposits formed are in many ways similar to the one in which the features described here are currently forming.     

Heat flow and thickness of the oceanic lithosphere

Existing thermal models of the oceanic lithosphere predict too sharp an increase of heat flow towards the ridge axis. A new mathematical model of a thickening lithosphere is presented. The temperature distribution is computed by the use of observed surface heat flow as a boundary condition. If observed heat flow values represent flow of heat from the mantle, the model predicts a rather rapid growth of the lithosphere within the first 30 m.y. and a nearly steady state after 100 m.y. Heat flow from the asthenosphere to the lithosphere shows a minimum near the ridge axis, suggesting a down-going convective flow in the asthenosphere at both sides of a spreading center.     

Origin of metalliferous sediments from the East Pacific Rise

The distribution of several metals in East Pacific Rise sediments, when normalized to Al₂O₃, exhibit stronger maxima near the rise crest than when simply plotted on a carbonate-free basis. The similarity (1) between the distribution of metals in ridge sediments and previously measured mean heat flow values and (2) between the composition of crestal sediments and terrestrial ore bodies associated with greenstone belts, strongly supports a hydrothermal origin for rise crest sediments.     

Fluxes of fluid and heat from the oceanic crustal reservoir

Recent discoveries define a global scale fluid reservoir residing within the uppermost igneous oceanic crust, a region of seafloor that is both warm and may harbor a substantial biosphere. This hydrothermal fluid reservoir formed initially within volcanic rocks newly erupted at mid-ocean ridges, but extends to the vastly larger and older ridge flanks. Upper oceanic crust is porous and permeable due to the presence of lava drainbacks, fissuring, and inter-unit voids, and this porosity and permeability allows active fluid circulation to advect measurable quantities of lithospheric heat from the crust to an average age of 65 Myr. A compilation of crustal porosities shows that this fluid reservoir contains nearly 2% of the total volume of global seawater. Heat flow and sediment thickness data allow calculation of reservoir temperatures, predicting 40°C mean temperatures in Cretaceous crust. Utilizing these temperature estimates, heat flow measurements and models for the thermal structure and evolution of the oceanic lithosphere, we have computed mean hydrothermal fluxes into the deep ocean as a function of plate age. The total hydrothermal volume flux into the oceans approaches 20% of the total riverine input and may contribute to the global seawater mass balance.     

Seawater transport and reaction in upper oceanic basaltic basement: chemical data from continuous monitoring of sealed boreholes in a ridge flank environment

Osmotically pumped fluid samplers were deployed in four deep-sea boreholes that were drilled during Ocean Drilling Program (ODP) Leg 168 on the eastern flank of the Juan de Fuca Ridge. Samplers were recovered from ODP Sites 1024 and 1027 and aliquots were analyzed for a variety of dissolved ions. Results from both of the samplers show a drastic change in the major ion composition within the first 20-40 days after the borehole was sealed at the seafloor followed by a more gradual change in composition. This gradual change ceased after 820 days at Site 1024 but continued throughout the 3-year deployment at Site 1027. We modeled this change in composition to estimate the flux of formation fluid through the open borehole. The rapid early change requires a flow of ∼1500 kg of formation fluid per day. The more gradual later change requires flow rates of 38 kg/day at Site 1024 and 17.5 kg/day at Site 1027. The latter fluxes require a minimum average specific discharge of meters to hundreds of meters per year through the surrounding basaltic matrix. Trace element data show surprisingly little contamination given the presence of steel casing, Li-organic-rich grease at each joint, cement, and drilling muds. Observed changes in trace element concentrations relative to those of bottom seawater provide a measure for the global significance of cool (23°C; ODP Site 1024) ridge flank hydrothermal systems relative to warm (64°C; Baby Bare and ODP Site 1027) hydrothermal systems and illustrate the importance of these cooler systems to global geochemical budgets.     

Hyperbolic echo zones in the eastern Atlantic and the structure of the southern Madeira Rise

The Madeira Rise is a 450 km northeast-trending structural-sedimentary feature which lies west of the Madeira Islands. Its northern half is controlled by a basement ridge, but its southern section consists of an apparent current-controlled sediment deposit. Its maximum sediment thickness is about 1 km over a relatively level basement. There are two reflectors which can be traced within the sediment pile. A shallow reflector (R1) may mark the termination of a rapid constructional phase of the drift. Sediment cores taken on the southern Madeira Rise have recovered brown marls and chalks with sedimentation rates of only about 1.5 cm/1000 years over the last 225,000 years.A striking zone of hyperbolic reflectors mapped around the flank of the southern Madeira Rise in the 4000–4800 m depth range may be the expression of bedforms created by contour-following currents. Another zone of hyperbolic echoes is found on the continental rise at about 24–26°N in a similar depth range. Trend determinations in this area suggest that the bedforms which give rise to the hyperbolae are oriented north-northeast and may be similar to the abyssal furrows discovered by DEEP-TOW observations on the Blake-Bahama Outer Ridge.The bottom photographs available from both hyperbolic echo zones show a tranquil bottom. This suggests that the hyperbolic bedforms are relict.     

Two heat flow profiles across the Atlantic Ocean

Measurements of the terrestrial heat flux at 76 localities along 2 profiles across the Mid-Atlantic Ridge at 19.5°N and 8.5°S latitude are presented. Two high heat-flow values were measured 800 to 1000 km east of the ridge crest at 8.5°S, but no high values were found at the ridge crest at this latitude. Detailed surveys and heat-flow measurements near the ridge crests on both profiles indicate that bottom topography influences the heat-flow variability. The average heat flow on both profiles, about 1.4 to 1.5 × 10⁻⁶ cal/cm² sec, is close to the average for other ocean basins, in contrast to previous studies indicating lower heat flow for the Atlantic.     

Progress in understanding the age and origin of the alpha ridge, arctic ocean

Multiparameter geophysical measurements and geological samples from CESAR suggest that close to Canada the Alpha Ridge is oceanic in nature and was built in part by volcanic activity. It is unclear whether this part of the ridge formed in an intraplate or a plate margin environment. Estimates from paleontological, heat flow and magnetic data place the construction of a volcanic ridge within the Cretaceous period between about 120 and 80 Ma ago, the interval in which the Canada Basin seafloor formed but clearly before the creation of the Lomonosov Ridge. The place of the Makarov Basin in this chronology remains unclear.     

BSR热流的三维地貌校正和流体汇聚探测(英文)

从天然气水合物稳定区底界的地震似海底反射BSR(Bottom Simulating Reflector)深度计算得到的BSR热流包含了海底地貌(热流在凹地型会聚,在凸地形发散)和增生楔内部流体活动的影响。从BSR热流中移除地貌效应的贡献就能揭示出流体是否发生了汇聚。在难以使用解析方法计算地貌效应的复杂海底区域,三维有限元方法可以高精度的模拟地貌对背景热流的影响,从而可以对BSR热流进行地貌效应校正,得到平坦地形条件下的BSR热流,并进一步通过与背景热流值的对比,识别目前仪器所不能探测的流体汇聚区。在北卡斯卡底(Cascadia)俯冲边缘陆坡中部的研究区应用该方法,显示黄瓜岭(Cucumber Ridge)高地及其周围的海底热流正异常显著(高出背景热流值10-20%),同时这些区域在地震成像上与海底的裂隙系统相对应,指示了流体沿着这些高渗透率通道进行汇聚,并且很可能导致较高的水合物富集度。     

Correlated Nd,Sr and Pb isotope variation in Walvis Ridge basalts and implications for the evolution of their mantle source

Basement intersected in DSDP holes 525A, 528 and 527 on the Walvis Ridge consists of submarine basalt flows and pillows with minor intercalated sediments. These holes are situated on the crest and mid and lower northwest flank of a NNW-SSE-trending ridge block which would have closely paralleled the paleo mid-ocean ridge [13, 14]. The basalts were erupted approximately 70 m.y. ago, an age equivalent to that of immediately adjacent oceanic crust in the Angola Basin and consistent with formation at the paleo mid-ocean ridge [14]. The basalt types vary from aphyric quartz tholeiites on the ridge crest to highly plagioclase phyric olivine tholeiites on the ridge flank. These show systematic differences in incompatible trace element and isotopic composition. Many element and isotope ratio pairs form systematic trends with the ridge crest basalts at one end and the highly phyric ridge flank basalts at the other.The low ¹⁴³Nd/¹⁴⁴Nd (0.51238), ²⁰⁶Pb/²⁰⁴Pb (17.54), ²⁰⁸Pb/²⁰⁴Pb (15.47), ²⁰⁸Pb/²⁰⁴Pb (38.14) and high⁸⁷Sr/⁸⁶Sr (0.70512) ratios of the ridge crest basalts suggest derivation from an old Nd/Sm-, Rb/Sr- and Pb/U-enriched mantle source. This isotopic signature is similar to that of alkaline basalts on Tristan de Cunha but offset to significantly lower Nd and Pb isotopic ratios. The isotopic ratio trends may be extrapolated beyond the ridge flank basalts with higher¹⁴³Nd/¹⁴⁴Nd (0.51270), ²⁰⁶Pb/²⁰⁴Pb (18.32), ²⁰⁷Pb/²⁰⁴Pb (15.52), ²⁰⁸Pb/²⁰⁴Pb (38.77) and lower ⁸⁷Sr/⁸⁶Sr (0.70417) ratios in the direction of increasingly Nd/Sm-, Rb/Sr- and Pb/U-depleted source compositions. These isotopic correlations are equally consistent with mixing od depleted and enriched end member melts or partial melting of an inhomogenous, variably enriched mantle source. However, observe ZrBaNbY interelement relationships are inconsistent with any simple two-component model of magma mixing, as might result from the rise of a lower mantle plume through the upper mantle. Incompatible element and Pb isotopic systematics also preclude extensive involvement of depleted (N-type) MORB material or its mantle sources. In our preferred petrogenetic model the Walvis Ridge basalts were derived by partial melting of mantle similar to an enriched (E-type) MORB source which had become heterogeneous on a small scale due to the introduction of small-volume melts and metasomatic fluids.     

Numerical and approximate solutions for lithospheric thickening and thinning

Lithospheric thickness, surface heat flow and subsidence are calculated numerically as functions of time since the formation of oceanic lithosphere at a ridge crest. These numerical results confirm the validity of a previously derived approximate solution to the same problem. A numerical solution is also calculated for lithospheric thinning due to a sudden increase of heat flux into the base of the lithosphere. With this solution as a standard, an approximate solution is derived which can accurately predict thickness, heat flow, and uplift as functions of time. Lithospheric thinning and surface uplift begin immediately with the increase of heat flux at depth, but the increase in surface heat flow begins 20–40 m.y. later.     

Geothermal regime and geodynamics of the North Pacific Ocean

The distribution of heat flow in the North Pacific Ocean has been examined, and a map of geothermal and geomagnetic fields for the Bering Sea as it is known today has been made. Reliable data are lacking regarding the time of origin for features of oceanic and continental genesis in the Bering Sea, which is an obstacle to the study of geodynamic processes in the North Pacific. Heat flow data were used to yield numerical estimates for the age of seafloor features in the Bering Sea: the Kamchatka Basin (21 Ma), Shirshov Ridge (95 Ma for the northern part and 33 Ma for the southern), the Aleutian Basin (70 Ma), Vitus Rise (44 Ma), Bowers Ridge (30 Ma), and Bowers Basin (40 Ma). These age estimates are corroborated by combined geological, geophysical, and plate kinematic data. A thermochemical model of global mantle convection has been developed in order to perform a numerical simulation of the thermal process involved in the generation of extended regional features in the North Pacific (the Emperor Fracture Zone, Chinook Trough, etc.). The modeling suggests a plume-tectonic origin for these features, yielding the optimal model for the tectonic evolution of the North Pacific. An integrated geological and geothermal analysis leads to the conclusion that the northern and southern parts of the Shirshov Ridge are different, not only in geologic age, but also in tectonic structure. The northern part is of imbricated-thrust terrane origin, while the southern part is of ensimatic island-arc origin, similar to that of Bowers Ridge. The seafloor of the Aleutian Basin is an outlier of the Upper Cretaceous Kula plate where, in the Vitus Rise area, backarc spreading processes originated during Eocene time. The terminating phase of activity in the Bering Sea began about 21 Ma by spreading in the older seafloor of the Kamchatka Basin. We developed plate-tectonic reconstructions of evolution for the North Pacific for the times 21, 33, 40, and 70 Ma in the hotspot system based on age estimates for the seafloor features derived from heat flow data and modeling of the thermal generation of regional faults, as well as on an analysis of geomagnetic, tectonic, and geological data.     

南海西南次海盆的地热流特征与分析

为系统地了解南海西南次海盆的地热流特征,本文通过对研究区及邻域地热流数据的补充采集、收集整理和统计分析,获得了87个有效的地热流数据、一批热导率和生热率的地热参数资料,地热流测点在空间上覆盖了整个区域.研究区的地热流数据分布结果表明,西南次海盆热流密度的平均值为98.1±14.8 mW·m^-2,洋陆过渡带为103.6±19.4 mW·m^-2,南沙岛礁区和西部陆缘分别为79.0±15.5 mW·m^-2和78.3±15.6 mW·m^-2.研究区表层沉积物热导率的平均值0.86±0.06 W·mK^-1,生热率的平均值1.11±0.17 μW·m^-3,海底温度的平均值为2.43±0.01℃.综合海底地形地貌、地质与地球物理资料,认为研究区的热流特征在空间上具有一定的分布规律,表现为：（1）洋盆区测点的热流密度平均值高于两侧陆缘;（2）东南缘洋陆过渡带上测点的地热流密度值高于邻近海盆和南沙岛礁区的测点,而西北缘这种特征不明显;（3）西北翼的热流密度值总体比东南翼高;（4）沿着古扩张中心方向,西南次海盆热流值具有自东北向西南端方向逐步增大的趋势,表明海盆区同时存在着洋中脊与大陆裂谷两种不同的热状态,西南段裂谷热流值比东北段洋中脊高.对西南次海盆沉积物的热导率和生热率值参数的测量及数据空间分析可见,这两种热参数的空间分布无明显规律性,可能与海盆形成之后复杂的沉积环境相关.根据热流-洋壳年龄之间的关系,在西南次海盆东北段26个测站数据中,发现靠近古扩张中心的数据与理论值呈负偏移,而远离古扩张中心的数据呈正偏移,此现象是海盆内地热流数据受不同类型的地下流体影响所致.     

Emplacement processes of proto-arc basalt in the Izu–Bonin–Mariana arc system

Ascertaining the emplacement mechanism of oceanic basaltic lavas is important in understanding how ocean floor topography is produced and oceanic plates evolve, particularly during the early stages of crustal development of a supra-subduction zone. A detailed study of the volcanic stratigraphy at International Ocean Discovery Program (IODP) Site U1438 in the Amami Sankaku Basin, west of the Kyushu–Palau Ridge, has revealed the development of lava accretion and ridge topography on the Philippine Sea plate at about 49 Ma. Igneous basement rocks penetrated at Site U1438 are the uppermost 150 m of ~6 km-thick oceanic crust, and comprise, in a downhole direction, sheet flows (12.6 m), lobate sheet flows (61.3 m), pillow lavas (50.7 m), and thin sheet flows (25.3 m). The lowermost sheet flows are intercalated with layers of limestone and epiclastic tuff. Lithofacies analysis reveals that the lowermost sheet flows, limestone, and tuff formed on an axial rise, the pillow lavas were emplaced on a ridge slope, and the lobate sheet flows formed off ridge on an abyssal plain. The lithofacies of the basement basalt corresponds to the upper portions of fast-spreading oceanic crust, suggesting that subduction initiation was associated with intermediate to fast rates of seafloor spreading. The surface sheet flows are olivine–clinopyroxene-phyric basalt and differ from the lower basalt flows that contain phenocrysts of olivine and plagioclase, with or without clinopyroxene. The depleted chrome-spinel composition and olivine–clinopyroxene phenocryst assemblage in the surface sheet flows suggests a slight contribution of water for magma generation not present for the lower basalt flows. Considering the lithological difference between the backarc and forearc oceanic crust in the Izu–Bonin–Mariana arc, with sheet flow dominant in the former, seafloor spreading occurred faster in the later stage of subduction initiation.     

Earthquake source mechanisms from body-waveform inversion and intraplate tectonics in the northern Indian Ocean

Double-couple point-source parameters for 11 of the largest intraplate earthquakes in the northern Indian Ocean during the last 20 y were determined from a formal inversion of long-period P and SH waveforms. Nine of the events have centroid depths at least 17 km below the seafloor, well into the upper mantle; two have centroid depths as great as 39 km. Using the source mechanisms of these earthquakes, we distinguish two major intraplate tectonic provinces in the northern Indian Ocean. To the west of the Ninetyeast Ridge, in the southern Bay of Bengal, intraplate earthquakes have thrust-faulting mechanisms with P axes oriented N-S. The centroid depths of these earthquakes range from 27 to 39 km below the seafloor. Lithospheric shortening in this region is thus accomplished by thrust faulting in the strong core of the oceanic upper mantle, while other geophysical evidence suggests that shallow sedimentary and crustal layers apparently deform predominantly by folding. In the immediate vicinity of the Ninetyeast Ridge, earthquakes display strike-slip mechanisms with left-lateral motion on planes parallel to the ridge. This type of faulting occurs from at least 10°S to the northern end of the Ninetyeast Ridge near 10°N, where the ridge meets the Sunda Arc. Seismic activity diminishes to the east of the Ninetyeast Ridge, but is also characterized by strike-slip faulting. Despite these variations in deformational style, the inferred orientation of greatest compressive stress in the northern Indian Ocean displays a consistent long-wavelength pattern over a large portion of the Indian plate, varying smoothly from nearly N-S in the Bay of Bengal to NW-SE in the northeastern Indian Ocean. This plate-wide stress pattern and the high level of intraplate seismicity in the northern Indian Ocean are likely the results of substantial resistance, along the Himalayan continental collision zone, to the continued northward motion of the western portion of the Indian plate. Oceanic intraplate earthquakes in other regions, where the level of deviatoric stress associated with the long-wavelength part of the stress field is likely to be smaller, need not be comparably reliable indicators of the plate-wide stress field.     

Neotectonic modelling of the western part of the Africa–Eurasia plate boundary: from the Mid-Atlantic ridge to Algeria

In this work we use the thin-shell approximation to model the neotectonics of the western part of the Africa–Eurasia plate boundary, extending from the Mid-Atlantic ridge to Tell Atlas (northern Algeria). Models assume a nonlinear rheology and include laterally variable heat flow, elevation, and crust and lithospheric mantle thickness. Including the Mid-Atlantic ridge permits us to evaluate the effects of ridge push and to analyse the influence of the North America motion on the area of the Africa–Eurasia plate boundary. Ridge push forces were included in a self-consistent manner and have been shown to exert negligible effects in the neotectonics of the Iberian Peninsula and northwestern Africa. Different models were computed with systematic variation of the fault friction coefficient. Model quality was scored by comparing predictions of anelastic strain rates, vertically integrated stresses and velocity fields to data on seismic strain rate computed from earthquake magnitude, most compressive horizontal principal stress direction, and seafloor spreading rates on the Mid-Atlantic ridge. The best model scores were obtained with fault friction coefficients as low as 0.06–0.1. The velocity boundary condition representing spreading on the Mid-Atlantic ridge is shown to produce concentrated deformation along the ridge and to have negligible effect in the interior of the plates. However, this condition is shown to be necessary to properly reproduce the observed directions of maximum horizontal compression on the Mid-Atlantic ridge. The maximum fault slip rates predicted by the model are obtained along the Mid-Atlantic ridge, Terceira ridge and Tell Atlas front. Relatively high slip rates are also obtained in the area between the Gloria fault and the Gulf of Cadiz. We infer from our modelling a significant long-term seismic hazard for the Gloria fault, and interpret the absence of seismicity on this fault as possibly due to transient elastic strain accumulation. The present study has also permitted better understanding of the geometry of the Africa–Eurasia plate boundary from the Azores triple junction to the Algerian Basin. The different deformational styles seem to be related to the different types of lithosphere, oceanic or continental, in contact at the plate boundary.