Basaltic reservoirs in the Earth's mantle transition zone.

The formation and preservation of compositional heterogeneities inside the Earth affect mantle convection patterns globally and control the long-term evolution of geochemical reservoirs. However, the distribution, nature, and size of reservoirs in the Earth's mantle are poorly constrained. Here, we invert measurements of travel times and amplitudes of seismic waves interacting with mineralogical phase transitions at 400-700-km depth to obtain global probabilistic maps of temperature and bulk composition. We find large basalt-rich pools (up to 60% basalt fraction) surrounding the Pacific Ocean, which we relate to the segregation of oceanic crust from slabs that have been subducted since the Mesozoic. Segregation of oceanic crust from initially cold and stiff slabs may be facilitated by the presence of a weak hydrated layer in the slab or by weakening upon mineralogical transition due to grain-size reduction.

Limited underthrusting of India below Tibet: 3He/4He analysis of thermal springs locates the mantle suture in continental collision.

During continent­continent collision, does the downgoing continental plate underplate far inboard of the collisional boundary or does it subduct steeply into the mantle, and how is this geometry manifested in the mantle flow field? We test conflicting models for these questions for Earth's archetypal continental collision forming the Himalaya and Tibetan Plateau. Air-corrected helium isotope data (3He/4He) from 225 geothermal springs (196 from our group, 29 from the literature) delineate a boundary separating a Himalayan domain of only crustal helium from a Tibetan domain with significant mantle helium. This 1,000-km-long boundary is located close to the Yarlung-Zangbo Suture (YZS) in southern Tibet from 80 to 92°E and is interpreted to overlie the "mantle suture" where cold underplated Indian lithosphere is juxtaposed at >80 km depth against a sub-Tibetan incipiently molten asthenospheric mantle wedge. In southeastern Tibet, the mantle suture lies 100 km south of the YZS, implying delamination of the mantle lithosphere from the Indian crust. This helium-isotopic boundary helps resolve multiple, mutually conflicting seismological interpretations. Our synthesis of the combined data locates the northern limit of Indian underplating beneath Tibet, where the Indian plate bends to steeper dips or breaks off beneath a (likely thin) asthenospheric wedge below Tibetan crust, thereby defining limited underthrusting for the Tibetan continental collision.

Growth and stabilization of induced seismicity rates during long-term, low-pressure fluid injection.

We examine the temporal evolution of sequences of induced seismicity caused by long-term fluid injection using a compilation of over 20 case studies where moderate magnitude (M > 3.0) induced events have been recorded. We compare rates of seismicity with injection rates via the seismogenic index and seismic efficiency parameters, computing both cumulative and time-windowed values. We find that cumulative values tend to accelerate steeply as each seismicity sequence initiates-most cases reach a value that is within 0.5 units of their maximum value within 1-3 years. Time-windowed values tend to increase to maximum values within 25%-35% of the overall sequence, before decreasing as levels of seismicity stabilize. We interpret these observations with respect to the pore pressure changes that will be generated in highly porous, high permeability reservoirs. In such situations, the rate of pore pressure change is highest during the early phases of injection and decreases with time. If induced seismicity scales with the rate of deformation, which in turn is controlled by the rate of pore pressure change, then it is to be expected that induced seismicity is highest during the early phases of injection, and then decreases with time. This article is part of the theme issue 'Induced seismicity in coupled subsurface systems'.

Are maximum magnitudes of induced earthquakes controlled by pressure diffusion?

There is an ongoing discussion about how to forecast the maximum magnitudes of induced earthquakes based on operational parameters, subsurface conditions and physical process understanding. Although the occurrence of damage caused by induced earthquakes is rare, some cases have caused significant economic loss, injuries and even loss of life. We analysed a global compilation of earthquakes induced by hydraulic fracturing, geothermal reservoir stimulation, water disposal, gas storage and reservoir impoundment. Our analysis showed that maximum magnitudes scale with the characteristic length of pressure diffusion in the brittle Earth's crust. We observed an increase in the nucleation potential of larger-magnitude earthquakes with time and explained it by diffusion-controlled growth of the pressure-perturbed part of faults. Numerical and analytical fault size modelling supported our findings. Finally, we derived magnitude scaling laws to manage induced seismic hazard of upcoming energy projects prior to operation. This article is part of the theme issue 'Induced seismicity in coupled subsurface systems'.