Annals 47, 1, 2004, 01/07def 107 ANNALS OF GEOPHYSICS, VOL. 47, N. 1, February 2004 Key words crustal exploration – magnetotellurics – geothermal areas – Southern Tuscany 1. Introduction The magnetotelluric (MT) method is used to estimate the electrical characteristics of Earth structures using naturally occurring electromag- netic fields. It is particularly suited to the task of exploring areas of high heat flow since these ar- eas are commonly associated with dynamic ac- tivities such as magma emplacement, crustal fracturing, and the circulation of hot, electri- cally conductive fluids. Southern Tuscany, which is characterized by a high heat flow and the presence of two of the most important geot- hermal areas of the world, is an excellent test site for resistivity characterization. Indeed, Larderello has been studied in detail by MT soundings since 1973 (Fiordelisi et al., 1998; and references therein). For a long time MT sur- veys in Southern Tuscany region have encoun- tered significant difficulties that mainly stem from three factors: i) the high level of electro- magnetic signals from industrial and cultural sources, which interfere with the natural fields used in the method; ii) the presence of a very conductive shallow formation, and hence the need to acquire data for long periods; and iii) the naturally-occurring structural complexity commonly found in these areas, which requires two- or three-dimensional MT modelling for in- terpretation purposes, thus significantly increas- ing the amount of effort required by modelling. Apart from the uncorrelated noise that mainly affects the high frequencies linked to the wide variety of cultural noise common to all densely populated countries, an extremely coherent noise was found to be pervasive in the period Resistivity and heterogeneity of Earth crust in an active tectonic region, Southern Tuscany (Italy) Adele Manzella CNR-Istituto di Geoscienze e Georisorse, Pisa, Italy Abstract Southern Tuscany, belonging to the inner zone of the Northern Apennines, Italy, is an ideal laboratory for ob- serving the physical features of the lithosphere and their evolution in tectonically active regions. Here the crust is very thin, with a thickness of less than 25 km, and heat flow is very high, hence only very shallow depths of exploration are needed to investigate many of the middle-deep crustal features that are common to many other parts of the world. The magnetotelluric (MT) surveys performed in this region have provided information on the resistivity structure, which is related to the extent and distribution of free fluids and to the partial melts in the crust. The picture emerging from these MT surveys is that of a resistivity structure that is only partly related to the heat flow regime of the area. A very low resistivity was found below the vapour-dominated geothermal sys- tem of Larderello and below areas that have no clear connection to any geothermal system, whereas this reduc- tion of resistivity is less conspicuous below the water-dominated geothermal system of Mt. Amiata. Mailing address: Dr. Adele Manzella, CNR-Istituto di Geoscienze e Georisorse, Via Moruzzi 1, 56124 Pisa, Italy; e-mail: manzella@igg.cnr.it 108 Adele Manzella range of 1-1000 s. This noise, produced by the currents discharged by the direct-current electri- fied railways at the passage of trains, required considerable study and data processing (Fiorde- lisi et al., 1995; Larsen et al., 1996). The most recent MT surveys, however, have shown that in most areas this noise can be effectively re- moved and the data interpreted without distor- tion by using appropriate data acquisition and processing techniques. The continuous profil- ing technique was also used in the most recent surveys to increase the resolution of the shal- lower structure. The MT surveys in Southern Tuscany were undertaken primarily for geothermal and deep crustal exploration (CROP, Deep Crust Explo- ration Projects) (De Angelis et al., 1998; Fiordelisi et al., 1998; Manzella et al., 1999). The study area described in this paper is crossed by the CROP 03 and CROP 18 profiles (fig. 1). A total of 232 MT sites were acquired over the last 10 years, 108 of which were acquired for geothermal purposes and were referenced to a remote site on the island of Capraia in order to remove the extremely coherent electromagnetic noise. The remainder were single-site or locally referenced and many are affected by noise. This paper describes the overall results ob- tained by modelling the entire set of data recorded after 1991. This gives a broad view of the deep electrical structure of the region, show- ing a conductive medium-low crust of regional scale but with some important exceptions that are not simply related to the geothermal sys- tems. The interpretation of the MT data could provide some indication of the role played by Fig. 1. Location of the magnetotelluric sites in a schematic geological map. The CROP profiles and the main electrified rail tracks are shown, and the main survey profiles are indicated by letters for reference in the text. 109 Resistivity and heterogeneity of Earth crust in an active tectonic region, Southern Tuscany (Italy) electrical conductivity in the middle-lower crustal regimes. 2. Geological background and geophysical outline Southern Tuscany belongs to the inner part of the Northern Apennines, and originated from the collision between the Corsica-Sardinia and Apulia microplates (Gianelli et al., 1997; and references therein). The region is characterised by NNW-SSE faults linked to the Alpine oroge- ny (fig. 2). After the post-Tortonian, prevalent- ly extensional tectonics, the region underwent a general uplift, so that older units crop out in some places. Its extension was accompanied by anatectic magmatism whose ages decrease in an eastward direction. Shallow intrusive bodies are the heat sources of the two most important ge- othermal areas, Larderello and Mt. Amiata, which have undergone intensive exploratory and production drilling for almost half a centu- ry. Direct exploration by means of hundreds of shallow and deep wells has permitted to achieve a detailed reconnaissance of the first 4.5 km of the crustal structure. At Larderello drill-holes within the first 4.5 km of depth encountered se- quences of sedimentary, metamorphic and ig- neous rocks. Below the Jurassic-to-Cenozoic ophiolite and flysch units (Ligurids) and a Tri- assic-to-Tertiary allochthonous unit (Tuscan Nappe), there is a zone of tectonic slices, in- volving Mesozoic rocks and Paleozoic forma- tions. The deepest metamorphic units are mi- caschist and gneiss. Buried dykes, subvolcanic bodies of granite and thermo-metamorphic minerals were frequently crossed by deep drillings. The Mt. Amiata area, with respect to Larderello area, is characterised by the pres- ence of lava flows covering the sedimentary se- quence and by the particular nature of its base- ment, which is mostly made up of graphite- bearing metasedimentary rocks. The vapour-dominated geothermal system of Larderello and the water-dominated geother- mal system of Mt. Amiata both consist of two reservoirs: a shallow reservoir within the car- bonate units at 0.5-1 km depth, and a deeper, more extensive reservoir within the metamor- phic basement, at depths of more than 2 km be- neath the known geothermal fields, at pressures of up to 70 bar and temperatures between 300 and 350°C (Barelli et al., 2000; Fiordelisi et al., Fig. 2. Location map of Southern Tuscany, showing the geothermal areas and the main faults. 110 Adele Manzella 2000). The geothermal fluids circulate through the larger fractures and faults, whose main ori- entation is the regional Apenninic one. Perme- ability is usually low and heterogeneously dis- tributed. The results of many studies on contact metamorphic and hydrothermal minerals sug- gest that there were two main stages of hy- drothermal activity: an early, high-temperature stage related to the intrusion of the granites, which produced contact metamorphic and metasomatic processes in the basement rocks, and a second stage characterised by precipita- tion of lower-temperature mineral assemblages filling veins at shallow-intermediate depths or, in certain places, replacing early-stage contact metamorphic or igneous minerals (Magro et al., 2003; and references therein). Since the struc- ture is not widely fractured and the alteration minerals are localized around the faults and fractures, the alteration is not widespread. Seismological studies have shown that the area extending from NE of Corsica to the inner part of the Northern Apennines is characterised by a thinned continental lithosphere and a large anomalous structure within the upper mantle (Gianelli et al., 1997; and references therein). The thin, elongated Tuscan Moho (depth ∼ 25 km) lies directly over a soft asthenosphere that has risen by 40-50 km. The effects of this thin- ning and the subsequent magma intrusions at shallow depth are evident in the present heat flow regime, which shows a regional value of 120 mW/m2 and higher values in the geother- mal areas (fig. 3a). The thermal regime leads to a lower seismic velocity and density in the crust and in the upper mantle, as revealed by seismic tomography and gravity data (Batini et al., 1995; Bernabini et al., 1995). Minimum Bouguer values in the geothermal areas (fig. 3b), filtered for shallow structures and mantle effects, have been modelled as a roughly drop- shaped low-density body between 7-8 and 20 km depth, and interpreted as shallow intrusive and still molten bodies (Bernabini et al., 1995; Fiordelisi et al., 1995). Seismic data provided valuable indications as to the shape and depth of the anomalous body in the Larderello region. Inversion of local earthquake arrival times and teleseismic travel-time residuals revealed the presence of a low velocity anomaly, 20 km Fig. 3a-c. Maps showing (a) heat flow regime, b) Bouguer anomaly and (c) K-horizon depth. a b c 111 Resistivity and heterogeneity of Earth crust in an active tectonic region, Southern Tuscany (Italy) wide at its top and possibly wider as it ap- proaches the Moho, at depths of more than 8 km below ground level (see fig. 4a,b) (Batini et al., 1995). A low velocity body was also identi- fied below the Amiata area, but resolution here is much lower than at Larderello (Chiarabba et al., 1995). Another clue is provided by a seis- mic horizon, K, which was identified in the Amiata area and given the same name as the seismic horizon first observed at Larderello. Seismic reflectors having the same seismic sig- nature as the K-horizon were recognised not on- ly in the geothermal fields of Larderello and Mt. Amiata, but also along the deep seismic re- flection profiles CROP-03 and CROP-18 (e.g., Cameli et al., 1998). These horizons are often referred to as a single horizon, but this is not strictly true. The K-horizon locally shows «bright spot» features and is located at 3-4 km depth below the central parts of the geothermal fields of Larderello and Mt. Amiata, at 5-6 km depth along the margins of these same fields, and at 10-12 km depth outside the geothermal areas, where it becomes relatively flat. The Fig. 4a,b. a) 2D inversion results of the TM-mode MT data along AA′ profile in the Larderello area; b) local earthquake and teleseismic inversion results along the same profile (from Batini et al., 1995, modified). a b 112 Adele Manzella form of the K-horizon shows many similarities with the heat flow and gravity maps (fig. 3c), rising sharply above the deep anomalous bodies defined by seismic tomography. The most com- mon interpretation of this horizon is that it rep- resents the brittle-ductile transition of the Tus- can crust, where fluid-filled fractures determine the bright-spot appearance in geothermal areas (Cameli et al., 1998; and references therein). 3. Resistivity data and discussion Resistivity data have never achieved a deci- sive importance in Southern Tuscany. Due to the low resistivity of the shallower formations, the thousands of geoelectrical data acquired in the geothermal areas have never proved to be of real interest since they were able at most to de- fine the resistivity distribution down to the re- sistive basement, whereas the main geothermal target for the past 30 years has been the deep reservoirs within the metamorphic basement. Moreover, due to the geological characteristics of the formations hosting the geothermal reser- voirs in Tuscany and to their low permeability, geothermal fluid circulation does not produce a widespread alteration zone and the contribution to resistivity of very conductive clayey minerals can be neglected. Hence, the decrease in resis- tivity caused by thermal fluid circulation was not expected to be high, especially in the Larde- rello geothermal area, where the dominant flu- id phase is steam. The geoelectrical data acquired in the Larderello and Mt. Amiata areas, together with the results of a few resistivity logs in the Larde- rello area, and laboratory data (Losito, 1991), have been used to determine the resistivity val- ues of the shallow geological units and to for- mulate hypotheses with regard to the deepest units. The shallow Neogenic deposits and Fly- sch units are electrically conductive (10-20 Ω⋅m). The carbonate formations are only mod- erately resistive (usually around 100 Ω⋅m), whereas the anhydritic formations, interbedded with the carbonates, are occasionally very resis- tive, reaching values of 1000-1500 Ω⋅m. The resistivity of the metamorphic units measured in a few deep wells in the Larderello area are quite low (about 100 Ω⋅m) compared to those of the anhydrites and to the values measured on laboratory samples. This resistivity is probably representative of the upper, most fractured part of the metamorphic basement and a bulk value of many thousands of Ω⋅m is expected on the basis of the laboratory data, even accounting for the very high temperature and pressure condi- tions (Losito, 1991). It is difficult to define the resistivity values in geothermal areas accurately because of the presence of geothermal fluids; by definition warm and salty, these fluids decrease the resis- tivity of the formations affected by fluid circu- lation. Formations that are usually resistive, such as metamorphic formations, in these con- ditions can assume values that are not so dif- ferent from those of carbonatic formations. The presence of geothermal fluids also triggers dispersion phenomena over the frequency- dependent MT sounding data (Stoyer, 1976; Patella, 1987). This is not a problem when many electrical sounding data are available, and comparison of the high-frequency MT da- ta and the electrical data defining the shallow- er structures gives an idea of the dispersion rate (Patella et al., 1991). In our case, many electri- cal data were available in the Larderello and Amiata geothermal areas and the dispersion ef- fect was found to be negligible. Outside the ge- othermal area the dispersion effect produced by hydrothermal fluid circulation is also prob- ably negligible, as shown in Giammetti et al. (1996). Looking at the distribution of sites in fig. 1 it is evident that most of the MT sites are dis- tributed over profiles with a roughly WSW- ENE trend. This is because the main regional strike is defined by the Apenninic structure and the coast, i.e. ESE-WNW, and the profiles were designed orthogonal to the main strike. The profiles and the sparse sites are too widely sep- arated to allow a 3D modelling of the MT data. Therefore, 2D interpretation of the data have been attempted. Depending on inductive and geological length scales of the target, 2D inter- pretation may provide useful information of the features of the subsurface, provided that near surface galvanic distortion are removed (Ledo et al., 2002). The MT data along the profiles 113 Resistivity and heterogeneity of Earth crust in an active tectonic region, Southern Tuscany (Italy) shown in this paper have been modelled using a 2D inversion code, which use a regularized in- version algorithm that minimizes the sum of the chi-squared measure of data misfit and the squared norm of the Laplacian of the model function (Rodi and Mackie, 2001). Since TM-mode data are less affected by complex effects caused by near 3D structures (Wannamaker et al., 1984; and references there- in), we refer here to models obtained by inver- sion of TM-mode data only. Ledo et al. (2002) have discussed the applicability of this proce- dure and found results consistent with previous work, which suggest that the models obtained by inversion of TM-mode data only are more robust than those obtained using both polariza- tion. However, the resistivity of the conductive 3D bodies is overestimated and the dipping structures are not well defined. The MT data span a frequency range of about 300-0.001 Hz, with some slight changes due to the different equipment used during field- work. The MT data considered in this study were acquired in the region after 1991. Before 1991 the data were acquired in the Larderello area as single-site or were referenced to a local site at a distance of, at most, 3 km. The strong electromagnetic noise affecting the data for fre- quencies lower than 1 Hz cannot be separated from the MT signal in this case, and only shal- low information can be gathered for these data (Fiordelisi et al., 1998). The first MT data acquired in remote-refer- ence mode and processed with robust algo- rithms were from Larderello, in 1992 (Fiorde- lisi et al., 1995). Two profiles, for a total of 34 sites at an average distance of 2 km, crossed the main geothermal fields and the surrounding ar- eas (AA′ and BB′ in fig. 1). These data were re- cently reprocessed (Manzella et al., 2002) and the results obtained by 2D inversion of TM da- ta on profile AA′ are shown in fig. 4a. The rms misfit of the model is 10.5 with a floor error of 5% and 2% for resistivity and phase, respec- tively. The modelling shows, on both profiles AA′ and BB′ (the latter is not shown here) that the area is characterized by many low resistivi- ty anomalies. There is a good correlation be- tween the most productive and exploited geot- hermal fields and the low resistivity anomalies. The main anomalies show a very large vertical extension, and their location is comparable to those of the low seismic velocity anomalies de- fined by seismic tomography and described in Section 2 (fig. 4b). Since only TM-mode data were inverted, the resistivity of the anomalies is probably overestimated. Unfortunately, the lateral and vertical exten- sion of the anomalies, and hence of the reser- voirs and their connection to deeper features, does not emerge in great detail since the dis- tance between sites was quite large and the structure is very complex. The distance be- tween MT sites was found to be of critical im- portance for defining the extension of these low resistivity areas. As an example, we consider two models obtained inverting the MT data ac- quired in the Mt. Amiata water-dominated geot- hermal area on 1994 (fig. 5a,b), along a profile named CC′ in fig. 1 (Manzella et al., 1999). The model on the top was obtained by inverting the TM-mode MT data from just 26 sites at an av- erage distance of 1 km, and the one on the bot- tom by inverting the TM-mode MT data from 126 sites at an average distance of 200 m on the same profile. The rms misfit of the two models is 1.65 and 1.87, respectively. The two models define very different features of the subsurface, in particular at large depth. The one on the bot- tom is probably a better and more detailed rep- resentation of the subsurface, taking into ac- count, however, the limitation of the 2D inter- pretation of 3D structures. A spacing of just 1 km between MT sites is unable to resolve the lateral changes in resistivity occurring in an area with a sharp transition to lower resistivity, since this is caused by an abrupt increase in per- meability through fractures and faults. The da- ta, for the moment, are unable to predict the ex- act location and direction of the fractures at such great depths, which would be of consider- able interest in geothermal exploration, but they identify the areas that are affected by deep cir- culation of geothermal fluids. The correspon- dence between areas of low resistivity inside the resistive basement and geothermal reser- voirs is very evident in the Mt. Amiata water- dominated system, which was investigated by another continuous profiling survey, named DD′ in fig. 1 (Volpi et al., 2003). It is notewor- 114 Adele Manzella Fig. 5a,b. 2D inversion results in the Amiata area along profile CC′: a) TM mode inversion of 26 MT sites along a profile; b) TM mode inversion of 126 sites along the same profile. Fig. 6. 2D inversion results in the Amiata area along profile DD′. a b 115 Resistivity and heterogeneity of Earth crust in an active tectonic region, Southern Tuscany (Italy) thy that, below the Mt. Amiata exploited geot- hermal area, the resistivity is never as low as in Larderello at a depth of more than 7-8 km (figs. 5a,b and 6). The main decrease in resistivity, comparable to that of some sites in the Larde- rello area, is found outside the geothermal area, about 15 km SW of Mt. Amiata. Apart from the surveys in geothermal areas, MT data were also acquired within the frame- work of deep crustal exploration projects. Al- though many of these data are too affected by noise to be of any use, some proved to be good and were considered in this study. The profile EE′ in fig. 1 included many sites where noise was absent or very low (De Angelis et al., 1998). The model obtained by 2D forward modelling is shown in fig. 7, and depicts a con- tinuous increase of resistivity with depth. De- spite the clear limitation of the 2D model where the sites have an average distance of 10 km, it is evident from the figure that the most conductive formations, with resistivity lower than 50 Ω⋅m, are confined to depths of less than 2-3 km. The only exceptions to this general trend occur near the eastern border of the model, with an uplift of the Tuscan Nappe that appears reversed at the surface, and below the outcropping vol- canics (site 11 in fig. 7), where the resistive lay- er rises to a depth of about 8-9 km. This shallow resistive body was interpreted as an intrusion of hot material that rose in the past and cooled completely. Other data, acquired on sites located 20 km west of Mt. Amiata, show a very low, deep re- sistivity that is comparable only to that of some sites in the Larderello area and SW of Mt. Ami- ata (Manzella and Fiordelisi, 1995). In general we observe that MT data can be divided into three main categories, which are named only for representation purposes: 1) nor- mal sites, located outside the geothermal area, whose apparent resistivity at frequencies lower than 10−2 Hz (i.e. at greater depths) is higher than 100 Ω⋅m; 2) deep, low resistivity sites, characterized by an apparent resistivity lower than 10 Ω⋅m at frequencies lower than 10−2 Hz; and 3) geothermal sites, located in intensely ex- ploited geothermal areas, whose apparent resis- tivity at lower frequencies is lower than in nor- mal sites but higher than in low resistivity sites. The geographical distribution of the three types of sites is shown in fig. 8, as well as a sample curve for each of the three categories. Compar- ing fig. 8 with fig. 3a, we observe that the areas where the heat flow values are highest do not always correspond to the areas with the lowest resistivity below 8-10 km depth. The low resis- tivity data in the Larderello area show a good correspondence both with the maximum heat flow values and the location of the low velocity anomaly defined by seismic tomography. This evidence, together with other geophysical and geological information described in Section 2, can be interpreted as due to the presence of a still hot and partially molten magmatic body, representing the heat source of the Larderello geothermal area. With the actual data it is im- possible to say whether there is one, large mag- matic body only or many different bodies, pos- sibly interconnected. Geochemistry has shown that the helium values at Larderello are anom- alous, suggesting a mantle origin for part of the geothermal fluids, and hence depicting a direct connection of the magmatic body/-ies with the upper mantle (Magro et al., 2003). With the da- ta available at present, however, we are unable to verify this hypothesis. Fig. 7. 2D forward modelling results along profile EE′. 116 Adele Manzella Although its reservoir temperatures can be as high as 320°C, the Amiata area appears quite different from Larderello. Both the MT and the seismic data are unable to identify any clear re- duction in resistivity and seismic velocity be- low this area. This difference with the situation at Larderello could be explained by assuming that the magmatic body is deeper and maybe smaller than at Larderello. Nor does the geo- chemistry show any evidence of a mantle origin in the geothermal fluids, so the difference could also be due to a different connection with deep- er crustal and mantle structures. Fifteen kilometres SW of Amiata we identi- fied an area with deep, low resistivity and no ap- parent relation to the Amiata area itself (Man- zella et al., 1999). This area also shows an in- crease in heat flow values, but no seismic to- mography data are available. In the other area with low resistivity at depth, located 20 km west of Mt. Amiata, neither heat flow nor seismic to- mography data are available for comparison. Fig. 8. Distribution of the three main categories of MT sites (see text): dotted ovals indicate areas of very low deep resistivity (few Ω·m); slashed ovals correspond to areas of average deep resistivity (100-200 Ω·m) and densely lined ovals indicates areas of high deep resistivity (1000-2000 Ω·m). An example curve of the three main categories of sites is given on the bottom. 117 Resistivity and heterogeneity of Earth crust in an active tectonic region, Southern Tuscany (Italy) Magmatic bodies could be present below these areas, at depths similar to those at Larderello, but have never been identified before because of a lack of surface manifestations. The general behaviour of resistivity in the medium crust is another matter of some interest here. Many papers report a decrease of resistiv- ity in the medium-deep crust, which was inter- preted as being linked to the transition from brittle, shallow to ductile, deep conditions, and, possibly, the presence of fluids (see, e.g., Simp- son, 1999; and references therein). A region such as Southern Tuscany, so richly endowed in geothermal fluids, some of which are of deep origin, should be a perfect example of how deep crustal fluids can affect resistivity. As already mentioned in Section 2, the brittle-ductile tran- sition in this region is possibly evidenced by the seismic marker K and it probably follows the same shape of the K-horizon «surface» shown in fig. 3c at a slightly different depth. Hence, the transition was definitely explored by the MT data described in this study, since it is quite shallow. Indeed, no evidence of a decrease of resistivity can be seen at depths corresponding to the ductile condition of the crust. 4. Conclusions The MT data show evidence of an extreme- ly heterogeneous distribution of crustal resistiv- ity in Southern Tuscany. MT penetration depths of more than 4 km were achieved only where electromagnetic noise was absent, or could be eliminated. In general, the MT data show that below this depth the resistivity reaches values of only a few thousands of Ω⋅m, which is far lower than in metamorphic formations in other parts of the world. This is not unusual, howev- er, for a tectonically active area such as South- ern Tuscany, where the high reflectivity and the high heat flow could indicate the presence of fluids and/or melts on a fairly wide scale. Be- sides the resistivity anomalies created by local conditions at depths of 2-5 km, the MT data al- so reveal that the resistivity distribution is not homogeneous below 8-10 km depth, and that low resistivity subsurface does not show a sim- ple relation to the heat flow distribution. The deep, low resistivity anomalies were found below the geothermal area of Larderello and in areas located west and south-west of the Amiata area but outside the main geothermal area, whereas deep resistivity in the central Amiata area has been observed to increase. The MT data seem to indicate that the intrusions lie at greater depth in the Amiata area than in the Larderello area, and that intrusions comparable to those at Larderello may be present in areas that have still to be explored. The connections at depth between these many magmatic bodies has still to be investigated. A clear transition to low resistivity in the medium and deep crust is ab- sent in Southern Tuscany, suggesting that the ef- fects of temperature and fluid distribution do not depend on the transition to ductile conditions. More MT data are needed to obtain a clear- er picture in the region. Future experiments will include the acquisition of data at lower fre- quency than have been acquired till now, and longer recordings in order to obtain a better quality. Acknowledgements The author thanks Enel and the COFIN2000 Project for providing data and for partly fund- ing this research, and Marnell Dickson for proof-reading the original manuscript. REFERENCES BARELLI, A., G. BERTINI, G. BUONASORTE, G. CAPPETTI and A. FIORDELISI (2000): Recent deep exploration results at the margins of the Larderello Travale geothermal system, in Proceedings of the World Geothermal Congress 2000, Kyushu-Tohoku, Japan, 965-970. BATINI, F., A. FIORDELISI, F. GRAZIANO and M.N. TOKSOZ (1995): Earthquake tomography in the Larderello geot- hermal area, in Proceedings of the World Geothermal Congress 1995, Firenze, Italy, 817-820. BERNABINI, M., G. BERTINI, G.M. CAMELI, I. DINI and L. OR- LANDO (1995): Gravity interpretation of Mt. Amiata geot- hermal area (Central Italy), in Proceedings of the World Geothermal Congress 1995, Firenze, Italy, 859-862. CAMELI, G.M., I. DINI and D. LIOTTA (1998): Brittle/ductile boundary from seismic reflection lines of Southern Tuscany (Northern Apennines, Italy), Mem. Soc. Geol. It., 52, 153-163. CHIARABBA, C, A. AMATO and A. FIORDELISI (1995): Seis- micity and velocity images of the Roman Magmatic 118 Adele Manzella Province, in Proceedings of the World Geothermal Congress 1995, Firenze, Italy, 827-832. DE ANGELIS, M., A. FIORDELISI, A. MANZELLA and A. ZAJA (1998): Two-dimensional analysis of a magnetotelluric profile in the CROP 03 area in Southern Tuscany, Mem. Soc. Geol. It., 52, 295-304. FIORDELISI, A., R. MACKIE, T. MADDEN, A. MANZELLA and S. RIEVEN (1995): Application of the magnetotelluric method using a remote-remote reference system for characterizing deep geothermal system, in Proceedings of the World Geothermal Congress 1995, Firenze, Italy, 893-898. FIORDELISI, A., R. MACKIE, A. MANZELLA and A. ZAJA (1998): Electrical features of deep structures in South- ern Tuscany (Italy), Ann. Geofis., 41 (3), 333-341. FIORDELISI, A., A. MANZELLA, G. BUONASORTE, J. LARSEN and R. MACKIE (2000): MT methodology in the detec- tion of deep, water-dominated geothermal systems, in Proceedings of the World Geothermal Congress 2000, Kyushu-Tohoku, Japan, 1121-1126. GIAMMETTI, S., D. PATELLA, A. SINISCALCHI and A. TRA- MACERE (1996): The Siena Graben: combined interpre- tation of DES and MT soundings, Ann. Geofis., XXXIX (1), 189-200. GIANELLI, G., A. MANZELLA and M. PUXEDDU (1997): Crustal models of the geothermal areas of Southern Tuscany (Italy), Tectonophysics, 281, 221-239. LARSEN, J.C., R.L. MACKIE, A. MANZELLA, A. FIORDELISI and S. RIEVEN (1996): Robust smooth magnetotelluric transfer functions, Geophys. J. Int., 124, 801-819. LEDO, J., P. QUERALT, A. MARTÍ and G. JONES (2002): Two- dimensional interpretation of three-dimensional mag- netotelluric data: an example of limitations and resolu- tion, Geophys. J. Int., 150, 127-139. LOSITO, G. (1991): Proprietà elettriche delle rocce lungo il profilo CROP 03, Studi Geologici Camerti, 1, 33-40. MAGRO, G., G. RUGGIERI, G. GIANELLI and S. BELLANI (2003): Helium isotopes in paleofluids and present-day fluids of the Larderello geothermal field: constraints on the heat source, J. Geophys. Res., 108 (B1), doi: 10.1029/2001JB001590. MANZELLA, A. and A. FIORDELISI (1995): Resistivity struc- ture on the western side of the Mt. Amiata region, in Proceedings of the World Geothermal Congress 1995, Firenze, Italy, 881-886. MANZELLA, A., R. MACKIE and A. FIORDELISI (1999): MT survey in the Amiata volcanic area: a combined methodology for defining shallow and deep structures, Phys. Chem. Earth, A, 24 (9), 837-840. MANZELLA, A., G. VOLPI, F. CORSI and A. FIORDELISI (2002): Distribution of resistivity in the upper crust of the Larderello geothermal area, in Studi per l’Interpre- tazione del Profilo Sismico CROP18, II Stato di Avan- zamento, 16-24. PATELLA, D. (1987): Tutorial: interpretation of magnetotel- luric measurements over an electrically dispersionale one-dimensional Earth, Geophys. Prospect., 35, 1-11. PATELLA, D., A. TRAMACERE, R. DI MAIO and A. SINIS- CALCHI (1991): Experimental evidence of resistivity frequency-dispersion in magnetotellurics in Newbury (Oregon), Snake River Plain (Idaho) and Campi Flegrei (Italy) Volcano-geothermal areas, J. Volcanol. Geot- herm. Res., 48 (1-2), 61-75. RODI, W. and R.L. MACKIE (2001): Nonlinear conjugate gradients algorithm for 2D magnetotelluric inversion, Geophysics, 66, 174-187. SIMPSON, F. (1999): Stress and seismicity in the lower con- tinental crust: a challenge to simple ductility and im- plications for electrical conductivity mechanisms, Surv. Geophys., 20 (3-4), 2001-227. STOYER, C.H. (1976): Consequences of induced polariza- tion in magnetotelluric investigation, Pure Appl. Geo- phys., 114, 435-449. VOLPI, G., A. MANZELLA and A. FIORDELISI (2003). Investi- gation of geothermal structures by magnetotelluric (MT): an example from the Mt. Amiata area, Italy, Ge- othermics, 32 (2), 131-145. WANNAMAKER, P., G. HOHMANN and S. WARD (1984): Mag- netotelluric responses of three-dimensional bodies in layered earths, Geophysics, 49, 1517-1533.