In this paper, the Ground Penetrating Radar (GPR) method was used to characterize concrete tubes and steel/plastic tanks buried in IAG/USP test site. The microwave tomography was used to improve the GPR images, aiming...In this paper, the Ground Penetrating Radar (GPR) method was used to characterize concrete tubes and steel/plastic tanks buried in IAG/USP test site. The microwave tomography was used to improve the GPR images, aiming to retrieve the geometry of the targets. The numerical modeling studies also were done in order to predict the GPR results of the buried targets and to give more reliability to the results interpretation. The targets were installed in the first shallow geophysical test site of the Brazil located at Institute of Astronomy, Geophysics, and Atmospheric Science (IAG) of the University of S?o Paulo (USP). GPR profiles of 200 MHz (shielded bistatic antennas) were acquired along three lines containing concrete tubes and steel/plastic tanks buried in subsoil. The concrete tubes show a hyperbolic reflector for the top, and the vertical tube also presented a reflection on its bottom. The horizontal steel tanks were characterized by a strong GPR reflection on their top. The empty plastic tank shows a strong reflector for the top with normal polarity. On the other hand, the plastic tank filled with water shows a weaker reflector for its top characterized by the inverted polarity of GPR signal when compared with empty plastic tank. The plastic tank filled with water also went characterized by the strong reflection to its bottom, being a good indicative to interpret GPR data on target in subsoil with some types of fluid inside of tank. The results of polarity difference for the top of tank can be used as guide pattern to identify buried tank empty or filled with water. The application of microwave tomography to the GPR data permitted to determine the position and get a good identification of the edges of the targets studied. The numeric modeling presented a good accordance with real data reducing the ambiguities in interpretation of results. These results can be used as a reference, and they can be extrapolated for areas where there is no subsurface information.展开更多
The calculation of apparent resistivity in electrical methods is relatively simple and consolidated in the literature. Furthermore, in the case of DC methods, the apparent resistivity values are for the most part, int...The calculation of apparent resistivity in electrical methods is relatively simple and consolidated in the literature. Furthermore, in the case of DC methods, the apparent resistivity values are for the most part, intuitive and enable an initial interpretation of the results. On the other hand, in TEM method the apparent resistivity values are not very intuitive and interpretation based only on the apparent resistivity values cannot be done reasonably. In this way, this paper presents a discussion about the reasons and effects of the negatives values of apparent resistivity in TEM soundings. The main objective of this paper is to clarify the meanings of these negatives values and the possible effects in TEM inversion.展开更多
Identifying underground utilities and predicting their depth are fundamental when it comes to civil engineering excavations, for example, to install or repair water, sewer, gas, electric systems and others. The accide...Identifying underground utilities and predicting their depth are fundamental when it comes to civil engineering excavations, for example, to install or repair water, sewer, gas, electric systems and others. The accidental rupture of these systems can lead to unplanned repair costs, delays in completing the service, and risk injury or death of workers. One way to detect underground utilities is using the GPR-Ground Penetrating Radar geophysical method. To estimate depth, the travel time (two-way travel time) information provided by a radargram is used in conjunction with ground wave velocity, which depends on the dielectric constant of materials, where it is usually assumed to be constant for the area under investigation. This procedure provides satisfactory results in most cases. However, wrong depth estimates can result in damage to public utilities, rupturing pipes, cutting lines and so on. These cases occur mainly in areas that have a marked variation of water content and/or soil lithology, thus greater care is required to determine the depth of the targets. The present work demonstrates how the interval velocity of Dix (1955) can be applied in radargram to estimate the depth of underground utilities compared to the conventional technique of constant velocity applied to the same data set. To accomplish this, synthetic and real GPR data were used to verify the applicability of the interval velocity technique and to determine the accuracy of the depth estimates obtained. The studies were carried out at the IAG/USP test site, a controlled environment, where metallic drums are buried in known positions and depths allowing the comparison of real to estimated depths. Numerical studies were also carried out aiming to simulate the real environment with variation of dielectric constant in depth and to validate the results with real data. The results showed that the depths of the targets were estimated more accurately by means of the interval velocity technique in contrast to the constant velocity technique, minimizing the risks of accidents during excavation.展开更多
We present here a series of Ground Penetrating Radar (GPR) survey carried out in different areas of the Dardanelos 1 archaeological site in order to generate information about subsurface anomalies associated with arch...We present here a series of Ground Penetrating Radar (GPR) survey carried out in different areas of the Dardanelos 1 archaeological site in order to generate information about subsurface anomalies associated with archaeological material that could be use in decision making within the environmental licensing process of the Dardanelos Hydroelectric Power Plant, located near of the Aripuan? city, Mato Grosso State, northwest region of Brazil. GPR surveys with 200 MHz antenna were carried out in two blocks aiming to locate archaeological resources and features. The analysis of GPR 2D and 3D results allowed detecting anomalous regions characterized by hyperbolic reflections, shallow elongated continuous targets with high amplitudes, as well as sub-horizontal reflectors. Microwave tomography allowed estimating the geometry of the GPR anomalies sources. Excavations were done by archaeologists at the locations where hyperbolic anomalies were found, revealing interesting structures related to urns in the middle of a rich in organic matter consisting of black and ceramic materials up to about 1 m deep. The first sub-horizontal reflector at approximately 1 m depth is related to the base of the black soil layer rich in organic matter and the second sub-horizontal reflector between 2 and 3 m deep suggests a lithological change or may be related to presence of the water table. The continuous elongated shallow targets observed in the depth slices are related to tree roots in the middle of the archaeological strata. The GPR results guided archaeological excavations, reduced the time and costs involved in research, and contributed to the preservation of Brazilian historical heritage.展开更多
Despite representing one of the largest cratons on Earth,the early geological evolution of the Amazonia Craton remains poorly known due to relatively poor exposure and because younger metamorphic and tectonic events h...Despite representing one of the largest cratons on Earth,the early geological evolution of the Amazonia Craton remains poorly known due to relatively poor exposure and because younger metamorphic and tectonic events have obscured initial information.In this study,we investigated the sedimentary archives of the Carajás Basin to unravel the early geological evolution of the southeastern Amazonia Craton.The Carajás Basin contains sedimentary rocks that were deposited throughout a long period spanning more than one billion years from the Mesoarchean to the Paleoproterozoic.The oldest archives preserved in this basin consist of a few ca.3.6 Ga detrital zircon grains showing that the geological roots of the Amazonia Craton were already formed by the Eoarchean.During the Paleoarchean or the early Mesoarchean(<3.1 Ga),the Carajás Basin was large and rigid enough to sustain the formation and preservation of the Rio Novo Group greenstone belt.Later,during the Neoarchean,at ca.2.7 Ga,the southeastern Amazonia Craton witnessed the emplacement of the Parauapebas Large Igneous Province(LIP)that probably covered a large part of the craton and was associated with the deposition of some of the world largest iron formations.The emplacement of this LIP immediately preceded a period of continental extension that formed a rift infilled first by iron formations followed by terrigenous sediments.This major change of sedimentary regime might have been controlled by the regional tectonic evolution of the Amazonia Craton and its emergence above sea-level.During the Paleoproterozoic,at ca.2.1 Ga,the Rio Fresco Group,consisting of terrigenous sediments from the interior of the Amazonia Craton,was deposited in the Carajás Basin.At that time,the Amazonian lithosphere could have either underwent thermal subsidence forming a large intracratonic basin or could have been deformed by long wavelength flexures that induced the formation of basins and swells throughout the craton under the influence of the growing Transamazonian mountain belt.展开更多
In general, micropulsation is generated by the free energy within the magnetosphere, which is mainly caused by the magnetosphere-solar wind interactions and/or magnetic field reconnections in the magnetospheric tail. ...In general, micropulsation is generated by the free energy within the magnetosphere, which is mainly caused by the magnetosphere-solar wind interactions and/or magnetic field reconnections in the magnetospheric tail. In this article, we studied the types of low-frequency magnetospheric micropulsations Pc5 and Pc6 (1 - 6 mHz), during three magnetic storms, using ground magnetic observations. One of the main motivations of this study was to produce scientific knowledge on the subject in order to fill knowledge gaps in this region, collaborating with existing bibliographies. The data were recorded in a series of six Brazilian stations close to the Dip Equator (dip = 0), with using data from the Ascension Island observatory. Pulse trains were found to occur simultaneously with almost identical waveforms, with significant amplitude increase near magnetic local noon. Amplitude spectra showed a clear similarity in all equatorial seasons. The results obtained were compared with studies performed in regions with similar low latitudes. The possible mechanisms that induced these micropulsations and the peculiarities presented by these data were also discussed in this article.展开更多
This research applies Electrical Resistivity Tomography (ERT) and Time Domain Electromagnetic Method (TDEM) to study the hydrogeology of the Taubaté basin, which is characterized by half-grabens with about 850 m ...This research applies Electrical Resistivity Tomography (ERT) and Time Domain Electromagnetic Method (TDEM) to study the hydrogeology of the Taubaté basin, which is characterized by half-grabens with about 850 m of maximum sediments thickness. The study area is in Taubaté city, São Paulo State, Brazil, where the Taubaté aquifer is an important water source. The Taubaté Group is the main sedimentary package of the basin;it is formed mainly by shales that form aquicludes, and thin layers of sandstones that form the aquifer. There are 40 groundwater exploration wells in Taubaté city that provide important information. The study purpose is to characterize the geoelectrical stratigraphy of the subsurface to locate the contact between the Quaternary and Tertiary sediments and to identify the Taubaté aquifer. The ERT is used for shallow investigations (tens of meters) and the TDEM can reach a great investigation depth (hundreds of meters). Therefore, these geophysical methods are complementary. The ERT data were acquired with the pole-dipole array with 20 m of electrodes spacing and 400 m length, and the TDEM data with the central-loop array with a 200 × 200 m transmitter loop. The results permit to define the contact between the Quaternary and Tertiary sediments around 15 m depth, the Pindamonhangaba Formation between 15 m and 30 m depth and the Taubate Group between 30 m and 300 m depth. The TDEM method defined the Taubaté Group as a single geoelectric layer because the shale and the sandstone layers are all very conductive. The basement is formed by gneiss, which is a very resistive rock. The TDEM method is not able to identify a high conductor/resistor contrast. Overall, the results are consistent with the known geology and the wells information.展开更多
This work uses 2D TEM (Transient Electromagnetic) modeling for a hydrogeological study in the Paraná sedimentary basin. The study area is located at the northern region of the state of S?o Paulo, Brazil, where gr...This work uses 2D TEM (Transient Electromagnetic) modeling for a hydrogeological study in the Paraná sedimentary basin. The study area is located at the northern region of the state of S?o Paulo, Brazil, where groundwater is exploited from two aquifer systems: one sedimentary, shallow, and the other crystalline, deep. The interest in applying the TEM method in this area owes to the high exploitation rates of groundwater from the crystalline aquifer system for irrigation, which is triggering considerable seismic activity locally. This aquifer system is composed of fractured basalt within the Serra Geral Formation and is about 120 m deep. Eighty-six TEM soundings were acquired at this location, but in nine cases the data did not fit the modelled curve for 1D geoelectrical models due to the geological complexity of the area. This paper shows 2D geoelectrical modeling results based on the FDTD (Finite Differences in Time Domain) method to explain the lateral resistivity variation within the geological setting. A 2D model was generated for each sounding and compared with 1D inversion models as well as with direct information from wells. The results show some vertical variations of about 10 to 30 meters on the upper interface of the basalt layer from Serra Geral Formation. They are located at approximately 60 meters from the center of the soundings. The existence of these 2D structures in the subsurface can be related to the drainage system in the study area. The presence of these structures may indicate a connection between the shallow and deep aquifer systems, acting like a conduit that may contribute to the seismic activity reported.展开更多
Double seismic zones are commonly observed in the subducting slabs in a global scale,serving as ideal examples for studying the seismogenetic mechanism of the intermediate-depth earthquakes.In this study,we relocate e...Double seismic zones are commonly observed in the subducting slabs in a global scale,serving as ideal examples for studying the seismogenetic mechanism of the intermediate-depth earthquakes.In this study,we relocate earthquakes and determine seismic velocity models using the double-difference seismic tomography method in the northern Chile subduction zone where a double seismic zone exists.The results suggest that the double seismic zone in northern Chile is located at about 50-140 km depth,with an interval of approximately 20 km between the two zones.The upper seismic zone is characterized by relatively low Vp(~7.8-8.0 km/s),low Vs(~4.4-4.5 km/s)and high Vp/Vs(~1.85)above the depth of~90 km,while the region below~90 km is distinguished by relatively high Vp(~8.2 km/s),high Vs(~4.8 km/s)and slightly high Vp/Vs(~1.75),which may be related to a series of dehydration reactions of hydrous minerals in the subducted oceanic crust.In comparison,the lower seismic zone is featured by the anomaly of low Vp/Vs(~1.7),although some local areas may consist of relatively high Vp/Vs values(~1.8),possibly due to the dehydration reaction of serpentine.Based on the Vp,Vs,Vp/Vs anomalies combined with previous petrological experiments and thermodynamic models,it can be derived that intermediate-depth earthquakes are mainly related to the dehydration of various hydrous minerals in the subducting slab.The dehydration process of hydrous minerals releases water into the subducting slab and subsequently leads to the increase of pore fluid pressure and the decrease of effective normal stress,thus causing the occurrence of brittle failure and intermediate-depth earthquakes in subduction zones.The imaging results of the northern Chile subduction zone further indicate that the existence of the double seismic zone is related to the dehydration process of different hydrous minerals.展开更多
文摘In this paper, the Ground Penetrating Radar (GPR) method was used to characterize concrete tubes and steel/plastic tanks buried in IAG/USP test site. The microwave tomography was used to improve the GPR images, aiming to retrieve the geometry of the targets. The numerical modeling studies also were done in order to predict the GPR results of the buried targets and to give more reliability to the results interpretation. The targets were installed in the first shallow geophysical test site of the Brazil located at Institute of Astronomy, Geophysics, and Atmospheric Science (IAG) of the University of S?o Paulo (USP). GPR profiles of 200 MHz (shielded bistatic antennas) were acquired along three lines containing concrete tubes and steel/plastic tanks buried in subsoil. The concrete tubes show a hyperbolic reflector for the top, and the vertical tube also presented a reflection on its bottom. The horizontal steel tanks were characterized by a strong GPR reflection on their top. The empty plastic tank shows a strong reflector for the top with normal polarity. On the other hand, the plastic tank filled with water shows a weaker reflector for its top characterized by the inverted polarity of GPR signal when compared with empty plastic tank. The plastic tank filled with water also went characterized by the strong reflection to its bottom, being a good indicative to interpret GPR data on target in subsoil with some types of fluid inside of tank. The results of polarity difference for the top of tank can be used as guide pattern to identify buried tank empty or filled with water. The application of microwave tomography to the GPR data permitted to determine the position and get a good identification of the edges of the targets studied. The numeric modeling presented a good accordance with real data reducing the ambiguities in interpretation of results. These results can be used as a reference, and they can be extrapolated for areas where there is no subsurface information.
文摘The calculation of apparent resistivity in electrical methods is relatively simple and consolidated in the literature. Furthermore, in the case of DC methods, the apparent resistivity values are for the most part, intuitive and enable an initial interpretation of the results. On the other hand, in TEM method the apparent resistivity values are not very intuitive and interpretation based only on the apparent resistivity values cannot be done reasonably. In this way, this paper presents a discussion about the reasons and effects of the negatives values of apparent resistivity in TEM soundings. The main objective of this paper is to clarify the meanings of these negatives values and the possible effects in TEM inversion.
文摘Identifying underground utilities and predicting their depth are fundamental when it comes to civil engineering excavations, for example, to install or repair water, sewer, gas, electric systems and others. The accidental rupture of these systems can lead to unplanned repair costs, delays in completing the service, and risk injury or death of workers. One way to detect underground utilities is using the GPR-Ground Penetrating Radar geophysical method. To estimate depth, the travel time (two-way travel time) information provided by a radargram is used in conjunction with ground wave velocity, which depends on the dielectric constant of materials, where it is usually assumed to be constant for the area under investigation. This procedure provides satisfactory results in most cases. However, wrong depth estimates can result in damage to public utilities, rupturing pipes, cutting lines and so on. These cases occur mainly in areas that have a marked variation of water content and/or soil lithology, thus greater care is required to determine the depth of the targets. The present work demonstrates how the interval velocity of Dix (1955) can be applied in radargram to estimate the depth of underground utilities compared to the conventional technique of constant velocity applied to the same data set. To accomplish this, synthetic and real GPR data were used to verify the applicability of the interval velocity technique and to determine the accuracy of the depth estimates obtained. The studies were carried out at the IAG/USP test site, a controlled environment, where metallic drums are buried in known positions and depths allowing the comparison of real to estimated depths. Numerical studies were also carried out aiming to simulate the real environment with variation of dielectric constant in depth and to validate the results with real data. The results showed that the depths of the targets were estimated more accurately by means of the interval velocity technique in contrast to the constant velocity technique, minimizing the risks of accidents during excavation.
文摘We present here a series of Ground Penetrating Radar (GPR) survey carried out in different areas of the Dardanelos 1 archaeological site in order to generate information about subsurface anomalies associated with archaeological material that could be use in decision making within the environmental licensing process of the Dardanelos Hydroelectric Power Plant, located near of the Aripuan? city, Mato Grosso State, northwest region of Brazil. GPR surveys with 200 MHz antenna were carried out in two blocks aiming to locate archaeological resources and features. The analysis of GPR 2D and 3D results allowed detecting anomalous regions characterized by hyperbolic reflections, shallow elongated continuous targets with high amplitudes, as well as sub-horizontal reflectors. Microwave tomography allowed estimating the geometry of the GPR anomalies sources. Excavations were done by archaeologists at the locations where hyperbolic anomalies were found, revealing interesting structures related to urns in the middle of a rich in organic matter consisting of black and ceramic materials up to about 1 m deep. The first sub-horizontal reflector at approximately 1 m depth is related to the base of the black soil layer rich in organic matter and the second sub-horizontal reflector between 2 and 3 m deep suggests a lithological change or may be related to presence of the water table. The continuous elongated shallow targets observed in the depth slices are related to tree roots in the middle of the archaeological strata. The GPR results guided archaeological excavations, reduced the time and costs involved in research, and contributed to the preservation of Brazilian historical heritage.
基金funded by grants of the Fundação AmparoàPesquisa do Estado de São Paulo(FAPESP2015/16235-2,2017/18840-6,2018/02645-2,2018/14617-3,2018/05892-0,2019/17732-0,2019/16066-7 and 2019/12132-5)+2 种基金the Conselho Nacional de Desenvolvimento Científico e Tecnológico(CNPq308045/2013-0 and 307353/2019-2)the Fundação AmparoàPesquisa do Minas Gerais(FAPEMIG project APQ-03793-16)。
文摘Despite representing one of the largest cratons on Earth,the early geological evolution of the Amazonia Craton remains poorly known due to relatively poor exposure and because younger metamorphic and tectonic events have obscured initial information.In this study,we investigated the sedimentary archives of the Carajás Basin to unravel the early geological evolution of the southeastern Amazonia Craton.The Carajás Basin contains sedimentary rocks that were deposited throughout a long period spanning more than one billion years from the Mesoarchean to the Paleoproterozoic.The oldest archives preserved in this basin consist of a few ca.3.6 Ga detrital zircon grains showing that the geological roots of the Amazonia Craton were already formed by the Eoarchean.During the Paleoarchean or the early Mesoarchean(<3.1 Ga),the Carajás Basin was large and rigid enough to sustain the formation and preservation of the Rio Novo Group greenstone belt.Later,during the Neoarchean,at ca.2.7 Ga,the southeastern Amazonia Craton witnessed the emplacement of the Parauapebas Large Igneous Province(LIP)that probably covered a large part of the craton and was associated with the deposition of some of the world largest iron formations.The emplacement of this LIP immediately preceded a period of continental extension that formed a rift infilled first by iron formations followed by terrigenous sediments.This major change of sedimentary regime might have been controlled by the regional tectonic evolution of the Amazonia Craton and its emergence above sea-level.During the Paleoproterozoic,at ca.2.1 Ga,the Rio Fresco Group,consisting of terrigenous sediments from the interior of the Amazonia Craton,was deposited in the Carajás Basin.At that time,the Amazonian lithosphere could have either underwent thermal subsidence forming a large intracratonic basin or could have been deformed by long wavelength flexures that induced the formation of basins and swells throughout the craton under the influence of the growing Transamazonian mountain belt.
文摘In general, micropulsation is generated by the free energy within the magnetosphere, which is mainly caused by the magnetosphere-solar wind interactions and/or magnetic field reconnections in the magnetospheric tail. In this article, we studied the types of low-frequency magnetospheric micropulsations Pc5 and Pc6 (1 - 6 mHz), during three magnetic storms, using ground magnetic observations. One of the main motivations of this study was to produce scientific knowledge on the subject in order to fill knowledge gaps in this region, collaborating with existing bibliographies. The data were recorded in a series of six Brazilian stations close to the Dip Equator (dip = 0), with using data from the Ascension Island observatory. Pulse trains were found to occur simultaneously with almost identical waveforms, with significant amplitude increase near magnetic local noon. Amplitude spectra showed a clear similarity in all equatorial seasons. The results obtained were compared with studies performed in regions with similar low latitudes. The possible mechanisms that induced these micropulsations and the peculiarities presented by these data were also discussed in this article.
文摘This research applies Electrical Resistivity Tomography (ERT) and Time Domain Electromagnetic Method (TDEM) to study the hydrogeology of the Taubaté basin, which is characterized by half-grabens with about 850 m of maximum sediments thickness. The study area is in Taubaté city, São Paulo State, Brazil, where the Taubaté aquifer is an important water source. The Taubaté Group is the main sedimentary package of the basin;it is formed mainly by shales that form aquicludes, and thin layers of sandstones that form the aquifer. There are 40 groundwater exploration wells in Taubaté city that provide important information. The study purpose is to characterize the geoelectrical stratigraphy of the subsurface to locate the contact between the Quaternary and Tertiary sediments and to identify the Taubaté aquifer. The ERT is used for shallow investigations (tens of meters) and the TDEM can reach a great investigation depth (hundreds of meters). Therefore, these geophysical methods are complementary. The ERT data were acquired with the pole-dipole array with 20 m of electrodes spacing and 400 m length, and the TDEM data with the central-loop array with a 200 × 200 m transmitter loop. The results permit to define the contact between the Quaternary and Tertiary sediments around 15 m depth, the Pindamonhangaba Formation between 15 m and 30 m depth and the Taubate Group between 30 m and 300 m depth. The TDEM method defined the Taubaté Group as a single geoelectric layer because the shale and the sandstone layers are all very conductive. The basement is formed by gneiss, which is a very resistive rock. The TDEM method is not able to identify a high conductor/resistor contrast. Overall, the results are consistent with the known geology and the wells information.
文摘This work uses 2D TEM (Transient Electromagnetic) modeling for a hydrogeological study in the Paraná sedimentary basin. The study area is located at the northern region of the state of S?o Paulo, Brazil, where groundwater is exploited from two aquifer systems: one sedimentary, shallow, and the other crystalline, deep. The interest in applying the TEM method in this area owes to the high exploitation rates of groundwater from the crystalline aquifer system for irrigation, which is triggering considerable seismic activity locally. This aquifer system is composed of fractured basalt within the Serra Geral Formation and is about 120 m deep. Eighty-six TEM soundings were acquired at this location, but in nine cases the data did not fit the modelled curve for 1D geoelectrical models due to the geological complexity of the area. This paper shows 2D geoelectrical modeling results based on the FDTD (Finite Differences in Time Domain) method to explain the lateral resistivity variation within the geological setting. A 2D model was generated for each sounding and compared with 1D inversion models as well as with direct information from wells. The results show some vertical variations of about 10 to 30 meters on the upper interface of the basalt layer from Serra Geral Formation. They are located at approximately 60 meters from the center of the soundings. The existence of these 2D structures in the subsurface can be related to the drainage system in the study area. The presence of these structures may indicate a connection between the shallow and deep aquifer systems, acting like a conduit that may contribute to the seismic activity reported.
基金We are grateful for two anonymous reviewers for their constructive comments.This research is supported by National Natural Science Foundation of China under grant(41861134009)the PIA ANID grant(AFB180004)the ANID grant(PCI PII-180003).
文摘Double seismic zones are commonly observed in the subducting slabs in a global scale,serving as ideal examples for studying the seismogenetic mechanism of the intermediate-depth earthquakes.In this study,we relocate earthquakes and determine seismic velocity models using the double-difference seismic tomography method in the northern Chile subduction zone where a double seismic zone exists.The results suggest that the double seismic zone in northern Chile is located at about 50-140 km depth,with an interval of approximately 20 km between the two zones.The upper seismic zone is characterized by relatively low Vp(~7.8-8.0 km/s),low Vs(~4.4-4.5 km/s)and high Vp/Vs(~1.85)above the depth of~90 km,while the region below~90 km is distinguished by relatively high Vp(~8.2 km/s),high Vs(~4.8 km/s)and slightly high Vp/Vs(~1.75),which may be related to a series of dehydration reactions of hydrous minerals in the subducted oceanic crust.In comparison,the lower seismic zone is featured by the anomaly of low Vp/Vs(~1.7),although some local areas may consist of relatively high Vp/Vs values(~1.8),possibly due to the dehydration reaction of serpentine.Based on the Vp,Vs,Vp/Vs anomalies combined with previous petrological experiments and thermodynamic models,it can be derived that intermediate-depth earthquakes are mainly related to the dehydration of various hydrous minerals in the subducting slab.The dehydration process of hydrous minerals releases water into the subducting slab and subsequently leads to the increase of pore fluid pressure and the decrease of effective normal stress,thus causing the occurrence of brittle failure and intermediate-depth earthquakes in subduction zones.The imaging results of the northern Chile subduction zone further indicate that the existence of the double seismic zone is related to the dehydration process of different hydrous minerals.