Source document: https://docs-lawep.s3.us-east-2.amazonaws.com/1697809424599.pdf (Mauritius)
Solid Earth, 11, 2557–2568, 2020 https://doi.org/10.5194/se-11-2557-2020 © Author(s) 2020.

This work is distributed under the Creative Commons Attribution 4.0 License.
Seismic gaps and intraplate seismicity around Rodrigues Ridge (Indian Ocean) from time domain array analysis Manvendra Singh1,2 and Georg Rümpker2 1Mauritius Oceanography Institute, Avenue des Anchois, Morcellement de Chazal, Albion, Mauritius 2Institute of Geosciences, Goethe University Frankfurt, Frankfurt am Main, Germany Correspondence: Manvendra Singh (msingh@moi.intnet.mu) Received: 17 April 2020 – Discussion started: 27 April 2020 Revised: 26 October 2020 – Accepted: 10 November 2020 – Published: 23 December 2020 Abstract.
Rodrigues Ridge connects the Réunion hotspot track with the Central Indian Ridge (CIR) and has been suggested to represent the surface expression of a sub- lithospheric ﬂow channel.

From global earthquake cata- logues, the seismicity in the region has been associated mainly with events related to the fracture zones at the CIR.

However, some segments of the CIR appear void of seismic events.
Here, we report on the seismicity recorded at a tem- porary array of 10 seismic stations operating on Rodrigues Island from September 2014 to June 2016.

The array analy- sis was performed in the time domain by time shifting and stacking the complete waveforms.

Event distances were es- timated based on a 1-D velocity model and the travel time differences between S and P wave arrivals.

We detected and located 63 new events that were not reported by the global networks.
Most of the events (51) are located off the CIR and can be classiﬁed as intraplate earthquakes.

Local magni- tudes varied between 1.6 and 3.7.

Four seismic clusters were observed that occurred to the west of the spreading segment of the CIR.

The Rodrigues Ridge appeared to be aseismic during the period of operation.
The lack of seismic activity along both Rodrigues Ridge and the sections of the CIR to the east of Rodrigues may be explained by partially molten upper-mantle material, possibly in relation to the proposed material ﬂow between the Réunion plume and the CIR.

1 Introduction The theory of plate tectonics is largely supported by the oc- currence of seismicity along (oceanic) ridge systems, sub- duction zones, and mountain belts.
However, a substantial part of seismicity occurs off the ridge axes within the oceanic basins.

Such events are usually denoted as intraplate earth- quakes (Krishna et al., 1998; Okal, 1983).

Gutenberg and Richter (1941) made the ﬁrst attempt to document these in- traplate events.

A more systematic approach was adopted by Sykes and Sbar (1973) and later by Bergman and Solomon (1980), where they compiled a catalog of 159 oceanic in- traplate events between 1939 and 1979.
Assessing the seis- micity inside the oceanic environment has always been a challenge as most of the permanent stations are located on the islands where the attenuation of the seismic signal along the oceanic path further reduces the detection capability.

An- other factor contributing to this problem is swell-generated noise in the band of 1–5 s (Okal, 1983).
Rodrigues Island is found at the eastern extremity of the Rodrigues Ridge near 63◦25 E and 19◦42 S, approximately 650 km east of Mauritius.

The east–west trending Rodrigues Ridge was formed between 7 and 10 Ma (Dyment et al., 2007; Saddul et al., 2002).
It is suggested that it repre- sents the surface expression of the interaction of the Réu- nion Hotspot with the Central Indian Ridge (CIR) through a sublithospheric ﬂow channel (Duncan and Hargraves, 1990; Dyment et al., 2007; Bredow et al., 2017) though Conrad and Behn (2010) and Becker and Faccenna (2011) suggest deeper mantle circulation as the possible mechanism of the for- mation of Rodrigues Ridge.
Samples collected from several sites along the ridge suggest that it is composed of basaltic Published by Copernicus Publications on behalf of the European Geosciences Union.

2558 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge Figure 1.

Tectonic setting of Rodrigues Ridge–Central Indian Ridge (CIR) region.

The ridge axis is shown by a solid line (red).

Event locations from USGS catalogue (January 2000–August 2020) are shown in yellow.
MCFZ stands for Marie-Celeste fracture zone, and EFZ stands for Egeria fracture zone.

The gaps in seismicity (GAP 1 and GAP 2) along the section of CIR between MCFZ and EFZ have been marked by white ellipses.

rocks (Dyment et al., 2001).

As Rodrigues Island is approxi- mately 250 km from the active southern part of CIR, seismic- ity around it is generally characterized by events recorded along this ridge system (Fig.

1).
1).

The largest magnitudes, 6.3 and 6.7, were recorded on 16 August 2010 and 26 July 2012, respectively, at the CIR near Rodrigues.

It was also observed that earthquakes with magnitudes 6 or greater are concen- trated along Marie-Celeste fracture zone (MCFZ) only.

Kr- ishna et al.

(1998) reported six events between October and November 1984 near the ridge, on the east of the spreading segment, between the Egeria fracture zone (EFZ) and MCFZ.

Similarly, Bergman et al.


(1984) have reported a large num- ber of “off-ridge” earthquakes in the region of the Southeast Indian Ridge.

Interestingly, there are two segments along the CIR, between MCFZ and EFZ, for which the global cata- logues are void of any seismicity (denoted GAP 1 and GAP 2 in Fig.

1).
1).

In this study we use seismological array techniques (Har- jes and Henger, 1973; Husebye and Ruud, 1989; Rost and Thomas, 2002) to characterize the seismic activity around the region of the Rodrigues Ridge and to conﬁrm the seismic gaps and provide possible explanations for them.

The data for this study were collected from temporary deployment of a seismic array on Rodrigues Island between September 2014 and June 2016 (Fig.

2).

Figure 2.
2).

Figure 2.

Locations of individual stations in Rodrigues array (yel- low triangles).

Red triangles denote permanent stations ROCAM and RODM.

Station RODM was operational between 10 November 2010 and 7 September 2014.

2 Methods and data analysis 2.1 Array conﬁguration In order to study the seismicity around Rodrigues Ridge, we deployed an array of 10 seismic stations on the island of Rodrigues located at about 19◦42 S and 63◦25 E approxi- mately 250 km west of the CIR.
Rodrigues array design is based on classical nine-element arrays that use three and ﬁve seismic sensors located along two concentric rings, respec- tively, with an additional sensor placed in the center.

The ben- eﬁt of this conﬁguration is that with irregular sensor spacing it provides a relatively sharp maximum of the array response function (Haubrich, 1968).
For Rodrigues array, we deployed 10 sensors in a similar conﬁguration with a ∼ 1.5 km ra- Solid Earth, 11, 2557–2568, 2020 https://doi.org/10.5194/se-11-2557-2020 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge 2559 dius of the inner ring and an outer ring radius of ∼ 2.5 km (Fig.

3).

The ﬁnal locations of the individual stations were then chosen according to the local conditions on the island (such as accessibility by roads, etc.).
Each station consisted of a MARK L-4C-3D (1 Hz) sensor and an Omnirecs CUBE data logger recording at a sampling rate of 100 Hz.

The relatively large aperture of the array (∼ 5 km) was chosen based on events listed in the USGS catalogue, which are located at the CIR near Rodrigues and were also recorded at the permanent station RODM.

This GEOSCOPE station was relocated during the course of the array deployment.

The dominant frequency of these events is close to about 2 Hz.
However, it became clear later that the dominant frequency of most earthquakes recorded by the newly installed array is ap- proximately 5 Hz.

Therefore, we decided to perform the array analysis in the time domain, or equivalently for a wide fre- quency range to reduce possible ambiguities resulting from side lobes in the array response function (Fig.

3), as will be discussed further below.

A similar approach was also used by Leva et al.

(2019).
(2019).

2.2 Epicentral distance and origin time For regional earthquakes, the slowness (apparent velocity) cannot be used to determine the epicentral distance of an event, as the ray path is mainly conﬁned within the upper- most mantle and the depth variations of velocity are not well constrained.

We, therefore, use the arrival-time difference be- tween the S and P waves in conjunction with a simpliﬁed 1-D velocity model of the crust and upper mantle to approx- imately determine the epicentral distance.
In this model, we keep the hypocentral depth ﬁxed at 6 km, and the ray path corresponds to a head wave (Fig.

4).

From receiver function analysis (Fontaine et al., 2015), a Moho depth of 10 km has been determined beneath Rodrigues Is- land.

We ﬁx the crustal thickness in our laterally homoge- neous model to this value, such that the thickened oceanic crust is accounted for on the receiver side leg of the ray path.
For the P wave velocity in the crust, we assume V C = P 6.1 km s−1 and in the mantle we assume V M = 7.9 km s−1 P with a VP /VS ratio of 1.80, as suggested by results of Chris- tensen (2004), Kong et al.

(1992), Wolfe et al.

(1995), and Grevemeyer et al.

(2013).

We explored the inﬂuence of these parameters on the determination of the epicentral distances for the events that occur at the distances of ∼ 120 and ∼ 265 km from the array.

As shown in Figs.
As shown in Figs.

S1 and S2 in the Supplement, results are most sensitive to variations of the P wave velocity in the mantle and VP /VS, owing to its con- tribution to the ray path.

The ﬁnal value for the epicentral distance of an event (with respect to the center of the array) is obtained by taking the mean of the distance values determined at all stations of the array.

Stations for which the recordings do not exhibit a clear onset for either the P or S phase are discarded from this cal- culation.
The standard deviation (SD) is used to deﬁne the error of the distance calculation.

To determine the event origin time, torigin, the travel time, ttt, was calculated from the distance obtained from the S and P wave travel time difference for the central (reference) station.

By subtracting this value from the manually picked P wave arrival time, tP , we obtain torigin = tP − ttt, where it is assumed that tP is given as an absolute time value.
2.3 Magnitude determination In order to calculate the magnitude of an event, we ac- count for the sensitivity of the recording system (CUBE data logger and 1 Hz MARK sensor) and integrate the ve- locity seismogram, which is then convolved with the Wood– Anderson transfer function to obtain the ground displace- ment in nanometers.
Magnitudes are determined based on the maximum (absolute) amplitude of the horizontal components using all available stations for which the recordings show a clear (dominant) S phase.

Again, the mean and the SD of the magnitude value is calculated.
As most of the events are lo- cated well within a 1000 km radius, we use the relation given by Havskov and Ottemöller (1999) in SEISAN package to determine the local magnitude ML = log10(A) + 1.1 × log10() + 0.00189 × − 2.09, (1) where A is the amplitude (in nm) and is the epicentral distance (in km).

2.4 Beamforming and array transfer function In seismic array analysis, the slowness and the back azimuth of an event are determined by beamforming.
Assuming a plane wave front of horizontal slowness s0 moving across the array with an apparent velocity(va = 1/|s0|), the waveform at station j is given by wj (t) = w (t − r j · s0), (2) where r j deﬁnes the position of the station with respect to a suitable coordinate system.
For an array consisting of M stations, the beam energy is calculated from the trace amplitudes within a suitable time window deﬁned by t1 and t2 according to (e.g., Harjes and Henger, 1973; Rost and Thomas, 2002) E = t2 t1 y2(t) dt = t2 t1 1 M M j =1 2 wj (t + r j · s) dt, (3) where s denotes the (trial) slowness for the current beam.

The beam energy reaches a maximum, if s = s0.


The back azimuth is then obtained from slowness components accord- ing to tan−1(s0x/s0y), where x and y correspond to the east- ern and northern components, respectively (see details given below).

https://doi.org/10.5194/se-11-2557-2020 Solid Earth, 11, 2557–2568, 2020 2560 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge Figure 3.

Conﬁguration of the seismic array deployed on Rodrigues Island (a).
Station locations are indicated by blue triangles.

The green circles represent inner and outer rings of the array with a radius of ∼ 1.5 and ∼ 2.5 km, respectively, from the central (reference) station.

The array transfer function of the Rodrigues array is shown at 2 Hz (b) and 5 Hz (c).

The inner and the outer red rings correspond to an apparent velocity of 8 and 5 km s−1, respectively.
In the real-data analysis, the apparent velocity at which the maximum occurs is used as a ﬁrst indication to discriminate between crustal and upper-mantle ray paths.

In view of Eq.
In view of Eq.

(1) and using Parseval’s theorem in combi- nation with the shift theorem of Fourier transforms, we can write the energy in the form E = 1 2π ω2 ω1 ¯y2(ω) dω = 1 2π ω2 ω1 ¯w2(ω) C (ω, s − s0) dω, (4) where the bar denotes Fourier-transformed functions and the array transfer function, C (ω, s − s0), in the frequency– slowness domain is given by (e.g., Schweitzer, 2002) C (ω, s − s0) = 1 M M j =1 eiω r j · (s−s0) 2 .
(5) It is further assumed that ¯w = 0 is outside of the range of in- tegration used in Eq.

(4).

The array transfer function deﬁnes the sensitivity and resolution of the array for seismic signals with frequency ω (Fig.

3).

2.5 Data example All the events for this study were detected by manual inspec- tion of hourly traces from all stations of the Rodrigues array using the SEISAN package (Havskov and Ottemöller, 1999).
A time window of 120 s was cut around each event using the GIPPtool software (http://www.gfz-potsdam.de, last access: August 2019).

As shown in Fig.

5, we estimate the epicentral distance from the arrival-time differences of the S and P waves.

The STA / LTA ratio of the Z-component trace is calculated to aid in the manual picking of the two arrivals.

All S picks were made independently and are based on visual inspection of horizontal components as well as vertical components where necessary.
The mean of the P phase picking time, calculated from all the picks, is used to determine the time window to calculate the beam energy later in the array analysis.
For this example, the arrival-time difference is 13.3 s, which corre- Solid Earth, 11, 2557–2568, 2020 https://doi.org/10.5194/se-11-2557-2020 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge 2561 analysis corresponds to a broadband energy stack and sup- presses the effects of unwanted side lobes, but there is an ad- ditional beneﬁt: the frequency domain approach usually re- quires selection of a common time window for all traces be- fore the Fourier transform is applied.
In cases of signiﬁcantly different arrival times of the phase to be analyzed (e.g., due to a large aperture of the array), a relatively wide common time window has to be selected such that the cut waveforms of individual traces may be signiﬁcantly different.

In the time domain, however, we can time shift the traces (with respect to the trial slowness value) before cutting and stacking, which may then be performed within a much narrower time win- dow.
This approach ensures that only the relevant waveform is contained within the stack, provided that the correct time shift has been applied.

Our analysis involves applying time shifts to all the traces with respect to reference trace (no.

1) for different values of slowness, deﬁned by a grid, and calculating the energy of the resulting stacked trace or beam within the time window of the ﬁrst arrival, which is selected using mean of the P wave picks on all the usable traces for a given earthquake.
Only vertical (Z) component traces from all stations were used for the beamforming.

Figure 6a shows the Z component traces for an earthquake on 6 April 2015 originating at 20:25:42 UTC.

A band pass between 2 and 10 Hz is applied.

The trace number on the ver- tical axis corresponds to the respective station number.

The amplitudes of trace 3 are set to zero as it did not show a clear event.

At ﬁrst, a relatively wide time window (∼ 2 s) is automat- ically chosen to zoom-in on the P phase (Fig.

6b).
6b).

Then, a narrower time window (∼ 0.6 s or about three dominant periods wide, as indicated by the red vertical lines) is au- tomatically selected with respect to reference trace 1 based on the P phase picking on all the traces.

Time shifts are ap- plied to all other (normalized) traces in correspondence to the slowness values of the grid which is used in the energy stack (Fig.

6c).

The stack, however, is only calculated from the trace segments within the narrow time window.
A predeﬁned grid with slowness from −0.3 to 0.3 s km−1, equally spaced over 248 × 248 = 61 504 points is used for calculating the appropriate time shifts for each trace and the resulting energy.

The red circle in Fig.

6c denotes the maxi- mum beam energy at s0 = (0.06, 0.07) s km−1, which corre- sponds to a back azimuth =∼ 38◦ and va =∼ 10.6 km s−1.
We assume that this relatively large value for the apparent ve- locity results from steepening of the ray path at shallow depth beneath Rodrigues Island (where the propagation velocity is smaller than 6.1 km).

The corresponding shifted traces and the beam (by summing the amplitudes of all traces) are shown in Fig.

6d.

The values of slowness components thus obtained, were used to calculate the corresponding apparent velocity, va and Figure 4.
Cartoon depicting the model employed to estimate epicentral distances from travel-time differences between S and P wave arrival.

The simple 1-D velocity model of the oceanic crust and uppermost mantle consists of a uniform layer overlying a uni- form half-space.

The parameters are set as follows: thickness of the crust, d = 10 km; hypocentral depth, h = 6 km; P wave veloc- = 6.1 km s−1; P wave velocity in the mantle, ity in the crust, V C P = 7.9 km s−1; VP /VS = 1.80 in both the crust and mantle.
Fig- V M P ures S1 and S2 in the Supplement show the inﬂuence of these pa- rameters on the determination of the epicentral distance.

Figure 5.

Estimating the epicentral distance from the arrival-time difference between S and P waves at the reference station 1.

The vertical, north, and east components are shown by traces 1, 2, and 3, respectively (a).

The vertical red lines show manually picked P and S phases.

The picking is aided by calculating the STA / LTA of the Z component (b).
sponds to a distance of approximately 119 km (Table S1).

Fi- nally, the magnitude of the event is calculated from Eq.

(1).

Here, we obtain a magnitude of 3.0.

In the following, we apply this methodology to all events detected that exhibit a clear P and S wave onset.

Conventionally, array analyses are performed in the fre- quency domain, which is computationally advantageous as the energy stacking can be limited to the dominant frequency or a narrow frequency band.
As explained above, however, we perform the array analysis in the time domain to include the complete waveform of the ﬁrst arrival.

The time domain https://doi.org/10.5194/se-11-2557-2020 Solid Earth, 11, 2557–2568, 2020 2562 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge Figure 6.

Array analysis of a regional event near Rodrigues.

(a) Vertical-component traces.

Station 3 was not recording at the time of the event.

(b) Zoom-in of the P waveforms.


After applying the time shift corresponding to a trial slowness s, the beam energy is calculated from squared amplitudes within the time window enclosed by the red circle.

(c) Beam energy as function of horizontal slowness components sx and sy .

The maximum is marked by a red line and corresponds to the slowness s0.

White circles denote apparent velocities, va, of 4, 6, and 8 km s−1.

(d) Shifted traces and beam (red trace) corresponding to the maximum indicated in (c).
The back azimuth of this event is determined at about 38◦ with an apparent velocity of 10.57 km s−1.

the back azimuth, of the event according to 3 Comparison with event locations from global networks va = 1/ s2 0x + s2 0y and = (180◦/π) tan−1(s0x/s0y), (6) where s0x and s0y are the values of slowness vector obtained for the beam with maximum energy.
To estimate the error in back azimuth, δ, we deﬁne the area enclosed by the contour at a level of 95% of the maxi- mum beam energy as a conﬁdence region.

The error in back azimuth in terms of kilometers (δkm) is then given by δkm = (π/180◦)δ (7) where is the epicentral distance (in km) and δ is given in degrees.
The larger earthquakes along the CIR picked up by the global networks are listed in the earthquake database provided by USGS (https://earthquake.usgs.gov/earthquakes/search/, last access: August 2016).

Thirty events reported in the catalogue were also recorded by Rodrigues array.

Six of these events exhibit sufﬁcient signal-to-noise ratio to perform the array analysis for a comparison with the USGS-reported locations.
Geographical coordinates, origin times, and magnitudes for the six events obtained from the array analysis are shown in Table 1.

They can be compared with the results provided by the USGS catalogue (Table 2).

Magnitudes obtained from the array analysis are slightly lower than those reported by USGS.
This may be in part due to different magnitude scales Solid Earth, 11, 2557–2568, 2020 https://doi.org/10.5194/se-11-2557-2020 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge 2563 Figure 7.

Estimating the epicentral distance from the arrival-time difference between S and P waves at reference station 1.

The vertical, north, and east components are shown by traces 1, 2, and 3, respectively (a).

The vertical red lines show manually picked P and S phases.
The picking is aided by calculating the STA / LTA of the Z component (b).

used (local magnitude versus body wave magnitude).

Ris- tau (2009) compared different magnitude scales for New Zealand earthquakes, where mb was lower than ML for deep focus (> 33 km) events but the results were fairly consistent for shallow (≤33 km) earthquakes.

In the current study, we consider that the data are not sufﬁcient to derive a relation- ship between mb and ML.
Amplitude variations related to the different radiation directions relevant for regional (at Ro- drigues Island) and teleseismic recordings may also play a role.

In addition, the amplitudes for the dominantly horizon- tal Pn ray paths from the recordings at Rodrigues may be more affected by regional attenuation processes, as described further in the next section.

The array analysis for the six events was performed as de- scribed above.

An example is given in Figs.


7 and 8 for an event of 14 February 2015 at 07:08:59 UTC.

P and S phases are clearly visible and the dominant frequency of the P wave arrival for this event is about 3.5 Hz.

From the array analy- sis, we obtain a clear maximum for s0 = (0.07, 0.08) km s−1, which corresponds to a back azimuth ∼ 39◦ and apparent velocity va ∼ 9.5 km s−1.
The time shifts derived from the obtained slowness value according to δtj = δr j · s0 (where δr j corresponds the position of station j with respect to the reference station) yield a good alignment of the P wave sig- nal.

The mean distance of the event was calculated using dis- tances from all the stations that provided clear P and S phases.

Figure 7 shows the three-component seismogram and the STA / LTA trigger function to identify the onset times.
The arrival-time difference of 24.63 s corresponds to a dis- tance of 231 km (based on the model described above).

For the magnitude, we derive ML = 3.5.

For three events, we obtain small variations in back az- imuth, which are well within the error estimations.

Events on 24 November 2014 at 22:23:23 UTC, 2 April 2016 at 18:01:49 UTC, and 2 June 2016 at 21:18:10 UTC show dif- ferences of about 11, 9, and 7◦, respectively.
Array-derived distances from the reference station are generally smaller than those given in USGS catalogue, except for the event of 2 April 2016 at 17:53:21 UTC, where the distance obtained from array analysis is well within the error range.

We attribute these differences to local inhomogeneities not accounted for in the array analysis and to the simple 1-D model used for the distance estimates (in addition to possible errors in the global locations).
Figure S3 in the Supplement shows a comparison of the results obtained by the array anal- ysis with those provided by the USGS catalogue.

Generally, the results agree well.

On average, location differences are about 17 km, which is a reasonable value considering the un- certainties of the approach.

https://doi.org/10.5194/se-11-2557-2020 Solid Earth, 11, 2557–2568, 2020 2564 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge Figure 8.
Array analysis of an event of 14 February 2015 at 07:08:55 UTC that occurred close to the CIR, east of Rodrigues Island, which is listed in the USGS catalogue.

(a) Vertical-component traces.

(b) Zoom-in of the P waveforms.

(c) Beam energy as a function of horizontal slowness components sx and sy .

(d) Shifted traces and beam (red trace) corresponding to the maximum indicated in (c).

The back azimuth of this event is determined at about 39◦ with an apparent velocity of 9.5 km s−1.

Table 1.
Table 1.

Events from USGS catalogue located using the Rodrigues array.

Date, time (yy/mm/dd, hh:mm:ss) Lat.

(◦) Long.

ML ± SD (◦) Dist.
± SD (km) 14/11/24, 22:23:23 15/02/14, 07:08:59 15/02/14, 11:13:33 16/04/02, 17:53:21 16/04/02, 18:01:49 16/06/02, 21:18:10 −18.404 −18.100 −18.210 −18.432 −18.472 −17.846 65.043 64.837 64.927 65.108 65.057 65.118 3.8 ± 0.1 3.5 ± 0.1 3.7 ± 0.1 3.4 ± 0.1 3.4 ± 0.1 3.7 ± 0.1 222.02 ± 3.23 231.27 ± 5.01 228.26 ± 2.22 225.21 ± 5.66 218.24 ± 3.89 271.93 ± 5.75 BAZ (◦) 49.39 39.48 43.27 51.14 51.14 40.61 BAZ error Min (◦) Max (◦) 46.73 35.55 38.67 46.73 47.64 36.13 53.87 44.06 45.00 54.60 55.61 43.27 4 Results and discussion Using the array technique, we were able to detect and locate 63 earthquakes in the Rodrigues–CIR region that are not re- ported by the USGS catalogue (Fig.
9).

The details of all the events, such as event location, origin time, and magnitude are summarized in Table S1 (Supplement).

The magnitudes of these events range from 1.6 to 3.7 and are spread out in a re- gion of radius up to 600 km from the array.

The nearest event to Rodrigues Island (ML 1.6) occurred on 12 February 2016 at 19:21:45 UTC to the north of the island at a distance of about 36 km.
Solid Earth, 11, 2557–2568, 2020 https://doi.org/10.5194/se-11-2557-2020 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge 2565 Table 2.

Details of events from USGS catalogue for comparison.

Date, time (yy/mm/dd, hh:mm:ss) Lat.

(◦) 14/11/24, 22:23:19 15/02/14, 07:08:55 15/02/14, 11:13:29 16/04/02, 17:53:18 16/04/02, 18:01:45 16/06/02, 21:18:07 −18.6335 −18.1178 −18.0478 −18.5407 −18.6435 −17.9694 Long.
(◦) 65.4567 64.9817 65.1849 65.0872 65.4247 65.4337 mb 4.8 4.7 4.8 4.5 4.8 4.7 Dist.

(km) 242.74 239.68 260.07 215.96 239.27 284.98 BAZ (◦) 60.78 42.62 44.99 53.24 60.61 47.54 Figure 9.

Locations of new events (and clusters) detected and located using array methods are shown by red circles.

The purple circles repre- sent event locations from the array analysis in comparison to the events (orange circles) from the USGS catalogue between September 2014 and June 2016.
The black bars pair the respective events.

Yellow circles denote events from the USGS catalogue for the same recording period that were not (clearly) detected at the Rodrigues array and that were therefore not used for a comparison.

The empty GAP 1 and GAP 2 have been demarcated by white ellipses.

The dashed white line marks the possible region of partially molten material below the eastern extension of Rodrigues Ridge.
CIR, Central Indian Ridge (solid red line); MCFZ, Marie-Celeste fracture zone; and EFZ, Egeria fracture zone.

Of the 63 events, 51 are located off the ridge axis and can therefore be classiﬁed as intraplate earthquakes.

A total of 12 events were located very close to the ridge axis.

A total of 24 events occur between back azimuths of 29◦ and 42◦ at a distance of about 120 km and exhibit magnitudes between 1.6 and 3.0 (Fig.

9).
9).

Almost all events in this back azimuth (20) occur from March to April 2015 (Fig.

10).

Magnitudes for this cluster of earthquakes are variable and do not follow a certain (main shock–aftershock) pattern (Cluster 1, Table S2 in the Supplement).

Detailed bathymetry indicates a step-like non-transform discontinuity at the southeastern end of this chain of events.
Three more clusters occur to the east, south, and north- west of Rodrigues Island at distances of about 265, 140, and 220 km with 6, 3, and 7 events, respectively.

This is further corroborated by Fig.

10, where events’ longitude and latitude are plotted as functions of time.

Also shown is the monthly number of events, which exhibits a clear maximum for March 2015 in relation to the activity of Cluster 1.
As Cluster 4 (∼ 220 km east of Rodrigues Island) is close to the ridge axis, events of this cluster are not considered as intraplate events.

As shown in the sensitivity test provided in the supple- mentary information (Figs.

S1 and S2), the location errors for Cluster 1 are much smaller as compared to Cluster 3, par- tially owing to their distance from the Rodrigues array.
Some inﬂuence of anisotropic velocity variation may also be pos- sible, as studies by Barruol and Fontaine (2013) and Scholz https://doi.org/10.5194/se-11-2557-2020 Solid Earth, 11, 2557–2568, 2020 2566 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge stresses and hence earthquakes.

Similarly, Collette (1974) and Turcotte (1974) suggested thermal contraction due to dif- ferent ages of the oceanic crust (Müller et al., 2008) as an- other possible mechanism.
This may be substantiated by the seismicity observed in clusters 2 and 3 as they occur along possible prolongations of fracture zones.

Mantle-derived car- bon dioxide discharge (e.g., Bräuer et al., 2003; Gold and Soter, 1984; Irwin and Barnes, 1980) can also provide an explanation for the intraplate seismic activity around the Rodrigues–CIR region, as most of the detected events occur in clusters 1 and 4, which are not related to fracture zones.
In continental regions, similar clustering has been associated with carbon dioxide discharge (e.g., Lindenfeld et al., 2012).

It is interesting to note that the two seismic gaps along the slow-spreading CIR, as indicated above, still show no seis- mic activity based on the new data.

Thus, these seismic gaps may indeed represent an anomalous section of the CIR that is deforming aseismically.
The anomalous nature of this re- gion is further corroborated by the fact that globally recorded events from the farther southern section of CIR are not de- tected by the array (see events south of −20◦ in Fig.

S3 in the Supplement).

This could be explained by a more extended region of partially molten material in the upper mantle that causes signiﬁcant attenuation of wave amplitudes for the cor- responding ray paths (Fig.

9), as also suggested by Mazzullo et al.

(2017).


In combination with the absence of any seis- mic activity along Rodrigues Ridge, this may be taken as ev- idence for the explanation of Rodrigues Ridge as a surface expression of the interaction of the Réunion Hotspot with the CIR through a sublithospheric ﬂow channel (Morgan, 1978; Dyment et al., 2007; Bredow et al., 2017).
Another explana- tion for the lack of seismicity in GAP 1 and GAP 2 could be excessive magmatism and a resulting relatively thin litho- sphere that does not support large earthquakes, as suggested by Cannat (1996) and Grevemeyer et al.

(2013).
(2013).

Based on hydro-acoustic data, similar seismic gaps were also observed along the intermediate-spreading Southeast Indian Ridge (Tsang-Hin-Sun et al., 2016) and the slow- spreading Mid-Atlantic Ridge (Escartin et al., 2008; Simao et al., 2010), where events were focused on ridge segment ends.
This uneven distribution of seismicity along CIR, thus suggests segment-scale variation in lithosphere structure, as also observed in Southeast Indian Ridge (Tsang-Hin-Sun et al., 2016) and Mid-Atlantic Ridge (Goslin et al., 2012).

5 Conclusions Figure 10.

Variation in longitude (a) and latitude (b) of the events detected and located using Rodrigues array.

Monthly distribution of the events is shown in (c).
The solid red line in (b) marks the events for Cluster 1 (∼ 120 km northeast of Rodrigues Island), Cluster 2 (∼ 140 km northwest of Rodrigues Island), and Cluster 3 (∼ 265 km south of Rodrigues Island).

et al.

(2018) suggest fast-axis directions trending east–west around the Rodrigues–CIR region.

Various mechanisms providing explanation to the cause of intraplate seismicity have been proposed previously.

De Long et al.
De Long et al.

(1977) suggested that different lithospheric ages across the fracture zones, as also observed in the Rodrigues– CIR region, are related to differential subsidence causing We installed a 10-station seismological array on Rodrigues Island to study the seismicity along a remote section of the CIR and nearby areas including Rodrigues Ridge.
The re- sults show that array analysis provides a valuable tool to study earthquake activity in oceanic regions which are rel- atively inaccessible otherwise.

The region around Rodrigues Island clearly shows evidence of intraplate seismicity.
Of the Solid Earth, 11, 2557–2568, 2020 https://doi.org/10.5194/se-11-2557-2020 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge 2567 63 events detected by Rodrigues array, the majority are lo- cated within Cluster 1 at a distance of ∼ 120 km from the island with back azimuths between 29 and 42◦.

The local magnitude (ML) of the events detected ranges between 1.6 and 3.7.

Three additional event clusters were identiﬁed.


Pos- sible explanations for the off-axis seismic activities are CO2 degassing from the mantle and differential thermal contrac- tion, whereas events in Cluster 4, being close to ridge axis, are not considered as intraplate events.
The lack of seismic activity along both Rodrigues Ridge and a section of the CIR to the east of Rodrigues Island (GAP 2) may be explained by partially molten upper-mantle material, possibly in rela- tion to the proposed material ﬂow from the Réunion plume and the CIR (Morgan, 1978; Dyment et al., 2007; Bredow et al., 2017).
This explanation is further supported by the ob- servation that relatively strong seismic events from the CIR, east to southeast of Rodrigues Island (which are listed in the USGS catalogue) are not detected by the array.

However, a detailed geodynamic model for the ridge–plume interaction in this region is still needed.
We anticipate that longer-term deployments of seismic arrays on Rodrigues and other re- mote islands of Mauritius, such as Agalega and St. Brandon, will provide further constraints on the seismic gaps along the CIR and the intraplate seismicity of the region.

Dedicated de- ployment of Ocean Bottom Seismometers or Hydrophones (OBS/H) at or near these targets is another option for future studies.

Data availability.

The data used in this research are currently re- stricted.
The data will be publicly available from GEOFON data archive of Deutsches Geoforschungszentrum Potsdam as from July 2021.

Supplement.

The supplement related to this article is available on- line at: https://doi.org/10.5194/se-11-2557-2020-supplement.

Author contributions.

MS and GR participated in ﬁeld work and data analyses, provided their interpretations, and wrote the article.

Competing interests.

The authors declare that they have no conﬂict of interest.

Disclaimer.
Disclaimer.

The study presented here was funded by the Mauri- tius Oceanography Institute, Mauritius, through the Ministry of Fi- nance and by the Deutsche Forschungsgemeinschaft (DFG) through a grant to GR.

Acknowledgements.

The instruments for the temporary array were provided by the Geophysical Instrument Pool Potsdam at the Deutsches Geoforschungszentrum, Potsdam.
We wish to thank the Government of Mauritius for allowing us to carry out the study, as well as the MOI and Rodrigues Regional Assembly (RRA) for the support provided during the ﬁeldwork.

We thank Jérôme Dyment of IPGP for his suggestions and comments regarding the bathymetry around the CIR region.

We also thank Frederik Link, Corrado Sur- manowicz, Olivier Pasnin, and Shane Sunassee for providing help during ﬁeld work at Rodrigues Island.
We thank the reviewers for the thorough reviews that helped to improve the paper.

Financial support.

This research has been supported by the German Research Foundation (grant no.

RU 886/7-1).

This open-access publication was funded by the Goethe University Frankfurt.

Review statement.

This paper was edited by Irene Bianchi and re- viewed by two anonymous referees.

References Barruol, G. and Fontaine, F. R.: Mantle ﬂow beneath La Réunion hotspot track from SKS splitting, Earth Planet.

Sci.
Sci.

Lett., 362, 108–121, https://doi.org/10.1016/j.epsl.2012.11.017, 2013.

Becker, T. W., and Faccenna, C.: Mantle conveyor beneath the Tethyan collisional belt, Earth Planet.

Sci.

Lett., 310, 453–461, https://doi.org/10.1016/j.epsl.2011.08.021, 2011.

Bergman, E. A. and Solomon, S. C.: Oceanic intraplate earthquakes: Implications for local and regional intraplate stress, J. Geophy.

Res.

Solid Earth, 85, 5389–5410, 1980.


Bergman, E. A., Nábˇelek, J. L., and Solomon, S. C.: An extensive region of off-ridge normal-faulting earthquakes in the southern Indian Ocean, J. Geophy.

Res.-Sol.

Ea., 89, 2425–2443, 1984.

Bräuer, K., Kampf, H., Sträuch, G., and Weise, S. M.: Isotopic trig- gered intraplate seismicity, J. Geophys.

Res., 108, 2070, https://doi.org/10.1029/2002JB002077, 2003.
(3He/4He, 13CCO2) of ﬂuid evidence Bredow, E., Steinberger, B., Gassmöller, R., and Dannberg, J.: How plume-ridge interaction shapes the crustal thickness pattern of the R éunion hotspot track, Geochem.

Geophy.

Geosy., 18, 2930– 2948, https://doi.org/10.1002/2017GC006875, 2017.

Cannat, M.: How thick is the magmatic crust at slow spreading oceanic ridges?, J. Geophys.

Res.-Sol.

Ea., 101, 2847–2857, https://doi.org/10.1029/95JB03116, 1996.
Christensen, N. I.: Serpentinites, peridotites, and seismology, Int.

Geol.

Rev., 46, 795–816, https://doi.org/10.2747/0020- 6814.46.9.795, 2004.

Collette, B. J.: Thermal contraction joints in a spreading seaﬂoor as origin of fracture zones, Nature, 251, 299–300, 1974.

Conrad, C. P. and Behn, M. D.: Constraints on lithosphere net rotation and asthenospheric viscosity from global mantle ﬂow models and seismic anisotropy, Geochem.

Geophy.

Geosy., 11, Q05W05, https://doi.org/10.1029/2009GC002970, 2010.
De Long, S. E., Dewey, J. F., and Fox, P. J.: Displacement history of oceanic fracture zones, Geology, 5, 199–202, 1977.

Duncan, R. A. and Hargraves, R. B.: 40Ar/39Ar Geochronology of basement rocks from the Mascarene plateau, the Chagos bank https://doi.org/10.5194/se-11-2557-2020 Solid Earth, 11, 2557–2568, 2020 2568 M. Singh and G. Rümpker: Seismic gaps and intraplate seismicity around Rodrigues Ridge and the Maldives ridge, in: Proceedings of Ocean Drilling Pro- gram, Scientiﬁc Results, Vol.
115, College Station, TX: Ocean Drilling Program, 43–51, 1990.

Dyment, J., Hémond, C., Guillou, H., Maia, M., Briais, A., and Gente, P.: Central Indian Ridge and Reunion hotspot in Ro- drigues area: Another type of ridge-hotspot interaction?, Eos Trans.

AGU, 82(47), Fall Meet.

Suppl., Abstract T31D–05, 2001.

Dyment, J., Lin, J., and Baker, E. T.: Ridge-hotspot interactions: What mid-ocean ridges tell us about deep earth processes, Oceanography, 20, 102–115, 2007.
Escartín, J., Smith, D. K., Cann, J., Schouten, H., Langmuir, C. H., and Escrig, S.: Central role of detachment faults in accretion of slow-spreading oceanic lithosphere, Nature, 455, 790–794, 2008.

Fontaine, F. R., Barruol, G., Tkalˇci´c, H., Wölbern, I., Rümpker, G., Bodin, T., and Haugmard, M.: Crustal and uppermost mantle structure variation beneath la Réunion hotspot track, Geophys.

J.

Int., 203, 107–126, 2015.


Grevemeyer, I., Reston, T. J., and Moeller, S.: Microseismicity of the Mid-Atlantic Ridge at 7◦ S–8◦15 S and at the Logatchev Massif oceanic core complex at 14◦40 N–14◦50 N, Geochem.

Geophy.

Geosy., 14, 3532–3554, 2013.

Gutenberg, B. and Richter, C. F.: Seismicity of the Earth, Geol.

Soc.

Am.

Spec.

Pap., 34, 1–126, 1941.

Gold, T. and Soter, S.: Fluid ascent through the solid lithosphere and its relation to Earthquakes, Pure Appl.
Geophys., 122, 492–530, 1984.

Goslin, J., Perrot, J., Royer, J. Y., Martin, C., Lourenço, N., Luis, J., Dziak, R. P., Matsumoto, H., Haxel, J., Fowler, M. J., and Fox, C. G.: Spatiotemporal distribution of the seis- micity along the Mid-Atlantic Ridge north of the Azores from hydroacoustic data: Insights into seismogenic processes in a ridge–hot spot context, Geochem.

Geophy.

Geosy., 13, https://doi.org/10.1029/2011GC003828, 2012.
Lindenfeld, M., Rümpker, G., Link, K., Koehn, D., and Batte, A.: Fluid-triggered earthquake swarms in the Rwenzori region, East African Rift – Evidence for rift initiation, Tectonophysics, 566, 95–104, https://doi.org/10.1016/j.tecto.2012.07.010, 2012.

Mazzullo, A., Stutzmann, E., Montagner, J. P., Kiselev, S., Maurya, S., Barruol, G., and Sigloch, K.: Anisotropic tomography around La Réunion island from Rayleigh waves, J. Geophys.

Res.-Sol.
Res.-Sol.

Ea., 122, 9132–9148, https://doi.org/10.1002/2017JB014354, 2017.

Müller, R. D., Sdrolias, M., Gaina, C., Morgan, W. J.: Rodriguez, Darwin, Amsterdam, ..., a second type of hotspot island, J. Geophys.

Res.-Sol.

Ea., 83, 5355–5360, 1978. and Roest, W. R.: Age, spreading rates, and spreading asymmetry of the world’s ocean crust, Geochem.

Geophy.

Geosy., 9, https://doi.org/10.1029/2007GC001743, 2008.

Okal, E. A.: Oceanic intraplate seismicity, Annu.

Rev.

Earth Plan.

Sci., 11, 195–214, 1983.


Ristau, J.: Comparison of magnitude estimates for New Zealand earthquakes: Moment magnitude, local magnitude, and teleseis- mic body-wave magnitude, B. Seismol.

Soc.

Am., 99, 1841– 1852, https://doi.org/10.1785/0120080237, 2009.

Thomas, applications, S. Methods https://doi.org/10.1029/2000RG000100, 2002.

Seismology: 40, Geophys., and and Array Rost, Rev.

C.: Saddul, P., Nelson, B. W., and Jahangeer-Chojoo, A.: Mauritius: A Geomorphological Analysis, vol.
3, Mahatma Gandhi Institute, Mauritius, 2002.

Scholz, J. R., Barruol, G., Fontaine, F. R., Mazzullo, A., Montagner, J. P., Stutzmann, E., Michon, L., and Sigloch, K.: SKS splitting in the Western Indian Ocean from land and seaﬂoor seismometers: Plume, plate and ridge signatures, Earth Planet.

Sci.

Lett., 498, 169–184, https://doi.org/10.1016/j.epsl.2018.06.033, 2018.

Schweitzer, J., Fyen, J., Mykkeltveit, S., and Kværna, T.: Manual of Harjes, H. and Henger, M.: Array-Seismologie, Z.
Geophys., 39, seismological observatory practice, 2002.

865–905, 1973.

Haubrich, R. A.: Array design, B. Seismol.

Soc.

Am., 58, 977–991, 1968.

Havskov, and Ottemoller, J. software, Seismol.

Res.

Lett., analysis https://doi.org/10.1785/gssrl.70.5.532, 1999.

L.: SeisAn Earthquake 532–534, 70, Husebye, E. and Ruud, B.: Array seismology – past, present and future developments, Observ.

Seismol., 23–153, 1989.


Irwin, W. P. and Barnes, I.: Tectonic relations of carbon dioxide discharges and earthquakes, J. Geophys.

Res., 85, 3115–3121, 1980.

Kong, L. S., Solomon, S. C., and Purdy, G. M.: Mi- croearthquake characteristics of a mid-ocean ridge along- J. Geophys.

Res.-Sol.

Ea., 97, 1659–1685, axis high, https://doi.org/10.1029/91JB02566, 1992.

Krishna, M. R., Verma, R., and Arora, S.: Near-ridge intraplate earthquakes in the Indian Ocean, Mar.

Geol., 147, 109–122, 1998.


Leva, C., Rümpker, G., Link, F., and Wölbern, beneath Fogo subcrustal earthquakes dence for injection, https://doi.org/10.1016/j.jvolgeores.2019.106672, 2019. induced J. Volcanol.

Geotherm.

Res., I.: Mantle volcano, Cape Verde: Evi- by magmatic 106672, 386, fracturing Simao, N., Escartin, J., Goslin, J., Haxel, J., Cannat, M., and Dziak, R.: Regional seismicity of the Mid-Atlantic Ridge: observa- tions from autonomous hydrophone arrays, Geophys.

J.

Int., 183, 1559–1578, 2010.


Sykes, L. R. and Sbar, M. L.: Intraplate earthquakes, lithospheric stresses and the driving mechanism of plate tectonics, Nature, 245, 298–302, 1973.

Tsang-Hin-Sun, E., Royer, J. Y., and Perrot, J.: Seismicity and ac- tive accretion processes at the ultraslow-spreading Southwest and intermediate-spreading Southeast Indian ridges from hydroa- coustic data, Geophys.

J.

Internat., 206, 1232–1245, 2016.


Turcotte, D. L.: Are transform faults thermal contraction cracks?, J. Geophys.

Res., 79, 2573–2577, 1974.

Wolfe, C. J., Purdy, G. M., Toomey, D. R., and Solomon, S. C.: Mi- croearthquake characteristics and crustal velocity structure at 29◦ N on the Mid-Atlantic Ridge: The architecture of a slow spread- ing segment, J. Geophys.

Res.-Sol.

Ea., 100(B12), 24449–24472, https://doi.org/10.1029/95JB02399, 1995.
Solid Earth, 11, 2557–2568, 2020 https://doi.org/10.5194/se-11-2557-2020