Introduction

Distributed acoustic sensing (DAS) is a new geophysical imaging technique that measures the ground strain (or strain rate) along a fiber-optic cable, using variations in the phase-coherent Rayleigh backscattering that are generated by a sequence of laser pulses (e.g.,Hartog 2017; Miah and Potter 2017; Lu et al. 2019; Lindsey et al. 2020). Although DAS was originally developed for the oil and gas industry, it has since been utilized to address a wide range of seismological research topics (e.g., Zhan 2020). It is a relatively inexpensive way to obtain spatially dense and temporally continuous measurements, when we can use fiber-optic cables that have already been installed for communication purposes. We can simultaneously measure the strain at many locations along the cable (meter-scale intervals) at a high-sampling rate (~ 1000 samples per second (sps)) by simply connecting an interrogator, which is a laser transmitter/receiver system, to one end of the cable. It has been reported that the noise levels of the strain measurements are sufficiently small to detect natural earthquakes (Lindsey et al. 2017). Wang et al. (2018) have demonstrated the similarity between DAS and seismograph strain measurements, and confirmed that the observed signal is axial strain along a fiber-optic cable in the ordinary frequency range for seismic observations (i.e., several Hertz). At lower frequencies, Ajo-Franklin et al. (2019) observed the surface waves from large teleseismic earthquakes, and Lindsey et al. (2020) demonstrated how to convert DAS records into broadband seismograms via a comparative analysis of DAS and seismometer signals in 1–120 s range. Further applications for subsurface structural investigations have been proposed, such as receiver function (Yu et al. 2019) and H/V spectral ratio (Spica et al. 2020a) analyses.

Ocean-bottom strain measurements are a promising research direction for future DAS observations since the installation of conventional geophysical sensors is difficult, and several studies have successfully installed DAS interrogators along submarine communication cables. Williams et al. (2019) used a telecommunication cable at shallow (< 40 m) depths offshore Belgium, and analyzed the spectral features of the ambient DAS signals. These signals mostly consisted of ocean surface gravity waves and Scholte waves, which reflect solid earth–ocean coupling. They also observed the seismic waves from a large (Mw 8.2) teleseismic (hypocentral distance of ~ 150°) earthquake, although no regional earthquakes were recorded due to the low seismicity. Sladen et al. (2019) also utilized a telecommunication cable offshore France for DAS observations to detect a small (M1.9) earthquake that occurred ~ 100 km from the cable. Shinohara et al. (2019) and Spica et al. (2020b) used a fiber-optic cable that was originally installed for ocean-bottom seismic sensors to monitor active seismicity offshore Tohoku region, Japan. In addition to observation of many small earthquakes, they showed that DAS observations could potentially resolve spatial heterogeneities in the very shallow seismic structure along the cable. Some ocean sediments possess shear-wave (S-wave) velocities as low as 300 m/s.

Here we used a fiber-optic cable that was originally used for ocean-bottom seismic observations off Cape Muroto, Nankai subduction zone, where the Philippine Sea Plate is subducting beneath western Japan at 3–5 cm/years. Many large (M8 +) historical earthquakes have occurred in this region at 100–200-year intervals; more than 70 years have now passed since the last large event, the 1946 Nankai Earthquake; therefore, various observation instruments have been installed in this region with the hope of obtaining comprehensive observations of the next expected large earthquake. Our fiber-optic cable was a part of such observation systems. Although the ambient seismicity in this region is low compared with other regions in Japan, various types of slow earthquake (e.g., Beroza and Ide 2011; Obara and Kato 2016), such as slow slip events and tectonic tremors, occur frequently in the region. The relationship between slow and large megathrust earthquakes is of great interests and the effectiveness of DAS observations in detecting such phenomena is important to evaluate. Our temporary observations confirm the feasibility of acquiring DAS observations for earthquake detection and analysis. Furthermore, we also demonstrate the potential to investigate deep-ocean water movement via DAS observations, since these observations mainly record temperature changes in the quasi-static, sub-seismic frequency range.

DAS observations and noise levels

The Japan Agency for Marine Earth Science and Technology (JAMSTEC) developed the Comprehensive Seafloor Monitoring System off Cape Muroto in 1997 (Momma et al. 1997). A submarine cable with six single-mode fibers connects the land station and ocean-bottom instruments (Fig. 1). The total cable length is 120 km, with two ocean-bottom seismometers deployed 70 and 100 km from the coast. The system was operated for more than 20 years before being ceased in 2019. The cable was not covered by sediments, except for a 2-km shallow section, at the time of installation; however, the current status of the seafloor cable is unknown. We used a 50-km section of the cable that extended from the coast for the DAS observation. No instruments are attached to this section. Figure 1 shows the regional bathymetry and estimated cable route. The bathymetry is complex in this region, and the cable extends to the south, passing across the edge of Muroto Basin and over Tosabae Ridge.

Fig. 1
figure 1

Locations of the seafloor fiber-optic cable and hypocenter distribution during the observation period. The bold line represents the cable route of the Comprehensive Seafloor Monitoring System off Cape Muroto, with the red section indicating the 50-km-long DAS observation section. The two small open triangles along the cable show the ocean-bottom seismometer locations. The triangles with different colors are the stations that were used for a comparison with the DAS observations: MURH and KTGH (Hi-net), UMJF (F-net), and MRA04 and MRG29 (DONET2). Red circles indicate the epicenters of earthquakes that occurred during the January 27–31, 2020 (JST), period, as determined by the Japan Meteorological Agency

We acquired temporary DAS observations during the January 27–31, 2020 (JST) period along this cable. The DAS system was a commercial product of AP Sensing (model N5200A) with the maximum measurement fiber length of 50 km. The DAS system is based on phase-based C-OTDR in 1550 nm wavelength class 1 laser, together with telecommunication grade low attenuation single mode fiber used in this experiment, allowing measurement of Rayleigh scattered waves from such a long distance. We used several different parameter settings. The most used parameter settings for our observations were a 500-sps temporal sampling rate, 5.10-m spatial sampling interval, and 40.4-m gauge length. There were 9788 channels in the 50-km section of cable. The size of the uncompressed records was ~ 1 TB/day, since the interrogator measured a 16-bit recording of the differential strain (strain rate) at each channel.

The sensitivity of the 9788 channels to ground strain varied along the cable. Some channels recorded almost nothing, even when clear seismic signals were recorded at other channels. These insensitive channels recorded low levels of ambient noise that corresponded to microseisms, which are discussed later in the text. We computed the standard deviation of the strain rate using a 5-min interval for the entire observation period to visualize the spatial heterogeneity of the channel sensitivity. Figure 2 shows the median values of the standard deviations for the entire observation period. The insensitive channels are those with very small standard deviations, which are supposed to be “uncoupled” from the ocean bottom. The observed amplitudes at channels 0–500 (first 2.5 km of cable from the coast) were often near the maximum of the dynamic range (\(-\pi\) to \(\pi\) in phase), and were likely saturated due to the skipping of some cycles by large noise, which made the analysis difficult. At channels 500–2200, the standard deviations were generally small with the exception of some channels, suggesting weak coupling. There were two sections with large standard deviations: channels 2200–3000 (~ 4 km length) and channels 3500–4000 (~ 2.5 km length). Above channel 4000, there were alternating sections of coupled and uncoupled channels. For the channels that were far from the coast, there was a smooth increase in the standard deviation with channel number. This increase likely arose from the high noise levels due to the attenuation of the laser pulses.

Fig. 2
figure 2

Standard deviation of observation at each channel. Median values of standard deviation of strain rate, measured for each 5-minute at each channels

We selected ten channels and computed their spectrograms, which consist of periodograms of the strain rate at 5-min intervals for the entire observation period (Additional file 1: Figure S1). Two of the ten channels (6000 and 7510) are examples of uncoupled channels, which only recorded small amplitudes during the entire observation period. Figure 3 summarizes the median of the spectral amplitude (5-min interval) for the entire observation period at these ten channels (Fig. 3a), together with the spectral amplitudes that were estimated in a similar manner using the accelerations at several nearby seismic stations (and their associated networks): MUTH and KTGH (Hi-net, a high-sensitivity seismic network with short period (1 Hz) sensors); UMJF (F-net, a broadband seismic network); and MRA04 and MRG29 (DONET2, an ocean-bottom network with broadband sensors) (Fig. 3b). All these seismic stations are maintained by the National Research Institute for Earth Science and Disaster Resilience, and the data are publicly available.

Fig. 3
figure 3

Comparison of the noise levels between the DAS and nearby seismometers. a Strain rate spectra at ten selected channels. Each line shows the median value of the periodograms calculated at 5-min intervals for the entire observation period (January 27–31, 2020 (JST)), which were computed at 5-min intervals during the entire observation period (Figure S1). b Acceleration spectra at five seismic stations. Solid and dashed lines are for the vertical- and horizontal-component seismometers, respectively. Bold gray lines show a standard noise model (Peterson 1993). MRA04 and MRG29 are DONET2 stations, and UMJF is a F-net station; each of these stations has a broadband sensor. We corrected the instrumental response for the MURH and KTGH spectra, which are Hi-net stations

The amplitudes of the acceleration spectra at the two DONET2 stations are ~ 10−6 m/s2/Hz1/2 above 2 Hz, whereas those of the DAS strain rate spectra are ~ 1 nanostrain/s/Hz1/2. These two levels are comparable if we assume plane waves with a velocity of ~ 1 km/s are producing this background noise (e.g., Benioff 1935). A peak at 0.3–0.5 Hz is visible in both the DONET2 and DAS observations, suggesting that similar ambient signals are being recorded above 0.1 Hz. This peak at 0.3–0.5 Hz is common in broadband ocean-bottom seismic observations (e.g., Araki et al. 2004; Suetsugu and Shuobara 2014), and is also visible at a coastal Hi-net station (MURH). This peak is small at inland stations (KTGH and UMJF), and shifted to 0.2 Hz, which is the secondary microseism peak that is commonly observed at many stations and in a standard noise model (Peterson 1993). Below 0.1 Hz, the DAS noise levels appear to be worse than the DONET2 noise levels. At ~ 0.05 Hz, the DAS noise levels also appear to be worse than the Hi-net noise levels after the instrument response correction, although they are lower at ~ 0.01 Hz. Moreover, DAS recorded similar peaks responsible to infragravity waves at ~ 0.02 Hz. We therefore expect our DAS observations to be qualitatively similar to the DONET2 observations above 0.01 Hz.

Example earthquake observations

Our DAS observations recorded several tens of earthquakes, despite relatively inactive regional seismicity (Fig. 1). Figure 4 shows three examples of observed earthquake signals. One of the most clearly recorded events was a M2.9 earthquake that occurred ~ 100 km from the cable. Figure 4a shows a comparison of the strain records at three channels, with each consisting of the integrated strain rate after bandpass filtering at 2–16 Hz. Both tiny P waves and large S waves are visible in most of the channels (Fig. 4d), with the exceptions being the uncoupled channels with small signal levels, as suggested from the differences in the ambient standard deviation of the strain rate (Fig. 2). The maximum signal amplitudes were ~ 10 to 50 nanostrain, with large variations among the channels. This large variation highlights the need to characterize each channel response when conducting a future quantitative analysis of these signals. The maximum amplitudes for the same event that were recorded at two DONET2 stations, MRA04 and MRG29, were ~ 50 and ~ 10 μm/s, which were located at epicentral distances of 56 and 131 km, respectively. A relation between the 10–50 nanostrain DAS observations and 10–50 μm/s ocean-bottom seismometer observation is explicable by assuming plane waves at a velocity of the order of 1 km/s, which is expected from Fig. 3.

Fig. 4
figure 4

Example earthquake signal observations. a A moderate-sized (M2.9) event on January 27, 2020, at 16:23:31.8 (UT). The observed strain at three DAS channels (top) and the horizontal velocity seismograms at two DONET2 stations (bottom) are shown; all of the waveforms have been bandpass filtered at 2–16 Hz. b A very small (M1.1) event on January 30, 2020, at 00:20:25.7 (UT). The observed strain at three DAS channels (top) and the horizontal velocity seismograms at two Hi-net stations (bottom) are shown; all of the waveforms have been bandpass filtered at 2–16 Hz. c A teleseismic Mw7.7 event on January 28, 2020, at 19:10:24 (UT) near Cuba. The observed strain at three DAS channels (top), and the horizontal velocity seismograms at a DONET station (MRG29) and F-net station (UMJF) (bottom) are shown; all of the waveforms have been bandpass filtered at 0.01–0.1 Hz. All of the velocity data are shown relative to the origin time of the event. d DAS records for all of the channels and the earthquake in a. The amplitudes are presented as the logarithm of the normalized absolute amplitude of the strain rate, and have been bandpass filtered at 2–8 Hz

The smallest earthquake recorded during these temporary observations was a M1.1 event that occurred 30–50 km from the cable (Fig. 4b). The signal was only visible at the channels with good coupling close to the source. Nevertheless, the signal-to-noise ratio was comparable with that of a nearby Hi-net station (MURH) at some channels. This similarity likely occurred because the dominant frequency of this event was 2–4 Hz, affected by the attenuation effects along the propagation path, whereas the coastal Hi-net station had relatively high noise levels in this frequency range.

The largest earthquake we observed was a teleseismic earthquake of Mw7.7 that occurred near Cuba on January 28, 2020 (Fig. 4c). The signal was not visible at high frequencies, but was clearly visible in the spectrograms below the microseism level at most of the channels (Additional file 1: Figure S1). The surface waves were detectable for ~ 1 h; however, the noise level was not good compared with the broadband seismometers, and even worse than the response-corrected Hi-net records. The amplitude differences among the DAS channels were large; channel 1290, which often observed large amplitudes at high frequencies, recorded small amplitudes at low frequencies. This again highlights the need to estimate each channel response as a function of frequency.

Quasi-static observations: periodic variations and migrating events

What can we observe using a DAS system at much lower frequencies? One may expect that the DAS observations would detect the tidal strain due to solid earth and ocean tides, since the DAS measures strain. Figure 5 shows the extremely low-frequency DAS records, with the mean strain rate measured at a 5-min interval over the entire observation period for all of the channels along the cable. We observed large spatial and temporal variations in the average strain rate above channel 1000. The observed variation at channels 1000–4000 exhibited a nearly daily cycle with an amplitude on the order of 0.1 nanostrain/s (10−10 strain/s); this is very large since the tidal deformation of the Earth’s surface is on the order of 10−12 strain/s. Noticeable variations were observed at channels 4000–8000, where the cable passes the Tosabae Ridge, and similar but smaller amplitude changes occurred above channel 8000. Temporal integration yields strain variations of > 20 microstrain, which cannot be explained by tidal deformation. Nevertheless, the spectrum of these changes shows clear peaks that correspond to diurnal and semidiurnal tidal periods (Figure S2). The very large change at around channel 6000 has a clear 12.4-h peak, which corresponds to the M2 constituent of the tidal signal. Therefore, these variations are obviously controlled by ocean tides.

Fig. 5
figure 5

Five-minute averages of the spatial and temporal changes in strain rate for the entire observation period. a Bathymetry profile along the cable route. b Space–time plot of the average strain rate (5-min interval) for the entire observation period. c Temporal changes in the average strain rate (5-min interval) at three selected channels. The traces have been shifted vertically for comparison. d Integration of the strain rate records shown in c. Note that the scale of 1 K includes a large uncertainty as explained in the text

The rapid change in strain rate migrated a long distance along the cable above channel 4000. Some large migrations of a negative anomaly propagated ~ 10 km (2000 channels) in ~ 6 h, and many small-scale migrations occurred repeatedly across many of the channels. The migration direction appears to be related to the geometry of the cable. Each of the negative anomalies propagated from the lower to higher elevations on both sides of Tosabae Ridge. Some large migrations, such as the events during the afternoon of January 28th and on the morning of January 29th, propagated over the summit of the ridge and continued further southward. The positive anomaly also migrated, but the migration speed was variable and behaved more randomly than those of the positive anomalies, which suggests that antisymmetric mechanisms influence the strain rate.

Although we have explained this change as a strain rate, the observed deformation is likely not an actual physical variation in the length of the fiber-optic cable, but rather a change in the refractive index of the fiber-optic cable due to variations in the ocean-bottom water temperature. The observed strain \(\in\) depends on the temperature \(T\), thermal expansion coefficient \(\alpha\), and change in the refractive index \(dn/dT\) as:

$$\frac{d \in }{{dT}} = n\alpha + \frac{dn}{{dT}}.$$

Typical value of \(n\), \({\alpha }\), and \(dn/dT\) are \(1.4\), \(4 \times 10^{ - 7}\), and \(10^{ - 5}\), respectively, which means that \(dn/dT\) is the dominant term. \(d \in /dT\) has been measured \(\sim 10^{ - 5}\) in laboratory experiments (Hocker 1979; Zumberge et al. 2018). We note that the real coefficient may be small in the ocean-bottom environment, as estimated by Zumberge et al. (2018), such that the degree of the temperature change described below may up to 3–5 times larger. Nevertheless, we assume \(dn/dT = 10^{ - 5}\) in our calculations, which converts the apparent strain rate changes of 1 and 10 microstrain to 0.1 and 1 K temperature changes for the basin (channels 1000–4000) and ridge (channels 4000–8000) sections, respectively.

Therefore, we can regard our DAS record as distributed temperature sensing (DTS) record, though the absolute values are unknown and the relative values are only approximately measured. DTS often uses Brillouin and Raman scattering and has been applied in various fields (e.g.,Hartog 2017; Miah and Potter 2017; Lu et al. 2019). Sinnett et al. (2020) demonstrated the usefulness of DTS in oceanographic researches, by measuring temperature fields along several ocean bottom cables at shallow (< 100 m) depth. They detected tidally modulated spatio-temporal temperature variations up to about 6 degrees, which migrate along the fiber. Rapid increase and relatively slow and more random decrease of temperature are characteristics commonly observed in their DTS record and our DAS record.

Our record suggests temporal temperature variation up to several degrees, given the abovementioned conversion is correct. Temperature variations of ~ 0.2 K were measured by the temperature sensors of the ocean-bottom pressure gauges along the Muroto cable (Inazu and Hino 2011), which is comparable with our observations in the basin section. However, we are unaware of any direct measurements of ocean-bottom temperature changes that are as large as several degrees. For example, Petruncio et al. (1998) reported large diurnal temperature change at a depth of 370 km in Monterey Canyon, USA, but the amplitude is usually less than 1 K. Conductivity-temperature-depth measurements have suggested large standard deviations in the water temperature of up to several degrees at several locations near our cable in the Nankai subduction zone (Watanabe and Uchida 2016). The variations at < 800 m depth are particularly large and likely due to the effect of the Kuroshio Current; our observations are likely related to such large variations. Nevertheless, the large uncertainty remains in our amplitude conversion, which may be modified in future calibration studies.

These temperature variations must reflect the movement of ocean water, with the rapid changes also suggesting movement at a sharp boundary, or interface, between two water masses with different temperatures. Internal surface waves are excited by tidal current and propagate on the interface much faster than the current itself. Complicated shape of the wave front especially during the increase of temperature suggests the breakage of the wave fronts and local water mixing. Although tidal modulation likely has a major influence on this movement, the effect of the Kuroshio Current, a strong surface ocean current that flows from west to east at > 1 m/s across the study region, may not be negligible. The overall orientation of the cable is at a large angle relative to the mean direction of the Kuroshio Current, such that the influence of this current on the DAS observations may be reduced. Nevertheless, the current direction near the ocean bottom might be different due to the complex bathymetry of the region, including Muroto Basin and Tosabae Ridge, which forms a narrow saddleback structure along the cable. Therefore, the ocean-bottom water mass likely moves as a result of complex interactions between relative steady ocean currents and quasi-periodic oceanic tides that are bounded by an irregular bathymetry. It should be noted that the migration velocity of 10 km/6 h (~ 0.5 m/s) is only the apparent velocity along the cable orientation, such that the speed of two- or three-dimensional boundaries can be much slower and yield a similar apparent migration velocity. Nevertheless, this rapid, long-distance migration suggests the existence of relatively large-scale movement near the deep (~ 1000 m) ocean bottom at a velocity comparable with the Kuroshio Current at the surface. This observation may prove useful in advancing our understanding of the mechanisms that drive the tidal mixing of deep ocean water.

Discussion and conclusion

We demonstrate that DAS can observe a suite of geophysical signals via DAS observations from an unused seafloor fiber-optic cable offshore Cape Muroto over a 5-day period. Our original interest was broadband seismic signals since the measured quantity is supposed to be the axial strain of the cable. We show that the performance of high-frequency (> 1 Hz) seismic observations is almost comparable with those of ocean-bottom permanent seismic stations, as long as the cable is well-coupled with the ocean bottom. Although only one component of the strain is measurable, high-density DAS observations at many channels could prove useful for quantifying different aspects of the local seismic structure and seismic source parameters.

The Nankai subduction zone is well known for various slow earthquake activity. Although our temporary observations did not detect any tectonic tremor or low-frequency earthquakes, which generally yields signals above 1 Hz, the noise level suggests that the cable has the ability to detect these phenomena if they occur within a few tens of kilometers of the cable. However, in the frequency range for very low frequency earthquakes, i.e., 0.02–0.05 Hz, the DAS observations are noisier than the short-period sensors, which makes such event detection impossible unless they occur very close to the cable.

The major noise source for the DAS measurements is the noise due to temperature-dependent changes in the refractive index; this effect is particularly large at very low frequencies. The amplitude of the temperature-dependent noise is as much as 10 microstrain for quasi-static (> 5 min) observations, which is about two orders of magnitude larger than expected for earth and ocean tide signals. However, this temperature-dependent noise can be used as a signal to monitor the temperature field, i.e., as DTS observation. Although DAS observation lacks the information of absolute temperature, it is relatively easy to detect systematic variations in space and time, as observed in our experiment. This strain-to-temperature conversion is only possible via continuous spatiotemporal monitoring, such as DAS observations. We can identify several migration sequences of rapid temperature changes that are up to several degrees, some of which propagated more than 10 km over several hours. We believe that this is the first report of such ocean-bottom temperature migration events. These events likely arise from interactions between the periodic tides, complex bathymetry, and strong Kuroshio Current, and may contribute to deep-water mixing processes near the ocean bottom.