Hydration- and dehydration-induced rheological heterogeneities on the deep subduction interface, and possible relationships to episodic tremor and slow slip

Author(s):  
Whitney Behr ◽  
Carolyn Tewksbury-Christle ◽  
Alissa Kotowski ◽  
Claudio Cannizzaro ◽  
Robert Blass ◽  
...  

<p>Episodic tremor and slow slip (ETS) is observed in several subduction zones down-dip of the locked megathrust, and may provide clues for preparatory processes before megathrust rupture. Exhumed rocks provide a unique opportunity to evaluate the sources of rheological heterogeneity on the subduction interface and their potential role in generating ETS-like behavior. We present data from two subduction interface shear zones representative of the down-dip extent of the megathrust: the Condrey Mountain Schist (CMS) in northern CA (greenschist to blueschist facies conditions) and the Cycladic Blueschist Unit (CBU) on Syros Island, Greece (blueschist to eclogite facies). Both complexes highlight the propensity for fluid-mediated metamorphic reactions to produce strong rheological heterogeneities:</p><p>In the CMW, hydration reactions led to progressive serpentinization of peridotite bodies that were entrained from the overriding plate and underplated along with oceanic-affinity sediments. The margins of each peridotite-serpentinite lens show extreme strain localization accommodated by dislocation glide and minor pressure solution in antigorite, whereas lens interiors show evidence for more distributed, alternating, frictional-viscous deformation, with abundant crack-seal veins occupied by antigorite, brucite and oxides that are in some places also ductilely sheared. Deformation in the surrounding metasedimentary matrix was purely viscous.</p><p>In the CBU on Syros Island, dehydration reactions in MORB-affinity basalts, subducted and underplated with oceanic and continental-affinity sediments, led to progressive development of strong eclogitic lenses within a weaker blueschist and metasedimentary matrix. The eclogite lenses are commonly coarse-grained and massive and show brittle deformation in the form of dilational and shear fractures/veins filled with quartz, white mica, glaucophane and/or chlorite. Brittle deformation in the eclogites is coeval with ductile deformation in the surrounding blueschist and metasedimentary matrix, indicating concurrent frictional-viscous flow.</p><p>Although we cannot easily distinguish transient deformation processes in exhumed rocks, we can use the following three approaches to assess whether these heterogeneities could have generated deformation behaviors similar to deep ETS: 1) We measure displacements within, and dimensions of the heterogeneities in outcrop/map-scale to estimate the maximum possible seismic moment that would be released when the frictional heterogeneities slip;  2) We compare deformation mechanisms inferred from field and microstructural observations to their expected mechanical behavior from rock deformation experiments; and 3) We use seismo-thermo-mechanical modeling to examine expected slip velocities and moment-duration ratios for frictional-viscous shear zones that are scaled to observations from nature and the lab.  </p><p>All three approaches suggest that frictional-viscous heterogeneities of the types and length-scales we observe in the exhumed rock record are compatible with ETS as documented in modern subduction zones.</p>

Geosphere ◽  
2019 ◽  
Vol 15 (4) ◽  
pp. 1038-1065 ◽  
Author(s):  
Alissa J. Kotowski ◽  
Whitney M. Behr

Abstract We use structural and microstructural observations from exhumed subduction-related rocks exposed on Syros Island (Cyclades, Greece) to provide constraints on the length scales and types of heterogeneities that occupy the deep subduction interface, with possible implications for episodic tremor and slow slip. We selected three Syros localities that represent different oceanic protoliths and deformation conditions within a subduction interface shear zone, including: (1) prograde subduction of oceanic crust to eclogite facies; (2) exhumation of oceanic crust from eclogite through blueschist-greenschist facies; and (3) exhumation of mixed mafic crust and sediments from eclogite through blueschist-greenschist facies. All three localities preserve rheological heterogeneities that reflect metamorphism of primary lithological, geochemical, and/or textural variations in the subducted protoliths and that take the form of brittle pods and lenses within a viscous matrix. Microstructural observations indicate that the matrix lithologies (blueschists and quartz-rich metasediments) deformed by distributed power-law viscous flow accommodated by dislocation creep in multiple mineral phases. We estimate bulk shear zone viscosities ranging from ∼1018 to 1020 Pa-s, depending on the relative proportion of sediments to (partially eclogitized) oceanic crust. Eclogite and coarse-grained blueschist heterogeneities within the matrix preserve multiple generations of dilational shear fractures and veins formed under high-pressure conditions. The veins commonly show coeval or overprinting viscous shear, suggesting repeated cycles of frictional and viscous strain. These geologic observations are consistent with a mechanical model of episodic tremor and slow slip (ETS), in which the deep subduction interface is a rheologically heterogeneous distributed shear zone comprising transiently brittle (potentially tremor-genic) sub-patches within a larger, viscously creeping interface patch. Based on our observations of outcrop and map areas of heterogeneous patches and the sizes, distributions, and amounts of brittle offset recorded by heterogeneities, we estimate that simultaneous brittle failure of heterogeneities could produce tremor bursts with equivalent seismic moments of 4.5 × 109–4.7 × 1014 N m, consistent with seismic moments estimated from geophysical data at active subduction zones.


2021 ◽  
Author(s):  
Carolyn Tewksbury-Christle ◽  
Alissa Kotowski ◽  
Whitney Behr

<p>The strength, or viscosity, of the subduction interface is a key parameter in subduction dynamics, influencing both long-term subduction plate speeds and short-term transient deformation styles. Fossil subduction interfaces exhumed from downdip of the megathrust record ductile deformation accommodated by diverse lithologies, including metasedimentary and metamafic rocks. Existing flow laws for quartz-rich rocks predict relatively low viscosities, in contrast to high viscosities predicted for basalt and eclogite, but the rheological properties of blueschists representative of metamorphosed oceanic crust of the down-going slab are poorly constrained. Two key questions remain: 1) are there significant viscosity contrasts between blueschists and quartz- or mica-rich metasedimentary rocks, and 2) what are the microscale mechanisms for creep in naturally deformed blueschists and how do they vary with pressure and temperature? To address these questions, we characterized deformation in natural samples from the Condrey Mountain Schist (CMS) in northern California, USA, and the Cycladic Blueschist Unit (CBU) on Syros Island, Cyclades, Greece, using outcrop-scale structural observations, optical microscopy, and Electron Backscatter Diffraction. The CMS and CBU record pressure-temperature conditions of 0.8-1.1 GPa, 350-450°C and 1.4-1.8 GPa, 450-550°C, respectively. </p><p>In the field, blueschists form m- to km-scale lenses that are interfolded with quartz schists, ultramafics, and, in the CBU, eclogites and marbles. At the outcrop scale in both localities, quartz-rich schists and blueschists each exhibit strong foliations and lineations and planar contacts at lithological boundaries. At the thin section scale, the prograde foliation and mineral lineation in blueschists are commonly defined by Na-amphiboles elongated in the lineation direction. Crystallographic preferred orientations in Na-amphibole in all samples have c-axes parallel to lineation and a-axes predominantly defining point-maxima perpendicular to the foliation, suggesting some component of dislocation activity for all temperature conditions in our sample suite. Microtextures in lower temperature CMS samples suggest strain accommodation primarily by dislocation glide and kinking in Na-amphibole, with extremely high-aspect-ratio grains and limited evidence for climb-controlled dynamic recrystallization. Some higher temperature CBU samples show large porphyroclasts with apparent ‘core-and-mantle’-type recrystallization textures and subgrain orientation analyses consistent with the (hk0)[001] slip systems. In contrast, epidote grains accommodate less strain than Na-amphibole, via some combination of rigid rotation, brittle boudinage, and minor intracrystalline plasticity.</p><p>Observations of evenly-distributed strain, despite lithological heterogeneity, suggest low viscosity contrasts and comparable bulk strengths of quartz schists and blueschists. Our microstructural observations suggest that Na-amphibole was the weakest phase and accommodated the majority of strain in mafic blueschists. Dislocation activity, and not just rigid-body-rotation or diffusional processes, accommodated some component of strain and possibly transitioned with increasing temperature from glide- to climb-controlled. Although effective viscosities appear to be similar, subduction interface shear zones dominated by blueschists may exhibit a power-law rheology consistent with dislocation activity, in contrast to the common inference of Newtonian creep in metasediments. Complementary experimental work on CMS and CBU rocks will also be presented at this meeting (see Tokle et al. and Hufford et al.).</p>


1993 ◽  
Vol 57 (386) ◽  
pp. 55-66 ◽  
Author(s):  
D. Brown ◽  
K. R. McClay

AbstractThe Vangorda Pb-Zn-Ag orebody is a 7.1 M tonne, polydeformed stratiform massive sulphide deposit in the Anvil mining district, Yukon, Canada. Five sulphide lithofacies have been identified within the desposit with a typical mineralogy of pyrite, sphalerite, galena, and barite. Pyrrhotite-sphaleritemagnetite assembalges are locally developed. Etched polished sections of massive pyrite ores display relict primary depositional pyrite textures such as colloform growth zoning and spheroidal/framboidal features. A wide variety of brittle deformation, ductile deformation, and annealing textures have been identified. Brittle deformation textures include thin zones of intense cataclasis, grain indentation and axial cracking, and grain boundary sliding features. Ductile deformation textures include strong preferred grain shape orientations, dislocation textures, grain boundary migration, dynamic recrystallisation and pressure solution textures. Post deformational annealing has produced grain growth with lobate grain boundaries, 120° triple junctions and idioblastic pyrite porphyroblasts. The distribution of deformation textures within the Vangorda orebody suggests strong strain partitioning along fold limbs and fault/shear zones, it is postulated that focussed fluid flow in these zones had significant effects on the deformation of these pyritic ores.


2020 ◽  
Author(s):  
Clément Herviou ◽  
Anne Verlaguet ◽  
Philippe Agard ◽  
Hugues Raimbourg ◽  
Michele Locatelli ◽  
...  

<p>Important amounts of fluids are released in subduction zones by successive dehydration reactions occurring both in the previously hydrated oceanic crust (and mantle) and overlying sedimentary cover. The release and circulation of such fluids in rocks have major consequences on both their mechanical and chemical behavior. Indeed, the presence of a free fluid phase strongly modifies the rock rheology, fracturing properties, and could be implicated in both intermediate-depth earthquake and slow slip events nucleation. Moreover, the scale of mass transfer, associated chemical changes in infiltrated rocks and element recycling in subduction zones are controlled by both the rock permeability and the amount and composition of such fluids. Thus, there is a crucial need to identify the major fluid sources, amounts and pathways to better constrain their impact on subduction dynamics.</p><p>Metamorphic veins, as well as mineralized fractures and shear zones in exhumed fossil subduction zones are the best witnesses of fluid-rock interactions and fluid circulation pathways. However, their interpretation in terms of fluid sources, residence time, scale of circulation requires a good knowledge of the composition of potential fluid sources. In order to determine the composition of the fluid released by both oceanic crust and sediments at various depth along their subduction, we analyzed the composition of fluid inclusions contained in vein minerals formed at peak P-T conditions, in rock units buried at various depths in the Alpine subduction zone.</p><p>The Schistes Lustrés complex is a slice-stack representing the deep, underplated part of the former Alpine accretionary wedge. These Alpine Tethys rocks are mainly composed of oceanic calcschists with fewer mafic and ultramafic rocks, buried to various depths before exhumation. From West to East, the juxtaposed Schistes Lustrés units show increasing peak P-T conditions from blueschist (300-350°C - 1.2-1.3 GPa) to eclogite facies (580°C - 2.8 GPa). This study focuses on the Schistes Lustrés - Monviso transect, which shows an almost continuous increase in metamorphic grade.</p><p>In the Schistes Lustrés blueschist-facies sediments, fluid inclusions were analyzed in quartz from high-pressure veins, i.e. quartz that co-crystallized with prograde to peak metamorphic minerals such as lawsonite and Fe-Mg carpholite. In the metamorphosed mafic rocks, we analyzed fluid inclusions from the peak metamorphic assemblages, i.e. glaucophane +/- omphacite in blueschist facies rocks, omphacite in eclogite-facies slices. Raman spectroscopy data on these fluid inclusions suggest that fluids released during dehydration of calcschists in blueschist-facies conditions are aqueous fluids with low-salinity and small amounts of CO<sub>2</sub> and CH<sub>4</sub>. In contrast, eclogitic fluids released from metagabbros are highly saline brines with low N<sub>2 </sub>content. These results, which will be associated with LA-ICP-MS analysis of fluid inclusions in metasedimentary quartz veins, will contribute to better constrain the evolution of composition of the fluids liberated by dehydration reactions with depth and protolith composition along the subduction interface.</p>


2021 ◽  
Author(s):  
Francesco Giuntoli ◽  
Giulio Viola

<p>The geological record of deep seismic activity in subduction zones is generally limited due to common rock overprinting during exhumation and only a few regions allow studying well-preserved exhumed deep structures. The Northern Apennines (Italy) are one such area, granting access to continental units (Tuscan Metamorphic Units) that were subducted to high-pressure conditions, were affected by brittle-ductile deformation while accommodating deep tremor and slip and then exhumed back to surface, with only minor retrogression.</p><p>Our approach is based on detailed fieldwork, microstructural and petrological investigations. Field observations reveal a metamorphosed broken formation composed of boudinaged metaconglomerate levels enveloped by metapelite displaying a pervasive mylonitic foliation. Shear veins occur in both lithologies, but are more common and laterally continuous in the metapelite. They are mostly parallel to the foliation and composed of iso-oriented stretched quartz and Mg-carpholite (XMg>0.5) fibres, which are single-grains up to several centimetres long. These fibres define a stretching direction coherent with that observed in the metaconglomerate and metapelite, which is marked by K-white mica and quartz. Thermodynamic modeling constrains the formation of the high-pressure veins and the mylonitic foliation to ~ 1 GPa and 350°C, corresponding to c. 30-40 km depth in the subduction channel.</p><p>Shear veins developed in subducted (meta)sediments are a key indicator of episodic tremor and slip (e.g. <sup>1</sup>). We propose that these structures reflect the repeated alternation of localised brittle failure, with shear veins development, and more diffuse viscous deformation. These cycles were probably related to the fluctuation of pore pressure that repeatedly reached lithostatic values. Concluding, these structures can be considered the geological record of episodic tremors and slip occurring at >30 km of depth in the Apenninic subduction channel.</p><p>1. Fagereng, Å., Remitti, F. & Sibson, R. H. Incrementally developed slickenfibers — Geological record of repeating low stress-drop seismic events? Tectonophysics <strong>510</strong>, 381–386 (2011).</p><p>This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 839779.</p>


2020 ◽  
Author(s):  
Carolyn Tewksbury-Christle ◽  
Whitney Behr ◽  
Mark Helper

<p>The low velocity layer (LVL) in modern subduction zones is a 3-5 km thick region that parallels the top of the downgoing slab and is characterized by anomalously high V<sub>p</sub>/V<sub>s</sub> ratios (1.8-2.5) consistent with 2.5-4% fracture porosity at near-lithostatic pore fluid pressures. The LVL has been previously interpreted as partially hydrated, relatively undeformed oceanic crust at the top of the downgoing slab, but collocation of the LVL with episodic tremor and slow slip events (ETS) in modern subduction zones suggests that the LVL may alternatively represent the seismic signature of a subduction interface shear zone. </p><p>To test this hypothesis, we use field & structural observations, geochronology, and seismic velocity calculations to compare and contrast the bulk seismic properties of a fossil subduction interface shear zone (Condrey Mountain Schist, CMS, northern CA) to properties of modern LVLs. Specifically, we 1) determined thicknesses of underplated packages (interpreted to represent the maximum thickness of the actively deforming interface) using depositional age discontinuities and high resolution structural mapping, 2) averaged the bulk rock seismic velocities weighted by mapped lithologic proportions and corrected for pressure-temperature effects, and 3) used field evidence of modifying factors (e.g., microcracks, fluid-filled veins, mineral anisotropy) to further refine the possible range of seismic velocities and effects on V<sub>p</sub>/V<sub>s</sub> ratio.</p><p>The CMS greenschist- to blueschist-facies units were subducted to ~25-35 km (450°C, 0.8-1.0 GPa) with limited retrogression or exhumational overprint. These rocks were underplated episodically at depth in three packages individually up to 4.5 km thick from 155-135 Ma, based on detrital zircon data. Each package is dominantly composed of metasedimentary rocks with m- to km-scale metamafic and serpentinized ultramafic lenses. Strain localization to ~1 km thick ductile shear zones between underplating episodes is collocated with km-scale serpentinized ultramafic lenses at the base of each package. Deformation was distributed and ductile with rare macro- or micro-scale prograde brittle failure in the metasedimentary or metamafic units. In the serpentinized ultramafics, ductile shear zones wrap massive blocks with prograde brittle fracture. Maximum fracture porosity estimated from relict veins is ~10%. Average V<sub>p</sub>/V<sub>s</sub> for the CMS is ~1.6 (lithology alone) but up to 3.0 (accounting for maximum fracture porosity).</p><p>The fossil subduction interface shear zone preserved in the CMS is consistent in both thickness and seismic signature with the LVL in modern subduction zones. Estimated V<sub>p</sub>/V<sub>s</sub> is higher than the LVL but assumes that all fractures are simultaneously open. The total thickness of the CMS (10+ km) is greater than the LVL, however, so previously underplated material must lose its anomalous seismic signature during underplating (e.g., due to fluid loss and transport up the slab during or after underplating). Our results support the hypothesis that LVLs in modern subduction zones represent the seismic signature of the subduction interface shear zone.</p>


2019 ◽  
Author(s):  
Mathieu Soret ◽  
Philippe Agard ◽  
Benoît Ildefonse ◽  
Benoît Dubacq ◽  
Cécile Prigent ◽  
...  

Abstract. This study sheds light on the deformation mechanisms of subducted mafic rocks metamorphosed at amphibolite and granulite facies conditions, and on their importance for strain accommodation and localization at the top of the slab during subduction infancy. These rocks, namely metamorphic soles, are oceanic slivers stripped from the downgoing slab and plastered below the upper plate mantle wedge during the first million years of intra-oceanic subduction, when the subduction interface is still warm. Their formation and intense deformation (i.e. shear strain ≥ 5) attest to a systematic and transient coupling between the plates over a restricted time span of ~1 My and specific rheological conditions. Combining micro-structural analyses with mineral chemistry constrains grain-scale deformation mechanisms and the rheology of amphibole and amphibolites along the plate interface during early subduction dynamics, as well as the interplay between brittle and ductile deformation, water activity, mineral change, grain size reduction and phase mixing. Results indicate, in particular, that increasing pressure-temperature conditions and slab dehydration (from amphibolite to granulite facies) lead to the crystallization of mechanically strong phases (garnet, clinopyroxene and high-grade amphibole) and rock hardening. In contrast, during early exhumation and cooling (from ~850 down to ~700 °C – 0.7 GPa), the garnet-clinopyroxene-bearing amphibolite experiences pervasive retrogression (and fluid ingression) and significant strain weakening essentially accommodated by dissolution-precipitation and grain boundary sliding processes. Observations also indicate cyclic brittle deformation near peak conditions and throughout the early exhumation, which contributed to fluid channelization within the amphibolites, and possibly strain localization accompanying detachment from the slab. These mechanical transitions, coeval with detachment and early exhumation of the HT metamorphic soles, controlled mechanical coupling across the plate interface during subduction infancy, between the top of the slab and the peridotites above. Our findings may thus apply to other geodynamic environments where similar temperatures, lithologies, fluid circulation and mechanical coupling between mafic rocks and peridotites prevail, such as in mature warm subduction zones (e.g., Nankai, Cascapedia), in lower continental crust shear zones and oceanic detachments.


Solid Earth ◽  
2016 ◽  
Vol 7 (2) ◽  
pp. 579-598 ◽  
Author(s):  
Francisco José Fernández ◽  
Sergio Llana-Fúnez ◽  
Pablo Valverde-Vaquero ◽  
Alberto Marcos ◽  
Pedro Castiñeiras

Abstract. High-grade, highly deformed gneisses crop out continuously along the Masanteo peninsula and constitute the upper part of the lower crustal section in the Cabo Ortegal nappe (NW Spain). The rock sequence formed by migmatitic quartzo-feldspathic (qz-fsp) gneisses and mafic rocks records the early Ordovician (ca. 480–488 Ma) injection of felsic dioritic/granodioritic dykes at the base of the qz-fsp gneisses, and Devonian eclogitization (ca. 390.4 ± 1.2 Ma), prior to its exhumation. A SE-vergent ductile thrust constitutes the base of quartzo-feldspathic gneissic unit, incorporating mafic eclogite blocks within migmatitic gneisses. A NW-vergent detachment displaced metasedimentary qz-fsp gneisses over the migmatites. A difference in metamorphic pressure of ca. 0.5 GPa is estimated between both gneissic units. The tectono-metamorphic relationships of the basal ductile thrust and the normal detachment bounding the top of the migmatites indicate that both discrete mechanical contacts were active before the recumbent folding affecting the sequence of gneisses during their final emplacement. The progressive tectonic exhumation from eclogite to greenschist facies conditions occurred over ca. 10 Ma and involved bulk thinning of the high-grade rock sequence in the high pressure and high temperature (HP–HT) Cabo Ortegal nappe. The necessary strain was accommodated by the development of a widespread main foliation, dominated by flattening, that subsequently localized to a network of anastomosing shear bands that evolved to planar shear zones. Qz-fsp gneisses and neighbouring mafic granulites were exhumed at > 3 mm yr−1, and the exhumation path involved a cooling of  ∼  20 °C/100 MPa, These figures are comparable to currently active subduction zones, although exhumation P–T trajectory and ascent rates are at the hotter and slower end in comparison with currently active similar settings, suggesting an extremely ductile deformation environment during the exhumation of qz-fsp gneisses within a coherent Cabo Ortegal nappe.


2020 ◽  
Author(s):  
Emily Warren-Smith ◽  
Bill Fry ◽  
Laura Wallace ◽  
Enrique Chon ◽  
Stuart Henrys ◽  
...  

<p><span>The occurrence of slow slip events (SSEs) in subduction zones has been proposed to be linked to the presence of, and fluctuations in near-lithostatic fluid pressures (P</span><sub><span>f</span></sub><span>) within the megathrust shear zone and subducting oceanic crust. In particular, the 'fault-valve' model is commonly used to describe occasional, repeated breaching of a low-permeability interface shear zone barrier, which caps an overpressured hydrothermal fluid reservoir. In this model, a precursory increase in fluid pressure may therefore be anticipated to precede megathrust rupture. Resulting activation of fractures during slip opens permeable pathways for fluid migration and fluid pressure decreases once more, until the system becomes sealed and overpressure can re-accumulate. While the priming conditions for cyclical valving behaviour have been observed at subduction zones globally, and evidence for post-megathrust rupture drainage exists, physical observations of precursory fluid pressure increases, and subsequent decreases, particularly within the subducting slab where hydrothermal fluids are sourced, remain elusive. </span></p><p><span>Here we use earthquake focal mechanisms recorded on an ocean-bottom seismic network to identify changes in the stress tensor within subducting oceanic crust during four SSEs in New Zealand’s Northern Hikurangi subduction zone. We show that the stress, or shape ratio, which describes the relative magnitudes of the principal compressive stress axes, shows repeated decreases prior to, and rapid increases during the occurrence of geodetically documented SSEs. We propose that these changes represent precursory accumulation and subsequent release of fluid pressure within overpressured subducting oceanic crust via a ‘valving’ model for megathrust slip behaviour. Our observations indicate that the timing of slow slip events on subduction megathrusts may be controlled by cyclical accumulation of fluid pressure within subducting oceanic crust.</span></p><p><span>Our model is further supported by observations of seismicity preceding a large SSE in the northern Hikurangi Margin in 2019, captured by ocean-bottom seismometer</span><span>s</span><span> and </span><span>absolute </span><span>pressure </span><span>recorders.</span> <span>O</span><span>bservations of microseismicity </span><span>during this period </span><span>indicate that a stress state conducive to vertical fluid flow was present in the downgoing plate prior to SSE initiation, before subsequently returning to a</span><span> down-dip</span><span> extensional state following the SSE. We propose this precursory seismicity is indicative of fluid migration towards the interface shear zone from the lower plate fluid reservoir, which may have helped triggering slip on the megathrust. </span></p><p><span>We also present preliminary results of a moment tensor study to investigate spatial and temporal patterns in earthquake source properties in SSE regions along the Hikurangi Margin. In particular, earthquakes near Porangahau – a region susceptible to dynamic triggering of tremor and where </span><span>shallow </span><span>SSEs occur every 5 years or so – exhibit distinctly lower double couple components than elsewhere along the margin. We </span><span>attribute this to elevated fluid pressures within the crust here, which is consistent with recent observations of high seismic reflectivity from an autocorrelation study. Such high fluid pressure may control the broad range of seismic and aseismic phenomena observed at Porangahau. </span></p>


Sign in / Sign up

Export Citation Format

Share Document