1 Introduction
The assessment of seismic hazard is generally based on one of two approaches. In stochastic studies, the aim is to estimate the probability of occurrence of earthquakes of different magnitude at the location in question during a specified period. Deterministic assessment deals primarily with the largest credible earthquake at the site during the current tectonic regime. Seismic zoning, being concerned with the effects rather than the causes of the seismicity, will therefore tend to draw on both schemes. In the USSR, for example, zoning combined seismological, geological and terrain factors, and some local building codes in the USA take account of probabilistic studies of peak ground velocity. Other procedures are more firmly rooted in structural geology albeit with a historical component. The seismic gap model, which is largely confined to subduction settings, identifies parts of the potential seismogenic zone that has not experienced a significant event during the period of record. The asperity model is effectively its converse, as it seeks locations along such a zone where earthquakes have occurred repeatedly.
A third approach is to define belts of heightened deformation within which seismicity is likely. It has been applied to peninsular India on the basis of intraplate buckling [16]. An earlier version of this proposal was advanced for Taiwan; it identified structural belts where fault geometry was favourable to strain release in major earthquakes [15]. One of the structural units that was analysed in order to assess variations in potential seismicity in Taiwan was the Hengchun Peninsula, the southern extremity of Taiwan (Fig. 1). The primary aim of this paper is to establish how far that assessment needs to be modified by new information on the Holocene uplift of the peninsula.
2 The Hengchun Peninsula
South of Taiwan the Eurasian plate is being overridden by the Philippine Sea plate at the Manila Trench [1,5]. The resulting accretionary wedge is marked by the Hengchun Ridge, which extends onto Taiwan as the Central Range and the fold-and-thrust zone west of it [11] (Fig. 1). In the southern part of the island, these two belts are separated by the Chaouchou thrust fault [7], which extends into the Hengchun Peninsula as the Hengchun fault (Fig. 2). The Miocene-Pleistocene strata that make up the peninsula display other reverse faults and recumbent folds.
Convergence of ∼86 mm yr−1 between the Philippine Plate and Eurasia is oblique at about 307° [11], and the suture, which is marked by the Longitudinal Valley, propagates southwards. The GPS evidence for 1990–1994 [19] indicates rapid shortening at the Longitudinal Valley (29.6 mm yr−1) and the fold-and-thrust belt (locally as high as 20 mm yr−1), and extension both at the Coastal Range (0.6 mm yr−1) and at the Central Range (8.6 mm yr−1). Levelling indicates uplift of the Central Range by 36–42 mm yr−1, compared with long-term rates of 7.15–15 mm yr−1 derived from fission-track and K/Ar dating [9].
Holocene palaeoshores on the east coast of Taiwan, in conjunction with microseismic records, point to distributed deformation on imbricate reverse faults within the Coastal Range [17], which could well account for geodetic evidence of extension by giving rise to anticlinal folding above blind faults [6]. Age/height data for the islands of Lanhsu and Lutao point to slow (1.3 mm yr−1) uplift of the former until about ∼5000 yr ago and continuing, more rapid (2.6 mm yr−1) uplift of the latter; the difference in tectonic history is consistent with a reduction in the rate of uplift as the plate suture closes progressively southward [15]. In contrast, the Longitudinal valley combines reverse faulting with aseismic strike slip, and in the fold-and-thrust western foothills relatively low angle structures favour strain release in large earthquakes of which the 1999 Chi-Chi earthquake is a prime example.
To judge from the instrumental evidence, the Hengchun area at the southern extremity of the Central Range experiences few significant earthquakes, although two substantial events have been located south of it: the Hengchun earthquake of 1959.08.15 (21°42′N, 121°18′E, ), and the event of 1955.4.4 (21°48′N, 121°E, ). Both are associated with the belt of active seismicity along the Luzon arc [10], that is the volcanic arc at the rear of the accretionary complex.
3 Uplift history
Shoreline data for the Hengchun Peninsula reported in the classic study by Liew and Lin [8] indicated a rate of uplift of about 4.7 mm yr−1 [15]. As in the Coastal Range, geodetic evidence for extension could be explained by near-surface folding above a series of reverse faults. Indeed, Liew and Lin [8] describe a folded raised beach that is cut near its culmination by a thrust fault (‘A1’) striking N20°–30°W and dipping steeply eastward (Fig. 2). They took their A1 fault to be part of the Hengchun fault system; the 1/50 000 geological map of the Hengchun Peninsula [2] shows the Hengchun fault as a high-angle thrust dipping east, parallel to and about 1.5 km west of the proposed location of A1 (Fig. 1). The Hengchun fault juxtaposes the Mio-Pliocene Kenting Formation with gently folded Pleistocene and Holocene alluvium, coral reef, siltstone and limestone. A second major reverse fault, this time with a westward dip (the Tameihsi Fault), crosses the coast at Yakoulu; to judge from published sections [2] the fault acted as backstop to a series of north–south folds in the Miocene clastics that make up much of the peninsula (Fig. 3). Both faults can be traced north to the Tulungyan Fault, which forms the western boundary of the Central Range [5].
In order to evaluate the extent of Holocene deformation, the existing twenty age/height values [8,15], including three which had been unaccountably omitted in the earlier account, were supplemented by a further eight on samples collected for the present study. The ages were calibrated using the 1999 (Rev. 4.1.2) version of the program by Stuiver and Reimer [13] after allowing for the small marine reservoir correction value of −17 yr proposed by P. Reimer for the western Pacific (http://radiocarbon.pa.qub.ac.uk/marine/). The elevations were adjusted for the glacio-hydroisostatic component (as before) using data derived by interpolation between sea-level curves for 22°N 121°E and 25°N 122°E by Prof. W.R. Peltier (July 1997, pers. commun.). As most of the samples consisted of redeposited shells and coral fragments, the ages must be considered to be maximum values for the sample sites.
Table 1 lists the various age/height values. The enlarged dataset supports the earlier estimate of Holocene uplift, as the least squares line indicates an emergence rate of 3.8 mm yr−1 () (Fig. 4a). In an attempt to distinguish between net uplift and the effect of fault-related folding the age/height values for the four sites furthest from the fault, namely Yakoulu, Fantzuliao, KT1 and Shanhaili, were plotted separately. The resulting graph () (Fig. 4b) indicates general uplift at about 2.5 mm yr−1, so that the remaining 1.3 mm yr−1 of uplift is attributable to deformation associated with the Hengchun Fault: the ‘upbuckling’ described by Liew and Lin [8].
14C ages for the Hengchun Peninsula
Âges carbone 14 relatifs à la péninsule de Hengchun
Site | Elev. (m) | Corr. elev. (m)1 | 14C age (yr BP) | Cal. age2 (yr BP) | Source3 |
Yakoulu | 10 | 9.5 | 3700 ± 250 | 3630 ± 680 | a |
Kenting | 1 | 1 | 1580 ± 20 | 1165 ± 250 | a |
12 | 11 | 4310 ± 120 | 4430 ± 375 | a | |
14 | 13.5 | 3900 ± 125 | 3885 ± 355 | a | |
14 | 13.5 | 4040 ± 120 | 4080 ± 650 | a | |
14 | 13.5 | 4600 ± 400 | 4820 ± 1100 | a | |
20 | 20 | 5200 ± 95 | 5580 ± 250 | a | |
20 | 20.5 | 5730 ± 480 | 6170 ± 1140 | a | |
20 | 37 | 8410 ± 155 | 8920 ± 475 | a | |
20 | 38.5 | 8660 ± 155 | 9055 ± 715 | a | |
22 | 32 | 7810 ± 155 | 8310 ± 290 | a | |
Shihniuchiao | 1.5 | 1.5 | 1710 ± 110 | 1275 ± 230 | a |
1.5 | 1.5 | 1780 ± 115 | 1325 ± 260 | a | |
15 | 14.5 | 5000 ± 300 | 5325 ± 770 | a | |
15 | 15.5 | 5670 ± 140 | 6100 ± 350 | a | |
15 | 20 | 7050 ± 145 | 7555 ± 260 | a | |
15 | 14 | 5000 ± 400 | 5325 ± 1020 | a | |
15 | 40.5 | 9600 ± 1000 | 10300 ± 2840 | a | |
KT5 | 17 | 38 | 8760 ± 60 | 9225 ± 535 | b (NTU-2636) |
KT6 | 17 | 23.5 | 8300 ± 60 | 8855 ± 410 | b (NTU-2637) |
KT | 1.5 | 7 | 7270 ± 60 | 7730 ± 110 | b (NTU-2643) |
KT8 | 0 | −0.5 | 4980 ± 50 | 5315 ± 150 | b (NTU-2644) |
Shanhaili | 12 | 11.5 | 4950 ± 125 | 5300 ± 360 | a |
Fantzuliao | 1 | 1 | 1300 ± 120 | 885 ± 240 | a |
KT1 | 3 | 17 | 6670 ± 50 | 7225 ± 90 | b (NTU-2640) |
KT2 | 6 | 5.5 | 5070 ± 40 | 5450 ± 140 | b (NTU-2641) |
KT3 | 6 | 6.5 | 5660 ± 50 | 6090 ± 155 | b (NTU-2642) |
KT4 | 2 | 2 | 1820 ± 40 | 1370 ± 110 | b (NTU-2649) |
1 Elevations rounded to nearest 0.5 m; calibrated ages expressed as ± maximum error at 2 s.d. Glacioisostatic correction after W.R. Peltier (see text).
2 Age calibration using CALIB Rev. 4.1.2 (1999): see [13]. No correction for fractionation as δ13C not determined. For marine reservoir effects, see text. The regional average for the western sub-tropical Pacific is [13].
3 a: [8] and sources therein; b: National Taiwan University radiocarbon laboratory.
4 Seismic hazard
An earlier attempt to identify zones of high potential seismicity on the basis of the neotectonic record proposed that the western fold-and-thrust belt was most at risk primarily because the geometry of its major fault structures was propitious both for strain accumulation and its seismic release [15]. The reasoning was that folding was driven by blind faults until their progressive rotation and steepening at the rear of the trench ruled out further slip. Support for the model comes from nascent folding in the coastal plain west of the fold-and-thrust belt [4,12].
The new data refine the scheme by showing that deformation in the Central Range, or at any rate in its southern part, is partitioned between general uplift and localised folding. A plausible interpretation is that the former represents energy storage and the latter its release. To judge from the published cross-sections [2], the Hengchun fault in places dips E but rather than being inactive the structure serves as the locus of folding. Some seismicity is doubtless promoted by bedding slip but any residual stored energy is dissipated seismically in low angle reverse faults to the west.
This interpretation conflicts with the outcome of probabilistic analyses. Tsai et al. [14], for example, had concluded that the highest seismic risk was in the Northeast and East of Taiwan and that the lowest level of risk was displayed by the central range, the fold-and-thrust belt, and the North and South, with an intermediate level near the active faults in the western coastal plain and in the Southwest of the island. Needless to say the severity of risk within each zone depends on many additional factors. Thus Wu and Rau [18] have argued that blind thrusts in the fold-and-thrust hills and the coastal plain are especially dangerous because they occur under populated urban areas; uplands prone to mass movement will evidently be more sensitive to uplift than lowlands.
Reliance on neotectonics in no way minimises the value of seismic analysis or geodesy; for, by identifying active zones, subsurface structures and fault polarity, these sources of data are an essential component of the assessment. What is at issue is the governing time scale. In many areas the recurrence interval between major events is measured in millennia, so that the historical evidence may lead one to underestimate the potential level of seismicity; on the other hand, geodetic measurements made over a few years give values for deformation which, though geometrically correct, can exaggerate the prevailing strain rate. The Hengchun Peninsula illustrates a third benefit of the long-term perspective: strain energy [3] may be stored across the entire accretionary prism but seismically released only where the available slip surfaces are favourably oriented and inclined. In sum, the proposed approach enables zones most susceptible to damaging seismicity to be identified and provides a rationale for the zoning of building regulations and emergency provision which is not compromised by inadequate instrumental or geodetic data.
Acknowledgements
C.V.-F. is am grateful to the Royal Society of London and the National Science Council, Taiwan, for support, Dr. P.J. Vita-Finzi for assistance in the field, Prof. T.K. Liu, National Taiwan University, for radiocarbon assays, and Prof. J. Angelier for comments on a draft of the paper.