Volume Of Crust Research Articles

The September 6, 2021 MW 5.4 Tofalaria earthquake at a weakly active segment of the Main Sayan fault (Eastern Siberia)

We study in detail the Mw 5.4 September 6, 2021 Tofalaria earthquake occurred in a mountain area of the Eastern Sayan which is characterized by a low level of seismic activity. An interest in the seismic event is caused, on the one hand, by poor knowledge about stress-strain field of the crust in the considered region, and, on the other hand, by its relation to the NW segment of the ancient Main Sayan fault – a structural boundary between the Sayan-Baikal fold belt and the tectonically stable Siberian platform. Seismic moment tensors and hypocentral depths of the mainshock and its largest aftershock (Mw 4.6) are inverted from intermediate-period surface wave amplitude spectra calculated at the stations located at teleseismic distances. Integral source parameters of the mainshock, characterizing its spatio-temporal development, are also estimated and the fault plane is determined. Epicenters of 31 aftershocks with M ≥ 1.8, occurred up to the end of 2021, are constrained from body waves recorded at regional seismic stations. The obtained results show that the Tofalaria earthquake occurred under the influence of the SW-NE compression, which is observed in Western Mongolia. Focal mechanism of the largest aftershock and the NE elongation of the aftershock epicentral field (22 km) indicate stress redistribution after the mainshock in a local crustal volume, bordered by small-scale faults.

Frequency-dependent shear wave attenuation across the Central Anatolia region, Türkiye

Abstract. The Central Anatolian Plateau with its volcanic provinces represents a broad transition zone between the compressional deformation in the east and the extensional regime in the west. The Central Anatolian Fault Zone separates the Kırşehir Block in the north and the Anatolide–Tauride Block in the south within the plateau. A proper understanding of physical properties such as seismic attenuation in the crustal volume of this region can provide hints toward the possible source for the geodynamic events in the past and present that likely lead to the observed deformation. In order to model intrinsic and scattering attenuation separately, we perform a nonempirical coda-wave modeling approach in which a fitting process between observed and synthetic coda-wave envelopes is performed for each earthquake in multiple frequency bands. Here, the acoustic radiative transfer theory, assuming multiple isotropic scattering, was utilized for the forward modeling of the synthetic coda-wave envelopes of local earthquakes. Our findings generally highlight the prominent nature of intrinsic attenuation over scattering attenuation, implying the presence of thick volcanic rocks with relatively high attenuation values beneath Central Anatolia. Overall, the spatial distribution of the attenuation at varying frequencies marks the Kırşehir Massif distinctively with its considerable high-attenuating character. Our findings, combined with early seismological and geo-electrical models, suggest a possible partial melt beneath most of the Central Anatolian Volcanic Province, and the resultant zones of elevated fluid-rich content exhibit high and dominant intrinsic attenuation. To the southeast, a gradual decrease in the observed attenuation coincides with the Central Taurus Mountains where high altitude is considered to be evolved following the slab break-off and resulting mantle upwelling.

A slowly cooled deep crust on asteroid 4 Vesta and the recent impact history of rubble pile vestoids recorded by diogenites

In this study, we investigate the 40Ar/39Ar systematics of nineteen diogenites thought to come from deep crustal levels of asteroid 4 Vesta. We applied both Electron Backscattered Diffraction (EBSD) and 40Ar/39Ar and methods to the unbrecciated diogenite LAP 031381. We obtained three plateau ages resulting in a combined weighted mean age of 4441 ± 15 Ma (P = 0.16). The EBSD analyses suggest that LAP 031381 displays minimal evidence of shock and, when combined with petrography observations, diffusion modelling and 40Ar/39Ar data, these results suggest that the crustal volume that initially contained this diogenite, reached a temperature of ca. 630 °C at ∼ 4.44 Ga. This corresponds to a linear cooling rate of ∼ 5 °C / Ma for a crystallization age of 4550 Ma. Independent thermal models suggest that these conditions were present at a depth of 60 to 65 km at 4.44 Ga.The other eighteen diogenites yielded 40Ar/39Ar results that indicate that they have been variously shocked by impact events and seven of them yielded plateau ages ranging from 2413 ± 189 Ma to 84 ± 162 Ma. We combined these results with 40Ar/39Ar ages from eucrites and howardites and propose that the HED (Howardite, Eucrite, Diogenite) meteorites recorded impact events at the surface of Vesta until ∼ 3.4 Ga when they were then ejected during a large collision. The eucrites, diogenites and howardites were then recombined into small rubble pile asteroids which probably make up a large part of the Vestoid family. After ejection, the K/Ar system in plagioclase crystals ceased in most cases to be fully reset by impact events as the temperature spikes reached during small impacts lack enough energy to trigger significant 40Ar* diffusion. On the other hand, ultra-transient and high-temperature – sensitive pyroxene crystals kept a more systematic record of small impacts until recent time. 38Arc cosmochron cosmogenic exposure ages on diogenites mostly range from 51 ± 7 Ma to 0 ± 1 Ma and when combined with other HED cosmochron ages, suggest that almost all the HED meteorites were continuously ejected from secondary rubble pile asteroids mostly between 50 Ma and present.

Investigating the Role of Fluids in the Source Parameters of the 2013–2014 Mw 5 Matese Seismic Sequence, Southern Italy

Abstract We investigate the variability of Brune stress drop (Δσ), apparent stress (τa), and Savage–Wood radiation efficiency (ηsw=τa/Δσ), in the 2013–2014 Mw 5.0 earthquake sequence that struck the Matese area in the southern Apennines range of Italy. The sequence is clustered in a relatively small crustal volume in the 13–22 km depth range, which is greater than that of background seismicity and normal-faulting sequences that occurred under the range axis, usually located in the first 15 km of the crust. We find high Savage–Wood radiation efficiency values for most of the analyzed earthquakes located in a narrow crustal volume, with values ranging from well above the self-similarity value to very high values as high as 0.55. In addition, a large variability in radiation efficiency (up to 90%) is observed for two similar magnitude events at different depths. Previous studies reported seismic evidence of fluid involvement in the nucleation process of the Matese earthquakes. By integrating our results with crustal geophysical data published recently, we propose that most of the earthquakes characterized by high values of ηsw are nucleated within high pore pressure zones located in the crystalline midcrust of Adria. We reckon that high pore pressure fluids of deep origin played a role in the rupture process and were responsible for the mixed shear-tensile sources inferred from the analysis of the S-wave/P-wave spectral amplitude ratio for most of 2013–2014 earthquakes.

Crustal Imaging of Portugal Mainland Using Magnetotelluric Data

AbstractThe comprehensive mapping of resistivities in depth and its interpretation is of great importance to unravel the electrical properties of deep‐seated rocks, providing significant insights into the structure of the crust. The first 3D resistivity model for Portugal mainland is here reported using data from 31 broadband MT soundings spaced 50 × 50 km apart. The model shows large crustal volumes with contrasting resistivity values. The central and northern regions are characterized by a high resistivity crustal domain spreading in‐depth (103–105 Ω.m) that correlates well with voluminous granitoid bodies and high‐grade metamorphic rocks. To the west and the south, roughly corresponding to the Portuguese Western Shore and the South Portuguese Zone (SPZ) respectively, the crustal resistivity tends to decrease unevenly to values between 1 and 103 Ω.m. Interconnected graphite in structural discontinuities, possibly forming a regional mid‐crustal décollement at ≈13–15 km depth, are highlighted by multiple low resistivity crustal domains (1–100 Ω.m). Similarly, the model reveals a very large E‐W low resistivity crustal domain in SPZ possibly representing a deep‐seated major décollement. The present study aims to contribute to a better understanding of the physical properties characterizing the concealed crust in Portugal mainland, providing also indirect information on its composition and structure.

Pressure increase at the magma-hydrothermal interface at Krafla caldera, North-Iceland, 2018–2020: Magmatic processes or hydrothermal changes?

Following volcano deflation since 1989, the deformation pattern in Krafla caldera, North Iceland, changed in 2018. Geodetic measurements reveal a difference in surface velocity fields for 2015–2018 and 2018–2020 periods, reflecting a change in the deformation pattern. The difference velocity field broadly fits deformation caused by a spherical pressure source within a uniform elastic half-space, with a volume change of 2.6–3.8 × 105 m3/yr and centre depth of 2.1–2.5 km, which is close to the brittle-ductile boundary in the area, at a depth of 1.8 to 2.2 km. Potential processes causing the deformation change are evaluated: magmatic processes such as magma inflow or accumulation of volcanic gas, changes in the geothermal area because of change in geothermal production, or a combination of these. In particular, we evaluate if the change in deformation may relate to about 0.1 MPa/yr pressure increase in the geothermal system as measurements in monitoring well KG-10 indicate, eventually due to changes in the geothermal exploitation strategy at the Krafla power plant. Modelling shows that inferred volume change may be due to a spherical source with 1.4 km radius with 0.1 MPa/yr pressure change if the surrounding crust has a Young's modulus E of about 7 GPa. However, the average regional Young's modulus for the upper crust in Iceland has been estimated to be 30 GPa. We use the Finite Element Method (FEM) to assess the influence on the displacement due to the presence of a local crustal volume 5 × 5 × 4 km (horizontal dimensions × depth), which envelops the source (within the Krafla caldera), with E = 7 GPa in the central area and 30 GPa in the far field, in a three-dimensional model. Such a model can reproduce significant features of the observed deformation. There are no changes in seismicity in 2018. In late 2019, the earthquake rate increases following a slight decrease over few months. The seismic moment release is relatively steady until the end of 2019, when the rate increases and is relatively constant in 2020. Gravity measurements in late 2019, when compared to limited measurements in 2018, are inconclusive regarding the nature of the deformation, but useful for further monitoring. No noticeable changes have been observed in the chemical composition of fumarole discharge in the Krafla field that relate to new intrusions.

Contemporaneous Thick- and Thin-Skinned Seismotectonics in the External Zagros: The Case of the 2021 Fin Doublet, Iran

In this work, we propose a geodetic model for the seismic sequence, with doublet earthquakes, that occurred in Bandar Abbas, Iran, in November 2021. A dataset of Sentinel-1 images, processed using the InSAR (Interferometric Synthetic Aperture Radar) technique, was employed to identify the surface deformation caused by the major events of the sequence and to constrain their geometry and kinematics using seismological constraints. A Coulomb stress transfer analysis was also applied to investigate the sequence’s structural evolution in space and time. A linear inversion of the InSAR data provided a non-uniform distribution of slip over the fault planes. We also performed an accurate relocation of foreshocks and aftershocks recorded by locally established seismographs, thereby allowing us to determine the compressional tectonic stress regime affecting the crustal volume. Despite the very short time span of the sequence, our results clearly suggest that distinct blind structures that were previously unknown or only suspected were the causative faults. The first Mw 6.0 earthquake occurred on an NNE-dipping, intermediate-angle, reverse-oblique plane, while the Mw 6.4 earthquake occurred on almost horizontal or very low-angle (SSE-dipping) reverse segments with top-to-the-south kinematics. The former, which cut through and displaced the Pan-African pre-Palaeozoic basement, indicates a thick-skinned tectonic style, while the latter rupture(s), which occurred within the Palaeozoic–Cenozoic sedimentary succession and likely exploited the stratigraphic mechanical discontinuities, clearly depicts a thin-skinned style.

Evolution of the preserved European continental crust, constrained by U-Pb, O and Hf isotopic analyses of river detrital zircons

Integrated U-Pb, O and Lu-Hf isotopic analyses of detrital zircons, collected from fourteen European river sands and one beach sand from northern Germany, were carried out to study the evolution of the preserved European continental crust. The drainage area of these sand samples covers large areas in eastern and western Europe, which evolved independently until their amalgamation during the Paleozoic. U-Pb ages of detrital zircons from eastern and western Europe generally match the ages of the major magmatic events in the two tectonic regions. The combined histogram of 2403 zircon U-Pb ages, obtained for this study, shows five major periods of zircon crystallisation: 2.80–2.65, 1.95–1.35, 1.25–0.85, 0.65–0.40, and 0.40–0.25 Ga, which correspond broadly with the timing of supercontinent assembly. The correlation could be due to enhanced granitoid magmatism and zircon growth during the supercontinent amalgamation, or the high preservation potential of igneous rocks generated during continental collisions. The O and Hf isotopic data of the analysed European river zircons show wide spreads in δ18O and εHf(t) values at ca. 1.9 Ga and 650–250 Ma, contrasting with the limited range for the Mesoproterozoic zircons. This indicates that extensive reworking of pre-existing crustal rocks occurred during the late Paleoproterozoic and Paleozoic in the evolution of the European continental crust, whereas the formation of new crust played an important role during the Mesoproterozoic.Arc mantle Hf model ages for the detrital zircons were calculated to estimate when European juvenile crust was extracted from the mantle. Model age uncertainties were estimated using the bootstrap method, considering three sources of uncertainty: (i) uncertainty in the Hf isotopic composition of the mantle reservoir, (ii) uncertainty in the 176Lu/177Hf for the crustal source region, and (iii) analytical uncertainties. Two types of model ages are calculated based on constant crustal 176Lu/177Hf (TAMC), and on variable176Lu/177Hf (TAMV), the latter being constrained by zircon δ18O values. Histograms of the calculated TAMC and TAMV model ages both show major juvenile extraction events between the late Paleo- and Mesoproterozoic. Although the principal periods in the European crustal growth contrast with those identified in the African, Russian and Mississippi river basins, the overall growth curves for the four continents show comparable results: (i) the formation of little or no continental crust before 4.0 Ga, and (ii) the formation of 50 % of the present crustal volume by ca. 2.0–1.6 Ga.

Shallow magmatic intrusion evolution below La Palma before and during the 2021 eruption

La Palma, Canary Islands, underwent volcanic unrest which culminated in its largest historical eruption. We study this unrest along 2021 using Interferometric Synthetic Aperture Radar (InSAR) and a new improved interpretation methodology, comparing achieved results with the crustal structure. We reproduce the final phase of La Palma volcanic unrest, highligthing a shallow magma accumulation which begins about 3.5 months before the eruption in a crustal volume charactherized by low density and fractured rocks. Our modeling, together with our improved pictures of the crustal structure, allows us to explain the location and characteristics of the eruption and to detect failed eruption paths. These can be used to explain post-eruptive phenomena and hazards to the local population, such as detected gases anomalies in La Bombilla and Puerto Naos. Our results have implications for understanding volcanic activity in the Canaries and volcano monitoring elsewhere, helping to support decision-making and providing significant insights into urban and infrastructure planning in volcanic areas.

The 2021 Greece Central Crete ML 5.8 Earthquake: An Example of Coalescent Fault Segments Reconstructed from InSAR and GNSS Data

The ML 5.8 earthquake that hit the island of Crete on 27 September 2021 is analysed with InSAR (Interferometry from Synthetic Aperture Radar) and GNSS (Global Navigation Satellite System) data. The purpose of this work is to create a model with sufficient detail for the geophysical processes that take place in several kilometres below the earth’s surface and improve our ability to observe active tectonic processes using geodetic and seismic data. InSAR coseismic displacements maps show negative values along the LOS of ~18 cm for the ascending orbit and ~20 cm for the descending one. Similarly, the GNSS data of three permanent stations were used in PPK (Post Processing Kinematic) mode to (i) estimate the coseismic shifts, highlighting the same range of values as the InSAR, (ii) model the deformation of the ground associated with the main shock, and (iii) validate InSAR results by combining GNSS and InSAR data. This allowed us to constrain the geometric characteristics of the seismogenic fault and the slip distribution on it. Our model, which stands on a joint inversion of the InSAR and GNSS data, highlights a major rupture surface striking 214°, dipping 50° NW and extending at depth from 2.5 km down to 12 km. The kinematics is almost dip-slip normal (rake −106°), while a maximum slip of ~1.0 m occurred at a depth of ca. 6 km. The crucial though indirect role of inherited tectonic structures affecting the seismogenic crustal volume is also discussed suggesting their influence on the surrounding stress field and their capacity to dynamically merge distinct fault segments.

Detection of Spatial and Temporal Stress Changes During the 2016 Central Italy Seismic Sequence by Monitoring the Evolution of the Energy Index

AbstractWe consider approximately 23,000 microearthquakes that occurred between 2005 and 2016 in central Italy to investigate the crustal strength before and after the three largest earthquakes of the 2016 seismic sequence (i.e., the Mw 6.2, 24 August 2016 Amatrice, the Mw 6.1, 26 October 2016 Visso, and the Mw 6.5, 30 October 2016 Norcia earthquakes). We monitor the spatiotemporal deviations of observed radiated energy, ES, with respect to theoretical values, ESt, derived from a scaling model between ES and M0 calibrated for background seismicity in central Italy. These deviations, defined here as Energy Index (EI), allow us to identify in the years following the Mw 6.1, 2009 L’Aquila earthquake a progressive evolution of the dynamic properties of microearthquakes and the existence of high EI patches close to the Amatrice earthquake hypocenter. We show the existence of a crustal volume with high EI even before the Mw 6.5 Norcia earthquake. Our results agree with the previously suggested hypothesis that the Norcia earthquake nucleated at the boundary of a large patch, highly stressed by the two previous mainshocks of the sequence. We highlight the mainshocks interaction both in terms of EI and of the mean loading shear stress associated to microearthquakes occurring within the crustal volumes comprising the mainshock hypocenters. Our study shows that the dynamic characteristics of microearthquakes can be exploited as beacons of stress change in the crust and thus be exploited to monitor the seismic hazard of a region and help to intercept the preparation phase of large earthquakes.

Reworking subducted sediments in arc magmas and the isotopic diversity of the continental crust: The case of the Ordovician Famatinian crustal section, Argentina

Since the onset of plate tectonics, continents have evolved through a balance between crustal growth, reworking, and recycling at convergent plate margins. The term “reworking” involves the re-insertion of crustal material into pre-existing crustal volumes, while crustal growth and recycling respectively represent gains from and losses to the mantle. Reworking that occurs in the mantle wedge (“source” contamination from slab material) or within the upper plate (“path” contamination), will have contrasting effects on crustal evolution. However, due to limited access to deep crustal and mantle rocks, quantifying source vs. path contamination remains challenging. Based on the 4-dimensional record of the fossil (Ordovician) Famatinian continental arc (Argentina), we demonstrate that source contamination plays a dominant role in imprinting mafic to granitic rocks with crustal oxygen-hafnium (O-Hf) isotopic compositions. We argue that source contamination at convergent plate margins significantly increased the diversity of O-Hf isotopic signatures of continents over geologic time. Our interpretation implies that crustal evolution models attributing this isotopic diversity dominantly to intra-crustal reworking may be over-simplistic and may underestimate continental growth in the last 2.5 billion years.

Multi-scale flow structure of a strike-slip tectonic setting: A self-similar model for the Liquiñe-Ofqui Fault System and the Andean Transverse Faults, Southern Andes (39–40°S)

The flow structure of a brittle crustal volume is defined by the multi-scale geometric and hydraulic properties of its fracture meshes. The length density distribution n(L,l) and the transmissivity distribution K(L,l) control the hydrologic scaling, where l is fracture length and L is the system size. The flow structure might display at most three key hydrologic scales: the connection scale, above which flow is focused in few critical paths; the channeling scale, above which flow is distributed in several paths; and the homogenization scale, above which permeability approaches a constant value. According to these scales, the hydrological structure could be distributed or clustered, thus having a clear impact in geothermal exploration campaigns and reservoir modeling. In this work, we determine the multi-scale flow structure for the Liquiñe-Ofqui Fault System (LOFS) and the Andean Transverse Faults (ATF) in the Southern Andes, by establishing the hydrologic scaling they follow. Using fractal statistics, we integrated geological data at the regional, meso‑ and micro-scale, including image analysis from X-ray microtomography. Our results suggest a self-similar, dense network with n(L,l)∼l−a and a = 2.6–2.9, from the regional scale where the LOFS and ATF interact to the meso‑ and micro-scale within highly fractured areas of the LOFS. Scaling models are constrained by the length distribution, and other power-law functions reflecting the geometric arrangement of fractures, as well as the spatial distribution of superficial geothermal occurrences. Thus, we expect the hydrologic scaling to depend on the transmissivity distribution. Lognormal transmissivity distribution yields a permeability increase with scale, from the connection to the homogenization scales; whereas power-law transmissivity distribution yields a permeability increase from the connection scale without a limiting value. Approximations of the connection scale are around 10−3–100 m; the channeling scale, around 100–104 m; and if the homogenization scale exists, it should be equal or greater than 103–104 m. Finally, the results presented here could to define the internal architecture of fracture meshes in fault-controlled fluid flow, and be used to select an appropriate hydrologic model according to the analyzed scale. Therefore, these findings must be taken into consideration in future geothermal prospecting, modeling and exploitation.

Inflation-Deflation Episodes in the Hengill and Hrómundartindur Volcanic Complexes, SW Iceland

Non-eruptive uplift and subsidence episodes remain a challenge for monitoring and hazard assessments in active volcanic systems worldwide. Sources of such deformation may relate to processes such as magma inflow and outflow, motion and phase changes of hydrothermal fluids or magma volatiles, heat transfer from magmatic bodies and heat-mining from geothermal extraction. The Hengill area, in southwest Iceland, hosts two active volcanic systems, Hengill and Hrómundartindur, and two high-temperature geothermal power plants, Hellisheiði and Nesjavellir. Using a combination of geodetic data sets (GNSS and InSAR; Global Navigation Satellite Systems and Interferometry Synthetic Aperture Radar, respectively) and a non-linear inversion scheme to estimate the optimal analytical model parameters, we investigate the ground deformation between 2017–2018. Due to other ongoing deformation processes in the area, such as plate motion, subsidence in the two geothermal production fields, and deep-seated source of contraction since 2006, we estimate 2017–2018 difference velocities by subtracting background deformation, determined from data spanning 2015–2017 (InSAR) or 2009–2017 (GNSS). This method highlights changes in ground deformation observed in 2017–2018 compared to prior years: uplift signal of ∼10 km diameter located in the eastern part of the Hengill area, and geothermal production-related temporal changes in deformation near Húsmúli, in the western part of the Hengill area. We find an inflation source located between the Hengill and Hrómundartindur volcanic complexes, lasting for ∼5 months, with a maximum uplift of ∼12 mm. Our model inversions give a source at depth of ∼6–7 km, located approximately in the same crustal volume as an inferred contracting source in 2006–2017, within the local brittle-ductile transition zone. No significant changes were observed in local seismicity, borehole temperatures and pressures during the uplift episode. These transient inflation and deflation sources are located ∼3 km NW from a source of non-eruptive uplift in the area (1993–1999). We consider possible magmatic and hydrothermal processes as the causes for these inflation-deflation episodes and conclude that further geophysical and geological studies are needed to better understand such episodes.

Dynamic earthquake triggering response tracks evolving unrest at Sierra Negra volcano, Galápagos Islands.

The propensity for dynamic earthquake triggering is thought to depend on the local stress state and amplitude of the stress perturbation. However, the nature of this dependency has not been confirmed within a single crustal volume. Here, we show that at Sierra Negra volcano, Galápagos Islands, the intensity of dynamically triggered earthquakes increased as inflation of a magma reservoir elevated the stress state. The perturbation of short-term seismicity within teleseismic surface waves also increased with peak dynamic strain. Following rapid coeruptive subsidence and reduction in stress and background seismicity rates, equivalent dynamic strains no longer triggered detectable seismicity. These findings offer direct constraints on the primary controls on dynamic triggering and suggest that the response to dynamic stresses may help constrain the evolution of volcanic unrest.

Coseismic vertical ground deformations vs. intensity measures: Examples from the Apennines

DInSAR data provide a powerful tool to recognize the Earth's surface where permanent deformations are concentrated and undergo the strongest ground motion during an earthquake, i.e., defining the epicentral area. We analyzed three recent seismic events in the Apennines belt in Italy related to extensional and contractional earthquakes and documented that the largest macroseismic intensity is recorded where the ground underwent the largest vertical component of motion and the stronger vertical component of the peak ground acceleration. Besides site amplification effects that may occur everywhere even outside the epicentral area, we infer that the vertical oscillations in the so-called near-field allow the horizontal shaking to be more effective, hence producing larger damage above the active crustal volume. The active volume is the one moving vertically and is contemporaneously crossed by seismic waves nucleated by the shear on the fault. The surrounding passive volume in the far field is only crossed by seismic waves. The near field areas are elliptical and cover areas of 300–600 km2 for a range of Mw 6–6.9, and they should be the areas where the highest seismic hazard must be expected. Therefore, their area is too large to be neglected for seismic hazard assessment.

The Campotosto linkage fault zone between the 2009 and 2016 seismic sequences of central Italy: Implications for seismic hazard analysis

In the last decade central Italy was struck by devastating seismic sequences resulting in hundreds of casualties (i.e., 2009-L′Aquila moment magnitude [Mw] = 6.3, and 2016-Amatrice-Visso-Norcia Mw max = 6.5). These seismic events were caused by two NW-SE−striking, SW-dipping, seismogenic normal faults that were modeled based on the available focal mechanisms and the seismic moment computed during the relative mainshocks. The seismogenic faults responsible for the 2009-L′Aquila Mw = 6.3 (Paganica Fault—PF) and 2016-Amatrice-Visso-Norcia Mw max = 6.5 (Monte Vettore Fault—MVF) are right-stepping with a negative overlap (i.e., underlap) located at the surface in the Campotosto area. This latter was affected by seismic swarms with magnitude ranging from 5.0 to 5.5 during the 2009 seismic sequence and then in 2017 (i.e., a few months later than the mainshocks related with the 2016 seismic sequence). In this paper, the seismogenic faults related to the main seismic events that occurred in the Campotosto Seismic Zone (CSZ) were modeled and interpreted as a linkage fault zone between the PF and MVF interacting seismogenic faults. Based on the underlap dimension, the seismogenic potential of the CSZ is in the order of Mw = 6.0, even in the case that all the faults belonging to the zone were activated simultaneously. This has important implications for seismic hazard assessment in an area dominated by the occurrence of a major NW-SE−striking extensional structure, i.e., the Monte Gorzano Fault (MGF). Mainly due to its geomorphologic expression, this fault has been considered as an active and silent structure (therefore representing a seismic gap) able to generate an earthquake of Mw max = 6.5−7.0. However, the geological evidence provided with this study suggests that the MGF is of early (i.e., pre- to syn-thrusting) origin. Therefore, the evaluation of the seismic hazard in the Campotosto area should not be based on the geometrical characteristics of the outcropping MGF. This also generates substantial issues with earthquake geological studies carried out prior to the recent seismic events in central Italy. More in general, the 4-D high-resolution image of a crustal volume hosting an active linkage zone between two large seismogenic structures provides new insights into the behavior of interacting faults in the incipient stages of connection.

The MS 6.9, 1980 Irpinia Earthquake from the Basement to the Surface: A Review of Tectonic Geomorphology and Geophysical Constraints, and New Data on Postseismic Deformation

The MS 6.9, 1980 Irpinia earthquake occurred in the southern Apennines, a fold and thrust belt that has been undergoing post-orogenic extension since ca. 400 kyr. The strongly anisotropic structure of fold and thrust belts like the Apennines, including late-orogenic low-angle normal faults and inherited Mesozoic extensional features besides gently dipping thrusts, result in a complex, overall layered architecture of the orogenic edifice. Effective decoupling between deep and shallow structural levels of this mountain belt is related to the strong rheological contrast produced by a fluid-saturated, shale-dominated mélange zone interposed between buried autochthonous carbonates—continuous with those exposed in the foreland to the east—and the allochthonous units. The presence of fluid reservoirs below the mélange zone is shown by a high VP/VS ratio—which is a proxy for densely fractured fluid-saturated crustal volumes—recorded by seismic tomography within the buried autochthonous carbonates and the top part of the underlying basement. These crustal volumes, in which background seismicity is remarkably concentrated, are fed by fluids migrating along the major active faults. High pore fluid pressures, decreasing the yield stress, are recorded by low stress-drop values associated with the earthquakes. On the other hand, the mountain belt is characterized by substantial gas flow to the surface, recorded as both distributed soil gas emissions and vigorous gas vents. The accumulation of CO2-brine within a reservoir located at hypocentral depths beneath the Irpinia region is not only interpreted to control a multiyear cyclic behavior of microseismicity, but could also play a role in ground motions detected by space-based geodetic measurements in the postseismic period. The analysis carried out in this study of persistent scatterer interferometry synthetic aperture radar (PS-InSAR) data, covering a timespan ranging from 12 to 30 years after the 1980 mainshock, points out that ground deformation has affected the Irpinia earthquake epicentral area in the last decades. These ground motions could be a result of postseismic afterslip, which is well known to occur over years or even decades after a large mainshock such as the 23 November 1980, MS 6.9 earthquake due to cycles of CO2-brine accumulation at depth and its subsequent release by Mw ≥ 3.5 earthquakes, or most likely by a combination of both. Postseismic afterslip controls geomorphology, topography, and surface deformation in seismically active areas such as that of the present study, characterized by ~M 7 earthquakes. Yet, this process has been largely overlooked in the case of the 1980 Irpinia earthquake, and one of the main aims of this study is to fill such the substantial gap of knowledge for the epicentral area of some of the most destructive earthquakes that have ever occurred in Italy.

Microcontinents and Continental Fragments Associated With Subduction Systems

AbstractMicrocontinents and continental fragments are small pieces of continental crust that are surrounded by oceanic lithosphere. Although classically associated with passive margin formation, here we present several preserved microcontinents and continental fragments associated with subduction systems. They are located in the Coral Sea, South China Sea, central Mediterranean and Scotia Sea regions, and a “proto‐microcontinent,” in the Gulf of California. Reviewing the tectonic history of each region and interpreting a variety of geophysical data allows us to identify parameters controlling the formation of microcontinents and continental fragments in subduction settings. All these tectonic blocks experienced long, complex tectonic histories with an important role for developing inherited structures. They tend to form in back‐arc locations and separate from their parent continent by oblique or rotational kinematics. The separated continental pieces and associated marginal basins are generally small and their formation is quick (<50 Myr). Microcontinents and continental fragments formed close to large continental masses tend to form faster than those created in systems bordered by large oceanic plates. A common triggering mechanism for their formation is difficult to identify, but seems to be linked with rapid changes of complex subduction dynamics. The young ages of all contemporary pieces found in situ suggest that microcontinents and continental fragments in these settings are short lived. Although presently the amount of in‐situ subduction‐related microcontinents is meager (an area of 0.56% and 0.28% of global, non‐cratonic, continental crustal area and crustal volume, respectively), through time microcontinents contributed to terrane amalgamation and larger continent formation.

Seismic Imaging of Convective Geothermal Flow Systems to Increase Well Productivity

A core feature of convective geothermal resource production is wellbore energy flow Q ~ ρC x T x V. E.g., for wellbore fluid of volume heat capacity ρC ~ 4.3MJ/m3∙oC, temperature T ~ 230oC, and volumetric flow V ~ 50L/s, wellbore heat energy production is Q~ 50MWth ~ 5MWe. Wellbore fluid flow V =2πr0φv0ℓ for open wellbore length ℓ is given in turn by the spatially variable product crustal porosity times crustal fluid velocity v ≡ φv0 at the wellbore radius r0. For a geothermal wellbore to be productive (nominal Q ~ 5MWe), locally variable bulk inflow rates v = φv0 across crustal volumes of dimension ℓ must be adequate to sustain high wellbore flows (nominal V ~ 50L/s). Wide-ranging crustal well productivity statistics show that few crustal wells flow at these rates. This is not surprising as local bulk flow ~ 10-2 m/s needed for production wellbores is decades greater than ambient bulk fluid flow ~ 10-8-10-7m/s characteristic of natural convective geothermal systems. Such rare high flow locales must be found. While existing crustal surveys generally fix resource temperatures T with nominal 500m spatial resolution, as yet no survey technique provides data to locate production-grade flow rates at nominal ℓ ~ 50m spatial resolution. In consequence, many costly unproductive wells are drilled before sites of productive local geothermal fluid bulk flow v = φv0 ~ 10-2 m/s are located. A seismic survey methodology sensitive to crustal flow v = φv0 at ℓ ~25m resolution has evolved from multi-channel seismic reflection technology applied to the production of shale hydrocarbons. The field-proven seismic flow-imaging methodology can be adapted from sedimentary terrains to volcanic terrains through appropriate seismic refraction means for generating effective seismic velocity models for convective geothermal flow volumes. Numerical simulation for characteristic ambient crustal heterogeneous flow distributions v = φv0 at ℓ ~50m resolution shows that travel-time data for seismic source energy refracted from deep geothermal wellbores to surface seismic sensor arrays can replicate the multi-channel seismic flow-imaging capability demonstrated for shale formations. Multi-channel seismic reservoir flow monitor detection and mapping of flow-induced seismic emissions in shale formations can be adapted to achieve similar mapping of convective geothermal flow system structures with v = φv0 ~ 10-2m/s at nominal ℓ ~ 50m resolution. Such seismic emission flow structure maps can greatly reduce the uncertainty and hence the cost of drilling effective production wells at brownfield sites and assessment/development wells at greenfield sites