Search

Intermediate-depth earthquakes (>40 km) have been observed beneath the central Himalayas over decades, with little known about their nature and characteristics. Here, we apply a state-of-the-art systematic processing routine, starting from continuous waveform data, to obtain the most comprehensive high-quality earthquake catalog with a focus on the intermediate-depth seismicity beneath the central Himalayas. We construct a catalog containing 414 robust earthquake locations with depths ranging from 40 to 110 km spanning from late 2001 till mid-2003. We calculate earthquake magnitudes in a consistent way and obtain values ranging between ML 0.8 and 4.5 with a magnitude of completeness of Mc 2.4. This information allows us to study the spatiotemporal characteristics of the seismicity in great detail. Earthquakes mainly take place in a cluster, consisting of two linear segments at ca. 35° azimuth difference, situated beneath the high Himalayas in NE Nepal and adjacent S. Tibet. Seismicity there does not feature any mainshock-aftershock patterns but presents a few sequences with potential seismicity migration rates compatible with linear or diffusive migration. This result, along with previous studies in the lower Indian crust, allows interpreting these events as related to metamorphic reactions involving dehydration processes. However, given the geodynamic context, a tectonic interpretation with a dextral basement fault zone propagating beneath the Himalaya and continuing as a westward propagating tear fault would also be possible. This represents a continuous fault zone from the deep crust in S. Tibet, across the Himalaya along the Dhubri-Chungthang fault zone (DCFZ) to the Shillong plateau, which could be an inherited tectonic feature.

The 1714 Bhutan earthquake was one of the largest in the Himalaya in the last millennium. We show that the surface rupture caused by this earthquake extended further to the east than previously known, it was at least 175 km long, with slip exceeding 11 m at our study site. The age of the surface rupture was constrained by a combination of radiocarbon and traditional optically stimulated luminescence dating of affected river sediments. Computations using empirical scaling relationships, fitting historical observations and paleoseismic data, yielded a plausible magnitude of Mw 8.1 ± 0.4 and placed the hypocentre of the 1714 Bhutan earthquake on the flat segment of the Main Himalayan Thrust (MHT), the basal décollement of the Himalayan orogen. Calculations of Coulomb stress transfer indicate that great earthquakes along the leading part of the MHT would cause surface rupture. In contrast, distal earthquakes may not immediately trigger surface rupture, although they would increase the stresses in the leading part of the MHT, facilitating future surface-rupturing earthquakes. Frontal earthquakes would also transfer stress into the modern foreland basin facilitating southward propagation of the MHT as a blind basal décollement. In conclusion, studies of surface-rupturing events alone likely underestimate the seismic slip along the Himalayan megathrust.

We present results from a joint inversion of new seismic and recently compiled gravity data to constrain the structure of a prominent geophysical anomaly in the European Alps: the Ivrea Geophysical Body (IGB). We investigate the IGB structure along the West-East oriented Val Sesia profile at higher resolution than previous studies. We deployed 10 broadband seismic stations at 5 km spacing for 27 months, producing a new database of ∼1000 high-quality seismic receiver functions (RFs). The compiled gravity data yields 1 gravity point every 1–2 km along the profile. We set up an inversion scheme, in which RFs and gravity anomalies jointly constrain the shape and the physical properties of the IGB. We model the IGB’s top surface as a single density and shear-wave velocity discontinuity, whose geometry is defined by four, spatially variable nodes between far-field constraints. An iterative algorithm was implemented to efficiently explore the model space, directing the search toward better fitting areas. For each new candidate model, we use the velocity-model structures for both ray-tracing and observed-RFs migration, and for computation and migration of synthetic RFs: the two migrated images are then compared via cross-correlation. Similarly, forward gravity modeling for a 2D density distribution is implemented. The joint inversion performance is the product of the seismic and gravity misfits. The inversion results show the IGB protruding at shallow depths with a horizontal width of ∼30 km in the western part of the profile. Its shallowest segment reaches either 3–7 or 1–3 km depth below sea-level. The latter location fits better the outcropping lower crustal rocks at the western edge of the Ivrea-Verbano Zone. A prominent, steep eastward-deepening feature near the middle of the profile, coincident with the Pogallo Fault Zone, is interpreted as inherited crustal thickness variation. The found density and velocity contrasts of the IGB agree with physical properties of the main rock units observed in the field. Finally, by frequency-dependent analysis of RFs, we constrain the sharpness of the shallowest portion of the IGB velocity discontinuity as a vertical gradient of thickness between 0.8 km and 0.4 km.

We use a novel technique named global-phase seismic interferometry (GloPSI) to image the lithospheric structure, and in particular the Moho, below two parallel north-south transects belonging to the GANSSER network (2013–2014). The profiles cross the Himalayan orogenic wedge in Bhutan, a tectonically important area within the largest continent-continent collision zone on Earth that is still undergoing crustal thickening and represents a challenging imaging target for the GloPSI approach. GloPSI makes use of direct waves from distant earthquakes and receiver-side reverberations with near vertical incidence. Reflections are isolated from earthquake recordings by solving a correlation integral and are turned into a reflectivity image of the lithosphere below the arrays. Our results compare favorably with first-order features observed from a previous receiver function (RF) study. We show that a combined interpretation of GloPSI and RF results allows for a more in-depth understanding of the lithospheric structure across the orogenic wedge in Bhutan.

Nepal, located above the convergent India-Eurasia plate boundary, has repeatedly experienced devastating earthquakes. During the 2015 magnitude 7.8 Gorkha earthquake, an often-reported experience was that people were not aware of the threatening seismic hazard and had an insufficient level of preparedness. An important source of the problem is that earthquake-related topics are not part of the school curriculum. Earthquake education reaching a broad group of the population early in their lives is therefore strongly needed. We established an initiative in Nepal to introduce seismology in schools, with a focus on education and citizen seismology. We have prepared educational materials adapted to the Nepali school system, which we distributed and also share on our program’s website: http://seismoschoolnp.org. In selected schools, we also installed a low-cost seismometer to record seismicity and to allow learning-by-doing classroom activities. Our approach was very well received and we hope it will help make earthquake-safe communities across Nepal. The seismic sensor which we installed in schools is a Raspberry Shake 1D (RS1D), this was selected based on its performance in laboratory tests and suitability for the field conditions. At a test site in Switzerland we were able to record magnitude 1.0 events up to 50 km distance with a RS1D. In Nepal, 22 such seismometers installed in schools create the Nepal School Seismology Network providing online data openly. The seismometer in each school allows students to be informed of earthquakes, visualize the respective waveforms, and estimate the distance and magnitude of the event. For significant local and regional events, we provide record sections and network instrumental intensity maps on our program’s website. In 6 months of network operation, more than 194 local and teleseismic earthquakes of M ≥ 4 have been recorded. From a local and a global catalog, complemented with our own visual identifications, we have provided an earthquake wave detectability graph in distance and magnitude domain. Based on our observations, we have calibrated a new magnitude equation for Nepal, related to the epicentral distance D [km] and to the observed peak vertical ground velocity PGVV [μm/s]. The calibration is done to best fit local catalog magnitudes, and yields the following equation: M = 1.05 × log10(PGVV) + 1.08 × log10(D) + 0.75.