Abstract
The recent dramatic increase in millimeter- to centimeter- resolution topographic datasets obtained via multi-view photogrammetry raises the possibility of mapping detailed offset geomorphology and constraining the spatial characteristics of active faults. Here, for the first time, we applied this new method to acquire high-resolution imagery and generate topographic data along the Altyn Tagh fault, which is located in a remote high elevation area and shows preserved ancient earthquake surface ruptures. A digital elevation model (DEM) with a resolution of 0.065âm and an orthophoto with a resolution of 0.016âm were generated from these images. We identified piercing markers and reconstructed offsets based on both the orthoimage and the topography. The high-resolution UAV data were used to accurately measure the recent seismic offset. We obtained the recent offset of 7â±â1âm. Combined with the high resolution satellite image, we measured cumulative offsets of 15â±â2âm, 20â±â2âm, 30â±â2âm, which may be due to multiple paleo-earthquakes. Therefore, UAV mapping can provide fine-scale data for the assessment of the seismic hazards.
Similar content being viewed by others
Introduction
High-resolution three dimensional (3D) data are essential for constraining the spatial characteristics of active faults as they provide a fundamental dataset for mapping the complexity of faulting and the details of the surficial features offset by the faults. The spatial characteristics of a fault can be addressed using datasets with different resolutions1, 2. Sub-meter-resolution satellite images (QuickBird, GeoEye, Pleiades, etc.) and airborne Lidar data are commonly used to map earthquake surface ruptures over hundreds of kilometers. However, the spatial and vertical accuracies of these datasets are not comparable to more precise measurement methods, such as surveys using total station, real-time kinematic Global Positioning System (GPS) or terrestrial Lidar data3. Still, the acquisition of such data requires long survey sessions and significant investments of time in data processing, and such acquisitions are not always possible due to logistics.
With the development of remote sensing devices and new mapping techniques, the Unmanned Aerial Vehicles (UAV) have been used to enhance the mapping efficiency and to obtain images and 3D data at an unprecedented resolution4, 5. UAV mapping, like other photogrammetry-based techniques (i.e., terrestrial digital photogrammetry), captures 3D information of features from two or more photographs of the same object, taken from different angles6,7,8. High-quality photos can be quickly and easily collected in the field by attaching modern digital cameras to a small flying platform. The platform is generally inexpensive and can be either a fixed-wing type or a multi-rotor UAV. Multiple investigations have used UAVs to map folds and fractures in low elevation regions9,10,11,12. However, using UAV techniques to map faults in high elevation areas is rare, mainly due to poor accessibility to the field sites and the difficulties associated with proper logistical support. In addition, operating UAVs at high elevated sites can be challenging due to issues such as lower air pressure13,14,15 and frequent signal disturbances, as the level of ionization increases with elevation.
In this study, we used a multi-rotor UAV that can fly at exceptionally low altitudes with high quality cameras. We aim to use the UAV to capture the evidence of paleo-seismic deformation along an active fault, the Altyn Tagh fault (Fig. 1). This fault is one of the largest strike-slip faults on Earth. We selected the Bangewa segment, where degradation of the landscape is minimal due to the dry climatic and where earthquake offsets are consequently well preserved. The well-preserved co-seismic and cumulative offsets, which characterize the spatial displacements along the fault, provide important information for assessing the regional seismic hazards. However, with the currently available satellite images (Fig. 2a,b and c), it is difficult to interpret the offset created by the most recent earthquake (MRE) or small cumulative offsets. Therefore, targeting this problem, our investigation aims to obtain high resolution data to provide accurate measurements of these recent offsets. Additionally, we validated the cumulative offsets based on SPOT satellite images obtained by Mériaux et al.16. In this study, we first introduce the UAV methodology used to obtain the images and topographic data. Then, we reconstruct the offsets based on the new data. Finally, we evaluate the results by comparing them with satellite images (Figs 2 and 3) and terrestrial laser scanning (TSL) results (Fig. 4). The methods used in this study can be applied to any active fault, and the results of the study reveal the best methods for efficiently mapping the details of fault ruptures in high resolution.
Overview of the northeast section of the Altyn Tagh fault (see red rectangle in the inset for location). NATF: North Altyn-Tagh fault, SATF: South Altyn-Tagh fault. The study site is located on the NATF, at the Bangewa site (BGW). Active faults were mapped onto a background image from Google (Image Landsat/Copernicus) based on the Adobe Illustrator mapublisher (http://www.avenza.com/mapublisher). The inset map background is based on a global digital elevation model (GTOPO30 DEM with a horizontal grid spacing of 30 arc seconds, courtesy of the U. S. Geological Survey).
(a) Pleiades satellite image (0.7âm resolution) of fans along the ATF near Bangewa. (b) Geomorphic interpretation of faulting across the fan apexes and river channels along the fault (revised after Meriaux et al.16). White boxes indicate the locations of Fig. 6. Black boxes indicate the location of Fig. 9c,d and eâh.
The Bangewa site (BGW) is located along the North Altyn Tagh fault (NATF) (Fig. 1). This is a fault where large earthquakes have not been recorded in the historical earthquake catalog or by modern instruments. This region lies in the rain-shadow of the Himalayan range, and its inhabitants are predominantly nomadic. Here, the effects of climatic or anthropogenic degradation on the landscape are minimal. As a result, it represents an extraordinary âmuseumâ of extremely well preserved historic and ancient strike-slip ruptures17. The fault system comprises two principal parallel strands bounding the Altun Shan, which represents a large push-up structure. The South Altyn Tagh fault (SATF) follows the western rim of the Qaidam basin. The North Altyn Tagh fault (NATF) cuts through the Altun Shan with a N84°E strike (Fig. 1). Systematic sinistral offsets of alluvial fans, fluvial terraces, risers and channels are observed all along the fault, whose principal trace is remarkably sharp and linear, indicating sustained active tectonic movement (Fig. 3). We conducted this UAV-based investigation to further constrain the spatial characteristics of the ATF, to validate the offset obtained by the previous study, i.e., approximately 20âm16, and to measure the recent small seismic offset.
Methodology
To obtain high resolution topographic data at the Bangewa site, we followed the typical procedures5 of using a UAV. Here, we first introduce the UAV platform and the workflow used to acquire high resolution topographic data. Then, based on the acquired data, we accurately reconstruct the offsets created by the recent earthquakes.
Platform and sensor
The self-assembled UAV platform used in this study is equipped with a Sony NEX 5âT digital camera with a Sony Alpha Wide-Angle E-Mount lens, which has a 16 mm focal length and a 16 megapixel resolution. The whole platform weighs approximately 3âkg, including the 3300 mAh battery. The total cost of the equipment is less than $4,000 USD. A single flight conducted at approximately 60âm above ground level can produce approximately 200â250 images using a standard operating configuration. Detailed specifications are listed in Table 1. The Hexa-copter has an onboard navigation system based on a navigation-grade GPS receiver (NAZA-M2). Imaging is triggered by a camera auto-trigger at a rate of approximately 1 photo/s, which provides ample image overlap.
Workflow for high-precision data collection
Multi-view photogrammetry works on the basis that two or more overlapping photographs can be used to 1) calculate the unique three-dimensional (3D) location of a set of homologous image points in both photographs, which in turn can be used to compute the position and orientation of the camera, and 2) obtain a mosaic orthoimage. Photos are typically collected with at least 50â60% overlap under near-parallel viewing conditions18, 19. In this study, we used the Structure from Motion (SfM) approach, an alternative to the classical digital photogrammetry approach based on the structured acquisition of images20. The process starts by considering the scaling and georeferencing requirements for the target surface2 (Fig. 5). Although the target surface can be fully reconstructed in 3D without any scale or position information, in order to extract oriented and scaled data, additional control data or direct georeferencing must be used. At least three ground control points (GCPs) are required to perform the scaling or referencing. However, to ensure accuracy in scale and full georeferencing, more GCPs are needed. These GCPs should be distributed widely across the target area21. To ensure visibility, 12 bright red plastic boards were used for the GCPs in this study, and the center point of each board was measured. We made two rounds of GPS observations (moving from GCP1 to GCP12 then repeating the process) using the same set of handheld GPS instruments. The measurements were completed within 40âminutes. We calibrated the two rounds of observed data for the final geo-referencing coordinates with differential correction to within 1âm accuracy. The GCPs were collected using a Trimble handheld GPS device which mounted with high-sensitivity GPS receiver. To minimize the horizontal error, each measurement was taken with at least 8 satellites. To quantify the geometric accuracy, the root mean square error (RMSE) was used22. The RMSE is the comparison of real world (ground truth) information to estimated (image-derived) measurements. The RMSE values of the GCPs were 0.38âm (X) and 0.26âm (Y) in the horizontal directions and 0.07âm (Z) in the vertical direction. The mean re-projection error was evaluated as approximately 1.1 pixel.
At the study site, wind speeds generally increase over the course of the day, and the UAV flights were therefore performed in the morning to maximize flight stability and image quality. The UAV was launched from a wooden board (1.2âmâÃâ2âm) and was manually flown in a rectangular path (Fig. 5). All of the nine flights performed were successful and yielded imagery usable for both the full coverage of the fault segment and the DEM generation. The images collected in the field were processed with the commercial software, Agisoft Photoscan Professional, following the detailed procedures introduced by Lucieer et al.23 (Fig. 5). The initial bundle adjustment formed a dense 3D point cloud of the terrain with 37 million points. A more detailed 3D model with 2 million facets was also reconstructed based on multi-view stereopsis. We obtained a high-resolution grid-based DEM with a vertical resolution of 0.065âm, and an orthophoto with a spatial resolution of 0.0016âm. The area of the mapping site is ~700âm long and ~180âm wide. The specific algorithms implemented in Photoscan and the commonly used parameters are described in Verhoeven24.
Reconstruction of the seismic offsets
A seismic offset is the displacement created by an earthquake12, 17, 25, 26. The main goal of the reconstruction of seismic offsets is to link the lateral displacements with individual paleoearthquakes. Researchers use imagery to observe and measure the earthquake produced surface offsets, which include the offset caused by the most recent earthquake (MRE) and the cumulative offsets27, 28. To document the distribution of offsets along the rupture trace of the fault, we reconstructed them based on high resolution orthoimages combined with the topographic profiles derived from the DEM to assist in the interpretation and measurement of the offset markers, which are mostly terrace risers and river channels. The offset of these features can only have accumulated since the abandonment of the terrace tread, assuming that the feature is contemporaneous with the terrace.
The fault trace is mapped and the preferred offset values are determined by realigning each marker (terrace riser, gully, and abandoned channel) that has been displaced by fault motion based on both image interpretation and the topographic data (Fig. 6a,b and d). Targeting the clearest offset geomorphic features, we determined the seismic offsets (Fig. 7) by linearly projecting the correlative markers into the fault trace by eye and measuring the range of horizontal offsets29. First, we interpreted/identified the geomorphic markers that have been displaced by the motion of the fault. Then, we cut the image along the fault lines and realigned the offset markers. Finally, we measured the length of offset caused by motion on the fault. To evaluate the uncertainty of the recent small offset, we carefully checked the image for offsets from 2âm to 9âm (Fig. 8). To evaluate the uncertainty of the cumulative offset, we referred to the method of Gold et al.29. We considered two projections for each marker: one based on the trend of the terrace riser/center of the river channel nearest the fault line and a second based on the average trend for the entire length of the terrace riser/center of the river. The error bars for these measurements were determined based on the maximum value between the best fit measurements and the maximum/minimum measurements26, 29.
(a) Orthoimage (0.016âm resolution) of the study area. (b) High resolution bare earth digital elevation model (0.065âm resolution) from the UAV mapping. (c) Digital elevation model (0.5âm resolution) calculated from a stero pair of Pleiades satellite images. (d) The fault and the six streams (R1âR6) used in the reconstruction of the seismic offsets.
Results
Recent small offset obtained from the UAV data
Along the fault at the Bangewa site, geomorphic piercing lines defined by stream channels and terrace risers were identified to measure the offsets. The degree of preservation of offset geomorphic markers was not constant from river channel (R) R1 to R6, and offset markers were not preserved in the most active out-wash. However, by using the high resolution data, especially the contour lines derived from the DEM, we were able to identify micro-scale geomorphic markers, which helped to add important constraints on the small offsets (Fig. 7). Restoring the terrace risers on the left bank of R2, the small offset of ~7.5âm was measured (Fig. 7). At R1, the terrace riser on the left bank is offset by ~7âm (Fig. 7a,b). We carefully checked the entire image for possible offsets of 2âm to 9âm. We found maximum markers for offset of ~7âm (the most clear matching markers are also shown in Fig. 7). We determined that the uncertainty was ±1âm, as we found almost same number of markers for offsets of 6âm and 8âm (Fig. 8).
Cumulative offset
We further restore the highest terrace risers of R2 and obtained a cumulative offset of 30â±â2âm (Fig. 7c). Additionally, we also restored the major active channels and the abandoned channels east of the study site based on the Pleiades satellite image (Fig. 9). Large offsets of 15â±â2âm, 20â±â2âm and 30â±â2âm were obtained. However, the observation of small offsets was difficult due to the relatively low resolution of the satellite image compare to that of the UAV data.
Reconstructions of cumulative offsets based on the Pleiades satellite image. Locations of a and b were indicated as white box in e. See locations of c, d and eâh in Fig. 3. Cumulative offset of 15â±â2âm, 21â±â2âm and 30â±â2âm were identified.
Discussion
Compared with other remote platforms such as satellites, UAVs can operate at low altitudes30. The UAV used in this investigation is capable of flying at low altitudes of approximately 60âm above the ground. Therefore, even on cloudy days or in the early morning when the wind is usually weak and light is not sufficient, UAVs can provide images with unparalleled spatial resolution compared to satellite images (Fig. 2). The high spatial resolution of the data enables researchers to study active tectonics to make more detailed interpretations and measurements of the faulted landscapes12.
The resolution and accuracy of the dataset are ideal for identifying micro-scale geomorphic features for offset reconstruction. The uncertainties in the horizontal offset measurements from high resolution satellite imagery restorations are usually greater than 10%31. Previous studies using satellite images obtained a smallest cumulative offset value of approximately 20âm16. We re-measured this offset and further constrained the uncertainty to be 20â±â2âm, which is consistent with the previous measurements.
The accuracy of the measurements derived from the imagery is strongly correlated with the acquisition height above the ground32, 33. The horizontal RMSE is less than 40âcm, and the vertical RMSE is at the millimeter scale. The errors are small regarding the measurements of the offsets of large earthquakes, which are commonly on the order of a few meters. The main error in the offset is mainly related to how well an offset geomorphic marker is preserved or its shape across the fault trace. The error tends to be larger for eroded markers than for straight markers, such as a linear incised stream or a riser, because the reconstructions of their initial geometry usually involves less ambiguity17. However, uncertainties can be reduced by using images in combination with the DEM data to identify the markers, as the relative chronology of the incised features can also be considered using elevation. In this study, we use both imagery and topographic data to reconstruct the offsets; therefore, the accuracies of the UAV-obtained orthoimage and the DEM are sufficiently high for the purpose of seismic offset measurements.
We selected 6 streams along the fault in the Bangewa site to reconstruct the local offsets. The offsets reported here are all from natural terrace risers and fluvial channels, and no man-made features were measured. Uncertainties related to the measurements are mostly due to the degree of preservation of the piercing line across the fault zone. Offsets are not regularly distributed along the fault due to the uneven preservation of the geomorphic markers. All 6 selected streams have been incised by active channels, which are progressively offset by the movement of the fault. Targeting the best-preserved markers and unambiguous offsets ensures that the measurement uncertainties are usually small17, 26. As offsets of stream channels can be interpreted in terms of individual paleoearthquakes only if the stream re-establishes its course across the fault between each earthquake25, we only use small abandoned streams that are unambiguously related to earthquakes.
The reconstruction uncertainties depend on the marker size and sharpness, the orientation relative to the fault trace and the width of the fault zone. Previous works have suggested that field measurements tend to provide minimum offset values34,35,36. In certain cases, researchers using field observations may mis-interpret the cumulative offset as a single event offset, thereby overestimating the latter37. Here, the relative simplicity and freshness of the rupture trace and the high quality of the images make the measurements more reliable. The offset reconstructions for this site were based on a combination of river channels and terrace risers. For reconstructions of successive fluvial risers, criteria associated with the use of either upper or lower terraces should be considered according to the actual geomorphic settings38. These criteria should be considered to determine the slip-rates when the ages of the different terraces (both lower and upper terraces) are available. Unfortunately, this study site contains few materials that can be dated with high temporal resolution.
To evaluate the DEM derived from the UAV, a survey using a Trimble VX Spatial Station was performed for the eastern part of the study site (Fig. 4)39. The Trimble VX Spatial Station was remotely operated by a TSCII Controller, scanning up to 15 data points persecond. Using infrared direct (ID) reflex technology, this instrument is capable of advanced optical surveying, metric imaging, and 3D scanning, all of which can be accessed in the field through the controller using the real-time video images captured by the station40. Since the Trimble VX Spatial Station has an estimated horizontal beam divergence of 4âcm per 100âm and a vertical divergence of 8âcm per 100âm, the scan resolution was typically set to 5âcm horizontal to 2.5âcm vertical, which is more than 10 times lower than the resolution of the UAV-obtained data. Compared to UAV mapping, this mapping technology requires more manpower and more time. In addition, the price of the equipment is also much higher than that of the UAV. The DEM (0.5 resolution) derived from the high resolution images revealed even less detail (Fig. 6c).
Compared with terrestrial Lidar and detailed topographic mapping using total-station instruments, the UAV survey is more efficient, as a minimum of only 1 person is needed to perform the task in the field for a similarly sized mapping area. The price of the terrestrial Lidar is also much higher than that of the UAV. Lidar collects a 3D point location with a single pulse of light and is thus capable of collecting ground surface points wherever the pulse of light is able to penetrate the vegetation, reflect off the ground surface, and return to the instrument. In densely vegetated areas, Lidar may yield higher quality DEMs.
To summarize, it is important to design a survey and choose the mapping tool that suits the resolution needed for the investigation. The methods initially developed to quantify seismic offsets in this study and in similarly little vegetated areas are repeatable and can be applied to the rest of the fault after further testing and improvements. Still, the UAV system is not fully automated and requires the user to make decisions. Future studies may offer an automated UAV approach that minimizes the required level of user attention. Detecting geomorphic markers and measuring seismic offsets could also be processed automatically for large areas28.
Based on the restoration of images and 3D topography, we obtained the recent small offset of approximately ~7âm in the UAV data. It is uncertain whether offsets smaller than 7âm have been eroded or not. Although in general, strike-slip faults typically produce less relief than dip-slip faults, they are less likely to destroy the evidence of slip via uplift or erosion, or to disrupt sedimentation, erosion, or other surficial processes that produce stratigraphic and geomorphic markers25. However, previous work shown that the offset of the most recent earthquake (MRE) for the south ATF was approximately ~5âm41. The slip-rates of the fault were also suggested to be decreasing eastward as the fault keep splaying to the east as thrusts42. Therefore, the possibility of that the record of the MRE may not be preserved. Based on the UAV and satellite images, the offset values of approximately 15âm, 20âm and 30âm imply that the offsets of past earthquakes were almost constant along this segment of the fault and might indicate characteristic slip behavior43.
There is no large earthquake (Mââ¥â7) has been recorded along this segment of the ATF44, 45 for a considerable amount of time. The latest paleoseimological investigations on this segment of the fault suggest that the timing of the two MREs occurred at 1200âyr B. P., and 700 yrs B. P., respectively39, 43. Nevertheless, the seismic hazard along this segment of the fault is considerable, and further detailed work is needed to date these paleo-events. This work would have implications for the seismic risk posed by a large seismic event for the cities of Dunhuang, Akesay and the Mogao cave, which are approximately 50âkm from the fault.
Conclusions
This study focused on an approximately 1-km-long active fault section in the northern Tibetan Plateau. UAV techniques have been shown to be a powerful tool for mapping active faults at high elevations due to the ability to collect large quantities of high-resolution 3D geospatial data. A significant achievement of this study is that it demonstrated that it is possible to generate accurate DEMs of the fault at high elevations using UAV. The method is fast, and the resolution of the data is high. Along the fault, we measured the most recent seismic offsets resulting from past large earthquakes and reached the following conclusions:
-
1)
Based on the high resolution data, we not only located the fault accurately but also reconstructed the offsets with topography, which enabled us to detect more geomorphic offset markers. The UAV-based approach is demonstrated to be straightforward.
-
2)
The UAV method provided data with unprecedented spatial resolution for the identification of recent small seismic offsets.
-
3)
In the future, UAV-based high resolution data acquisition can be extensively applied to map active faults and assess seismic hazards.
References
McCaffrey, K. J. W. et al. Unlocking the spatial dimension: digital technologies and the future of geosciences fieldwork. J. Geol. Soc. Lond. 162, 927â938 (2005).
Bemis, P. S. et al. Ground-based and UAV-Based photogrammetry: A multi-scale, high-resolution mapping tool for structural geology and paleoseismology. Journal of Structural Geology 69, 163â178 (2014).
Mancini, F. et al. Using Unmanned Aerial Vehicles (UAV) for High-Resolution Reconstruction of Topography: The Structure from Motion Approach on Coastal Environments. Remote Sens. 5, 6880â6898 (2013).
Eltner, A. & Schneider, D. Analysis of Different Methods for 3D Reconstruction of Natural Surfaces from Parallel-Axes UAV Images. The Photogrammetric Record 30(151), 279â299 (2015).
Eltner, A. et al. Image-based surface reconstruction in geomorphometry - merits, limits and developments of a promising tool for geoscientists. Earth Surface Dynamics Discussions 3, 1445â1508 (2016).
Yathunanthan, V., Eun-Jung, H., Peter, K. & Steven, M. An interactive image segmentation method for lithological boundary detection: A rapid mapping tool for geologists. Computer & Geosciences 100, 27â40 (2017).
Haneberg, W. C. Using close range terrestrial digital photogrammetry for 3-D rock slope modeling and discontinuity mapping in the United States. Bull. Eng. Geol. Environ. 67, 457â469 (2008).
Linder, W. Digital Photogrammetry: A Practical Course, 3rd ed. Springer. (2009).
Stefan, A. V. & Alexander, R. C. Mapping folds and fractures in basement and cover rocks using UAV photogrammetry, Cape Liptrap and Cape Paterson, Victoria, Australia. Journal of Structural Geology 85, 168â187 (2016).
Yathunanthan, V., Eun-Jung, H., Peter, K. & Steven, M. Semi-automatic mapping of geological Structures using UAV-based photogrammetric data: An image analysis approach. Computers & Geosciences 69, 22â32 (2014).
Tavani, S., Corradetti, A. & Billi, A. High precision analysis of an embryonic extensional fault-related fold using 3D orthorectified virtual outcrops: The viewpoint importance in structural geology. Journal of Structural Geology 86, 200â210 (2016).
Johnson, K. et al. Rapid mapping of ultrafine fault zone topography with structure from motion. Geosphere 10, 969â986 (2014).
Immerzeel, W. W. et al. High-resolution monitoring of Himalayan glacier dynamics using unmanned aerial vehicles. Remote Sens. Environ. 150, 93â103 (2014).
Ragettli, S. et al. Unraveling the hydrology of a Himalayan catchment through integration of high resolution in situ data and remote sensing with an advanced simulation model. Advances in Water Resources 78, 94â111 (2015).
DallâAsta, E. et al. Use of UAS in a High Mountain Landscape: the Case of Gran Sommetta Rock Glacier (AO). ISPRS â Int. Arch. Photogramm. XL-3/W3, 391â397 (2015).
Mériaux, A.-S. et al. The Aksay segment of the northern Altyn Tagh fault: tectonic geomorphology, landscape evolution, and Holocene slip rate. Journal of Geophysical Research 110, B04404 04410.01029/02004JB003210 (2005).
Li, H. et al. Co-seismic and cumulative offsets of the recent earthquakes along the Karakax left-lateral strike-slip fault in western Tibet. Gondwana Research 21, 64â87 (2012).
Kraus, K. Photogrammetry: Geometry from Images and Laser Scans, 2nd edition, Berlin (2007).
Abdullah, Q., Bethel, J., Hussain, M. & Munjy, R. Photogrammetric project and mission planning. In: McGlone, J.C. (Ed.), Manual of Photogrammetry. American Society for Photogrammetry and Remote Sensing 1187â1220 (2013).
James, M. R. & Robson, S. Straightforward reconstruction of 3D surfaces and topography with a camera: Accuracy and geoscience application. Journal of Geophysical Research 117, F03017 (2012).
Smith, M. W. & Vericat, D. From experimental plots to experimental landscapes: topography, erosion and deposition in sub-humid badlands from Structure-from-Motion photogrammetry. Earth Surf. Process. Landforms 40, 1656â1671 (2015).
Wang et al. Gully erosion mapping and monitoring at multiple scales based on multi-source reomote sensing data of the Sancha River catchment, northeast China. ISPRS Int. J. Geo-Inf. 5, 200 (2016).
Lucieer, D. J., Turner, D., King, D. H. & Robinson, S. A. Using an Unmanned Aerial Vehicle (UAV) to capture micro-topography of Antarctic moss beds. International Journal of Applied Earth Observation 27, 53â62 (2014).
Verhoeven, G. Taking computer vision aloftâArchaeological threeâdimensional reconstructions from aerial photographs with PhotoScan. Archaeol. Prospect 18, 67â73 (2011).
McCalpin, P. J., Rockwell, K. T. & Weldon, J. R. Chapter 6 Paleoseismology of StrikeâSlip Tectonic Environments. International Geophysics 95, 421â496 (2009).
Klinger, Y. et al. High-resolution satellite imagery mapping of the surface rupture and slip distribution of the Mw ~7.8, 14 Novembre 2001 Kokoxili earthquake, Kunlun fault, northern Tibet, China. Bulletin of the Seismological Society of America 95, 1970â1987 (2005).
Zielke, O., Klinger, Y. & Arrowsmith, R. J. Fault Slip and Earthquake Recurrence Along Strike-Slip FaultsâContributions of High-Resolution Geomorphic Data. Tectonophysics 638, 43â62 (2015).
Haddon, E. K., Amos, C. B., Zielke, O., Jayko, A. S. & Buergmann, R. Surface Slip During the Large Owens Valley Earthquake. G-cubed (2016).
Gold, R. D. et al. Faulted terrace risers place new constraints on the late Quaternary slip rate for the central Altyn Tagh fault, northwest Tibet. Geological Society of America Bulletin 123, 958â978 (2011).
Peña, J. M., Torres-Sánchez, J., Serrano-Pérez, A., de Castro, A. I. & López-Granados, F. Quantifying Efficacy and Limits of Unmanned Aerial Vehicle (UAV) Technology for Weed Seedling Detection as Affected by Sensor Resolution. Sensors 15, 5609â5626 (2015).
Klinger, Y. Relation between continental strike-slip earthquake segmentation and thickness of the crust. Journal of Geophysical Research. 115, no. B07306 (2010).
Hunt, E. R. J. et al. Acquisition of NIR-green-blue digital photographs from unmanned aircraft for crop monitoring. Remote Sens. 2, 290â305 (2010).
Scaioni, M. et al. Rc-Heli and Structure & Motion Techniques for the 3-D Reconstruction of a Milan Dome Spire. In Proceedings of the 3rd ISPRS International Workshop 3D-ARCH 2009: â3D Virtual Reconstruction and Visualization of Complex Architecturesâ, Trento, Italy, 25â28 February 82, (2009).
Lienkaemper, J. 1857 slip on the San Andreas fault southeast of Cholame, California. Bull. Seismol. Soc. Am 91, 1659â1672 (2001).
Zielke, O., Arrowsmith, J. R., Grant Ludwig, L. & Akçiz, S. O. Slip in the 1857 and earlier large earthquakes along the Carrizo Plain, San Andreas fault. Science 327, 1119â1122 (2010).
Xu, X. et al. Late Quaternary sinistral slip rate along the Altyn Tagh fault and its structural transformation model. Science in China (Series D) 48, 384â397 (2005).
Sieh, K. The repetition of large-earthquake ruptures. Proc. Natl. Acad. Sci. 93, 3764â3771 (1996).
Cowgill, E. Impact of riser reconstructions on estimation of secular variation in rates of strikeâslip faulting: Revisiting the Cherchen River site along the Altyn Tagh Fault, NW China. Earth and Planetary Science Letters 254, 239â255 (2007).
Li, K. et al. Paleoseismic events in Bangouba trench along Aksay segment of the Altyn Tagh fault zone. Seismology and Geology 38, 670â679 (2016).
Heng, F. H. 3D Terrestrial laser scanning for application in earthwork and topographical surveys. University of Southern Queensland, Dissertation. (2008).
Elliott, A. J., Oskin, M. E., Liu-Zeng, J. & Shao, Y. Rupture termination at restraining bends: The last great earthquake on the Altyn Tagh fault. Geophys. Res. Lett. 42, 2164â2170 (2015).
Xu, X. et al. Near-surface character of permanent geologic deformation across the mega-strike-slip faults in the northern Tibetan Plateau. Seismology and Geology 29, 201â217 (2007).
Klinger, Y., Etchebes, M., Tapponnier, P. & Narteau, C. Characteristic slip for five great earthquakes along the Fuyun fault in China. Nature Geoscience 4, 389â392 (2011).
Seismological Bureau of Xinjiang Uygur Autonomous Region (SBX). Catalog for Earthquakes in Xinjiang 1970â2000. (2008).
Feng, X. Palaeo-earthquakes of Xinjiang. Science, Technology and Health Publishing House, Wurümqi. 1â250 (in Chinese) (1997).
Acknowledgements
We thank Mr. Allan Derrien and Dr. Kyle Bradley for assistance in assembling the UAV. We appreciate the constructive and helpful comments from Editor Francesca Cigna and the three anonymous reviewers. We thank our colleague Dr. Kang Li, who conducted the Trimble VX scanning and provided the data. This is EOS, NTU contribution 129.
Author information
Authors and Affiliations
Contributions
M. Gao and P. Tapponnier conceived and designed the investigation; M. Gao performed the image acquisition and analyzed the data; Y. Klinger, J. Van der Woerd and X. Xu contributed to discussion/revision of the paper; and M. Gao wrote the paper.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Gao, M., Xu, X., Klinger, Y. et al. High-resolution mapping based on an Unmanned Aerial Vehicle (UAV) to capture paleoseismic offsets along the Altyn-Tagh fault, China. Sci Rep 7, 8281 (2017). https://doi.org/10.1038/s41598-017-08119-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-08119-2
This article is cited by
-
Experimental verification of seafloor crustal deformation observations by UAV-based GNSS-A
Scientific Reports (2023)
-
Detailed mapping of the surface rupture of the 12 February 2014 Yutian Ms7.3 earthquake, Altyn Tagh fault, Xinjiang, China
Science China Earth Sciences (2021)
-
Real worldâbased immersive Virtual Reality for research, teaching and communication in volcanology
Bulletin of Volcanology (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.