Research article 27 Aug 2019
Research article  27 Aug 2019
Visual analytics of aftershock point cloud data in complex fault systems
 ^{1}Guangdong Key Laboratory of Urban Informatics, School of Architecture & Urban Planning, Shenzhen University, Shenzhen, 518000, China
 ^{2}Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, 518000, China
 ^{3}Key Laboratory of Urban Land Resources Monitoring and Simulation, Ministry of Land and Resources, Shenzhen, 518000, China
 ^{4}Shenzhen Urban Public Safety and Technology Institute, Shenzhen, 518000, China
 ^{5}Department of Earth and Space Sciences, Academy for Advanced Interdisciplinary Studies, Southern University of Science and Technology, Shenzhen, 518000, China
 ^{1}Guangdong Key Laboratory of Urban Informatics, School of Architecture & Urban Planning, Shenzhen University, Shenzhen, 518000, China
 ^{2}Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, 518000, China
 ^{3}Key Laboratory of Urban Land Resources Monitoring and Simulation, Ministry of Land and Resources, Shenzhen, 518000, China
 ^{4}Shenzhen Urban Public Safety and Technology Institute, Shenzhen, 518000, China
 ^{5}Department of Earth and Space Sciences, Academy for Advanced Interdisciplinary Studies, Southern University of Science and Technology, Shenzhen, 518000, China
Correspondence: Jincheng Jiang (jiangjincheng0305@126.com)
Hide author detailsCorrespondence: Jincheng Jiang (jiangjincheng0305@126.com)
Aftershock point cloud data provide direct evidence for the characteristics of underground faults. However, there has been a dearth of studies using stateoftheart visual analytics methods to explore the data. In this paper, we present a novel interactive visual analysis approach for visualizing the aftershock point cloud. Our method employs a variety of interactive operations, rapid visual computing functions, flexible display modes, and various filtering approaches to present and explore the desired information for the fault geometry and aftershock dynamics. The case study conducted for the 2016 Central Italy earthquake sequence shows that the proposed approach can facilitate the discovery of the geometry of the four main fault segments and three secondary fault segments. It can also clearly reveal the spatiotemporal evolution of the aftershocks, helping to find the fluiddriven mechanism of this sequence. An opensource prototype system based on the approach is also developed and is freely available.
 Article
(3600 KB) 
Supplement
(630 KB)  BibTeX
 EndNote
Fault complexity of earthquakes is universal in all types of faults (strikeslip, normal, and reverse). Usually, the faults are divided into a number of subparallel segments along the strike by geometrical discontinuities (Manighetti et al., 2015). In some cases, the depth segmentation of a fault is also observed (Elliott et al., 2011) where the cascading earthquake occurs successively along the dip. The geometrical discontinuities in the complex fault structure can sometimes impede and stop the rupture propagation and limit the earthquake magnitude. In some cases these act as a highpermeability conduit to allow fluid migration and control the spatiotemporal sequence of the subsequent earthquakes (Walters et al., 2018). The elasticrebound theory based on the earthquake cycle concept is commonly applied for longterm earthquake prediction. However, the complex fault structure complicates the seismic hazard assessment because theoretical models based on simplistic assumptions (e.g., single dynamic rupture) were found to be inapplicable to natural faults (Barbot et al., 2012). The time intervals for cascading fault rupture due to fault segmentation range from seconds to years on a casebycase basis, while the control factors of the faults are still not clearly understood. Therefore, it is highly important to understand the complex fault geometry, and such understanding may provide useful insight into fault mechanics and help to constrain the existing theoretical models.
Due to the past lack of ability to conduct highresolution observations, it is always difficult to define the geometry of the active fault system. Recent developments of both geodetic and seismic observation techniques make it possible to reveal the fault geometry at the surface and depth. Advanced synthetic aperture radar (SAR) geodesy provides nearfield deformation and can be used to directly map the fault surface and constrain the deep fault through deformation inversion. The construction of an increasing number of seismic stations allows for the identification of the locations of a large number of aftershocks following the main shocks. The point cloud of the located aftershocks contains the information regarding both the threedimensional coordinates and precise time of occurrence that can directly reveal the fault geometry and temporal evolution of the earthquake sequence. Current interpretation methods for the aftershock cloud are mostly casespecific rather than general, and to date, no standard processing procedure for the visualization of these data has been available. While many studies plot the crosssection profile of the aftershock point cloud to reveal the dip of the fault plane (Irmak et al., 2012; Hengesh and Whitney, 2016; Li et al., 2010; Elliott et al., 2013), these simple visualization methods usually have the following limitations:

Lack of interactive operations: the visualization methods provide a static display without an interactive interface to assist the users in effectively performing analyses through various human–computer interactions.

Lack of visual analysis: the current methods simply present these data without computing the feature (e.g., plane fitting or local outlier factor value calculation) to enhance the structure.

Lack of noise filtering: while the aftershock locations normally contain many uncertainties, most studies do not consider either the identification or removal of these noisy points.
As was shown in many recent studies, aftershock data can reveal much more information if these data are presented well, and this is particularly true for a complex fault system (Walters et al., 2018). Therefore, it is necessary to develop an interactive visual system to effectively mine information from the aftershock data.
Visual analytics is a research direction of the field of data visualization (Wong and Thomas, 2004), which integrates new computational and theorybased tools with innovative interactive techniques and visual representations to enable human–information discourse. How to efficiently transform large earth observation data into knowledge is a tremendous challenge in the era of big data (Xia et al., 2018). Visual analytics is highly advantageous and has good potential in complex data analytics for intuitive information representation and dynamic data exploration (Ward et al., 2015). The importance of scientific visualization has been recognized for a long time (McCormick, 1988). Specifically, in geophysics, the visualization of seismic refraction data has already been widely applied to explore the distribution of petroleum and gas (Liu et al., 2018; Raposo et al., 2016). However, to date, studies using stateoftheart visual analytics methods to explore other kinds of seismic data (e.g., aftershock point cloud) to gain information about earthquakes have been lacking.
In this study, we propose a novel interactive visualization method for exploring the 3D aftershock point cloud data. This method can help researchers to better understand the complex fault geometry in three ways:

A set of interactive operations is provided to stimulate creative analysts and exploit the researchers' background knowledge.

Onthefly computation of the local outlier factor (LOF), 2D projection, and plane fitting are supported in the visual analysis so that the hidden features of the fault structure can be discovered dynamically.

Noisy or irrelevant aftershock points can be filtered to reduce the distractions from invalid information.
Figure 1 shows the flowchart of the proposed approach. Initially, we imported the aftershock point data from an earthquake catalog or the literature into our system. At the same time, the LOF values for all of the points are computed automatically. Then, the points are presented in a 3D view window, allowing free rotation, panning, and zooming to observe the point structure. The earthquake fault, with data either from an existing source or from inferred information, can be imported into the 3D view as well. After defining the parameters for the projection (i.e., strike and dip of the plane and azimuth of the projection direction), the 2D points are rapidly computed and presented. Then, some interactions, i.e., coloring and filtering, can be performed on the 3D point clouds and 2D plane projection to enhance the structure. From the 2D plane projection, two interactions are further provided for exploring the fault geometry and aftershock dynamics. The first is the selection and plane fitting of a point cluster on a potential plane. After interactively selecting the clustered points, a plane will be automatically fitted by a robust estimator to reveal the potential fault structure. The second is the observation of aftershock migration. Animation of the aftershocks is provided so that a preliminary recognition of the migration direction can be obtained. Then, the user can define a propagation direction to further explore the aftershock evolution. After assigning the direction, the distance traveled by the aftershocks along the specified direction as a function of time will be plotted, which is helpful in discovering the driving factor of aftershock migration. We note that when the preliminary exploration of the fault geometry or earthquake dynamics is finished, the users can return to the previous interactions to refine their interpretation of the data.
2.1 Visual design
As shown in Fig. 2, the user interface includes four main components. The first component is the operation tools panel with all the buttons, sliders, text boxes, and labels. It includes most interactions in the visual analytics. The second component is the 3D view panel for presenting cloud points and fault planes. The third component is the 2D view panel for presenting the points in the 2D projection plane. The fourth component is a temporal panel that switches between 3D fault plane fitting plotting and propagation distance–time plotting, corresponding to the explorations of the fault geometry and aftershock dynamics, respectively.
2.2 Interactive operations
The easytouse interactions represent a significant difference between the proposed method and conventional aftershock interpretation methods. The traditional static displays prevent the analyst from directly and rapidly exploring the aftershock data. However, the easytouse interactive displays can lead to deeper insights into 3D aftershock cloud point data and the discovery of new fault structures or spatiotemporal evolution patterns. These interactive operations are described as follows.
Threedimensional view manipulation. The point cloud data can be displayed in a 3D view control (Fig. 2b). Some typical 3D interactions are available including zooming, panning, rotating, data tips, and data brushing (Fig. 2a), allowing interactive exploration and editing of the plotted data to improve the visual display.
Interactive projection. The crosssection plot is the main visualization approach for exploring the complex fault geometry (Fig. 2c). The faults are normally represented by linear features when the aftershocks are projected onto a plane perpendicular to the fault. We offer an interface to input the strike and dip of the projection plane for rapid plotting of a desired plane projection profile.
Dynamic filtering. Noisy or irrelevant data can significantly affect the interpretation of the aftershock point cloud. The filtering or denoising interaction is an important step in removing noisy or irrelevant aftershock points while retaining the fault geometry details. The filtering is based on the four features of the aftershocks, namely the depth, LOF, time, and magnitude (Fig. 2a). LOF is an index for finding anomalous data points that will be introduced below (Breunig et al., 2000). Filtering based on LOF is used to denoise the data.
Color switching. The feature for coloring the aftershock points can be quickly switched between the depth, LOF, time, and magnitude, allowing the users to effectively understand the data from different dimensions (Fig. 2a).
Fault plane fitting. To test the probability of an identified fault, the users first brush the points that form a linear feature in the crosssection plot (Fig. 2c). A fault plane will be automatically fitted to these points based on a robust estimator (Fig. 2d). This interaction combines expert knowledge and visual computing, allowing the user to quantitatively locate the potential fault structure.
Fault data fusing. Users can import the faults to the 3D view in order to assist in the visualization and analysis of the aftershock point cloud (Fig. 2b). These faults can be hypothetical seismogenic faults or background fault systems.
Animation. The aftershocks can be plotted gradually to form an animation in both the 3D view and the 2D crosssection plot (Fig. 2b and c). An animated visualization has unique advantages and aids the viewers in interpreting the evolution of the aftershocks, making comparisons between different points in time, chunking the aftershock points into smaller datasets, focusing attention, and anticipating and understanding the spatiotemporal changes in the aftershocks.
Propagation distance–time plotting. To further quantify the migration process of the aftershocks, users can define a propagation direction (Fig. 2c) and plot the distance along this direction as a function of time (Fig. 2d). The temporal trend may reveal the mechanism (e.g., fluid diffusion) underlying the aftershock migration.
2.3 LOF calculation
Because of the limitations due to the resolution of seismic data inversion, there is some uncertainty in the location of the aftershocks, so that some earthquake points show large deviations from the fault plane. Usually, the original aftershock data do not contain any information about this uncertainty, and we cannot eliminate these outliers. In this study, we proposed to use the LOF, which is commonly used in the identification of the outliers in the lidar point cloud (Wang et al., 2019) to detect the anomalous points that deviate distinctly from the fault plane. The LOF (Breunig et al., 2000) is based on the concept of local density, which is estimated by the reachability distance of the object to its k nearest neighbors. Points with a substantially lower density than their neighbors are considered outliers. By comparing the local reachability density of the object to those of its neighbors, its LOF value can be obtained. The LOF is calculated as follows.
Let kd(p) be the k distance of a point p, defined as the distance from point p to the kth nearest neighbor. The reachability distance between points p and o is then given by
where d(p,o) is the actual distance between p and o.
The local reachability density of an object p is then defined by
where N_{k}(p) is the kth nearest neighborhood, defined as the set of points within the distance of the k distance (p).
The LOF for point p is then obtained by
Figure 3 shows an example of the calculated LOF values, where an LOF value significantly larger than 1 indicates a sparse region.
2.4 Twodimensional projection
Projecting the 3D aftershock point cloud to a 2D plane can significantly enhance the fault geometry features. As mentioned above, we offer an interactive way to rapidly compute the projection when a plane and projection direction are assigned. To make the interactive operation userfriendly for geologists, we use geological terms, i.e., strike and dip (to define the orientation of the projection plane), and azimuth angle (to define the projection angle). We firstly calculate 3D geographic coordinates of the projected points on the projection plane for aftershocks according to the line–plane intersection equation. The 3D geographic coordinates are then transformed to a 2D local coordinate system by two element rotations: a z axis rotation by an angle of strike − 90^{∘} and an x axis rotation by an angle of − dip, as shown in Fig. 4a.
The rotation matrices (R_{z} and R_{x}) are given by
and
The projected coordinates (P) from the aftershock point (u) is given by
where the x and y coordinates of P correspond to the coordinates in the projection plane, and the z coordinate represents the distance to the plane.
2.5 Plane fitting
When the analyst identifies a potential linear feature from the crosssection plot, the plane fitting can rapidly evaluate the possible feature and provides the precise geometric parameters. In this study, we used an algorithm based on the singular value decomposition (SVD) to fit a plane to the given 3D aftershock points. The calculation of the parameters of the plane ($ax+by+cz=d$) based on the singular value decomposition method is carried out as follows.
Given a set of points $\left({x}_{i},{y}_{i},{z}_{i},\phantom{\rule{0.125em}{0ex}}i=\mathrm{1}\mathrm{\dots}n\right)$ and their center point ($\overline{x},\phantom{\rule{0.125em}{0ex}}\overline{y},\phantom{\rule{0.125em}{0ex}}\overline{z}$), we form a matrix A as given by
We apply an SVD to matrix A:
The parameter vector $[a,b,c]$ is just the smallest singular vector corresponding to the least singular value, and the last parameter d can be obtained by
To avoid disturbance from noisy points, we include an iterative step to identify outliers and exclude them from plane fitting. This step is carried out using the following algorithm:

Fit a plane using the available points.

Calculate the distance (D) of the points to the plane.

Calculate the root mean square error (RMSE) of all of the distances (D).

Find points with distances 3 times larger than RMSE and define them as outliers.

Remove these points, and then return to step 1.

Stop when there is no outlier found, and return the plane parameter.
Given the normal vector of the fitted plane ($a,b,c$), the strike and dip of the fault plane are given by (Fig. 4b)
and
where atan2(Y,X) is the four quadrant arctangent of the elements X and Y, $\mathit{\pi}<=\mathrm{atan}\mathrm{2}(Y,X)<=\mathit{\pi}$, and sign(.) is the signum function that returns 1 for a positive element and −1 for a negative element.
2.6 Prototype system
The proposed visualization approach has been implemented in a software package (Aftershock Visualization, AFV1.0; https://github.com/caigenszu, last access: 22 August 2019) and is freely accessible to the scientific community. This prototype system was developed in MATLAB with a graphical user interface (GUI), providing pointandclick control for all operations. An analyst with no programming background does not need to learn a new language or type commands to run the application. Meanwhile, the source code of the package is also provided so that advanced users can easily read the data and results throughout the visualization processing, and they can modify the package to include their custombuilt modules. This package requires an input file consisting of the aftershock point cloud table in plain text with each row representing a point and the columns representing UTMx, UTMy, depth, magnitude, time (days from reference time), LOF, and an event flag (whether or not the shock is the main shock) in that order. An alternative input file is the fault plane used to assist data interpretation. Each fault plane is represented by one row, and the columns represent strike, dip, center UTMx, center UTMy, length, upper depth, and bottom depth.
On 24 August 2016 (01:36 UTC, local time 03:36), a destructive earthquake (M_{w} 6.2) occurred in Central Italy (the Amatrice earthquake). The epicenter was located at 42.70^{∘} N, 13.23^{∘} E, between the towns of Norcia and Amatrice. Two months later, two major earthquakes were triggered to the northwest of the Amatrice earthquake on 26 October (19:18 UTC) and 30 October (06:40 UTC) (Papadopoulos et al., 2017), with individual moment magnitudes of 6.1 (Visso earthquake) and 6.6 (Norcia earthquake). The earthquake sequence occurred across several fault systems (Fig. 5), including the Mt. Vettore–Mt. Bove Fault system (VBFS) and the Mt. Gorzano Fault system (GFS). The VBFS is a SWdipping fault system composed of different splays and segments. The fault system scarps are exposed along the SW foothills of Mt. Vettore, Mt. Porche, and Mt. Bove, with a length of ∼27 km. The GFS is a SWdipping fault with a ∼26 km long fault scarp developing along the foothills of Mt. Gorzano. These two faults are segmented by existing tectonic structures inherited from preQuaternary compressional tectonics (Pizzi and Galadini, 2009). A ∼3 km thick layer in which small events and some large extensional aftershocks occur is found below the seismogenic fault, limiting the seismicity to the first 8 km of the crust (Chiaraluce et al., 2017).
Many recent studies reported in the literature have revealed a very complex fault geometry for this earthquake sequence. In addition to the main rangefront normal faulting structures that were adopted in many early studies (Liu et al., 2017; Tinti et al., 2016), recent studies find increasing evidence for some additional minor antithetic and synthetic faults (Chiaraluce et al., 2017; Scognamiglio et al., 2018; Walters et al., 2018; Cheloni et al., 2017). The activation of this secondary fault structure has important implications for evaluating the seismic hazard in this section of the Central Apennines. The complex fault geometry also plays an important role in channeling fluid diffusion and affects the spatial order of the cascading earthquakes. The aftershock distribution provides direct evidence for the fault geometry and underground dynamics. By interactively visualizing these data, we can obtain useful information regarding both the fault geometry and earthquake dynamics. Additionally, a relocated aftershock dataset has been provided in the work of Chiaraluce et al. (2017), making this earthquake sequence an ideal case for applying our proposed interactive visualization method. Below, we demonstrate the benefits of our method in identifying fault geometry and aftershock migration.
3.1 Fault geometry identification
According to previous studies, several fault segments are involved in this earthquake sequence. Some of these are the main fault segments located in the GFS–VBFS, while others are secondary fault segments including a NEdipping normal fault antithetic to the GFS–VBFS and a preexisting compressional structure that is likely related to a segment of the Olevano–Antrodoco–Sibillini Thrust. We will show how our interactive visualization method can facilitate the recognition of these structures.
Following the flowchart shown in Fig. 1, we first import the relocated aftershock data from the results obtained by Chiaraluce et al. (2017) to the 3D view. Some of the basic interactions provided by this package, i.e., zooming, panning, rotating, data tips, and data brushing, can be used to explore the data. LOF values were then calculated for all points to identify the aftershock outliers. We can interactively select to show some of the points according to their location, depth, magnitude, LOF, and time, and we can project them to a desired plane. After we find a linear feature in the projection plane, we select the points comprising the linear feature and automatically fit a plane to the feature (Fig. 6).
First, we identify the four main fault planes by projecting the cloud point to a plane with a 70^{∘} strike and a 90^{∘} dip. This plane is approximately perpendicular to the GFS–VBFS, and it should reveal the cross sections of the fault plane well. To better illuminate the linear features of the faults, we filter the aftershock points according to the fault location and the spanning time of the aftershocks on the fault. To denoise the point cloud, the aftershocks with the LOF value larger than 2 are filtered. We then select the points around the linear feature for each fault plane, and we automatically obtain the suggested geometry parameters of the fitted plane as shown in Table 1, finding that the obtained parameters are in overall agreement with the main seismogenic fault planes in previous studies.
Second, we explore the point cloud to find other secondary faults. Since there is no hint of the existence of these faults, we use the basic interaction tools in the 3D view and the filtering function in our package to find the possible fault geometry. An obvious cluster of aftershock points is found on an antithetic NEdipping plane in the Norcia area. We therefore also project the cloud point to a plane with a 70^{∘} strike and a 90^{∘} dip. After filtering the points by location and LOF value, we obtain a very strong linear fault signal as shown in Fig. 6. Fault fitting gives a strike of −38^{∘} and a dip of 59.8^{∘} for the Norcia antithetic fault. However, other fault structures revealed in the literature are not very obvious in the aftershock cloud point. We then project the point cloud to a simplified seismogenic fault ($\text{strike}=\mathrm{160}{}^{\circ}$, $\text{dip}=\mathrm{50}{}^{\circ}$) along a horizontal vector with an azimuth of 15^{∘}, using the projection settings from Walters et al. (2018). We filter the points by the magnitude, location, timing, and LOF, and then, two linear structures are observed as shown in Fig. 6. The structure on the right is the NWdipping Pian Piccole fault that generates a significant surface displacement signal in the interferometric synthetic aperture radar (InSAR) measurements. The structure on the left is an eastdipping fault that was only inferred in Walters et al. (2018). We pick the points around these linear structures for plane fitting and obtain a strike of −114^{∘} and a dip of 46.5^{∘} for the Pian Piccole fault and a strike of 7.8^{∘} and a dip of 85.3^{∘} for the inferred fault.
3.2 Migration of aftershocks
Besides the static analysis of aftershocks to find the fault geometry, the timing information of aftershocks has the potential to reveal earthquake dynamics related to fluid diffusion, afterslip, and/or pressure transients. Fluid diffusion is a main factor in triggering seismicity as observed in many tectonic regimes (Vidale and Shearer, 2006). It was also proposed to have been present in the nearby 2009 L'Aquila earthquake (Malagnini et al., 2012). The migration of the aftershocks could provide some evidence for such diffusive processes. Walters et al. (2018) found a significant fluiddriven aftershock migration between the Amatrice earthquake and Visso earthquake controlled by the fault intersection in this earthquake sequence.
To observe the spatiotemporal process, we project the aftershocks of the seismogenic fault to be along a horizontal vector with an azimuth 15^{∘} as described above. The points are colored by their rupture time. As shown in Fig. 7a, we can observe a trend by the color change from green to red in the northward direction, implying the clear directionality of the aftershock migration. We then use the animation operation to plot the aftershock points step by step as a function of time in both the 2D projection plane and the 3D view. The directional propagation of the aftershocks can be observed much more clearly.
Next, we can define a starting point and an assumed direction in the projection plane. After drawing the starting point and the propagation line, the distances of the aftershocks to the starting point along this direction can be calculated automatically. They are then plotted as a function of time to further reveal the time evolution of the aftershocks. Following the work of Walters et al. (2018), we defined a starting point in the peak slip of the Amatrice earthquake and a direction perpendicular to the Pian Piccolo fault. A clear temporal trend of the aftershocks was then observed in the plot (Fig. 7b). Therefore, the observed spatiotemporal pattern of aftershock migration can be interpreted by the seismologist to understand the earthquake mechanism and infer the possible controlling factor for the fault rupture. In this case, the diffusivelike aftershock temporal trend can be further interpreted using a pore fluid source model. Through the interactive visualization operations, we can observe that aftershock propagation is consistent with a diffusive process and can also observe the spatial concentration of the aftershocks along the fault intersection. This information can be used to derive the conclusion, as described in the work of Walters et al. (2018), that the intersecting structures act to channel the fluids and control the timing and order of the subsequent earthquakes.
We present a novel interactive approach and develop a prototype system to illustrate 3D aftershock point clouds that can help the geophysicist and seismologist to better understand the geometry of a complex fault system and the spatiotemporal pattern of the aftershocks. Moreover, it can be applied in an educational scenario to facilitate student learning. In this approach, we design a set of interactive operations to facilitate the exploration of the data. Fast computation of LOF, 2D projection, and plane fitting are implemented in addition to the visualization to reveal additional, hidden information, which is not obvious in the original data. Additionally, the point cloud can be visualized in many ways (e.g., 3D view, 2D plane projection, and animation), aiding viewers in interpreting the aftershock point cloud. A wide range of data filtering options (e.g., by magnitude, depth, location, time, and LOF) enable the users to avoid interference from invalid data in order to focus on the useful information.
It is widely accepted that visual analytics with interactive visual interfaces can amplify human cognitive capabilities, and it is suitable for solving some complex problems relying on closely coupled human and machine analysis. When applied to the aftershock catalog, visual analytics can therefore help to discover some new fault planes or weak signals of aftershock migrations, which might be hard to observe directly from the data. Specifically, the designed visual analytics procedure can assist the knowledge discovery from the aftershock catalog in several ways:

by accelerating the fault segment discovery processing through 3D view manipulation (e.g., zooming, panning, and rotating);

by enhancing the recognition of fault segment patterns through rapid visual computing functions (e.g., interactive projection, plane fitting, and fault data fusing);

by reducing the influence from lowquality or irrelevant aftershock points through various filtering methods (e.g., depth, LOF, time, and magnitude filtering);

by enabling the aftershock migration exploration by providing more cognitive resources (e.g., animation and propagation distance–time plotting).
A case study of the 2016 Central Italy earthquake sequence shows that the proposed approach can complete a series of complex visualization operations in a userfriendly interactive manner. We can conveniently infer the fault geometry for both the main faults and the secondary faults, and we clearly observe the spatiotemporal pattern of the fluiddriven aftershock migration. In future work, we would like to examine the capability of this tool for more aftershock point data applications, such as analyzing tremor migration and tracking aseismic slip. Also, we plan to include more kinds of fault information (e.g., InSAR, GPS, source model, and fault scarp) in the visualization in order to further facilitate the analysis of earthquake fault behavior.
The research data are included in the Supplement.
The supplement related to this article is available online at: https://doi.org/10.5194/se1013972019supplement.
This work was conducted by CW and JJ under the supervision of QL. The paper was prepared by CW and JJ with contributions from all coauthors.
The authors declare that they have no conflict of interest.
This research has been supported by the Shenzhen Scientific Research and Development Funding Program (grant nos. KQJSCX20180328093453763, JCYJ20180305125101282, JCYJ20170817104236221, and JCYJ20170412142144518), the Open Fund of the State Key Laboratory of Earthquake Dynamics (grant no. LED2016B03), the NSFC (grant nos. 61802265 and 41974006), Guangdong Provincial Natural Science Foundation (grant no. 2018A030310426), and the Open Fund of Key Laboratory of Urban Land Resources Monitoring and Simulation, Ministry of Land and Resource (grant no. KF201803004).
This paper was edited by Caroline Beghein and reviewed by two anonymous referees.
Barbot, S., Lapusta, N., and Avouac, J. P.: Under the hood of the earthquake machine: toward predictive modeling of the seismic cycle, Science, 336, 707–710, 2012.
Breunig, M. M., Kriegel, H., Ng, R. T., and Sander, J.: LOF: identifying densitybased local outliers, International Conference on Management of Data, 29, 93–104, 2000.
Cheloni, D., De Novellis, V., Albano, M., Antonioli, A., Anzidei, M., Atzori, S., Avallone, A., Bignami, C., Bonano, M., and Calcaterra, S.: Geodetic model of the 2016 Central Italy earthquake sequence inferred from InSAR and GPS data, Geophys. Res. Lett., 44, 6778–6787, 2017.
Chiaraluce, L., Di Stefano, R., Tinti, E., Scognamiglio, L., Michele, M., Casarotti, E., Cattaneo, M., De Gori, P., Chiarabba, C., and Monachesi, G.: The 2016 Central Italy seismic sequence: a first look at the mainshocks, aftershocks, and source models, Seismol. Res. Lett., 88, 757–771, 2017.
Elliott, J., Copley, A., Holley, R., Scharer, K., and Parsons, B.: The 2011 M_{w} 7.1 van (eastern Turkey) earthquake, J. Geophys. Res.Sol. Ea., 118, 1619–1637, 2013.
Elliott, J. R., Parsons, B., Jackson, J. A., Shan, X., Sloan, R. A., and Walker, R. T.: Depth segmentation of the seismogenic continental crust: The 2008 and 2009 Qaidam earthquakes, Geophys. Res. Lett., 38, L06305, https://doi.org/10.1029/2011GL046897, 2011.
Hengesh, J. V. and Whitney, B. B.: Transcurrent reactivation of Australia's western passive margin: An example of intraplate deformation from the central IndoAustralian plate, Tectonics, 35, 1066–1089, 2016.
Irmak, T. S., Doğan, B., and Karakaş, A.: Source mechanism of the 23 October, 2011, Van (Turkey) earthquake (M_{w}=7.1) and aftershocks with its tectonic implications, Earth Planets Space, 64, 991–1003, 2012.
Li, Y., Jia, D., Shaw, J. H., Hubbard, J., Lin, A., and Wang, M.: Strutural interpretation of the coseismic faults of the Wenchuan earthquake: 3D modeling of the Longmen Shan foldandthrust belt, J. Geophys. Res., 390, 275–286, https://doi.org/10.1029/2009JB006824, 2010.
Liu, C., Zheng, Y., Xie, Z., and Xiong, X.: Rupture features of the 2016 M_{w} 6.2 Norcia earthquake and its possible relationship with strong seismic hazards, Geophys. Res. Lett., 44, 1320–1328, 2017.
Liu, R., Chen, S., Ji, G., Zhao, B., Li, Q., and Su, M.: Interactive stratigraphic structure visualization for seismic data, J. Visual Lang. Comput., 48, 81–90, 2018.
Malagnini, L., Lucente, F. P., Gori, P. D., Akinci, A., and Munafo, I.: Control of pore fluid pressure diffusion on fault failure mode: Insights from the 2009 L'Aquila seismic sequence, J. Geophys. Res.Sol. Ea., 117, B05302, https://doi.org/10.1029/2011JB008911, 2012.
Manighetti, I., Caulet, C., Barros, L., Perrin, C., Cappa, F., and Gaudemer, Y.: Generic alongstrike segmentation of Afar normal faults, East Africa: Implications on fault growth and stress heterogeneity on seismogenic fault planes, Geochem. Geophys. Geosyst., 16, 443–467, 2015.
McCormick, B. H.: Visualization in scientific computing, ACM SIGBIO Newsletter, 10, 15–21, 1988.
Papadopoulos, G., Ganas, A., Agalos, A., Papageorgiou, A., Triantafyllou, I., Kontoes, C., Papoutsis, I., and Diakogianni, G.: Earthquake triggering inferred from rupture histories, DInSAR ground deformation and stresstransfer modelling: the case of Central Italy during August 2016–January 2017, J. Pure Appl. Geophys., 174, 3689–3711, 2017.
Pizzi, A. and Galadini, F.: Preexisting crossstructures and active fault segmentation in the northerncentral Apennines (Italy), Tectonophysics, 476, 304–319, https://doi.org/10.1016/j.tecto.2009.03.018, 2009.
Raposo, A., Reis, L., Rodrigues, V., Pinheiro, R., Elias, P., and Cherullo, R.: A system for integrated visualization in oil exploration and production, IEEE Comput. Graph., 36, 10–16, 2016.
Scognamiglio, L., Tinti, E., Casarotti, E., Pucci, S., Villani, F., Cocco, M., Magnoni, F., Michelini, A., and Dreger, D.: Complex fault geometry and rupture dynamics of the M_{w} 6.5, 2016, October 30th central Italy earthquake, J. Geophys. Res.Sol. Ea., 123, 2943–2964, 2018.
Tinti, E., Scognamiglio, L., Michelini, A., and Cocco, M.: Slip heterogeneity and directivity of the ML 6.0, 2016, Amatrice earthquake estimated with rapid finitefault inversion, Geophys. Res. Lett., 43, 10745–10752, 2016.
Vidale, J. E. and Shearer, P. M.: A survey of 71 earthquake bursts across southern California: Exploring the role of pore fluid pressure fluctuations and aseismic slip as drivers, J. Geophys. Res.Sol. Ea., 111, B05312, https://doi.org/10.1029/2005JB004034, 2006.
Walters, R. J., Gregory, L. C., Wedmore, L. N. J., Craig, T. J., McCaffrey, K., Wilkinson, M., Chen, J., Li, Z., Elliott, J. R., Goodall, H., Iezzi, F., Livio, F., Michetti, A. M., Roberts, G., and Vittori, E.: Dual control of fault intersections on stopstart rupture in the 2016 Central Italy seismic sequence, Earth Planet. Sc. Lett., 500, 1–14, https://doi.org/10.1016/j.epsl.2018.07.043, 2018.
Wang, C., Shu, Q., Wang, X., Guo, B., Liu, P., and Li, Q.: A random forest classifier based on pixel comparison features for urban LiDAR data, ISPRS Journal of Photogrammetry Remote Sensing, 148, 75–86, 2019.
Ward, M. O., Grinstein, G., and Keim, D.: Interactive data visualization: foundations, techniques, and applications, AK Peters/CRC Press, New York, 2015.
Wong, P. C. and Thomas, J.: Visual analytics, IEEE Comput. Graph., 24, 20–21, 2004.
Xia, J., Yang, C., and Li, Q.: Building a spatiotemporal index for Earth Observation Big Data, Int. J. Appl. Earth Obs., 73, 245–252, 2018.