years BP. This allows increasing the basic geological knowledge necessary for determining the seismic hazard associated with cortical structures in the Aysén region in southern Chile. Fault cartography and field observations suggest dominant dextral-reverse strike slip along north-south and locally NNW-striking faults, dextral-normal strike slip along NE to NNE- striking faults, and sinistral strike slip along east-west faults. This kinematics is consistent with regional SW-NE shortening in the context of a major transpressional fault zone. Holocene and even historic monogenetic and sub-aquatic volcanism occurred in this tectonic setting in a close spatial relationship and probably favored by the activity and local architecture of faults.

Submarine fault scarps and deformed sediments observed at the bottom of the Aysén Fjord were associated with the destructive April 2007 Mw6.2 earthquake located along the LOFZ. Our observations show that this earthquake occurred along dextral 15-20 km long N-S structure named Punta Cola Fault (PCF). This fault system is located some kilometres to the east of the main N-S Río Cuervo Fault (RCF). Most of the epicentres of the seismic swarm during 2007 were located along or in between both structures. The study area is a transference zone between N-S regional branches of the LOFZ. The cartography of fault segments proposed here together with geophysical and geologic data suggest that large earthquakes Mw6.2-6.5 can be typically expected along most of the active faults.

Be- sides, seismic hazard assessment should also consider the possibility of earthquake magnitude in the order of 7.1 along the main fault systems like the RCF.

The close relationship between Quaternary faults and volcanic activity along this fault zone (Lara et al., 2008; Cembrano and Lara, 2009), suggest that the circulation of fluids with a possible magmatic component can plays a role in the generation of earthquakes along this fault zone (Legrand et al., 2011). Here, we assess the cartography and kinematics of Late Pleistocene- Holocene fault systems and spatially associated monogenetic volcanic centres in the Aysén Fjord area of southern Chile, with the aim to improve critic knowledge for seismic hazard assessment in the Patagonian Andes.

The Liquiñe-Ofqui fault zone (LOFZ) is a major trench parallel intra-arc fault system which is recog- nized for more than 1,200 km along southern Chile, between 38°S and the triple junction of the Nazca, South American and Antarctic plates at ca. 47°S (Hervé, 1976, 1994; Forsythe and Nelson, 1985; Cembrano et al., 1996, 2002; Rosenau et al., 2006; Fig. 1). The right lateral strike-slip, oblique slip, reverse slip ductile to brittle deformation charac- terizing a transpressional regime along this system has been active during the late Cenozoic (Cembrano et al., 1996, 2002; Lavenu and Cembrano, 1999; Arancibia et al., 1999; Thomson, 2002; Rosenau et al., 2006; Adriasola et al., 2006), associated to

143 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 the oblique northeasterly (N77ºE) subduction of the Nazca plate beneath the South American plate, with an average of 66 mm/yr at present (Angermann et al., 1999). South of 46-47ºS, the convergence margin is characterized by orthogonal subduction of the Antarctic plate beneath the South American plate at an average rate of 20 mm/yr (DeMets et al., 1994). Besides oblique subduction contact, strong inter- plate coupling, thermally weak continental crust, and FIG. 1. Regional tectonic context at the triple junction of the Nazca, Antarctic and SouthAmerican plates, and location of the 2007 seismic swarm that occurred in the Patagonian Andes.

White lines show the approximate maximum extend of the Pleistocene ice sheet during the Last Glacial Maximum (modified from Thomson et al., 2010). Right: Location of the April 2007 seismic swarm in the area of the Aysén Fjord (Legrand et al., 2011), and the relief shown by the means of a SPOT satellite image over a Digital Elevation Model (DEM). Circles represent seismic events of Mw ≤3.0, as well as some events of Ms>5. The red star indicates the location of the main shock (Mw 6.2), according to GUC (2007). The location of the Macá, Cay and Hudson volcanoes is also shown. Left bottom: Distribution of number of events versus time, and focal mechanisms for the main shocks according to GUC (2007) and HRV (2007).

While the Concepción and Maule earthquakes ruptured for almost 500 km along the subduction contact located to the north of the Arauco Peninsula (Fig. 1), the seismic rupture of the Valdivia earthquake propagated for about 1,000 km along the subduction margin from the southern part of this peninsula (38ºS) towards the triple plate junction, which coincides with the latitudinal extent of the LOFZ (Plafker and Savage, 1970; Barrientos and Ward, 1990; Melnick et al., 2009). Based on geological, seismological and geodetic data, the LOFZ has been considered the eastern boundary of a fore-arc sliver characterized by northward displacement and decoupled from stable South America (Wang et al., 2007; Melnick et al., 2009), with a northern limit that should be a buttress located at the latitude ofArauco Peninsula (37ºS, Fig.

The Quaternary volca- nism along the SVZ between 37° and 46°S is mainly basaltic to andesitic, with some explosive products of more silicic composition (López et al., 1995; Cembrano and Lara, 2009). This kind of volcanism is closely associated to NNE and NE-striking tension cracks with primitive minor Holocene monogenetic eruptive cones aligned along the master faults of the LOFZ (Lara et al., 2008; Cembrano and Lara, 2009), like the historically active Llaima volcano which developed along a NE fissural range with almost 40 minor centres aligned along a NE-striking arc structure (Moreno and Naranjo, 2003; Naranjo and Moreno, 2005).According to Lara et al.

(2006) and Cembrano and Lara (2009), dacitic to rhyolitic volcanism can result from more complex evolution in some cases associatedwithWNW-strikingbasementstructureslike the Puyehue-Cordón Caulle, which is a NW-striking central volcano and fissure system associated to the LOFZ that experienced large explosive eruptions on 1960 and since 2011. The recent tectonic history of the Patagonian Andes includes destructive crustal earthquakes like the 2007 seismic swarm occurred at 45.5°S, which was closely associated with master faults of the LOFZ (Barrientos et al., 2007; Comte et al., 2007; Cembrano et al., 2007; Mora et al., 2008; Legrand et al., 2011).

Since January 2007 and during subsequent months, thousands of shallow earthquakes with focus located up to 9-10 km depth were recorded by local seismic network deployed in the Aysén Fjord and surrounding region (GUC, 2007; Barrientos et al., 2007; Legrand et al., 2011; Fig. 1). After an earthquake Mw5.3 occurred on January 23rd 2007, the greatest magnitudes increased through time (5.3, 5.2, 5.3, 5.7, 6.1, 6.2) within the following three months (Legrand et al., 2011). The last Mw 6.2 earthquake triggered rock slides of ca. 17-22, 7 and 2 million cubic metres volume, which occurred in a ravine next to Punta Cola and in front of Isla Mentirosa on the northern shore of the fjord (Fig.

2), and in front of Punta Cola on the southern shore of the fjord, respectively (Clavero et al., 2008; Sepúlveda and Serey, 2009; Naranjo et al., 2009; Sepúlveda et al., 2010; Yugsi Molina et al., 2012). While the first reached the fjord as a large debris avalanche, the second and third discharged directly into the basin, locally generating large tsunami waves that produced casualties and damage to salmon farms (Sepúlveda and Serey, 2009; Naranjo et al., 2009). Given the timing, magnitude and energy distribu- tion of seismicity, together with the focal mechanism of the six largest earthquakes, mainly NS to NNE dextral strike-slip with one NE normal slip (GCMT catalogue, 2007), this seismic swarm was interpreted as the result of the movement of fluids that at the same time reactivated strike-slip and normal faults

145 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 of the LOFZ (Legrand et al., 2011). The inferred magmatic component as one of the possibilities for the origin of those fluids should be compatible with the presence of many monogenetic volcanic cones and complex Quaternary volcanoes in the area (Legrand et al., 2011), including the fjord itself (Lara, 2009). In this paper we show submarine and subaerial evidence of co-seismic surface rupture associated to the last seismic swarm and previous events in the Aysén Fjord, and we provide geomorphologic, structural and geochronological evidence for a close spatial relationship between Holocene volcanic and tectonic activity along the LOFZ.

Through detailed surface mapping of postglacial geomorphologic features affected by faults, together with field ob- servations of mesoscale outcrops of faults affecting rocks and sediments, we propose a cartography of active fault segments and a kinematic model for recent tectonics, providing new data for seismic hazard assessment in one of the most important urban and commercial areas in the Aysén region of southern Chile.

2. Methods and source of data Using a 3.5 kHz chirp-system sub-bottom profiler (Bathy, 2010), we obtained bathymetric and high-resolution seismic reflection data from the upper ~50-200 m of the sedimentary infill from the Aysén Fjord and surrounding lakes, near Puerto Aysén (Fig. 3). These techniques provides information about the stratigraphy of sediments that forms the subsurface of the sea bottom and about the structures that affect the sediments, at vertical resolution that, in terms of distance and in favourable cases such as laminated fine sediments, can be of 0.5 m and even as fine as 8 cm.

I94- I96, I105-I107, I117- I119, and I129-I131). DEMs obtained from the Shuttle Radar Topographic Mission dataset were also used to complete the observations concerning the main geomorphologic features in the area. Fieldwork and aerial recognition including the study in the field of the accessible outcrops complemented this analysis. The most important faults were determined through the cartography of recent tectonic forms, like lineaments and fault scarps affecting glacial or postglacial geomorphologic units, complemented with structural observations and local kinematic measurements, as well as with observations from sub-bottom profile data from the Aysén Fjord and FIG.

2. SPOT satellite image and photographs showing the area where the largest rockslides occurred during the main earthquake on April 21, 2007, as well as the impact of the tsunami waves generated during this event.

146 Submarine earthquake rupture, active faulting and volcanism along the major Liquiñe-... surrounding lakes. In particular, the timing for the last Pleistocene ice sheet retreat in this area of the Patagonian Andes has been poorly assessed. In this work we show evidence that provide a temporal framework for the assessment of postglacial features affected by faults. Because of the vegetation cover and due to the strong erosional and depositional dynamics associated with the humid climate of the western slope of the PatagonianAndes (3,000-4,000 mm of annual rainfall), the preservation and hence the observation of outcrops associated with recent fault scarps and other geomorphologic features affected by tectonics is generally difficult.

Therefore, the fault cartography proposed in this work provides a new and detailed representation of fault segments constituting the main structural systems, showing evidences of superficial ruptures during the latest Pleistocene and Holocene times in the area. The differentiation between observed and inferred fault segments that resulted from this geological- geomorphic mapping considers subjective criteria with respect to the degree and nature of observations (from imagery, field outcrops and indirect observations from subbottom profiles) that supports the proposed fault cartography at scale 1:50,000.

The geomorphology of monogenetic volcanic cones was analysed following the morphometric methods proposed by Tibaldi (1995), with the aim to test their association with recent tectonic featu- res in the area. Together with their alignment, the morphometry of pyroclastic cones can be used to infer the geometry of fracture feeding of magma to the cones (Tibaldi, 1995; Corazzato and Tibaldi, 2006). The base length and width together with axis orientations of the cones were mapped from FIG. 3. Location of the sub-bottom seismic profiles obtained from the Aysén Fjord and Riesco Lake, using a 3.5 KHz (Bathy, 2010) device in the area of the 2007 seismic swarm (GUC, 2007), and subaerial and submarine reliefs shown by the means of a Digital Elevation Model, based on 1:50,000 topographic sheets from the Chilean Army Geographic Institute (IGM) and bathymetric information from the Chilean Navy Hydrographic and Oceanographic Survey (SHOA).

The figure shows the location of the epicentres of the seismic events (yellow circles) and the epicentre of the Mw 6.2 earthquake (red star), according GUC (2007) and Legrand et al. (2011). Selected segments from profiles F1, F8 and R1 are shown in figures 4, 5 and 11, respectively.

147 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 SPOT satellite images and DEM, from which the elongation or degree of ellipticity using the ratio between the length of minor and major cone base axes was obtained. This was done for a total of 40 monogenetic cones, both as a whole and separated by geographically distributed groups. 3. Submarine evidence of seismic rupture in the Aysén Fjord Subbottom profiles evidenced different seismic facies, which can be related to different geomorpho- logic and stratigraphic context. Massive to poorly laminated seismic reflection facies constitute most of the sedimentary infill of marine and lacustrine basins in the area, while well-stratified seismic reflection facies were observed close to the main riverine dis- charges, like theAysén River and Cuervo River (Figs.

4 and 5). Chaotic seismic reflection facies were associated with moraine and landslide deposits, and stratified deformed sediments were observed near Punta Cola, to the east of the Cuervo River (Fig. 4). From the analysis of seismic reflection data we defined two units that forms the upper sedimentary infill of the fjord bottom (Figs. 4 and 5): the Lower Unit (of >50 m thick) constituted by massive or poorly stratified seismic facies and the Upper Unit (of ca. 25 m thickness) formed by well-stratified and locally deformedsediments(Fig.4).Considering that massive seismic facies can be associated with massive sandy sediments that constitute for example fluvial terraces of theAysén River at PuertoAysén, we interpret that the Lower Unit is most probably composed of sandy facies, while the Upper Unit of alternating sandy and FIG.

4. Selected segments from sub-bottom profiles F1 and F8 from the Aysén Fjord (location in Fig. 3), showing fresh submarine fault scarps and deformed sediments in the area of the Mw 6.2 earthquake, along the Punta Cola Fault (Fig. 6). Massive and stratified sediments correspond to the Lower and Upper units respectively, of the sedimentary infill of the Aysén Fjord (see text).

148 Submarine earthquake rupture, active faulting and volcanism along the major Liquiñe-... fine facies. Following sequential stratigraphy, we interpret that the sedimentation of the upper unit is probably associated to alternate supply of material from the marine basin and fluvial sediments from the Aysén River and other drainage systems in the area, most probably during the Holocene global high sea level stage since 7 cal. ka BP(Lambeck et al., 2002). Fresh submarine fault scarps disrupting the sea bottom with vertical offsets of around 6-8 m were observed in subbottom profiles in the area of the 2007 seismic swarm in the Aysén Fjord (Fig.

4). Seismic profiles obtained from that area revealed strong deformation affecting sediments of the Upper and Lower seismic units (Fig. 4). This deformation is characterized by folded strata and loss of stra- tification, suggesting soft-sediment deformation which is particularly evident in the case of the Upper Unit. The area is affected and limited by faults with reverse slip component, which together with the deformed units define a positive flower structure (Fig. 4: Profiles F1 and F8). The most prominent fault, which is located at the SE border of this area, together with secondary faults defines a local pop-up structure at the sea floor, affecting both the upper and the lower units (Fig.

4). At its western margin, relatively higher complexity seems to be related to westward propagating low-angle thrust faults along the contact between the massive-to poorly laminated- Lower Unit and the -well stratified- Upper Unit, as well as separating non-deformed and deformed sediments of this last unit (Fig. 4). Other minor scarps can also be observed in the middle of the rupture zone. The geometry of the entire submarine fault rupture zone defines a flower structure of ca. 2.5 km width measured along profile, with higher elevations at the SE and NW extremities (Fig. 4). The absence of any sedimentary layer covering fault scarps or disrupted sediments in this fjord, which is characterized by large supply of sediments from the Aysén River (Fig.

3), provides a strong argument to support that this submarine rupture zone is associated with very recent seismic events and possibly in part with the last seismic events, including the one of April 2007.

The sub-bottom profiles provide evidence of previous activities along faults affecting the sedi- mentary infill of the fjord (Fig. 5). In most of the cases, faults with reverse slip component affect the FIG. 5. Selected segment from sub-bottom profile F8 from the Aysén Fjord (location in Fig. 3), showing part of a submarine ridge located close to the Río Cuervo Fault (Fig. 6), which corresponds to the main structure of the Liquiñe-Ofqui Fault Zone in the area. This submarine ridge is associated with a bathymetric down step of ca. 100 m from the eastern part towards the western part of the fjord.

It is interpreted as a submarine barrier, which in part limits the sedimentation of the most recent and well- stratified sediments constituting the Upper Unit of the sedimentary infill of the Aysén Fjord to the east (see text).

149 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 underlying poorly stratified sediments of the lower unit, but not the entire overlying stratified upper unit, suggesting previous seismic ruptures in the area (Fig. 5). Together with complementary bathymetric information from SHOA, the sub-bottom profile data shows a submarine ridge affected by faults and elongated in a SW-NE direction across the fjord in the area of the 2007 seismic swarm (Figs. 3 and 5; Profile F8). This ridge is associated with a bathy- metric difference of about 100 m in the fjord bottom from 210 m to 310 m water depth from the eastern side towards its western side, respectively (Fig.

5). A submarine conic structure of about 130 m height that exhibits internal stratification and constitutes part of the ridge (Fig. 5) can be interpreted as a volcanic centre, following the inference of Lara (2009) from dark vesicular basalts dredged between fine sediments in the same area. The sub-bottom profile data show that the ridge constitutes the western limit for the upper stratified unit of the sedimentary infill in the fjord, and probably acts like a barrier for sediments supplied from the Aysén River.

4. Fault cartography and kinematics Recent fault scarps and lineaments are largely disposed along N-S, SW-NE to SSW- NNE and E-W directions, which coincide with the more frequent directions of the valleys in the area, as in the case of the main lacustrine basins of Yulton and Meullín lakes (Fig. 6). Valleys oriented NW-SE are more frequent in the eastern part of the study area, which are generally affected by fault scarps and lineaments disposed along N-S, SW-NE and E-W directions. Besides, they also exhibit high degree of glacial erosion like U-shaped valleys and glacial striations (Fig. 7).

Observations realized from outcrops located some kilometres to the east of PuertoAysén evidenced outwash and fluvial sediments directly overlaid by volcanic ash fall deposits including buried soils at the base, from which we obtained calibrated radiocarbon ages of 10,350 and 12,000 cal. years BP (Table 1, Fig. 7). These results provide maximum ages for latest Pleistocene-Holocene development of fluvial sedimentary infill in the Aysén River valley, after the last ice sheet retreat in the area, which occurred shortly before 12 cal. ka BP.

The cartography of faults in the area, together with local observations from outcrops and kinematic indicators (Table 2), suggests dextral strike slip as the dominant mechanism for recent deformation, which dominates along the N-S and SW-NE oriented faults (Table 2; Fig. 6, Río Cuervo Fault, Punta Cola Fault) and locally along NNW-SSE directions, especially in the case of the north-eastern part of the study area (Table 2; Fig. 6, Los Palos Fault). Some east-west oriented faults exhibit kinematics indicators of sinistral strike slip mechanism (Table 2; Fig. 6).As shown in TABLE 1. CONVENTIONAL AND CALIBRATED RADIOCARBON RESULTS (CAL.

150 Submarine earthquake rupture, active faulting and volcanism along the major Liquiñe-... the geological maps (Fig.6), the major structures of the different systems displace each other. In some cases, EW-trending faults displace NS faults, which in turn displace or curve the EW faults. This would suggest that all structures are active and responding to the same stress state, possibly with differences of slip rate.

In general, east-west oriented faults are deflected towards the south in the eastern blocks along north-south oriented fault systems, probably reflecting faster ate Quaternary slip rates along the later systems with respect to the first (Fig. 6).Addi- tionally, a reverse component is associated with the dominantly dextral slip along N-S and NNW-SSE faults, while dextral-normal slip can be inferred from the fault geometry and from field observations along NE-SW fault planes (Fig. 6).According to the cartography, major structures have the segments as much as 40 km long, while secondary faults have those in the order of 10-20 km (Fig.

6). The most prominent north-south fault is he Río Cuervo Fault (RCF), which can be characterized by fault segments up to 30-40 km long (Fig. 6). The main fault trace is located along the western side of Cuervo River valley, crossing theAysén Fjord close to the location of the submarine ridge (Figs. 1, 5 and 6). This fault system can be correlated with a major north-south branch of the regional LOFZ (Hervé, 1994; Cembrano et al., 1996, 2002; Thompson, 2002; Fig. 6). Intense brittle shear deformation affecting granitoids along the RCF defines hundred metres width damage zones, visible from the Aysén Fjord (Table 2).

Outcrops along this fault zone show subvertical faults dipping 60°-90°E, metric fault breccias, as well as mm to cm shear zones, with striations altogether with secondary minerals like iron oxides and incipiently chlorite and epidote in fault planes (Table 2).

The Punta Cola Fault (PCF) is a north-south striking 15-20 km long structure that crosses the Aysén Fjord about 5 km to the east of the RCF (Figs. 6 and 8). Rocks affected by intense brittle deformation including fault planes and fractures, breccia and centimetre shear zones were observed along the several hundred metres width damage zones at the northern side of the fjord (Fig. 8; Table 2). The location of the major Mw 6.2April 2007 earthquake given by the local seismic network (GUC, 2007; Legrand et al., 2011; Figs. 1 and 3), coincides with the location of this north-south fault (Fig. 6), and with the location of the submarine rupture observed in sub-bottom profiles (Fig.

4). Large landslides triggered by this main shock occurred close to this fault in the fjord on 2007 (Fig. 2). In particular, the largest landslide that occurred in the ravine next to Punta Cola (Clavero et al., 2008; Naranjo et al., 2009; Sepúlveda and Serey, 2009; Sepúlveda et al., 2010), was limited on its eastern side by north-south 80-85°W dipping faults associated with the PCF (Figs. 6 and 8).

Sub-bottom profiles and surface features indi- cate that several SW-NE striking secondary faults, with segments up to 10 km in length, are located between the RCF and the PCF (Fig. 6). Besides, deformed submarine sediments and/or surface manifestations of postglacial fault activity were observed in seismic profiles from the fjord bottom, that we associate with SW-NE structures traversing the fjord between the RCF and PCF (Fig. 6). This area corresponds to the zone where the epicentres of hundreds of seismic events (Mw ≤3.0) associated with the last seismic swarm were localized (GUC, 2007; Legrand et al., 2011; Fig.

3). In the area of the RCF deltaic sediments from the Cuervo River are interbedded with coastal marine sediments and form a terrace located up to 8-10 m above the fjord sea level. AMS radiocarbon results from levels with abundant tree remnants provided calibrated ages between 7750 and 8500 cal. years BP (Table 1; Fig. 9). This suggests several metres of surface uplift since the early Holocene (Fig. 9). These se- diments are affected by NE-striking faults and are overlain by gravel deposits, suggesting previous surface seismic rupture associated with secondary faults between the RCF and PCF (Table 2; Fig.

9). Two other north-south fault systems are the Punta Mano-Lago Yulton Fault (PMLYF) and Los Palos Fault, near Puerto Aysén (LPF; Fig. 6). Outcrops from LPF show fault breccias and mm to cm shear zones associated with subvertical N-S and NNW- SSE fault planes dipping 50-80°, with striations and secondary minerals evidencing dominant dextral strike slip mechanism (Table 2; Fig. 10). In the area of Puerto Aysén, a fluvial terrace of the Aysén River occurs up to 8-10 m over the thalweg of the river. In this area the altitude of the Puerto Aysén fluvial terrace increases downstream from the east towards the centre of the city, suggesting a tectonic control in its geometry that, besides fluvial dynamics, could be associated to recent activity of Los Palos Fault.

Sub-bottom profile data from the Riesco Lake, which is located 6 km to the south of Puerto Aysén, evidenced recent deformation

The number of observations from local measurements of kinematic indicators is shown together with each local fault plane solution obtained using Faultkin (Marret and Allmendinger, 1990; Allmendinger et al., 2012). The fault plane solution from the Mw 6.2 earthquake of 2007 is also shown (HRV, CMT, 2007). The cartography of the monogenetic centres, which are mostly clustered along the main NE to NNE-trending structures, is provided. Morphometric parameters from monogenetic centres (shown in Fig. 15), were characterized from areas A, B, C, D and E, which correspond to Cay, Macá, Bahía Pérez, Meullín and Pescado River zones, respectively.

B. Geomorphologic chart of the study area located to the north of the Aysén Fjord.

153 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 affecting postglacial massive -probably sandy flu- vial dominated- sediments of this lacustrine basin, dammed by Pleistocene moraines to the west. This deformation can be related with a reverse slip com- ponent along a north-south dextral fault segment of LPF (Fig. 11), which supports the active character of this fault system. A major SW-NE striking fault system is located in the southern part of the area, which is the Riesco Fault (RF, Fig. 6). This last fault (RF) coincides with the main NE-SW branch of the LOFZ, which has been named Mañihuales Fault to the northeast of the study area (Thompson, 2002; Cembrano et al., 2002).

Results from SPOT satellite imagery may be interpreted as showing systematic right- lateral displacement of recent drain-age systems along the RF (Fig. 12), which could be taken as preliminary evidence for Holocene strike slip deformation along the LOFZ in the area. Local outcrops along this fault system reveal subvertical fault planes with mm to cm shear zones, together with striations and secondary minerals in fault planes evidencing dominant dextral strike slip fault kinematics associated with structures oriented N-S to NNE-SSW (Table 2; Fig. 12).

Another NE-SW main branch of the LOFZ in the region is the Quitralco Fault (QF; Cembrano et al., 1996, 2002; Thompson, 2002). This fault system joints the RCF on the southern shore of the fjord 2-3 km to the SE of the submarine ridge observed in sub-bottom profiles, within an area where several monogenetic volcanic centres are located (Fig. 6). 5. Monogenetic volcanic centres The geological analysis revealed that dozens of well preserved monogenetic volcanic centres FIG. 7. A, B. Outcrops showing outwash and fluvial sediments directly overlaid by volcanic ash fall deposits, located some kilometres to the east of PuertoAysén, and calibrated radiocarbon results from buried soils (Table 1); C.

Glacial striations affecting bedrocks in the study region. These observations point to a last Pleistocene ice sheet retreat before 12,000 cal. years BP in the area.

156 Submarine earthquake rupture, active faulting and volcanism along the major Liquiñe-... FIG. 8. Photographs showing outcrops associated to Punta Cola and Río Cuervo faults, to the north of the Aysén Fjord (Fig.

6a). A, B. North-south fault planes and fault scarps limiting a large rockslide deposit formed during the last Mw 6.2 earthquake along the area of the Punta Cola Fault; C. Rockslide deposit associated to a NNE-striking fault in the area of the Río Cuervo Fault; D. Local kinematic observations from fault planes of the RCF and PCF (Table 2).

157 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 are clustered along main faults oriented following NE-SW, NNE-SSW and N-S directions (Fig. 6). A number of monogenetic volcanic centres are located close to the Pleistocene-Holocene basaltic to dacitic Macá and Cay stratovolcanoes (Stern et al., 1976; López et al., 1993; D’Orazio et al., 2003), some of them constituting possible parasite structures of the former volcano. Other monogenetic volcanic centres are located in the area of the Cuervo River close to the Meullín Lake and on the southern shore of the Aysén Fjord along the Pescado River valley (Demant et al., 1994; Lahsen et al., 1997; Fig.

13). Minor eruptive centres are characterized by hundreds of metres high conic structures with diameter up to ca. 2 km. Most of them are elongated and exhibit variable degree of erosion and vegetation cover, suggesting that eruptive activities occurred throughout the Holo- cene period. The most recent minor volcanic centres are located in the western part of the study area, at Bahía Pérez (Puerto Pérez area; Fig. 6), southwest of the Macá Volcano (Fig. 13a-d). Radiocarbon data from carbonized tree remnants under a basaltic lava flow from one of these centres yielded ages of 386±106 and 482±54 cal.

yr BP, which means that this centre erupted between the years AD 1450 and 1670, indicating its historical age (Table 1, Figs. 13c FIG. 9. Holocene littoral and deltaic sediments from the Cuervo River affected by faults on the northern shore of the Aysén Fjord. A, B and C. Views showing deltaic sediments (A, B), affected by faults (C), in the area of the Río Cuervo Fault (RCF), and cali- brated radiocarbon results (Table 1). D. Local kinematic data from faults affecting sediments in the area of the RCF (Table 2).

160 Submarine earthquake rupture, active faulting and volcanism along the major Liquiñe-... and d). Among the most vegetated features are the monogenetic centres located at Pescado River on the southern shore of the Aysén Fjord close to the RCF, and to the south of the submarine ridge. This area is characterized by major north-south oriented escarpments and lineaments that we interpret as the continuation of the dextral Río Cuervo Fault system towards the south of the fjord, with associated se- condary NE-SW faults inferred from SPOT satellite imagery (Figs. 13e and f).

The geomorphological analysis of the minor volcanic centres confirmed that the orientations of the major axes of the monogenetic cones coincide with those of the master faults of the LOFZ, with preferential strikes of NNE to NE (Figs.

14 and 15). Nonetheless, within each monogenetic group clus- tered in different areas it is possible to find several orientations, suggesting a local association with the faults. The degree of ellipticity of the monogenetic cones is moderate according to Tibaldi (1995), with an average ratio of minimum to maximum axis length of 0.80±0.13. The strongest ellipticity is detected in a monogenetic cone of the Meullín group close to the RCF system, with a ratio of 0.45, although in this zone there is a great dispersion ranging from circular to elongated cones (Fig. 14). As a group, the parasitic cones of Macá volcano have the hig- her average ellipticity with a degree of 0.71±0.10, while the others vary between 0.81 and 0.85 (Figs.

6 and 15).

The sub-bottom profiles, together with bathymetric information from theAysén Fjord, show that the sub- marine ridge is elongated in a NE-SW direction and that is located close to the junction of the NE-striking Quitralco Fault (QF) and the north-south-striking RCF (Figs. 5 and 6). The stratified conic structure forming part of this ridge is also elongated in a NE- SW direction, and is characterized by a diameter of more than 500 m and a height exceeding 100 m from its base at the eastern side of the ridge, reaching minimum water depths in the order of 60-70 m b.s.l. (Fig. 5). The internal stratification and the dimension of this submarine structure forming part of the ridge is similar to other volcanic monogenetic centres in the area, like those located at Pescado River (Figs.

5 and 13e and f), which supports the interpretation of its volcanic origin.

Sub-bottom profile and bathymetric data obtai- ned from the north-eastern side of the Yulton Lake revealed a conic structure that may be interpreted as a subaquatic volcanic centre located at the base of the Cay volcano (Fig. 16). In fact, this is a conic structure elongated in a SW-NE direction, which exhibits internal stratification that suggests a volcanic origin under shallow subaquatic conditions (Fig. 16). On its south-eastern border, this feature is affected by faults with normal slip component suggesting post- glacial deformation under local extensional regime, which can be related to the SW-NE-striking Lago Yulton Fault (LYF; Figs.

6 and 16), inferred from bathymetric,sub-bottom profile data and local field observations on the northern shore of theYulton Lake (Fig. 6; Table 2). Faults with reverse slip component were deduced from fault scarps disrupting the lake bottom on the north western border of this feature that can be associated to north-south striking faults interpreted from sub-bottom profiles and bathymetric information. Such an interpretation is consistent with fault kinematics deduced from the fault cartography, locally supported by kinematic indicators from the entire study area and especially from the northern shore of the Yulton Lake (Figs.

6 and 16). 6. Discussion The cartography of active faults indicates that recent deformation in the area occurs mainly along dextral N-S and NE-SW faults. E-W striking faults accommodate sinistral strike slip. The geological evidences show local reverse slip component asso- ciated with N-S and locally NNW-striking faults, which favours the generation of pop-up structures, as suggested by submarine evidences of the last seismic ruptures in the Aysén Fjord bottom, as well as by topographic and structural data (Fig. 6). According to the fault cartography and local field data, NE- SW secondary faults accommodate dextral-normal slip.

We interpret that the activity of the main N-S faults together with NE-SW faults (Fig. 16), would have favored the development of extensional duplex structures which together with strong glacial erosion have driven the formation of the Yulton and Meullín lake basins.

161 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 FIG. 13. SPOT satellite images and photographs of some of the monogenetic volcanoes in the study area. A-C. show the location of monogenetic centres in the area of Bahía Pérez, southwest from the Macá Volcano.

C and D show the location of carbonized tree remains under a basaltic lava flow, associated with a poorly vegetated monogenetic centre dated by radiocarbon analyses. Calibrated radiocarbon results evidence the historical age of this centre (Table 1). E and F show monogenetic centres located at Pescado River, near the southern shoreline of the Aysén Fjord, in the area of the Río Cuervo Fault (RCF) and close to a volcanic conic structure within a submarine ridge shown in sub-bottom profiles (see text). White arrows represent inferred kinematics along the faults.

162 Submarine earthquake rupture, active faulting and volcanism along the major Liquiñe-... tions obtained from faults measurements in rocks and sediments (Figs. 8, 9 and 10), together with the striations affecting secondary minerals in some fault planes, suggest that mesoscale fault outcrops are exhumed fault systems that continuous active under similar tectonic regime at the scale of Late Pleistocene-Holocene. The overall data is consis- tent with a dominant transpressional regime driven by dominant NE-SW shortening during the Late Pleistocene and Holocene, as suggested previously for the late Cenozoic in the region (Hervé, 1994; Cembrano et al., 1999, 2002; Rosenau et al., 2006), which includes the development of duplex basins, crustal-scale pop-up structures, denudation and ex- humation of granitoids along the LOFZ (Thompson, 2002).

According to Lara et al. (2008), the coeval transpression or compression associated both with a NE-striking σ1 along this major intra-arc fault zone can explain in part the vertical growth of the relief and associated volcanism. In the study area the different pattern that characterizes the positive relief along LPF, with respect to the geomorphology and associated volcanic centres along the RCF-LMF and LYF, could represent local differences in the context of a dominant transpression regime. Particularly, exhumation and fluid circulation along the main structures of the LOFZ seem to play a significant role in faulting mechanisms, as suggested by secondary minerals observed in fault planes in the area (Table 2).The focal mechanisms with N-S and NE-SWplane FIG.

163 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 solutions obtained from most of the large shocks (Ms >5) during the last seismic event in the area (HRV, 2007; Legrand et al., 2011; Fig. 1), including the Mw6.2 main shock (Fig. 6), and those deduced from previous seismic surveys performed north of the study area (Lange et al., 2008), can be explained in the context of this fault kinematics. Furthermore, Quaternary neotectonics associated to fault kinema- tics and regional denudation, landslides, fluvial and volcanic geomorphologic processes have overprinted conspicuous glacial landforms, controlling most of the latest Pleistocene-Holocene geomorphologic evolution of the valleys in the study area (Fig.

6). Most of the monogenetic volcanoes occurred clustered along NE to NNE-striking faults as shown by the geological mapping (Fig. 6), which coincides with the maximum stress direction deduced from the fault geometry in the area.Additionally, results from the morphometric analysis of the cones evidenced that most of them exhibit elongated cone bases along NE-striking directions as well as along north-south (Fig. 15). These observations support the idea that this kind of volcanic activity has been partly fa- voured by local extension along secondary faults or by inflections along master faults of the LOFZ, following Nakamura (1977); Lara et al.

(2008) and Cembrano and Lara (2009), as in the case of the Meullín area (Fig. 6 and 14), where, in spite of large dispersion, most of the cones are slightly elongated and aligned along NNE and NE-striking directions. The geochemistry of some of this type of volcanoes in the region, which suggest the unevolved or poorly evolved character of the associated magmatism, could support this interpretation (D’Orazio et al., 2003; Gutiérrez et al., 2005).Additional research must be done to better constrain the possibility of this or other mechanisms for the origin of monogenetic cones in the area, like for example decompression melting along the LOFZ as suggested by Lara et al.

(2008). The different degree of erosion and vegetation cover of the monogenetic centres together with radiocarbon results from carbonized trees that gave historic ages (Table 1, Fig. 13), suggest that this volcanism has been intense throughout the Holocene, after the retreat of the last Late Pleistocene ice sheet in the area. This includes subaquatic centres inferred from the Aysén Fjord and Yulton Lake bottoms. From geomorphologic and radiocarbon data from the western flank of the PatagonianAndes, Denton et al. (1999) estimated last large development of piedmont FIG. 15. Morphometric relationships between cone-length and cone-width from monogenetic centres in the study area, and directions of the main axis of the cone base from those centres.

164 Submarine earthquake rupture, active faulting and volcanism along the major Liquiñe-... glaciers at ca. 14.8 ka BP in the Llanquihue region at 41°S, and Heusser (1990) inferred a sudden retreat of the Chilotan piedmont glaciers at 42-43°S, at ca. 13 ka BP. In the study area, glacial retreat occurred before 12 cal. ka BP (10 ka BP), according to the new radiocarbon and geomorphologic observations concerning volcanicash falldeposits directly overlying fluvial and outwash sediments located to the east of Puerto Aysén, provided in this work. Additionally, subbottomprofilesprovidedindirectobservationsabout thetwounitsthatconstitutethesedimentaryinfillofthe Aysén Fjord (Figs.

4 and 5). We interpret that massive to poorly stratified facies that characterize the Lower Unit can be associated in partly with glacial sediments (till and outwash deposits), but mostly with fluvial postglacial massive sandy sediments produced in a terrestrial environment just following the last ice sheet retreat in the area. The well stratified seismic facies that characterize the Upper Unit can be interpreted as the result of an alternated deposition of deltaic and marine sediments under subaquatic conditions that favoured the development and preservation of that type of deposits in a marine environment.

Thus, since this Upper Unit is most probably associated with the supply of fluvial sediments from the main rivers in the area, mainly from the Aysén River, following classic sequential stratigraphy, it should represent the progradation of these sediments within the marine fjord basin during the Holocene global sea level highstand since 7 cal. ka BP (Lambeck et al., 2002). Considering that the submarine ridge in the Aysén Fjord acts as a barrier for sediments supplied by the Aysén River (Figs. 5 and 6), given that the strata that constitute the Upper Unit are onlapping the ridge (Fig. 5), and taking into account that the conic volcanic structure forming part of this ridge FIG.

16. Sub-bottom profile and bathymetric data from the Yulton Lake at the base of the SE flank of the Cay Volcano. The morphology and the internal structure of this feature are interpreted as the results of shallow subaquatic volcanism of possible fissure origin, flanked to the east by normal faults along the inferred Lago Yulton Fault system (LYF; Fig. 6).

165 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 have not been affected by glacial erosion (Fig. 5), we interpret a postglacial-pre sea level highstand age for the formation of this submarine volcanic feature, possibly between 12 cal. ka BP and 7 cal. ka BP.Atectonic control on its origin is a reasonable possibility considering the close spatial association of this feature and other postglacial monogenetic cones with master faults of the LOFZ (Lara et al., 2008). Additional research is needed to better constrain the shape, nature and origin of this submarine feature, which is critic for the assessment of seismic and volcanic hazard in the Aysén Fjord.

Several strata within deltaic sequences in the Cuervo River area forming a terrace located up to 8-10 m above the present-day sea level, with abundant remnants of trees and soils dated ca. 7750-8500 cal. years BP (Fig. 9), suggest episodes of marine flooding concomitantly with a rising shoreline close to the maximum global sea level high stand during the Holocene ca. 7 cal. ka BP (Lambeck et al., 2002). These observations allow us to estimate a minimum uplift rate of 1.1- 1.4 mm/yr during the last 7 ka. This rate is close to the average long-term denudation rates of

166 Submarine earthquake rupture, active faulting and volcanism along the major Liquiñe-... 2011), which are in the order of those deduced from previous seismic survey in the northern part of the LOFZ (Lange et al., 2008). These maximum focal depths are consistent with high geothermal gradient and thin upper crust (6-11 km), according to fission track-results from Miocene magmatic rocks (10 Ma) in the area of the Aysén Fjord (Thomson, 2002). Assuming 9 km depth for the brittle-seismogenic zone in the area, the Mw 6.2 earthquake occurred in April 2007 is coherent with a surface rupture length in the order of 15 km, which is within the scale of the surface trace-length of the PCF, dipping about 80ºW in the area of theAysén Fjord and to the north (Figs.

6 and 8).

Scale relationships between faults and earthquakes (Scholz, 1982;Wells and Coppersmith, 1994), suggest that the Mw 6.2 seismic event could have entrained surface displacements of only few decimetres. Such a dimension is consistent with the scarcity of evidence of surface rupture that can be clearly associated to this earthquake in the study area. In spite of that, fresh fault scarps and fault planes were observed along the PCF, particularly at the southeastern limit of the largest landslide occurred in the ravine next to Punta Cola. Together with fresh submarine fault scarps observed at the fjord bottom in front of that area, these observations can be taken as evidences for surface seismic rupture along this fault system (Figs.

4 and 8). Surface displacements in the order of some decimetres during the last large seismic event appears to be out of scale with respect to 6-8 m high submarine fault scarps observed in sub-bottom profiles in the fjord bottom in front of Punta Cola (Fig. 4).This suggest that those submarine fault scarps observed in the area of the epicentre of the 2007 Mw 6.2 earth- quake accumulates more than one seismic rupture, and thus could represents a cumulative process which includes the last seismic event and those previous that could have been recurrent during the last hundreds or thousands of years.

Some historic chronicles related to one seismic episode occurred in 1927 -shortly after the foundation of the city of Puerto Aysén during the late 1920s- which was probably located at Canal Moraleda located to the west of he study area and that entrained landslides and tsunamis in the area of the Aysén Fjord (Greve, 1964; Naranjo et al., 2009), provide additional information for historic large events along the LOFZ in the region. The fault cartography presented here suggests that major faults are composed of 20- 30 km length active segments that can reach up 40 km length like in the case of the Río Cuervo Fault and Riesco Fault, implicating similar or greater potential seismic hazard in the region.

Considering the possibility for a surface rupture length of 40 km, earthquake in the order of Mw 7.1 could be expected in the area according fault geometry and kinematics of the main structures together with scale relationships between active faults and earthquakes (Wells and Coppersmith, 1994), at period of time of hundreds of years, if we take into account Slemmons and Depolo (1996) and inferences about the order of magnitude of slip rates discussed previously. Besides, most of the faults are composed of segments in the order of 10-20 km length, implicating potential earthquakes in the order of Mw 6.2-6.5 at period of time of dozens of years (Wells and Coppersmith, 1994; Slemmons and Depolo, 1996).

This is particularly the case of Los Palos Fault defined in this work, which traverses the westernlimit of the city of PuertoAysén, from which recent activity can be inferred from sub- bottom profile data (Figs. 6 and 11). Further research focused on paleoseismology, seismic and geodetic survey as well as modelling of fault seismic ruptures is needed to provide more precise estimation of seismic hazard especially for the city of Puerto Aysén in this key area of southern Chile. at can reach up to 40 km length like in the case of the Río.

7. Concluding remarks Because of the steep mountains, dense vegetation cover, high precipitation rates and scarce access, the Aysén region of southern Chile is a difficult zone to develop neotectonic research and especially detailed mapping of active faults. From imagery, DEMs, field recognition, structural characterization, marine geology-techniques and radiocarbon geochrono- logy, we performed geomorphologic analysis that resulted in a proposed cartography of active fault segments, monogenetic volcanic cones and recent kinematics, improving basic knowledge for seismic hazard assessment in the area of the Aysén Fjord, which is crucial for the sustainable development of the Patagonian Andes.

Our observations reveal a strong association between Holocene surface manifestations of active faulting, destructive crustal earthquakes and monoge- netic volcanism along the Liquiñe-Ofqui Fault Zone (LOFZ) in the area of theAysén Fjord. Geomorpho-

167 Vargas et al. / Andean Geology 40 (1): 141-171, 2013 logic observations and radiocarbon results presented here indicate that the last Pleistocene ice sheet retreat in this region occurred before 12000 cal. years BP. Exhumed fault outcrops and postglacial geological evidences including radiocarbon-dated Holocene sediments affected by faults are in agreement with a transpressive regime associated to a NE-SW shor- tening, that defined deformation along right lateral strike-slip N-S and locally NNW-striking faults with reverse component, which is in turn consistent with the development of pop-up structures at surface.

Together with local extension along dominantly dextral SW-NE to NNE-striking faults, we inter- pret that this kinematics has controlled the origin of duplex structures favouring the development of pull-apart basins at Quaternary time scale, like the Yulton and Meullín Lakes, which geomorphologic evolution have also been driven by strong glacial erosion overprinted by fluvial, volcanic and lands- lide processes. In this context, the cartography of the Holocene monogenetic volcanic cones shows a close spatial association of most of these centres with the main fault systems in the area, especially with N-S and NE-striking faults.

Together with their morphometric characterization, this is consistent with a tectonic control for the origin of a part of them. Calibrated radiocarbon results from carbonized tree remnants under a lava flow associated to one of the most recent monogenetic volcanic centres yielded ages between the yearsAD 1450 and 1670, indicating that this volcanism occurred even during historic times in the area. From the analysis of geomorphologic features together with calibrated radiocarbon results we estimate a minimum Holocene uplift rate in the order of 1.1-1.4 mm/yr in the eastern block located along the dextral N-S Río Cuervo Fault system, which is one of the major structures in the area.

Some systematically deflected drainage networks observed from satellite imagery provided prelimi- nary evidence for recent right lateral displacement, suggesting significant Holocene dextral strike slip along the main NE-striking Riesco Fault system. Submarine fault scarps from sub-bottom profile data and surface neotectonics reveal that the large Mw 6.2 earthquake that affected the Aysén Fjord during 2007 was associated with a subvertical-dipping west 15-20 km long dextral strike-slip north-south structure, named Punta Cola Fault. This fault is located along an area that can be taken as a main north-south branch of the LOFZ, 6-8 km to the east of the major north-south Río Cuervo Fault system.

Most of the epicentres of the seismic swarm that occurred between January and April 2007 were located along or in between these structures. The cartography of active faults proposed in this work suggests that this area is a transference zone between main regional branches of the LOFZ, which implies potential seismic hazard in this portion of the Patagonian Andes, including the city of Puerto Aysén particularly associated to the active Los Palos Fault system proposed here. Our results, together with previous observations and scale relationships between active faults and earthquakes, suggest that large seismic events Mw 6.2-6.5 can be typically expected along most of the active structures, while earthquake magnitude in the order of 7.1 could be expected along the main fault systems like the Río Cuervo Fault.

Acknowledgements This work was supported by Energía Austral. The authors gratefully acknowledge additional support from Millennium Nucleus on Seismotectonics and Seismic Hazard, International Earthquake Research Centre, Montessus de Ballore (CIIT-MB), Grant P06- 064-F, and CEGA(Andean Geothermal Centre of Excelence) Fondap Project 15090013. 3.5 kHz sub-bottom profiler and ma- rine geology facilities were provided by the Department of Geology, University of Chile. Fieldwork and applied scientific research were greatly assisted by C. Nuñez, S. Moreira and A. Maurín from Energía Austral, J. Castillo and L.

Cares from Knight Piésold, together with crucial logistic support from J. Gangas, M. Contreras and field staff from Centro de Ecología Aplicada (CEA). We also thank the crew of the R/V Petrel of CONAF.A.Amigo, G. Orozco, J.Arriagada and T. Villaseñor provided additional support during fieldwork and/or data processing. We gratefully acknowledge reviews from L. Lara, one anonymous reviewer, and the editor M. Suárez. Also, we thank useful comments to a first version of this manus- cript from J. Cembrano. We thank M. Gardeweg and C. Lamperein who provided valuable comments with respect to the first results during the realization of the applied scientific research.

References Adriasola,A.;Thomson,S.;Brix,M.;Hervé,F.;Stöckhert, B. 2006. Postmagmatic cooling and late Cenozoic de- nudation of the North Patagonian Batholith in the Los

2003. The Quaternary calc-alkaline volcanism of thePatagonianAndesclosetotheChiletriplejunction: geochemistry and petrogenesis of volcanic rocks from theCayandMacavolcanoes,Chile(~45ºS).Journalof South American Earth Sciences 16: 219-242.