洋壳渗透率结构对其内部热液对流的影响研究——基于热力学数值模拟

洋底特殊环境的限制使数值模拟成为研究海底热液对流与成矿机制的有效方法.本文通过数值模拟的方法,研究洋壳渗透率单因素变化对热液对流系统的形态和输出参数(热液喷发温度和热流值)进行研究.模拟结果表明,洋壳平均渗透率分别与热液喷发温度和热量输出呈反相关和正相关关系,符合达西流体热对流的基本解析规律.另外,洋壳渗透率的垂向变化不会使洋壳内部的对流系统产生明显的横向偏移,只使热液羽规模在垂向上产生规律性变化,且渗透率越大,热液羽规模越小.渗透率在水平方向的变化则是引发热液羽和热液喷发位置横向偏移的主要因素,且只有在特定的渗透率水平变化率范围内,上涌热液羽才会发生向高渗区的明显偏移,这也从一定程度上解释了现今拆离断层相关的热液系统均未发育于断层终止线上的事实.以大西洋Trans-Atlantic Geotraverse(TAG)热液区为例,断裂带高渗区必然影响相邻洋壳内部热液的上涌路径,但受区内拆离断裂带周缘的渗透率水平变化规律的影响,上涌热液羽不至于完全偏移至断裂带,从而产生拆离断层上盘的热液活动.     

Excess3He in the deep water over the Mid-Atlantic Ridge at 26°N: evidence of hydrothermal activity

We have observed distinct and significant patterns of excess³He in bottom waters over the Mid-Atlantic Ridge in the “TAG hydrothermal field”. The lateral regional gradients are comparable in magnitude to gradients observed in the Galapagos rift, an area of confirmed hydrothermal activity. Together with this, the³He patterns, magnitudes and vertical gradients all indicate on-going and continuous hydrothermal activity.     

Tectonic pattern of the Azores spreading centre and triple junction

The major tectonic elements of the Azores triple junction have been mapped using long-range side-scan sonar. The data enable the Mid-Atlantic Ridge axis to be located with a precision of a few kilometres. Major faults and other tectonic and volcanic elements of the ridge maintain their regional trend of 010° to 020° past the triple junction area. There is no oblique spreading, and only minor transform offsets of the Mid-Atlantic Ridge occur here. The main effect of the triple junction or Azores hot spot is to diminish the amplitude of the median valley to 200 m or less. There is no axial high: a topographic high seen on several profiles is located to the east of the Mid-Atlantic Ridge spreading axis and does not appear to have any fundamental significance.The third arm of the triple junction includes the Azores srreading centre which appears to have developed as a series of en echelon rifted basins (the Terceira Rift) extending from Formigas Trough at 36.8°N, 24.5°W to a point near 39.3°N, 28.8°W. There are indications that recent activity in the spreading centre may be concentrated in a series of ridges which flank the older rifted basins. Until recently the northwest end of the Terceira Rift was connected to the Mid-Atlantic Ridge axis either directly at an RRR junction, or via a transform fault. The triple junction has probably moved south during the last 6 Ma to a positin on the Mid-Atlantic Ridge near 38.7°N.Initiation of the Azores spreading centre may have occurred during the 36 Ma B.P. rearrangement of poles, with an RFF triple junction north from the East Azores fracture zone to the North Azores fracture zone and transferring a wedge of European plate to the African plate.The tectonic elements revealed by this study are in good agreement with inferred earthquake mechanisms and with the RM2 plate tectonic model of Minster and Jordan, but east-west motion between North America and Africa does not seem to be compatible with the other motions at the triple junction unless it is of very recent (2>3 Ma) origin.     

The Kane fracture zone in the Central Atlantic Ocean

The Kane fracture zone has been traced as a distinct topographic trough from the Mid-Atlantic Ridge near 24°N to the 80-m.y. B.P. isochron (magnetic anomaly 34) on either side of the ridge axis for a total of approximately 2800 km. Major changes in trend of the fracture zone occur at approximately 72 m.y. B.P. (anomaly 31 time) and approximately 53–63 m.y. B.P. (anomaly 21–25 time) which are the result of major reorientations in spreading directions in the central Atlantic Ocean.     

The 3He distribution in deep water over the Mid-Atlantic Ridge

³He/⁴He ratios in dissolved helium at GEOSECS stations 115, 117, 30, and 120 provide an east-west section across the Mid-Atlantic Ridge at 30°N. Below 1500 m depth, the³He/⁴He profiles show little structure and have values within a few percent of the atmospheric ratio, indicating that the Mid-Atlantic Ridge is not a significant source of injected³He at this latitude. Mean³He/⁴He ratios calculated for the deep water at each station show that the³He/⁴He ratio in the western Atlantic at 30°N is 2–3% higher than in the eastern basin, probably due to mixing between a³He-rich boundary current in the western basin and low-³He deep water to the east.     

Geochemistry of basalts drilled in the North Atlantic by IPOD Leg 49: Implications for mantle heterogeneity

IPOD Leg 49 recovered basalts from 9 holes at 7 sites along 3 transects across the Mid-Atlantic Ridge: 63°N (Reykjanes), 45°N and 36°N (FAMOUS area). This has provided further information on the nature of mantle heterogeneity in the North Atlantic by enabling studies to be made of the variation of basalt composition with depth and with time near critical areas (Iceland and the Azores) where deep mantle plumes are thought to exist. Over 150 samples have been analysed for up to 40 major and trace elements and the results used to place constraints on the petrogenesis of the erupted basalts and hence on the geochemical nature of their source regions.It is apparent that few of the recovered basalts have the geochemical characteristics of typical “depleted” midocean ridge basalts (MORB). An unusually wide range of basalt compositions may be erupted at a single site: the range of rare earth patterns within the short section cored at Site 413, for instance, encompasses the total variation of REE patterns previously reported from the FAMOUS area. Nevertheless it is possible to account for most of the compositional variation at a single site by partial melting processes (including dynamic melting) and fractional crystallization. Partial melting mechanisms seem to be the dominant processes relating basalt compositions, particularly at 36°N and 45°N, suggesting that long-lived sub-axial magma chambers may not be a consistent feature of the slow-spreading Mid-Atlantic Ridge.Comparisons of basalts erupted at the same ridge segment for periods of the order of 35 m.y. (now lying along the same mantle flow line) do show some significant inter-site differences in Rb/Sr, Ce/Yb,⁸⁷Sr/⁸⁶Sr, etc., which cannot be accounted for by fractionation mechanisms and which must reflect heterogeneities in the mantle source. However when hygromagmatophile (HYG) trace element levels and ratios are considered, it is the constancy or consistency of these HYG ratios which is the more remarkable, implying that the mantle source feeding a particular ridge segment was uniform with respect to these elements for periods of the order of 35 m.y. and probably since the opening of the Atlantic. Yet these HYG element ratios at 63°N are very different from those at 45°N and 36°N and significantly different from the values at 22°N and in “MORB”.The observed variations are difficult to reconcile with current concepts of mantle plumes and binary mixing models. The mantle is certainly heterogeneous, but there is not simply an “enriched” and a “depleted” source, but rather a range of sources heterogeneous on different scales for different elements — to an extent and volume depending on previous depletion/enrichment events. HYG element ratios offer the best method of defining compositionally different mantle segments since they are little modified by the fractionation processes associated with basalt generation.     

The distribution of helium in oceanic basalt glasses

We have determined the concentrations and isotopic composition of helium in oceanic basaltic glass both by melting and by crushing in vacuo. A significant fraction of the helium is released by crushing, confirming that it resides within the vesicles. Comparison of volume percent vesicles to the fraction of helium contained in the vesicles gives qualitative agreement with experimental gas-melt partitioning studies. Measured concentrations are therefore, a function of original helium content, magmatic history, vesicle size and quantity, and grain size analyzed. Helium released by crushing is isotopically indistinguishable from that contained in the glass. Diffusion rates for helium in basaltic glass (in the temperature range 125–400°C) determined using the method of stepwise heating, yielded an activation energy of 19.9 ± 1 kcal/mole andlnD₀ = ?2.7 ± 0.6 (cgs units). Extrapolation of these results to ocean floor temperatures (0°C) gives a diffusivity of 1.0 ± 0.6 × 10^?17 cm²/s, indicating that diffusion is an insignificant mechanism for helium loss from fresh basaltic glasses. The³He/⁴He ratios are remarkably constant (at 1.10 ± 0.03 × 10^?5) for samples from the Mid-Atlantic Ridge (FAMOUS area and 23°N), the Juan de Fuca Ridge, the Galapagos spreading center, the Mid-Cayman Rise, and the Central Indian Ocean Ridge. This result is interpreted in terms of similar geochemical histories within the source regions for these samples.     

Asymmetric fracture zones and sea-floor spreading

An asymmetric pattern is observed in the orientation of minor fracture zones about the axis of the Mid-Atlantic Ridge at five sites where relatively detailed studies have been made between latitudes 22°N and 51°N. The minor fracture zones intersect the axis of the Mid-Atlantic Ridge in an asymmetric V-shaped configuration. The V's point south north of the Azores triple junction (38°N latitude) and point north south of that junction.The rates and directions of sea-floor spreading are related to the asymmetric pattern of minor fracture zones at the sites studied. Half-rates of sea-floor spreading averaged between about 0 and 10 m.y. are unequal measured perpendicular to the ridge axis. The unequal half-rates of spreading are faster to the west north of the Azores triple junction and faster to the east south of that junction. The half-rates of sea-floor spreading calculated in the directions of the asymmetric minor fracture zones are equal about the ridge axis within the uncertainty of the direction determinations.A discrepancy exists between minor fracture zones that form an asymmetric V about the axis of the Mid-Atlantic Ridge, and major fracture zones that follow small circles symmetric about the ridge axis. To reconcile this discrepancy it is proposed that minor fracture zones are preferentially reoriented under the influence of a stress field related to interplate and intraplate motions. Major fracture zones remain symmetric about the Mid-Atlantic Ridge under the same stress field due to differential stability between minor and major structures in oceanic lithosphere. This interpretation is supported by the systematic variation in the orientation of minor fracture zones and the equality of sea-floor spreading half-rates observed about lithospheric plate boundaries.     

Sulfur isotope systematics at the 21°N site,East Pacific Rise

Hydrogen sulfide in hydrothermal vent fluid at the 21°N site is enriched in ³⁴S relative to Mid-Ocean Ridge basalts, probably by addition of H₂S reduced from seawater sulfate by FeO-bearing basalt. Metalliferous sulfides are depleted in ³⁴S relative to the fluid from which they apparently precipitated, the degree of depletion reflecting the microenvironment in which each mineral crystallised and/or kinetic effects. Isotopic compositions of coexisting sulfides in a basal mound are consistent with equilibration at around 445°C, though heating to such a high temperature seems unlikely. Similar sulfides in a black smoker and in a dead chimney are out of isotopic equilibrium at any temperature, apparently reflecting a complex series of replacement mineralisations and post-depositional oxidation, respectively.     

Geothermal constraints on the hydrological regime of the TAG active hydrothermal mound, inferred from long-term monitoring

During August 1994 to March 1995, a period that included ODP Leg 158 drilling, bottom-water and sub-bottom temperatures were continuously logged by a long-term temperature monitoring system ‘Daibutsu’ at the base of the central black-smoker complex (CBC) and within the low heat flow zone at the TAG hydrothermal mound on the Mid-Atlantic Ridge. The temperature of hydrothermal fluid at CBC was also measured with a small high-temperature probe ‘Hobo’. Bottom-water temperature variations measured with Daibutsu at both sites have predominant semi-diurnal periods, causing the sub-bottom temperatures to fluctuate at these periods with reduced amplitudes and phase delays at sub-bottom depths. Seawater entrainment into the mound has been previously suggested at the low heat flow zone. We quantitatively evaluate the seawater entrainment rate at both sites from a one-dimensional numerical model, combined with a heat conduction model for the semi-diurnal variations. The entrainment rate of seawater at the base of CBC is estimated as 1.3±0.5×10⁻⁵ m/s, at least from August 17 to 30, 1994. On the other hand, the seawater entrainment rate at the low heat flow zone was undetected by long-term temperature monitoring at shallow sub-bottom depth. Nevertheless an increase in heat flow observed at the low heat flow zone during ODP drilling can be interpreted as a decrease in the entrainment rate of seawater. Before ODP Leg 158, Daibutsu measured three sub-bottom temperature anomalies at the base of CBC not derived from bottom-water temperature variations and Hobo also detected a CBC fluid temperature anomaly, indicating some natural changes in fluid flow within the mound. Daibutsu and Hobo also measured temperature anomalies during and after drilling at the ODP TAG-1 area. The Hobo temperature anomalies are inferred to have occurred when the cold fluid entrained through the drill holes at TAG-1 site reached or cooled the main fluid path to CBC. The entrained seawater through the drill holes appears to have contributed to dissolution and precipitation of anhydrite within the mound and perhaps affected the local permeability structure inside the mound. The temperature anomalies measured with Daibutsu at the base of CBC may have been induced by the change in the fluid flow pattern as a result of such permeability changes within the mound.     

The Mid-Atlantic Ridge at 33°N: the Hayes Fracture Zone

A geophysical study was conducted over the Mid-Atlantic Ridge between 32–39°N and 30–40°W. A particularly deep fracture was observed which offset the ridge crest 110 km in the vicinity of 33°N. A pole of relative motion between the North American and African plates was deduced from this fracture zone as being at 63.1°N, 17°W.     

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.     

Gypsum and halite from the Mid-Atlantic Ridge,DSDP Site 395

Gypsum and halite crystals, together with saponite and phillipsite, were found in a vein in a basalt sill 625 m below the sea floor at DSDP Site 395A, located 190 km west of the crest of the Mid-Atlantic Ridge. The δ³⁴S value of the gypsum (+19.4‰) indicates a seawater source for the sulfate. The δ¹⁸O values of the saponite (+19.9‰) and phillipsite (+18.1‰) indicate either formation from normal seawater at about 55°C or formation from¹⁸O-depleted seawater at a lower temperature.The gypsum (which could be secondary after anhydrite) was formed by reaction between Ca²⁺ released from basalt and SO₄^2? in circulating seawater. The halite could have formed when water was consumed by hydration of basalt under conditions of extremely restricted circulation. A more probable mechanism is that the gypsum was originally precipitated as anhydrite at temperatures above 60°C. As the temperature dropped the anhydrite converted to gypsum. The conversion would consume water, which could cause halite precipitation, and would cause an increase in the volume of solids, which would plug the vein and prevent subsequent dissolution of the halite.     

Zinc stable isotopes in seafloor hydrothermal vent fluids and chimneys

Many of the heaviest and lightest natural zinc (Zn) isotope ratios have been discovered in hydrothermal ore deposits. However, the processes responsible for fractionating Zn isotopes in hydrothermal systems are poorly understood. In order to better assess the total range of Zn isotopes in hydrothermal systems and to understand the factors which are responsible for this isotopic fractionation, we have measured Zn isotopes in seafloor hydrothermal fluids from numerous vents at 9–10°N and 21°N on the East Pacific Rise (EPR), the TAG hydrothermal field on the Mid-Atlantic Ridge, and in the Guaymas Basin. Fluid δ⁶⁶Zn values measured at these sites range from + 0.00‰ to + 1.04‰. Of the many physical and chemical parameters examined, only temperature was found to correlate with fluid δ⁶⁶Zn values. Lower temperature fluids (< 250 °C) had both heavier and more variable δ⁶⁶Zn values compared to higher temperature fluids from the same hydrothermal fields. We suggest that subsurface cooling of hydrothermal fluids leads to precipitation of isotopically light sphalerite (Zn sulfide), and that this process is a primary cause of Zn isotope variation in hydrothermal fluids. Thermodynamic calculations carried out to determine saturation state of sphalerite in the vent fluids support this hypothesis with isotopically heaviest Zn found in fluids that were calculated to be saturated with respect to sphalerite. We have also measured Zn isotopes in chimney sulfides recovered from a high-temperature (383 °C) and a low-temperature (203 °C) vent at 9–10°N on the EPR and, in both cases, found that the δ⁶⁶Zn of chimney minerals was lighter or similar to the fluid δ⁶⁶Zn. The first measurements of Zn isotopes in hydrothermal fluids have revealed large variations in hydrothermal fluid δ⁶⁶Zn, and suggest that subsurface Zn sulfide precipitation is a primary factor in causing variations in fluid δ⁶⁶Zn. By understanding how chemical processes that occur beneath the seafloor affect hydrothermal fluid δ⁶⁶Zn, Zn isotopes may be used as a tracer for studying hydrothermal processes.     

The structure, mass and interactions of the hydrothermal plumes at 26°N on the Mid-Atlantic Ridge

Water-column surveys by combined nephelometer/CTD (NCTD) tows contributed to the 1985 discovery of the first black smokers on the Mid-Atlantic Ridge. Subsequent regional water-column mapping has helped define the extent, mass and interactions of the suspended particulate matter phase (SPM) of the hydrothermal plumes emanating from the known and other nearby sources. The results of 29 NCTD cast/tows, covering 25–30 km² of ridge segment, indicate the presence of as many as two additional sources based on SPM concentration gradients and plume-top doming over source areas. Plume doming, documented here for the first time from field observations, conforms strikingly with laboratory experiments and can serve as a marker for source field location. A comparison of the plumes' SPM with potential temperature and salinity distributions indicates close correlation in water-column anomaly patterns for each, confirming modification of the regional potential temperature and salinity structure by hydrothermal plumes, which is expressed by wide separation and sloping of isotherms and isohalines.     

Sea-floor spreading: The central anomaly magnetization high

Indications of a narrow region of high magnetization within the central magnetic anomaly on some mid-ocean ridges are found on near-bottom and sea surface magnetic profiles. This zone, which probably represents the most recent extrusions onto the ocean floor, is similar to the narrow region of high magnetization found on the Mid-Atlantic Ridge at 45°N with a suite of dredge samples. This narrow region is probably the result of the initial high magnetization of pillow basalts when they are extruded onto the ocean floor and the subsequent rapid oxidation of the outer variolitic zone of the pillows. The large-amplitude, short-wavelength (<15 km) magnetic anomaly found within the central anomaly over both slow- and fast-spreading ridges is produced by this narrow magnetization high. This magnetic anomaly can be used to locate the region of most recent extrusions on most ridges. The absence of this short-wavelength anomaly on some ridges may reflect the episodicity with which basalts are extruded onto the ocean floor.     

Spinels in Mid-Atlantic Ridge basalts: Chemistry and occurrence

Three groups of spinels have been identified in dredged basalts from the Mid-Atlantic Ridge in the Azores region (30–40°N): (1) magnesiochromites with 0.4–0.5 Cr/(Cr + Al) are most common and characteristic of olivine tholeiites of the region; (2) titaniferous magnesiochromites are found in an olivine basalt with alkali affinities, of local occurrence and evolved in relatively high fugacity of oxygen; (3) chromian spinels with 0.23 Cr/(Cr + Al) occur in unusual high-Al picrites of local occurrence and possible high-pressure origin. Spinels are restricted in occurrence to the least fractionated lavas, with FeO^*/FeO^* + MgO ratio less than 0.575 and with Cr content greater than 350 ppm. A close relationship between Al content of liquidus spinel and Al content of magma has been observed for basaltic types. High-Al spinels deviating from this relationship, such as those found in picritic lavas from the Mid-Atlantic Ridge, may have crystallized at high pressure. The use of spinels as geobarometers in magmas of a restricted compositional range seems a promising prospect. There is no evidence of systematic variation in spinel chemistry of occurrence along the Mid-Atlantic Ridge, such as could be related to different mantle sources of the basalts, plume versus non-plume or binary mantle mixing.     

Crystal morphologies in pillow basalts: implications for mid-ocean ridge processes

We present the results of a detailed petrological study of a sparsely phyric basalt (MAPCO CH98-DR11) dredged along the Mid-Atlantic Ridge (30°41′N). The sample contains microphenocrysts of olivine that display four different rapid-growth morphologies. Comparison of these morphologies with those obtained in dynamic crystallization experiments allows us to constrain the thermal history of the sample. The dendritic morphology (swallowtail, chain and lattice olivine) is directly related to the final quenching during magma–seawater interaction. In contrast, the three other morphologies, namely the complex polyhedral crystal, the closed hopper and the complex swallowtail morphology result from several cycles of cooling–heating (corresponding to a maximum degree of undercooling of 20–25°C) during crystal growth. These thermal variations occurred before eruption and are interpreted to be the result of turbulent convection in a small magmatic body beneath the ridge. The results suggest that the Mid-Atlantic Ridge is underlain by a mush zone that releases batches of liquid during tectonic segregation. Aphyric basalts are emitted during eruptions controlled by the tectonic activity, whereas phyric basalts correspond to small fractions of magma from the mush zone mobilized by reinjections of primitive magmas.     

The Newfoundland Basin: ocean-continent boundary and Mesozoic seafloor spreading history

The ocean-continent boundary in the Newfoundland Basin is defined as the seaward limit of a continental margin magnetic smooth zone. East of the Grand Banks this boundary is marked by a prominent NNE-trending magnetic anomaly that is correlated with the J-Anomaly (115 m.y.). South of Flemish Cap the smooth zone boundary strikes approximately 060° and is approximately 15 m.y. younger. Magnetic anomaly trends suggest two directions of motion during separation of Iberia and North America. The first phase of motion, commencing at J-Anomaly time with a spreading center strike of 015°, produced a rifted margin along the Grand Banks south of the Newfoundland Seamounts. No spreading occurred north of the seamounts during this phase, implying a counter-clockwise rotation of Iberia and no Grand Banks-Galicia Bank separation. The second phase began at about 102 m.y. with a shift of the pole of rotation to a location near Paris, producing a ridge orientation of approximately 060°. This spreading center extended north and east into the northern Newfoundland Basin and Bay of Biscay, producing a rifted margin south of Flemish Cap and opening of Biscay. This ridge geometry produced a component of extension across the Newfoundland Fracture Zone and the southeastward migration of the resultant “leaky” transform fault between 102 m.y. and the next pole shift produced the volcanic edifice of the Southeast Newfoundland Ridge. Fracture zone trends during this phase also exerted strong control on volcanism within the Newfoundland Seamount province; this activity ceased at about 97 m.y. The date at which the second phase ended is not well defined by presently available data. A RRR triple-junction existed in the northeastern Newfoundland Basin-western Biscay region for a short time prior to anomaly33/34 (80 m.y.) which marks the inception of a continuous Mid-Atlantic Ridge spreading center between the Newfoundland and Charlie Gibbs Fracture Zones.     

Distributions of hydrothermal activity along the Mid-Atlantic Ridge: interplay of magmatic and tectonic controls

Hydrothermal activity has been investigated along three different sections of the slow-spreading Mid-Atlantic Ridge (MAR): 11°20′–30°N, 36–38°N and the Reykjanes Ridge, 57°45′–63°06′N. When considered in total, the incidence of venting along these three sections of the MAR compares well with the predictions of a model in which frequency of venting is linearly related to ridge-crest spreading-rate. At the scale of individual study areas, however, departure from the model is observed by up to an order of magnitude. Venting is anomalously rare along the Reykjanes Ridge but anomalously abundant along the MAR 36–38°N. Whilst such variability may be within the error of the linear spreading-rate model, we note that the interplay between magmatic and tectonic processes also differs between the three study areas. In the case of the Reykjanes Ridge we propose that the low incidence of venting reported may reflect a limitation of the sampling/investigative strategy because the style of venting which predominates may not give rise to conventional black-smoker hydrothermal plumes. Along the oblique and broadly segmented MAR 36–38°N, we propose that vigorous hydrothermal venting in broad segment-end non-transform discontinuities may be focussed along deeply penetrating active faults with the requisite heat supply being supported through some combination of along-axis magmatic intrusions and thermal release associated with the serpentinisation of crustal peridotites.