首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 890 毫秒
1.
We present 23 new ages from three volcanic complexes of the Lesser Antilles arc in Martinique Island (French West Indies). These ages obtained with the K–Ar Cassignol–Gillot technique are distributed within the whole Quaternary. They allowed us to reconstruct a detailed history of successive volcanic growth and flank collapse stages. Trois Ilets Volcanism has been active during at least 2 Ma, between 2.35 ± 0.03 Ma and 346 ± 27 ka, with monogenetic volcanoes of basaltic-andesite to andesitic compositions. We here propose that magma mixing, which characterizes this volcanism, could have been initiated between 617 and 346 ka by the activation of arc-parallel and arc-transverse fault systems. Meanwhile, the Carbet complex was active 25 km to the north from 998 ± 14 to 322 ± 6 ka, and was partially destroyed by a flank collapse after 602 ± 10 ka. Together with geochemical data, our ages show that Mount Conil and Mount Pelée volcanoes are parts of the same edifice sharing a single magmatic reservoir. Mount Conil started to emerge before 543 ± 8 ka, and andesites erupted until 127 ± 2 ka, when a flank collapse destroyed the western flank of the edifice, probably triggering the emplacement of Piton Marcel, the last eruption of this first stage. We note that this collapse occurred during the transition from oxygen stages 6 to 5, i.e. during glacial to interglacial change, when eustatic level rapidly increased. After that, and until present, Mount Pelée volcano was built with periods of cone growth intercalated by flank collapse events. We here show that a peak of activity occurred between 550 and 330 ka in western Martinique within the three complexes, which are spaced of 15–25 km. Since 330 ka volcanic activity is limited to the northernmost Mount Conil–Mount Pelée complex. Our data are in agreement with the regional scale observations that the whole recent Lesser Antilles arc was subject to a high volcanic activity since 600 ka, probably linked to an increase in magma production. This permanent establishment of rising magma in regularly spaced batches and tectonically controlled, could explain the individual chemical evolution of each edifice and the different eruptive dynamisms occurring at the same time along the recent arc.  相似文献   

2.
Magmas are transported through pre-existing fractures in many repeatedly erupting volcanoes. The study of this special process of magma transport is fundamentally important to understand the mechanisms and conditions of volcanic eruptions. In this paper, we numerically simulate the magma propagation process through a pre-existing vertical fracture in the crust by using the combined finite difference method (FDM), finite element method (FEM) and discontinuous deformation analysis (DDA) approach. FDM is used to analyze magma flow in the pre-existing fracture, FEM is used to calculate the opening of the fracture during magma intrusion, and DDA is used to deal with the contact of the closed fracture surfaces. Both two-dimensional (2D) and three-dimensional (3D) examples are presented. Parametric studies are carried out to investigate the influence of various physical and geometric parameters on the magma transport in the pre-existing fracture. We have considered magma chamber depth ranging from 7 km to 10 km under the crust surface, magma viscosity ranging from 2 × 10−2 to 2 × 10−7 MPa s, and the density difference between the magma and host rock ranging from 300 to 700 kg/m3. The numerical results indicate that (1) the fluid pressure p varies gradually along the depth, (2) the shape of the magma body during propagation is like a torch bar and its width ranges from 2 m to 4 m approximately in the 3D case and 10 m to 50 m in the 2D case for the same physical parameters used, (3) the crust surface around the pre-existing fracture begins to increase on both sides of the fracture, forms a trough between them, then gradually uplifts during the transport of the magma, and finally takes the shape of a crater when the magma reaches the surface. We have also examined the influence of physical and geometric parameters on the minimum overpressure for magma transport in the 3D case. The numerical results show that our numerical technique presented in this paper is an effective tool for simulating magma transport process through pre-existing fractures in the crust.  相似文献   

3.
Coloumbo submarine volcano lies 6.5 km offshore the NE part of the Santorini island complex and exhibits high seismicity along with vigorous hydrothermal activity. This study models the local stress field around Coloumbo's magma chamber and investigates its influence on intrusion emplacement and geometry. The two components of the stress field, hoop and radial stress, are calculated using analytical formulas that take into account the depth and radius of the magma chamber as these are determined from seismological and other observations. These calculations indicate that hoop stress at the chamber walls is maximum at an angle of 74° thus favouring flank intrusions, while the radial stress switches from tensile to compressive at a critical distance of 5.7 km from the center of the magma chamber. Such estimates agree well with neotectonic and seismological observations that describe the local/regional stress field in the area. We analyse in detail the case where a flank intrusion reaches the surface very near the NE coast of Thera as this is the worst-case eruption scenario. The geometrical features of such a feeder dyke point to an average volumetric flow rate of 9.93 m3 s−1 which corresponds to a Volcanic Explosivity Index of 3 if a future eruption lasts about 70 days. Hazards associated with such an eruption include ashfall, ballistic ejecta and base surges due to explosive mixing of magma with seawater. Previous studies have shown that areas near erupting vents are also foci of moderate to large earthquakes that precede or accompany an eruption. Our calculations show that a shallow event (3–5 km) of moment magnitude 5.9 near the eruptive vent may cause Peak Ground Acceleration in the range 122–177 cm s−2 at different locations around Santorini. These values indicate that seismic hazard even due to a moderate earthquake near Coloumbo, is not trivial and may have a significant impact especially on older buildings at Thera island.  相似文献   

4.
We invert measurements of coseismic displacements from 139 continuously recorded GPS sites from the 2010, Jiashian, Taiwan earthquake to solve for fault geometry and slip distribution using an elastic uniform stress drop inversion. The earthquake occurred at a depth of ~ 23 km in an area between the Western Foothills fold-and-thrust belt and the crystalline high mountains of the Central Range, providing an opportunity to examine the deep fault structure under Taiwan. The inferred rupture plane is oblique to the prominent orientation of thrust faults and parallel to several previously recognized NW-striking transfer zones that appear to connect stepping thrusts. We find that a fault striking 318°–344° with dip of 26°–41° fits the observations well with oblique reverse-sinistral slip under a low stress drop of about 0.5 MPa. The derived geodetic moment of 2.92 × 1018 N-m is equivalent to a Mw = 6.24 earthquake. Coseismic slip is largely concentrated within a circular patch with a 10-km radius at the depth between 10 and 24 km and maximum slip of 190 mm. We suggest this earthquake ruptured the NW-striking Chishan transfer fault zone, which we interpret as a listric NE-dipping lateral ramp with oblique slip connecting stepping thrust faults (ramps). The inferred slip on the lateral ramp is considerably deeper than the 7–15 km deep detachment identified in previous studies of western Taiwan. We infer an active basal detachment under western Taiwan at a depth of at least ~ 20–23 km based on these inversion results. The earthquake may have nucleated at the base of the lateral ramp near the intersection with the basal detachment. Coulomb stress change calculations suggest that this earthquake moved several NE-striking active thrust faults in western Taiwan nearer to failure.  相似文献   

5.
The South Sandwich volcanic arc is sited on a young oceanic crust, erupts low-K tholeiitic rocks, is characterized by unexotic pelagic and volcanogenic sediments on the down-going slab, and simple tectonic setting, and is ideal for assessing element transport through subduction zones. As a means of quantifying processes attending transfer of subduction-related fluids from the slab to the mantle wedge, boron concentrations and isotopic compositions were determined for representative lavas from along the arc. The samples show variable fluid-mobile/fluid-immobile element ratios and high enrichments of B/Nb (2.7 to 55) and B/Zr (0.12 to 0.57), similar to those observed in western Pacific arcs. δ11B values are among the highest so far reported for mantle-derived lavas; these are highest in the central part of the arc (+ 15 to + 18‰) and decrease toward the southern and northern ends (+ 12 to + 14‰). δ11B is roughly positively correlated with B concentrations and with 87Sr/86Sr ratios, but poorly coupled with other fluid-mobile elements such as Rb, Ba, Sr and U. Peridotites dredged from the forearc trench also have high δ11B (ca. + 10‰) and elevated B contents (38–140 ppm). Incoming pelagic sediments sampled at ODP Site 701 display a wide range in δ11B (+ 5 to ? 13‰; average = ? 4.1‰), with negative values most common. The unusually high δ11B values inferred for the South Sandwich mantle wedge cannot easily be attributed to direct incorporation of subducting slab materials or fluids derived directly therefrom. Rather, the heavy B isotopic signature of the magma sources is more plausibly explained by ingress of fluids derived from subduction erosion of altered frontal arc mantle wedge materials similar to those in the Marianas forearc. We propose that multi-stage recycling of high-δ11B and high-B serpentinite (possibly embellished by arc crust and volcaniclastic sediments) can produce extremely 11B-rich fluids at slab depths beneath the volcanic arc. Infiltration of such fluids into the mantle wedge likely accounts for the unusual magma sources inferred for this arc.  相似文献   

6.
Lengshuiqing is part of the late Proterozoic igneous province from the western margin of the Yangtze craton. The Lengshuiqing area comprises five ultramafic–mafic intrusions, emplaced in the late Proterozoic Yanbian Group. The intrusions from Lengshuiqing contain cumulate ultramafic zones (peridotite + olivine pyroxenite), with cumulus olivine and Cr-spinel, and intercumulus pyroxenes, hornblende, phlogopite and plagioclase. Ni–Cu ore (pyrrhotite + pentlandite + chalcopyrite) is hosted in the ultramafic zones. Olivine-free diorite–quartz diorite ± gabbro and granite zones commonly occur above the ultramafic rocks. The genesis of the intrusions (conduit-related accumulation or differentiation from stagnant magma) was investigated. The amount of sulphides in the intrusions from Lengshuiqing is one order of magnitude bigger than the sulphides that can be dissolved by a volume of mafic magma similar with the volume of the intrusions. Most intrusions from Lengshuiqing have bulk composition (peridotite ± diorite ± granite) more magnesian (MgO = 21–22%; Mg# > 78) than the deduced composition of their parental magma (MgO = 9–11%; Mg# = 64–67). This indicates the accumulation of sulphide and mafic silicates from a volume of magma much bigger than the volume of the intrusions, which can be explained by the fractionation from magma ascending through the intrusions to shallower depths. A continuous supply and vent of magma is consistent with the lack of chilled margins, the melting of the wall rocks and the generation of high-temperature mineral assemblages (K-feldspar, diopside, and sillimanite) in the Yanbian Group. The intrusions from Lengshuiqing are seen as microchambers on conduits draining olivine-, Cr-spinel-, and sulphide-bearing mafic magma from a larger staging chamber.  相似文献   

7.
The 1.0 Ma Kidnappers supereruption (~ 1200 km3 DRE) from Mangakino volcanic centre, Taupo Volcanic Zone, New Zealand, produced a large phreatomagmatic fall deposit followed by an exceptionally widespread ignimbrite. Detailed sampling and analysis of glass shards and mineral phases have been undertaken through a proximal 4.0 m section of the fall deposit, representing the first two-thirds of erupted extra-caldera material. Major and trace element chemistries of glass shards define three distinct populations (types A, B and C), which systematically change in proportion through the fall deposit and are inferred to represent three magma types. Type B glass and biotite first appear at the same level (~ 0.95 m above base) in the fall deposit suggesting later tapping of a biotite-bearing magma. Plagioclase and Fe–Ti oxide compositions show bimodal distributions, which are linked to types A and B glass compositions. Temperature and pressure (T–P) estimates from hornblende and Fe–Ti oxide equilibria from each magma type are similar and therefore the three magma bodies were adjacent, not vertically stacked, in the crust. Most hornblende model T–P estimates range from 770 to 840 °C and 90 to 170 MPa corresponding to storage depths of ~ 4.0–6.5 km. Hornblende model T–P estimates coupled with in situ trace element fingerprinting imply that the magma bodies were individually well mixed, and not stratified. Compositional gaps between the three glass compositional types imply that no mixing between these magmas occurred. We interpret these data, coupled with the systematic changes in shard compositional proportions through the fall deposit, to reflect that three independent melt-dominant bodies of magma contributed large (A, ~ 270 km3), medium (B, ~ 90 km3) and small (C, ~ 40 km3) volumes (as reflected in the fall deposits) and were systematically tapped during the eruption. We propose that the systematic evacuation of the three independent magma bodies implies that there was tectonic triggering and linkage of eruptions. Our results show that supereruptions can be generated by near simultaneous multiple eruptions from independent magma chambers rather than the evacuation of a large single unitary magma chamber.  相似文献   

8.
The relation of magma and crustal activity has been studied from spatial distribution of 3He/4He ratios of gas and/or water samples over the Izu Peninsula, where significant crustal deformation associated with seismic swarm activities has been observed since 1970s. The air-corrected values of 3He/4He ratios ranged from 3.5 to 8.2 RA, where RA is the atmospheric 3He/4He ratio = 1.4 × 10? 6, indicating that helium is mostly of magmatic origin. Among the three pressure sources proposed to explain the crustal deformation, two inflation sources beneath the inland of northeast and the mid east coast of the Izu Peninsula locate in the broad distribution of high 3He/4He ratios, which supports relation of magma to the crustal uplift. In contrast, the distribution of 3He/4He ratios around the tensile fault assumed in the area of seismic swarms appears not to indicate existence of significant amount of magma below the tensile fault. Alternatively, the results suggest magma below a point several kilometers south of the tensile fault. The seismic swarms are explained either by fluid pressurization of thermal water heated by this magma or by intrusion of magma to the tensile fault moved obliquely from the deep magma reservoir.  相似文献   

9.
Serpentinization of the mantle wedge is an important process that influences the seismic and mechanical properties in subduction zones. Seismic detection of serpentines relies on the knowledge of elastic properties of serpentinites, which thus far has not been possible in the absence of single-crystal elastic properties of antigorite. The elastic constants of antigorite, the dominant serpentine at high-pressure in subduction zones, were measured using Brillouin spectroscopy under ambient conditions. In addition, antigorite lattice preferred orientations (LPO) were determined using an electron back-scattering diffraction (EBSD) technique. Isotropic aggregate velocities are significantly lower than those of peridotites to allow seismic detection of serpentinites from tomography. The isotropic VP/VS ratio is 1.76 in the Voigt–Reuss–Hill average, not very different from that of 1.73 in peridotite, but may vary between 1.70 and 1.86 between the Voigt and Reuss bonds. Antigorite and deformed serpentinites have a very high seismic anisotropy and remarkably low velocities along particular directions. VP varies between 8.9 km s? 1 and 5.6 km s? 1 (46% anisotropy), and 8.3 km s? 1 and 5.8 km s? 1 (37%), and VS between 5.1 km s? 1 and 2.5 km s? 1 (66%), and 4.7 km s? 1 and 2.9 km s? 1 (50%) for the single-crystal and aggregate, respectively. The VP/VS ratio and shear wave splitting also vary with orientation between 1.2 and 3.4, and 1.3 and 2.8 for the single-crystal and aggregate, respectively. Thus deformed serpentinites can present seismic velocities similar to peridotites for wave propagation parallel to the foliation or lower than crustal rocks for wave propagation perpendicular to the foliation. These properties can be used to detect serpentinite, quantify the amount of serpentinization, and to discuss relationships between seismic anisotropy and deformation in the mantle wedge. Regions of high VP/VS ratios and extremely low velocities in the mantle wedge of subduction zones (down to about 6 and 3 km.s?1 for VP and VS, respectively) are difficult to explain without strong preferred orientation of serpentine. Local variations of anisotropy may result from kilometer-scale folding of serpentinites. Shear wave splittings up to 1–1.5 s can be explained with moderately thick (10–20 km) serpentinite bodies.  相似文献   

10.
We elucidate the ecology of Recent Ostracoda from a deep brackish lake, Tangra Yumco (30°45′—31°22′N and 86°23′—86°49′E, 4595 m a.s.l.) and adjacent waters on the southern Tibetan Plateau. Ostracod associations (living and empty valves) in sixty-six sediment samples collected from diverse aquatic habitats (lakes, estuary-like water and lagoon-like water waters, rivers, ponds and springs) were quantitatively assessed.Eleven Recent Ostracoda were found (nine living and two as empty valves only). Cluster analysis established two significant (p < 0.05) habitat specific associations; (i) Leucocytherella sinensis, Limnocythere inopinata, Leucocythere? dorsotuberosa, Fabaeformiscandona gyirongensis and Candona xizangensis are lacustrine fauna. (ii) Tonnacypris gyirongensis, Candona candida, Ilyocypris sp., Heterocypris incongruens and Heterocypris salina are temporary water species.Ostracod distribution and abundance are significantly (p < 0.05) correlated to physico-chemical variables. The first two axes of a canonical correspondence analysis (CCA) explain 30.9% of the variation in the species abundance data. Conductivity and habitat types are the most influential ecological factors explaining the presence and abundance of ostracods. Spearman correlation analysis reveals that: (i) Two species, L.? dorsotuberosa (r = 0.25) and L. inopinata (r = 0.36) have a significant positive correlation with conductivity while one species, T. gyirongensis (r = −0.68) displays a significant negative correlation with conductivity. Limnocythere inopinata correlates significantly positive (r = 0.37) with alkalinity. Fabaeformiscandona gyirongensis correlates significantly positive (r = 0.28) with water depth.Key indicator living assemblages are: (i) L. sinensis dominates Ca-depleted brackish waters although ubiquitously distributed; (ii) L.? dorsotuberosa dwells in fresh to brackish waters; (iii) L. inopinata predominates in mesohaline to polyhaline waters; (iv) F. gyirongensis inhabits exclusively brackish-lacustrine deeper waters; (v) C. candida populates freshwaters; (vi) T. gyirongensis and Ilyocypris sp. are restricted to shallow temporary waters; (vii) H. incongruens occurs in ponds.Water depth indicators are F. gyirongensis and L.? dorsotuberosa, useful in ostracod assemblages for palaeo-water depth reconstruction.Our results expand the knowledge of the ecological significance of Recent Tibetan Ostracoda ecology. This is a new insight on habitat chacteristics of both living assemblages and sub-Recent associations of ostracods in mountain aquatic ecosystems. The new modern ostracod dataset can be used for the quantitative reconstruction of past environmental variables (e.g., conductivity) and types of water environment. The key indicator ostracods are relevant in palaeolimnological and climate research on the Tibetan Plateau.  相似文献   

11.
Garnet geochronology was used to provide the first direct measurement of the timing of eclogitization in the central Himalaya. Lu–Hf dates from garnet separates in one relict eclogite from the Arun River Valley in eastern Nepal indicate an age of 20.7 ± 0.4 Ma, significantly younger than ultra-high pressure eclogites from the western Himalaya, reflecting either different origins or substantial time lags in tectonics along strike. Four proximal garnet amphibolites from structurally lower horizons are 14–15 Ma, similar to post-eclogitization ages published for rocks along strike in southern Tibet. P–T calculations indicate three metamorphic episodes for the eclogite: i) eclogite-facies metamorphism at ～ 670 °C and ≥ 15 kbar at 23–16 Ma; ii) a peak-T granulite event at ～ 780 °C and 12 kbar; and iii) late-stage amphibolite-facies metamorphism at ～ 675 °C and 6 kbar at ～ 14 Ma. The garnet amphibolites were metamorphosed at ～ 660 °C. Three models are considered to explain the observed P–T–t evolution. The first assumes that the Main Himalayan Thrust (basal thrust of the Himalayan thrust system) cuts deeper at Arun than elsewhere. While conceptually the simplest, this model has difficulty explaining both the granulite-facies overprint and the pulse of exhumation between 25 and 14 Ma. A second model assumes that (aborted) subduction, slab breakoff, and ascent of India's leading edge occurred diachronously: ～ 50 Ma in the western Himalaya, ～ 25 Ma in the central Himalaya of Nepal, and presumably later in the eastern Himalaya. This model explains the P–T–t path, particularly heating during initial exhumation, but implies significant along-strike diachroneity, which is generally lacking in other features of the Himalaya. A third model assumes repeated loss of mantle lithosphere, first by slab breakoff at ～ 50 Ma, and again by delamination at ～ 25 Ma; this model explains the P–T–t path, but requires geographically restricted tectonic behavior at Arun. The P–T–t history of the Arun eclogites may imply a change in the physical state of the Himalayan metamorphic wedge at 16–25 Ma, ultimately giving rise to the Main Central Thrust by 15–16 Ma.  相似文献   

12.
This study reconstructs the late stages in the exhumation history of a nappe derived from the Meliatic accretionary wedge in the Western Carpathians by means of zircon and apatite (U-Th)/He dating. The Meliatic accretionary wedge formed due to the closure of the Neotethyan Triassic–Jurassic Meliata–Hallstatt Ocean in the Late Jurassic. The studied fragments of the blueschist-bearing Meliatic Bôrka Nappe were metamorphosed at low-temperature and high- to medium-pressure conditions at ca. 160–150 Ma and included into the accretionary wedge. The time of the accretionary wedge formation constrains the beginning of the Bôrka Nappe northward thrusting over the Gemeric Unit of the evolving Central Western Carpathians (CWC) orogenic wedge. The zircon (U-Th)/He data on four samples recorded three evolutionary stages: (i) cooling through the ∼180 °C isotherm at 130–120 Ma related to starting collapse of the accretionary wedge, following exhumation of the high-pressure slices in the Meliatic accretionary wedge; (ii) postponed exhumation and cooling of some fragments through the ∼180 °C isotherm from 115 to 95 Ma due to ongoing collapse of this wedge; and (iii) cooling from 80 to 65 Ma, postdating the thrusting (∼100–80 Ma) of the Bôrka Nappe slices during the Late Cretaceous compression related to formation of the CWC orogenic wedge. The third stage already documents cooling of the Meliatic Bôrka Nappe slices in the CWC orogenic wedge. The apatite (U-Th)/He data may indicate cooling of a Bôrka Nappe slice to near-surface temperatures at ∼65 Ma. The younger AHe age clusters indicate that at least one, or possibly two, reheating events could have occurred in the longer interval from ∼40 to ∼10 Ma during the Oligocene–Miocene. These were related to sedimentary burial and/or the magmatism as documented in other parts of the CWC.  相似文献   

13.
Magnetic Resonance Sounding (MRS) is nowadays accepted as a new geophysical method that can be used for a reliable determination of the ground water content distribution in the top 150 m. A great effort has also been made in MRS development to deduce the hydraulic transmissivity, based on empiric relationships of the permeability with a factor F which is calculated with NMR parameters measured at laboratory scale. To use this relationship under field conditions a calibration coefficient CT = Tpt / F has to be previously established, which demands the knowledge of the transmissivity Tpt evaluated in the pumping test. The transmissivity can then be calculated at any other site of the same aquifer using the relation Tmrs = CTF. The CT values reported suggest a certain relationship with the lithology, but with a great dispersion and contradictory results. MRS surveys carried out in alluvial aquifers in Spain have shown that the value of CT evaluated at one site may not be valid at another place of the same aquifer, because of the great heterogeneity of this kind of geological environment. The demand of a pumping test at each site where a MRS is measured invalidates the method actually used for MRS transmissivity evaluation. More than 50 MRS have been used to propose a new methodology. The aquifers visited cover a great range of transmissivities (from 2 × 10− 6 to 9 × 10− 3 m2/s). The MRS signal amplitude varies between 20 and 1400 nV, the signal/noise ratio is in the range from 0.6 to 42, and the value of the decay time constant varies from 200 to 800 ms. It has been demonstrated that when the transmissivity increases, the value of F decreases, and CT increases, except for certain groups of MRS taken at the same aquifer or part of one aquifer, for which F increases with Tpt, keeping CT constant. A function CT(F) of the type CT = mF− n has been obtained that allows the transmissivity evaluation without the need of Tpt. Considering that both values of transmissivity, Tpt and Tmrs, are subjected to deviations due to the experimental errors as well as due to evaluation errors, the prediction achieved by the proposed equation is rather good. To perform a better evaluation of the values of the coefficients m and n it is necessary to have a greater number of MR soundings of good quality and with a trustworthy inversion at locations where a really comparable and good performed pumping test is available, covering a sufficient range of transmissivities. Though the data we have used do not always fulfil these conditions, the result is promising. Once a trustable function is available, the forecast of the transmissivity using MRS will not need the existence of any pumping test in the area. The general extension of this methodology demands the availability of MRS taken at all kinds of geological and hydrogeological environments, which is impossible without the existence of a universal MRS data base.  相似文献   

14.
Laboratory experiments on simulated faults in rocks clearly show the temperature dependence of dynamic rock friction. Since rocks surrounding faults are permeable, we have developed a numerical method to describe the thermo-mechanical evolution of the pre-seismic sliding phase which takes into account both the rate-, state- and temperature-dependent friction law and the heat advection term in the energy equation. We consider a laminar fluid motion perpendicular to a vertical fault plane and assume that fluids move away from the fault plane. A semi-analytical temperature solution which accounts for the variability of slip velocity and stress on the fault has been found. This solution has been generalized to the case of a time varying fluid velocity and then was used to include the thermal pressurization effect. After discretizing the temperature solution, the evolution of the system is obtained by the solution of a system of first order differential equations which allows us to determine the evolution of slip, slip rate, friction coefficient, effective normal stress, temperature and fluid velocity. The numerical solutions are found using a Runge-Kutta method with an adaptative stepsize control in time. When the thermal pressurization effects can be neglected, the heat advection effect gives rise to a delay, with respect to the purely conductive case, of the earthquake occurrence time. This delay increases with increasing permeability H of the system. When the thermal pressurization effects are taken into account the situation is opposite, i.e. the onset of instability tends to precede that of the purely conductive case. The advance in the time of occurrence of instability increases with increasing coefficient of thermal pressurization. In the small permeability range (H ≤ 10?18 m2), the seismic moment and nucleation length of the pre-seismic phase are significantly smaller than those predicted by the purely conductive model.  相似文献   

15.
We have measured magnetic hysteresis loops, zero-field-cooled (ZFC) and field-cooled (FC) remanence, and low-field AC susceptibility as a function of temperature between 2 and 40 K for a single crystal several mm in size and for two powders of manganese carbonate (mineral rhodochrosite, MnCO3), one ground from a natural precipitate (grainsize ∼100 μm) and another synthesized in the laboratory (grainsize ∼10 μm). For the single crystal, measurements carried out both in the basal (easy magnetization) plane and along the trigonal (hard magnetization) axis yielded, expectedly, grossly different magnetic properties. In the basal plane, hysteresis appears to be mostly controlled by domain wall movement at the two lowest temperatures studied, 5 and 15 K, as indicated by a fairly broad switching field distribution. At 25 K and above, however, magnetization reversal occurs at a single, well defined magnetic field, which we interpret as a characteristic field of the in-plane magnetic anisotropy. Hysteresis in the basal plane is observed up to 36 K which is above the nominal Néel temperature of rhodochrosite (34.3 K). In addition, a sharp coercivity peak occurs at 34.5 K. Rather unexpectedly, hysteresis is also observed for the magnetic field applied along the trigonal axis. It is very small at 5 K but develops gradually with increasing temperature, coercivity reaching maximum of 100 mT at 28 K and remanence peaking at slightly higher temperature (30–31 K). Hysteresis along the trigonal axis is observed up to 37 K. Hysteresis temperature dependence conforms with the AC susceptibility versus temperature curve which shows a maximum at 36.5 K. ZFC/FC remanence curves also closely match the temperature dependence of remanence extracted from hysteresis loops. We suggest that this behavior could be due to the presence of a minor, about 1 at.% amount of Fe2+ substituting for Mn in the crystalline lattice of rhodochrosite. Hysteresis measurements on powders have revealed a significant enhance in coercivity, up to 50 mT for the 100-μm powder and up to 150 mT for the 10-μm one. FC/ZFC ratio amounts to about 2 for the natural powder, while for the synthetic one, which is essentially pure material, it barely exceeds unity. FC/ZFC ratio can thus be viewed as a sensitive indicator of iron incorporation into rhodochrosite.  相似文献   

16.
Two mineralogically and chemically distinct rhyolite magmas (T1 and T3) were syn-erupted from the same conduit system during the 21.9 ka basalt intrusion-triggered Okareka eruption from Tarawera volcano, New Zealand. High spatial resolution U–Th disequilibrium dating of zircon crystals at the ~ 3–5 μm scale reveals a protracted yet discontinuous zircon crystallization history within the magmatic system. Both magma types contain zircon whose interiors predate the eruption by up to 200 ka. The dominant age peak in the T1 magma is ~ 30 ka with subordinate peaks at ~ 45, ~ 75, and ~ 100 ka, whereas the T3 magma has a dominant zircon interior age peak at ~ 90 ka with smaller modes at ~ 35 and ~ 150 ka. These patterns are consistent with isolated pockets of crystallization throughout the evolution of the system. Crystal rim analyses yield ages ranging from within error of the eruption age to at least ~ 90 ka prior to eruption, highlighting that zircon crystallization frequently stalled long before the eruption. Continuous depth profiling from crystal rims inward demonstrates protracted growth histories for individual crystals (up to ~ 100 ka) that were punctuated by asynchronous hiatuses of up to 30 ka in duration. Disparate zircon growth histories can result from localized thermal perturbations caused by mafic intrusions into a silicic reservoir. The crystal age heterogeneity at hand-sample scale requires considerable crystal transport and mixing. We propose that crystal mixing was achieved through buoyancy instabilities caused by mafic magma flow through crystal mush. A terminal pre-eruptive rejuvenation event was capable of mobilizing voluminous melts that erupted, but was too short (< 102–103 years) to result in extensive zircon growth. The contrasting, punctuated zircon histories argue against closed-system fractional crystallization models for silicic magmatism that require protracted cooling times following a mostly liquid starting condition.  相似文献   

17.
We report new high-precision laser fluorination three-isotope oxygen data for lunar materials. Terrestrial silicates with a range of δ18O values (− 0.5 to 22.9‰) were analyzed to independently determine the slope of the terrestrial fractionation line (TFL; λ = 0.5259 ± 0.0008; 95% confidence level). This new TFL determination allows direct comparison of lunar oxygen isotope systematics with those of Earth. Values of Δ17O for Apollo 12, 15, and 17 basalts and Luna 24 soil samples average 0.01‰ and are indistinguishable from the TFL. The δ18O values of high- and low-Ti lunar basalts are distinct. Average whole-rock δ18O values for low-Ti lunar basalts from the Apollo 12 (5.72 ± 0.06‰) and Apollo 15 landing sites (5.65 ± 0.12‰) are identical within error and are markedly higher than Apollo 17 high-Ti basalts (5.46 ± 0.11‰). Evolved low-Ti LaPaz mare-basalt meteorite δ18O values (5.67 ± 0.05‰) are in close agreement with more primitive low-Ti Apollo 12 and 15 mare basalts. Modeling of lunar mare-basalt source composition indicates that the high- and low-Ti mare-basalt mantle reservoirs were in oxygen isotope equilibrium and that variations in δ18O do not result from fractional crystallization. Instead, these differences are consistent with mineralogically heterogeneous mantle sources for mare basalts, and with lunar magma ocean differentiation models that result in a thick feldspathic crust, an olivine–pyroxene-rich mantle, and late-stage ilmenite-rich zones that were convectively mixed into deeper portions of the lunar mantle. Higher average δ18O (WR) values of low-Ti basalts compared to terrestrial mid ocean ridge basalts (Δ=0.18‰) suggest a possible oxygen isotopic difference between the terrestrial and lunar mantles. However, calculations of the δ18O of lunar mantle olivine in this study are only 0.05‰ higher than terrestrial mantle olivine. These observations may have important implications for understanding the formation of the Earth–Moon system.  相似文献   

18.
《Journal of Geodynamics》2007,43(2):320-329
A 2.5-month long gravity sequence, encompassing the starting period of the 2002–2003 Etna eruption and coming from a summit station only 1 km away from the new fractures, is presented and discussed. The sequence comprises four hours-long anomalies that have a great chance to reflect mass redistributions linked to the ensuing activity. In particular, the start of the eruptive activity on the northeastern flank was marked by a gravity decrease as strong as about 400 μGal, which reverted soon afterwards. This strong decrease/increase anomaly is interpreted as the opening, by tectonic forces, of a fracture system along the Northeastern Rift of Mt. Etna, followed by an intrusion of magma from the central conduit to the new fractures. They were used by the intruding magma as a path to the eruptive vents at lower elevations.Afterwards, on three occasions, in November and December 2002, 6–12 h-lasting gravity decreases, with amplitude ranging between 10 and 30 μGal, were observed simultaneously with increases in the amplitude of the volcanic tremor from four seismic stations. A correlation analysis, between the gravity signal and the overall spectral amplitude of each tremor sequence is performed over the 7 November–9 December period. A marked anti-correlation is found over each contemporaneous gravity decrease/tremor increase, while, over the rest of the investigated period, the correlation is negligible. Accordingly, a joint source is inferred to have acted during the occurrence of the three common anomalies. On the grounds of some volcanological observations spanning the period covered by our analysis, we propose the temporary accumulation of a gas cloud at some level within the plumbing system of the volcano to have acted as a joint source.The present work is a further evidence of the potential of continuous gravity observations as a tool to monitor and study active volcanoes and encourages their employment in spite of the difficulty of running spring gravimeters in a continuous fashion under the adverse conditions normally encountered on the summit zone of an active volcano.  相似文献   

19.
Understanding the processes at the origin of explosive events is crucial for volcanic hazard mitigation, especially during long-lasting eruptions at andesitic volcanoes. This work exposes the case of Tungurahua volcano, whose unrest occurred in 1999. Since this date, the eruptive activity was characterized by low-to moderate explosiveness, including phases with stronger canon-like explosions and regional ash fallout. However, in 2006, a sudden increase of the explosiveness led to pyroclastic flow-forming eruptions on July 14th (VEI 2) and August 16–17th (VEI 3). All magmas emitted from 1999 to 2005, as well as the samples from the 2006 eruptions, have homogeneous bulk-rock andesitic compositions (58–59 wt.% SiO2), and contain the same mineral assemblage consisting of pl + cpx + opx + mag ± ol. However, during the August 16–17th event, the erupted tephra comprise two types of magmas: a dominant, brown andesitic scoria; and scarce, light-grey pumice representing a subordinate, silica-rich juvenile component. For the andesitic magma, thermobarometric data point to magmatic temperatures ranging from 950 to 1015 °C and pressures in the range of 200 to 250 MPa, which corresponds to 7.5–9.5 km below the summit. Disequilibrium textures in plagioclase and pyroxene phenocrysts, particularly thin overgrowth rims, indicate the recharge of this magma body by mafic magma. Between 1999 and 2005, repeated injections from depth fed the intermittent eruptive activity observed while silica-rich melts were produced by in-situ crystallization in the peripheral parts of the reservoir. In April 2006, the recharge of a primitive magma produced strong convection and homogenisation in the reservoir, as well as pressure increase and higher magma ascent rate after seven years of only moderately explosive activity. This work emphasizes the importance of petrological studies in constraining the pre-eruptive magmatic conditions and processes, as a tool for understanding the fundamental causes of the changes in the eruptive dynamism, particularly the occurrence of paroxysmal phases in andesitic systems with open-vent behaviour.  相似文献   

20.
We present new analyses of volatile, major, and trace elements for a suite of glasses and melt inclusions from the 85°E segment of the ultra-slow spreading Gakkel Ridge. Samples from this segment include limu o pele and glass shards, proposed to result from CO2-driven explosive activity. The major element and volatile compositions of the melt inclusions are more variable and consistently more primitive than the glass data. CO2 contents in the melt inclusions extend to higher values (167–1596 ppm) than in the co-existing glasses (187–227 ppm), indicating that the melt inclusions were trapped at greater depths. These melt inclusions record the highest CO2 melt concentrations observed for a ridge environment. Based on a vapor saturation model, we estimate that the melt inclusions were trapped between seafloor depths (～ 4 km) and ～ 9 km below the seafloor. However, the glasses are all in equilibrium with their eruption depths, which is inconsistent with the rapid magma ascent rates expected for explosive activity. Melting conditions inferred from thermobarometry suggest relatively deep (25–40 km) and cold (1240°–1325 °C) melting conditions, consistent with a thermal structure calculated for the Gakkel Ridge. The water contents and trace element compositions of the melt inclusions and glasses are remarkably homogeneous; this is an unexpected result for ultra-slow spreading ridges, where magma mixing is generally thought to be less efficient based on the assumption that steady-state crustal magma chambers are absent in these environments. All melts can be described by a single liquid line of descent originating from a pooled melt composition that is consistent with the aggregate melt calculated from a geodynamic model for the Gakkel Ridge. These data suggest a model in which deep, low degree melts are efficiently pooled in the upper mantle (9–20 km depth), after which crystallization commences and continues during ascent and eruption. Based on our melting model and the assumption that CO2 is perfectly incompatible, we show that the highest CO2 concentrations of the melt inclusions (～ 1600 ppm) are consistent with the calculated CO2 concentrations of primary undegassed melts. The highest measured CO2/Nb ratio (443) of Gakkel Ridge melt inclusions predicts a mantle CO2 content of 134 ppm and would result in a global ridge flux of 2.0 × 1012 mol CO2/yr.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号