Active tectonics of the onshore Hengchun Fault using UAS DTM combined with ALOS PS-InSAR time series ( Southern Taiwan )

Characterizing active faults and quantify their activity are major concern in Taiwan, especially following the major Chichi earthquake of September 21 st , 1999. Among the targets that still remain poorly known in terms of active 15 tectonics, are the Hengchun and Kenting Faults (Southern Taiwan). From the geodynamic point of view, they affects the outcropping top of the Manila accretionary prism of the Manila subduction zone that runs from Luzon (N. Philippines) to Taiwan. In order to better settle the location and quantify the activity of the Hengchun Fault, we start from the existing geological maps, Chinese Petroleum Corporation (CPC, 1993) and Central Geological Survey (Sung, 1991) ones, which we update thanks to the use of two products derived from Unmanned Aircraft System's acquisitions : 1) the very high precision 20 (<50 cm) and resolution (<10 cm) Digital Terrain Model (DTM), and the 2) the georeferenced aerial photograph mosaic of the studied area. Moreover, the superimposition of the resulting structural sketch map with new Persistent Scatterers (PS) InSAR results obtained from PALSAR ALOS-1 images, validated by GPS and leveling data, allow characterizing and quantifying the surface displacements during the monitoring period (2007-2011). We confirm herein the geometry, the characterization and the quantification of the active Hengchun Fault deformation that act as an active left-lateral 25 transpressive fault. As the Hengchun ridge faces one of the last major earthquake of Taiwan (December 26 th 2006, depth: 44 km, ML=7.0), it is needed to better constrain the Hengchun Peninsula active tectonics in order if possible to prevent major destructions in the near future.


Introduction
The island of Taiwan is the location of the rapid NW-SEtrending oblique convergence of the Eurasian and Philippine sea plates (e.g., Ho, 1986Ho, , 1988. A southwest propagation of the deformation in Taiwan was proposed by Suppe (1981Suppe ( , 1984. It has been confirmed by continuous GPS measurements that highlight displacements in Southern Taiwan of about 60 mm yr −1 (e.g., Yu et al., 1997). This local shortening rate is the maximum range in the world. We focus on the Southern Taiwan area (Hengchun Peninsula) that is geodynamically interpreted as the Taiwan incipient collision zone (Lallemand et al., 2001, Fig. 1), which is the transition of the Manila subduction area toward onshore Taiwan collision. The Hengchun Peninsula is situated at the summit of the outcropping onshore Manila accretionary prism (Malavieille et al., 2002;Chang et al., 2003).  . (b) Western Hengchun Peninsula hill shade from the 5 m ground-resolution digital terrain model: the heavy red lines are Hengchun and Kenting faults (modified from Sung, 1991, andChang et al., 2009b). The black triangles are GPS stations; blue and yellow arrows represent horizontal and vertical GPS displacements, respectively; and the green quadrangle corresponds to the UAS-monitored area. A-A' is the location of the (c) W-Etrending synthetic geological cross section (modified from Chang et al., 2009b;Zhang et al., 2016) (4) white star is the GPS ground fixed station GS59 and base of this PS-InSAR data set. The fault displacements come from Chang et al. (2003). This is a geologically complex area, as its lithology is composed of both soft folded turbidites (Mutan formation) and the Kenting Mélange, situated close to the Hengchun Fault zone (Fig. 1) and composed of highly deformed shales and large blocks of various lithologies. Evidence of this structural complexity are the five different 1/50 000 scale geological maps (Huang et al., 1997(Huang et al., , 2006Chang et al., 2003Chang et al., , 2009bZhang et al., 2016) that have been drawn independently and that are geologically incompatible. Two different major NNW-SSE-trending fault zones (Hengchun and Kenting faults) are parallel to the eastern side of the 15 km long Hengchun Valley and have been recognized by previous authors as potentially active (e.g., Huang et al., 1997Huang et al., , 2006Chang, et al., 2003Chang, et al., , 2009a. In addition, these two faults affect soft, muddy and alluvial deposits and classical microtectonics; structural studies are impossible as there is no clear preserved fault plane and it is impossible to get any tectonic markers. As the Hengchun Peninsula is the location of sensitive structures, such as the Taipower Nuclear Power Plant No. 3 (Maanshan Nuclear Power Plant), it is a major concern to better decipher the neotectonics and the active tectonics of the area. Consequently, we address in this paper an indirect approach using an unmanned aircraft system (UAS), a digital elevation model (DEM), structural photointerpretation and quantitative measurements in order to locate, characterize and quantify the active deformation of the Hengchun Fault.
Below, we first focus on the inputs of a high-resolution DEM and the orthorectified images of the Hengchun Fault zone acquired from UAS in order to locate active faults. From a terminological point of view, we used a UAS, which encompasses all the aspects of deploying those aircrafts and is less restrictive than an unmanned aerial vehicle (UAV), which corresponds only to the platform itself. Then, following the classical photointerpretation techniques (e.g., Deffontaines et al., 1991), we carefully mapped and interpreted the new morphostructures in terms of geological mapping, revealing the active tectonic structures using a high-resolution and high-precision UAS digital surface model (DSM). In order to characterize and quantify their tectonic activity, we benefit from Persistent Scatterer Interferometry Synthetic Aperture Radar (PS-InSAR hereafter) results derived from ALOS radar images that give us the active interseismic displacement along the line of sight (LOS) of the radar through the 2007-2011 period. The latter is validated by two E-W-trending precise leveling lines across the Hengchun Fault and a few GPS stations (see Academia Sinica GPS, 2018, with a reference line in between Taipei and the Penghu Island) that give the local vertical and absolute 3-D displacements, respectively. Finally, the discussion leads us to propose a simple geodynamic and active tectonic model for the Hengchun Fault that fits with all the available data. The black hexagons, red circles and rectangles indicate cities, ground control points and check points, respectively. One may notice the contrasting relief on both sides of the Hengchun Valley: its western flank is a homogenous structural surface slightly dipping east, in contrast with the highly eroded eastern part. For both (a) and (b), the black and white hill-shade DTMs correspond to a part of the 5 m ground-resolution Taiwan DTM.

High-resolution DSM obtained from UAS
A UAS, commonly known as a drone, is an aircraft without human pilot on board that flies autonomously and is being used in many areas because of its convenience, high resolution, etc. (Huang and Chang, 2014;Fernandez Galarreta et al., 2015;Giordan et al., 2015;Tokarczyk et al., 2015;Bühler et al., 2016;Deffontaines et al., 2017). The UAS used in this study is a modified version of the already-available Skywalker X8 delta-wing aircraft reinforced by carbon fiber rods and covered partially by Kevlar fiber sheets. The drone is launched by hand, flies and takes photos autonomously, then glides back down to the ground by using a pre-programmed flight plans organized by a ground control system; it is controlled by the ground control station and remote controller. The autopilot system used in this study is composed of and modified from the open-source APM (ArduPilot Mega 2.6 autopilot) firmware and open-source software Mission Planner, transmitted by ground-air XBee telemetry. In order to generate a high-resolution DSM (defines the terrain data including buildings and tree canopy) and mosaic orthorectified images, a total of 3767 photos were gathered by Sony ILCE-QX1 camera mounted on the crafts, during five fly missions, for a total coverage area of 33 km 2 with 8 cm ground sam-pling distance. The coverage of the adjacent photos is kept at 85 % overlap and 45 % side lap. The raw images are processed by ContextCapture and Pix4Dmapper software mutually; the data sets generated in this study are orthoimages and DSMs, with a grid spacing of 10 cm. Prior to the morphotectonic analysis (e.g., Deffontaines and Chorowicz, 1991;Deffontaines et al., 1993Deffontaines et al., , 1994Deffontaines et al., , 1997Deffontaines et al., , 2017Pubellier et al., 1994), the quality of the data set (see whole data set in Fig. 2) needs to be evaluated. Eighteen ground control points are extracted with grid size of 2 m resolution of airborne light detection and ranging (lidar) data set and from the airborne lidar-associated 25 cm resolution orthorectified image. Most of the ground control points are situated on crossroads, targeted and georeferenced from orthorectified images, and the elevation is from the airborne lidar data. The comparison of the UAS DSM with airborne lidar data gives a root-meansquare deviation of 4.1 cm with maximum error of 42.5 cm from 26 sites of open, bare ground area, e.g., roads, school playgrounds, unvegetated terrains and parking lots. The elevation of the check points is averaged from an area of 4 m 2 , equal to the grid size of airborne lidar data. The distribution of the ground control points and checked points are located in Fig. 2b. Figure 3 highlights the precision and resolution difference on the available topographic data sets. The 40 m grid DEM was generated in 1986 by the Aerial Survey Office from aerial photos, the 5 m grid DEM was generated by several aerial photogrammetric agencies in 2004, whereas most of the aerial photos were acquired in 2003. The project to create a 5 m DEM data set was funded by the Taiwan Min-istry of the Interior. The airborne lidar data were acquired in 2011 and used to generated a 2 m grid DTM (defines the geomorphologic elevation after removing the buildings, trees and vehicles) and DSM, but this data set is not authorized to be published even though it was generated by the authors. . The first line of images highlights the active tectonic scarp of the Hengchun Fault, which differentiates two blocks: uplifted (+), locally eroded terraces to the east and a flat lowland sedimentation area to the west. The second line of images highlights the left-lateral offset geometry (red arrows) of tributaries (blue line) on a glacis which is uplifted (+) to the east of the Hengchun Fault. The third line of images shows a deformed N-S-trending elongated mud volcano (MV -Ma'anshan) previously evidenced by Giletycz (2015), situated close to NPP (see infrastructures on the lower left corner of the deformed mud volcano) and a flat, uplifted (+) terraces (paleosurface of Fig. 6) on the western side of the captions.
The ground control points used in this study were extracted from the airborne lidar data set we acquired some years ago. Figure 3 highlights visually the difference of quality, precision and resolution of each digital topographic data set. One can see the low quality of the 40 m DEM in contrast to the precision and resolution of our UAS DSM. Figure 4 illustrates some case examples of the quality of both the UAS DSM as well as the orthorectified image of the Hengchun Fault. The morphostructures analyses based on photo interpretation are conducted accordingly (Fig. 4c,d,g,h,k,and i).
3 Geology of the Hengchun area 3.1 Hengchun geological state of the art From a geologic and geodynamic point of view, the Hengchun Peninsula is interpreted as the northern tip of the Manila accretionary prism (Malavieille et al., 2002;Chang et al., 2003). It is mainly composed of several lithological formations that are described below from the older to the younger ones (see their location on the cross section of Fig. 1c).
The Mutan formation is the largest outcropping formation of the Hengchun Peninsula. Middle-to-late Miocene in age, it is composed of folded classical turbidites made of shales interbedded with sandstones and conglomerates (channels and levees) (Sung and Wang, 1986;Sung, 1991). The Mutan formation is a very thick pile (locally several thousands of meters) and highly deformed turbidites (locally even overturned: e.g., CGS geological map by Sung, 1991;. The soft complex Ma'anshan formation, late Miocene-Quaternary in age, is a mix of submarine erosion and depressions filled in with both muddy Mutan and Kenting formations (Cheng and Haung, 1975;Page and Lan, 1983;Lin and Wang, 2001;Huang, 2006). The Ma'anshan formation is folded and outcrops partly in the south of the Hengchun Basin (CGS geological map by Sung, 1991); see panorama of the Hengchun Basin in Fig. 5. We believe that the Ma'anshan formation composes partly the basement of the Hengchun Basin and the offshore Manila accretionary prism to the south (Lundberg et al., , 1997Reed, et al., 1992).
The Kenting Mélange, which only outcrops east of the Hengchun Basin, is known as a tectonic mélange (e.g., Tsan, 1974;Page and Lan, 1983;Huang, 1984;Pelletier and Stephan, 1986;Huang et al., 1997;Chang et al., 2003Chang et al., , 2009aDeffontaines et al., 2016;Zhang et al., 2016;Malavieille et al., 2016). It is interpreted as a chaotic tectonostratigraphic unit containing blocks of various origin (ophiolitic and sedimentary), size and lithology embedded in a sedimentary matrix (Malavieille et al., 2016). It is herein interpreted as a large tectonic breccia of composite ages that ranges from 1 to 10 Ma (see Malavieille et al., 2016) as it mixes various sizes of different blocks and different turbiditic and shale bodies (Page and Lan, 1983;Lin and Wang, 2001;Chang et al., 2003). In our opinion, this shaly, blocky tectonic mélange is injected (intruded) within the major fault zones south of Taiwan due to both the plate convergence and high pressure of fluids at depth (see also Deffontaines et al., 2016). Future scientific work is needed regarding the differing, incompatible assumptions and hypotheses in order to both characterize the processes that create this tectonic Kenting Mélange and to reconstruct (or inverse) its evolution through time (see synthesis in Malavieille et al., 2016).
The slightly tilted Hengchun limestone overlays the unconformable Ma'anshan formation on the western side of the Hengchun Valley and overlying the Kenting Mélange on the eastern side of the Hengchun Fault. As the Hengchun limestone is folded differently and continuously on both sides of the Hengchun Fault, it is unfortunately difficult to use the existed dating results to estimate a precise vertical uplift rate in the absence of the precise location of sampling (Chen and Liu, 1993).

Updated Hengchun geology and neotectonics
Geological mapping in the fields lead us to look for active tectonic cracks within anthropic concrete dykes. They reveal the clear lateral transpressive and thrusting component of the northern part of the Hengchun Fault (Fig. 6). As the concrete dykes are manmade and very recent, it is easy to reconstruct and inverse the deformation following the Angelier's methodology settled along the longitudinal valley fault (see also Lee et al., 1998Lee et al., , 2000Lee et al., , 2001Lee et al., , 2003Lee et al., , 2005Lee et al., , 2006. From a structural point of view, the Hengchun Fault geological map that we propose (Fig. 7) is the result of (1) field studies (Figs. 5 and 6 taken 3 January 2007), (2) the synthesis of previous geological mappings (Rokkaku and Makiyama, 1934;CPC, 1993;CGS in Sung, 1991;Chang et al., 2003;Chen et al., 2005;Giletycz et al., 2015) and (3) a detailed photointerpretation of both the orthorectified mosaic of the UAS aerial photograph as well as the hill shading of the highresolution UAS-derived DSM (Fig. 7). The detailed structural photointerpretation of the Hengchun Fault area checked in the field takes into account the basic morphostructural principles (Deffontaines et al., 1991(Deffontaines et al., , 1993(Deffontaines et al., , 1997Pubellier et al., 1994) such as the geometry of the drainage pattern (bayonet tracks, curves and alignments) and the alignments of small scarps present in the flat Hengchun Basin (Fig. 4e, f, g, h, and Fig. 7, index 5) that help to characterize the morphotectonic features (Fig. 7, faults that affect the topography, index 2 for instance). We also added on the Hengchun Fault geological mapping, geomorphological objects such as paleosurfaces (index 5 in Fig. 7), which are easily recognized by smooth textures, and structures above the Kenting Mélange. Those paleosurfaces are slightly tilted, folded and uplifted east of the Hengchun Fault and are witnesses of the vertical uplift of the transpressive tectonic activity of the Hengchun Fault. Unfortunately, due to the slight tilt, it remains difficult to extract fold axes from those paleosurfaces. It would be interesting to date these paleosurfaces in order to get the mid-to long-term activity of the deformation. We select also the chevrons from the morphostructural analysis of the high-resolution UAS DSM (index 6 in Fig. 7). Chevrons correspond to the top of harder rocks of the tilted structural surfaces within the Kenting Mélange east of the Hengchun Fault. They reveal opposite tilted strata, indicating the presence of folds (anticlines and synclines, index 7 in Fig. 7) close and parallel to the Hengchun Fault zone. Some tectonically offset and eroded mud volcanoes are also identified and mapped herein (index 3, in Fig. 6), especially in the South Hengchun Basin and close to Taipower Nuclear Power Plant No. 3. Those situated to the south of the study area were already evidenced by Giletycz (2015). They appear in the UAS DSM as characteristic "elongated conic volcanic" shape relief above the flat marine and alluvial lowland terrace. Their western side is locally affected by landslides (Fig. 7). We interpret these mud volcanoes as the outcropping traces of the shale intrusions within the Hengchun Fault zone, very probably due to fluids overpressure at depth in relation to the intense transpressive tectonic stress (see also Deffontaines et al., 2016).
The Hengchun Fault is recognized as active (Bonilla, 1976;Hsu and Chang, 1979;Lee, 1999;Lin et al., 2000;Lin et al., 2009;Deffontaines et al., 2001Deffontaines et al., , 2016. Effectively, geological mapping in the field reveals many active cracks within anthropic concrete dykes. These cracks highlight the clear transpressive and thrusting component of the Hengchun Fault (Figs. 5 and 6). We cannot recognize clearly in the topography the trace of the Hengchun Fault as a unique straight line thrusting (see Fig. 7) in contrast to the CGS geological mapping (Sung, 1991) or its sinusoidal shape on the CPC geological map (1993). We only highlight herein the Hengchun Fault zone, which is made of locally aligned segments of faults that affect the eastern part of the very recent Hengchun Basin deposits (see index 2 in Fig. 7). Note that regional geological maps aim to show both stratigraphy and structure, where the fault trace is at the maximum (or just approximately) stratigraphic separation. In contrast, the active fault trace is situated strictly at active lines where last rupture is located. Thus, these two facts may create large differences. Moreover, on land, the last fault rupture can only be identified by geomorphic features, not by drilling or geophysical prospecting.
To date the deformation, the vertical long-term slip rate has been deduced in marine terrace studies and estimated  Chen and Liu, 1993). Other marine terrace studies from Chen et al. (2005) show similar long-term slip rate.

Inputs of PS-InSAR interferometry on the Hengchun Fault
In order to deduce the regional tectonic activity of the Hengchun Fault, we processed 13 radar images (ALOS PAL-SAR images acquired in L-band; λ = 23.6 cm, average incidence angle = 34.3 • ) in ascending mode (flight direction 350 • E of N) between 17 January 2007 and 28 January 2011, thus providing 4 years of monitoring, through StaMPS software (Hooper et al., 2007). A series of 12 interferograms has been generated, all with respect to a unique "super master" as a reference image situated in the mid-time series that minimizes both spatial and temporal decorrelation (Pathier et al., 2003;Champenois, 2011, Champenois et al., 2011, 2012). Note that their perpendicular baselines situated between −1194 and 2417 m are significantly below the critical limit where all interferometric coherence is lost (Hooper, 2009). We used the 40 m ground-resolution Taiwan DTM in order to remove both the topographic and orbital components of the interferometric phase. All interferograms shows an excellent coherence in the Hengchun Plain as well as in the hilly area. The method identifies 20 133 PS pixels with an average density on the whole Hengchun Peninsula of 120 PS km −2 (see Fig. 9), characterized by their phase stability over the whole period. PS analyses lead us to extract a map of mean LOS surface velocities for each PS and also to catch their displacements through time series. Note that these LOS data are not projected onto the vertical component. The PS-InSAR base (fixed GPS station GS59 corresponds to the black and white star see its location close to Checheng and north of the Hengchun Valley in Figs. 1 and 9) is chosen as it presents a stable (to very low) continuous deformation during the monitoring time period. Nevertheless, this location might be subject to small continuous uplift or small continuous subsidence that may explain local discrepancies with the GPS average annual displacement. We compare our PS-InSAR results (Champenois, 2011;Champenois et al., 2011Champenois et al., , 2014 to three fixed GPS stations (HENC, GS57, GS59; see Fig. 9) and two leveling lines acquired in the same monitoring time period (see Fig. 10). For each of the three fixed GPS stations (HENC, GS57, GS59), an average displacement has been calculated, projected into the radar LOS by taking into consideration the various local incidence angles along the distance axis (Hanssen, 2001;Cigna et al., 2011). The results (Fig. 9) show an impressive density of PS (more than 140 PS km 2 in the surroundings of the Hengchun Fault) in this area, where both human activity and luxuriant vegetation prevail. The PS map reveals a clear continuous difference in LOS velocities across the Hengchun Fault where the Hengchun Valley (western side of the Hengchun Fault) subsides, with values ranging from −5 to −10 mm yr −1 , in contrast with the westernmost Hengchun ridge, which is uplifting. This is evidenced in Fig. 9 by a rapid transition from blue/green to orange/red PS dots. Taking into account the PALSAR ascending satellite orbit and the NNW-SSE Hengchun Fault, trending of both geometries is close to the best configuration to detect and measure the fault displacements (Champenois, 2011). Consequently, Fig. 9 led us to draw, with much more accurate detail than Shyu et al. (2005)   Fault is characterized by the rapid increasing uplift of the Kenting Mélange up to a maximum that seems to parallel the Hengchun Fault. Farther east the PS LOS velocities decrease with the distance to the Hengchun Fault so the eastern Hengchun Peninsula appears to be slightly tilted. We observe in the Hengchun Peninsula a similar dissymmetric active vertical behavior as in the coastal range (eastern Taiwan): the western part of the relief is highly actively uplifting and folding in contrast to the eastern part, which is characterized by a decreasing uplifting and which is slightly tilted to the east (see also Deffontaines et al., 2017). We note in Figs. 10b and 11 the coherence of the PS LOS velocities and the GPS measurements (colored large dots that we put into the same LOS radar geometry). Consequently, those PS-InSAR results offer a precise mapping of the active interseismic splays of the Hengchun Fault and allow the quantification of the tectonic displacements.
Regarding the Kenting Fault, with the data used in this study we see no gradient of PS LOS velocity nor noticeable differential displacements on the GPS stations situated around the Kenting Fault (see Fig. 9 and profiles of Fig. 10ab). Moreover, from the geological mapping analysis , the Kenting Fault has a sinusoidal planimetric geometry that follows more or less the topographic isocontour of the relief. So the Kenting Fault is a low-dipping thrust fault. Consequently, this fault may be either locked at the surface or an inherited inactive fault, or it may have a tectonic activity that we do not observe on this PS-InSAR data set during the monitoring time period.
We note a great coherence in between the PS and the available GPS measurements (colored large dots in Fig. 9) provided on two stations. For comparison purposes, the average annual GPS displacements acquired in the InSAR monitor- ing time period have been compared with the radar LOS. We compared also these PS results with leveling data acquired during the same period of time (between March 2008 and April 2011, during three campaigns of measurement). Figure 10 compares the leveling data, the mean LOS velocities of PS along the leveling lines and GPS data for the two stations situated close to the leveling lines. Note the agreement between all of them. The local heterogeneity of the displacements (e.g., small subsidence and small uplift) within the Hengchun Valley is highlighted by both the HENC GPS sta- tion that gives the absolute uplift displacement and the PS-InSAR data with various colored dots, ranging from blue to orange, situated on both sides of the O value.
Moreover, in order to better characterize the Hengchun Fault activity and quantify the LOS velocity variation on both sides of Hengchun Fault, we computed the LOS velocity offset across the fault, using 31 profiles perpendicular to the fault line. Each 10 km long PS profile is perpendicular to the Hengchun Fault and characterized by an equidistance of 400 m. In addition, each PS data point situated on both sides (200 m vicinity) is projected on each PS profile. Each data point is only represented on only one profile. The number of PS dots per profile varies between 400 and 1200 (Champenois, 2011). This offset characterizes the slip rate of the fault. We present two specific profiles in Fig. 11 (see location Fig. 9). The two PS profiles, which run transverse to the Hengchun Fault, one north of the Hengchun Valley close to the Checheng city and the second one southward (close to Hengchun city), show LOS offset less than 10 ± 2.5 mm yr −1 . Note the high number of PS reflectors situated in the cities of Checheng and Hengchun that confirm that on both profiles we evidence a clear jump in the PS magnitude that corresponds to the aseismic displacement of the Hengchun Fault (Champenois, 2011).
So thanks to a series of more than 30 profiles perpendicular to the fault, we computed the interseismic LOS velocity offset along the fault, from north to south, with LOS values (from min 4.9 mm yr −1 to max 10.2 mm yr −1 ), reaching an average 8 ± 2.5 mm yr −1 during the monitoring time pe-riod, which reveals, using two linear regressions calculated on both side of the Hengchun Fault, the clear tectonic activity of the Hengchun Fault (Fig. 12). Note the good coherence between slip rates derived from PS analysis and the GPS measurement (projected in LOS by taking into account each local incidence angle along the distance axis) represented by the yellow triangle in see Fig. 12. This good fit between the fixed GPS stations validates our PS results. This is a major result from this study as it demonstrates the along-strike interseismic activity of the fault (Champenois, 2011;Champenois et al., 2011Champenois et al., , 2012 this study) and allows us to obtain the average annual interseismic LOS displacement of the Hengchun Fault and its spatial variation along the fault. Variations of the along-strike velocity offset of Hengchun Fault need to be carefully monitored and analyzed in order to better characterize potential seismic hazards, especially the low interseismic creeping areas where stress may accumulate (see Fig. 12) and may correspond to a partially locked part of the fault source of future earthquakes (close to Hengchun Township, for instance).
StaMPS also gives access to time series of displacement, from 2007 to 2011 (see Champenois, 2011;Champenois et al., 2011Champenois et al., , 2012. We can thus derive at each image acquisition date the displacement related to the fault (a cumulated slip value) and then monitor its slip variation through time. The time series of slip values (Fig. 13) does not show any clear seasonal accelerations or decelerations, but only expresses a rather linear behavior of the fault during the monitoring time period. This time series should be updated with new SAR im-  ages in order to get a longer monitoring period (Champenois, 2011;Champenois et al., 2011). It will help to know whether the Hengchun Fault has a linear displacement through time (or not).
From all used geodetic techniques (GPS, leveling lines and interferometric), if we generalize the two leveling sections, which are coherent to both the three fixed GPS stations and the PS-InSAR results, the western Hengchun Peninsula is slowly uplifting (less than +10 mm yr −1 ), contrasting the slow subsidence of the eastern part of the Hengchun Peninsula (0 to −2.5 mm yr −1 ). So as we reveal a continuous active differential displacement across both sides of Hengchun Fault, which is deduced from leveling lines, GPS and PS-InSAR, Hengchun Fault shows the existence of a clear interseismic creeping component during the monitoring time pe-riod, with an average value of 8 mm yr −1 along the LOS, confirmed by the leveling monitoring (see leveling line 1: 8 mm offset - Fig. 10b). The variation of displacements on both sides of the Hengchun Fault zone also reveals the progressive folding of the Hengchun Fault zone: the eastern hanging wall is progressively shortening and uplifting (as revealed by GS57 and the two leveling lines) contrasting to the progressive subsidence of the Hengchun Basin (footwall situated west of Hengchun Fault -see Figs. 10b, c and 11, profiles 1 and 2). These displacement results (Figs. 9, 10 and 11) combined with the structural photointerpretation (Fig. 7) show that the fault is only locally outcropping and that active folding occurs, increasing progressively with time, stress and strain at depth.
Consequently, due to the rectilinear structural geometry, the Hengchun Fault appears to act as an almost vertically active left-lateral transpressive strike-slip fault with an active uplifting and folding (anticline) on its eastern hanging wall, in contrast to the active subsidence of the western footwall (which corresponds to the easternmost part of the dissymmetric Hengchun Basin, Fig. 14). Moreover, the Hengchun Fault is not outcropping along one small, straight, continuous north-to-south line but is instead a rather large "fault zone" made of little active fault segments (see location mapped in details, Fig. 7). One may note that the thrusting and left-lateral strike-slip motion have been deduced from both from the morphostructural analysis of the high-resolution DTM (see Fig. 4c,d,g and h) and the obliquity of the GPS arrows already evidenced in the bibliography (see blue arrows in Figs. 1 and 9, and Chang et al., 2003). Consequently, Fig. 14 proposes a simple model where the Hengchun Fault acts as a 160 • E N-trending transpressive left-lateral strikeslip fault zone. Its active hanging wall is overthrusting and is actively folding the topographic surface that tilts and folds the paleosurfaces: the Hengchun limestone and marine terraces (see Lin, 1987, andLiu, 1993). The Hengchun Basin (Hengchun Fault footwall and west of it) under the transpression is also folded and is actively subsiding with a parallel to the fault syncline axis that create the lagoonal/marine deposition of the Pleistocene-Holocene Hengchun Basin (Chen et al., 1991).
Thus, the Hengchun Fault, as many faults in central and southern Taiwan, presents both left-lateral strike-slip and thrusting components such as the Chelungpu Fault reactivated during the Chichi earthquake (21 September 1999) as well as any Taiwan foothill transfer fault zones (e.g., Deffontaines et al., 1997).
Regarding the amplitude and the components of the active Hengchun Fault deformation, the two leveling lines reveal the vertical offset of the Hengchun Fault, consistent with interferometric displacement toward the LOS, with a value of 8 ± 2 mm yr −1 , which corresponds to the vertical geodetic interseismic slip rate. Note that this short-term value is a bit different than the vertical long-term slip rate deduced from the marine terraces dating 6.3 mm yr −1 at Hȃikȏu and 3.8-6.1 mm yr −1 at Nanwan (Liew and Lin, 1987;Chen and Liu, 1993). From this slight difference in between the long-term and the vertical geodetic slip rates, we may deduce the fact that the Hengchun Fault may be not only a creeping fault. But due to the folding component of the eastern side of the Hengchun Fault, it is needed to locate much precisely the in situ dated samples (Chen and Liu, 1993) to reconstruct and inverse correctly the long-term deformation history. Chang et al. (2003) characterized the active horizontal slip rate of the western Hengchun Peninsula by comparing the directions and amplitude of the three fixed GPS stations. Note that these values take into account both Hengchun and Kenting faults without individualization.

Conclusions and perspectives
First, the flights of UAS with image acquisition allow us to produce through photogrammetry a high-resolution DSM of the Hengchun Fault area with a less than 7 cm ground resolution and less than 40 cm precision. However, for some places, the texture may change rate quicker and higher than the image capture acquisition, such as the lakes or water surfaces due to sunlight reflection, which causes noises in DEMs. In addition, concerning the DEM photointerpreta-tion, especially for the active geological structures, the terrain is easily modified by human activity in Taiwan: cities and farmland are rapidly growing in some areas, concealing morphotectonic structures. Fortunately, the orthomosaic images help to recognize those artifacts. Overall, the autonomous UAS and well-developed photogrammetry technique allow us to generate conveniently highly detailed topographic information and help to carry out large-scale, highresolution morphotectonic study. Integrating leveling lines, three fixed GPS stations and PS-InSAR processing that gives the short-term active movement of the Hengchun Fault and high-resolution DSM structural photointerpretation acquired from UAS demonstrates longer-term topographic deformations and anomalies, thus providing one new scope that may be able to locate, characterize and quantify the active tectonics features.
Then with these useful high-resolution topographic data, we undertook a classical geological and geomorphological photointerpretation that led us to update and refine the preexisting CPC (1993) and CGS (Sung, 1991) geological maps on both sides of the Hengchun Fault (see Figs. 4 and 7). The latter led us to characterize the structural geometry of the Hengchun Fault in the peninsula (Fig. 7). Moreover, the leveling lines and the PS-InSAR interferometric processing that we carried out over the Hengchun Peninsula from 2007 to 2011 with more than 12 coherent interferograms reveals clear interseismic displacements of the Hengchun Fault. Moreover, our interferometric results, which are in agreement with leveling lines and fixed GPS data, reveal the active sections and folding of the eastern hanging wall of the Hengchun Fault of about 8 ± 2 mm yr −1 . Furthermore, the footwall is also subject to a complementary active subsidence related to the development of an active syncline parallel to the Hengchun Fault. These displacements are a cumulative vector that takes into account both planimetric and vertical components of the deformation. Our PS-InSAR results are fully coherent with both E-W-trending leveling lines across the Hengchun Fault and three continuous GPS stations situated on both sides of the Hengchun Fault. Consequently, we expand our interpretation due to the PS density (more than 140 PS km −2 ), which offers a spatial density of measurements that is greatly higher than that offered by GPS, allowing us a global and whole coverage and mapping of the active Hengchun Fault deformations (Fig. 9).
From the active tectonics conclusions, we deduce that the Hengchun Fault, characterized by a high dip angle, acts as a large, creeping interseismic left-lateral strike-slip fault zone with a clear transpressive vertical component of 8 ± 2 mm yr −1 (leveling lines and LOS) associated with active folding. We show that the Hengchun Fault presents a different tectonic behavior from north to south: in its northern part, the Hengchun Fault is narrow and clearly outcrops below the marine terrace with a measurable offset (see Fig. 4d). That contrasts with the central and southern part of the Hengchun Fault, where it is wider or larger (see many parallel and collinear faults and cracks that offset the topography in Fig. 7) and submitted to active folding, and the fault is only locally partially locked (close to the Hengchun city, see Fig. 12). In the central and southern parts, as the PS interferogram spectrum increases rather continuously from west to east, transverse and above the whole Hengchun Fault zone, the Hengchun Fault actively folds at depth. Moreover, the lack of clear along-strike PS discontinuity across the Hengchun Fault zone proves that it is not clearly offsetting the topography (and outcropping) during the acquisition monitoring time period (2007)(2008)(2009)(2010)(2011). This should sign the slow continuous increase of tectonic stress at depth, which is a major indicator of a future earthquake in that area.
Unfortunately, the active differential displacement on both sides of the Kenting Fault is much more difficult to highlight with the present geodetic data sets and our PS-InSAR results due to (1) the low fault dip angle deduced from its mapping geometry, (2) an inherited inactive fault, (3) a locked fault at the surface, and possibly (4) the ductile rheology of the underlying formation (Mutan clay turbidites -MT). Consequently, better geologically mapping and monitoring of the Kenting Fault is needed in order to better understand its geometry, rheology and active tectonic behavior.
Nevertheless, the shallow seismicity (Deffontaines et al., 2016) confirms the high onshore shallow tectonic activity that hits the Hengchun Peninsula. The major, deep offshore Hengchun earthquake (26 December 2006, depth: 44 km, M L = 7.0) also participated in the active deformation of the Hengchun Peninsula.
Further dating ought to first be done with respect to the Hengchun Peninsula folding then undertaken in order to better characterize the ages of the terraces, paleosurfaces and the different deposits of the Hengchun Basin in order to inverse and reconstruct its tectonic history (Liew and Lin, 1987;Chen and Liu, 1993).
Finally, due to the average Hengchun Fault along-strike variations (around 8 ± 2 mm yr −1 displacement toward the LOS confirmed by leveling lines - Fig. 12), it is definitely important to better constrain the Hengchun Peninsula's active tectonics in order to prevent major destruction and major failure in the near future, such as damage to sensitive energy industries (e.g., Taiwan Nuclear Power Plant No. 3, situated along the southern tip of the Hengchun Fault).