• Keine Ergebnisse gefunden

Corresponding author. Tel.: +49 511 643 2149; fax: +49 511 643 3663.

E-mail address:kai.berglar@bgr.de(K. Berglar).

0025-3227/$see front matter © 2008 Elsevier B.V. All rights reserved.

doi:10.1016/j.margeo.2008.04.006

Contents lists available atScienceDirect

Marine Geology

j o u r n a l h o m e p a g e : w w w. e l s e v i e r. c o m / l o c a t e / m a r g e o

(DeMets et al., 1994), using geomagnetic time scale, gives values of 73 mm/y (N19°E) and 61 mm/y (N15°E) respectively.

Due to the curvature of the margin the plate convergence gradually changes from nearly perpendicular subduction off Java to highly oblique subduction off northern Sumatra (Moore et al., 1980). This results in slip partitioning along the northwestern Sunda island arc and the development of arc-parallel strike-slip faults. The most pro-nounced strike-slip shear zone is the Sumatran Fault Zone (McCarthy and Elders, 1997; Sieh and Natawidjaja, 2000) located on the Sumatran mainland along the volcanic arc which forms the Barisan Mountains (Fig. 1). The Sumatran Fault Zone accommodates most of the right-lateral stress of the relative plate motion. However, a smaller amount of the arc-parallel stress is taken up by right-lateral strike-slip fault systems along the western edges of the forearc basins, namely the Mentawai Fault Zone and West Andaman Fault (Diament et al., 1992;

Malod and Kemal, 1996; McCaffrey, 1991). The Mentawai Fault Zone extends from the Sunda Strait up to south of Nias Island where it is probably connected with the Sumatran Fault Zone along the Batee Fault (Milsom et al., 1995). The West Andaman Fault proceeds southwards from the Andaman Islands to the Simeulue basin (Curray, 2005). As pointed out byCurray et al. (1979)the Sumatran forearc acts as a sliver plate bounded to the west by the trench, below by the subducting plate, and to the east by the Sumatran Fault Zone. As a

consequence the forearc sliver consists of elongated strips moving to the northwest. This was further rened byMalod and Kemal (1996) proposing two forearc microplates between the outer arc high and the Mentawai Fault Zone, separated by the Batee Fault. The western border of the northern microplate is represented by the West Andaman Fault.

Besides strain partitioning, subducting ridges on the oceanic plate namely the Wharton Ridge and the Investigator Fracture Zonemay have had an effect on the deformation of the forearc and may have accounted for segmentation into several forearc basins off northern and central Sumatra.

3. Data acquisition and processing

Data used in this study were acquired during two cruises with the German research vessel SONNE in 2006 (SO186-2 SeaCause and SO189-1 SUMATRA). Approximately 1500 km of multi-channel seismic (MCS) data were obtained in the area of the Simeulue basin (Fig. 2). Shot distance was 50 m and we used a digital 240 channel streamer of 3 km length with a receiver spacing of 12.5 m, towed at a water depth of 6 m. The acoustic signal was generated by a tuned G-gun array of 16 units with a total volume of 50.8 l operating at air pressure of 14.5 MPa. Data were recorded with a sampling interval of 2 ms and 14 s length. Stacking velocities were picked at regular Fig. 1.Regional tectonic setting of the Sumatran subduction zone. Sumatran Fault Zone (SFZ), Mentawai Fault Zone (MFZ), Batee Fault (BF), West Andaman Fault (WAF) and deformation front are based onSieh and Natawidjaja (2000). Ages of the oceanic crust are afterMüller et al. (1997)andDeplus et al. (1998). Gray arrows indicate relative plate movements based on NUVEL-1A (DeMets et al., 1994), black arrows based on CGPS (Prawirodirdjo and Bock, 2004).

2 K. Berglar et al. / Marine Geology 253 (2008) 1–13

intervals of 3 km along every line. Pre-stack processing included resampling to 4 ms, trace editing, CMP-sort (nominal 30-fold cover-age, 6.25 m spacing), Ormsby bandpass lter (61260160 Hz), polygon fklter (window of 60 traces and 1 s length), zerophase spiking deconvolution (52 ms operator length, 1 s design window beginning shortly below seabottom reection), amplitude correction for spherical divergence based on stacking velocities (1 / (t×v2)), normal moveout correction (40% stretch mute), and Radon velocity lter for multiple suppression (rejecting velocities differing more than

± 20% of corresponding stacking velocity). After stack we applied a space and time variant Ormsby bandpasslter (upper window: 10 2060100 Hz, lower window: 61250100 Hz), a minimum phase predictive deconvolution and a post-stack Kirchhoff time migration with 90% of stacking velocities. From the interpreted 2-D proles grids were calculated using the Schlumberger-GeoQuest CPS-3 convergent gridding algorithm for the major stratigraphic boundaries (spacing:

200 m) and used for generation of TWT-isopach maps. Such maps may be of limited use for estimation of sedimentation rates as they do not consider compaction, erosion or density of the sediments. Also the accuracy of the grids is constricted by limited density of the data.

Nevertheless, they provide insights into the three-dimensional setting of the study area. In addition high resolution parametric echosounder recordings (difference frequency of 3.5 kHz) were evaluated. The data were resampled to 8 kHz, bandpassltered (1.752.13.84 kHz) and the envelope seismic attribute applied for visualization.

4. Chronostratigraphic calibration

Chronostratigraphic calibration of the data was accomplished using three wells located on the northwestern inner shelf of the

Simeulue basin (Fig. 2), drilled by Union Oil Company of California in the 1970s (Tuba-1A, Meulaboh-1 and Tripa-1). These wells were reviewed by several authors (Beaudry and Moore, 1985; Karig et al., 1979; Rose, 1983) with partially contradicting interpretations (see e.g.

Meulaboh-1 well; Fig. 3, upper section). The borehole had a total depth of 3072 m and the lowermost section was described by all authors of pre-Neogene age, ranging from Eocene to Oligocene.

Recrystallized belemnites may also indicate Jurassic strata, but it is more likely they have been reworked into the Paleogene (Rose, 1983).

The latter author described steeply dipping sediments in the Eocene to Oligocene which are overlain unconformably by the successive strata of Miocene age. The basal Miocene strata were described as a thin sequence of nearshore marine and non-marine clastics containing coal beds. Sparse fossil recovery led to a tentative age of Mid Miocene or possibly latest Early Miocene according toRose (1983). Above, a carbonate complex that is overlain by mudstone or shale was described by all three authors. This succession was reported to have an Early Miocene age byBeaudry and Moore (1985)and byKarig et al.

(1979), whereasRose (1983)interpreted a Mid Miocene age. These mudstones are unconformably (Beaudry and Moore, 1985) overlain by another carbonate complex of reportedly Late Miocene (Beaudry and Moore, 1985; Rose, 1983) or Pliocene age (Karig et al., 1979). Further differentiation was done byBeaudry and Moore (1985)describing an additional unconformity between the Pliocene and Pleistocene.

For a rough correlation of the seismic image with the well data the migrated time section was depth-stretched using interval velocities derived from seismic processing. The result shows a reasonable match between the marker horizons and the unconformities in the seismic image with those reported from the well (Fig. 3, upper section). In particular, we were able to correlate two distinct groups of high-Fig. 2.Bathymetric map of the Simeulue basin (Smith and Sandwell, 1997). The dotted gray line indicates the extent of the Simeulue basin. Dashed lines show the location of acquired seismic sections, solid bold lines the location ofFigs. 3 to 8 and 10. Crosses denote positions of wells used for chronostratigraphic calibration. Stars display the epicenter locations of the December 2004 (north of Simeulue) and March 2005 (Banyak Islands) mega-earthquakes. The sketch shows earthquake epicenter locations (mbN4) since 2000. All earthquake locations are from USGS NEIC.

3 K. Berglar et al. / Marine Geology 253 (2008) 1–13

amplitude and low-frequency reectors, interpreted as carbonate complexes, with the reported carbonates from the well. The interpretation of carbonates in the seismic image is further supported by the mound-like shape of the upper reectors and a velocity pull-up below, resulting from higher travel times in the carbonate rocks with respect to those of the surrounding sedimentary strata. The reectors above the upper carbonate complex are dipping slightly southeast-wards and may indicate an aggrading/prograding system. Atop, a group of parallel reectors give a good match to the Pleistocene sediments described byBeaudry and Moore (1985).

All available wells are located at the inner shelf of the basin in shallow water below 100 m depth. However, it was possible to unequivocally trace the unconformities and their correlating confor-mities into the deeper basin parts (Fig. 3, lower section). The boundary between pre-Neogene and Miocene rocks was identied in unison by all groups interpreting the well data (Beaudry and Moore, 1985; Karig

et al., 1979; Rose, 1983). It is obviously an unconformity of regional extent, as it was also found in forearc basins off central and southern Sumatra and south of Java (Karig et al., 1980; Schlüter et al., 2002;

Susilohadi et al., 2005; van der Werff, 1996). We correlated this unconformity with a distinct high-amplitude reector that also marks a signicant change in seismic facies and was found to be present in all our MCS sections across the forearc basin.

The Mid-Late Miocene boundary is equally distinct in our seismic data. At the eastern basin edge it is developed as an erosional surface with overlying build-up structures indicating subaerial exposure and shallow water environment respectively. In deep water regions a exured surface with thick onlapping strata above can be observed, denoting areas of subsidence and rapid sedimentaryll.

The boundary between the Pliocene and the Pleistocene is clearly visible in the area of the shallow inner shelf where it is characterized by parallel reectors overlying basinward-dipping Fig. 3.Chronostratigraphic calibration of seismic lines. The upper part of thefigure depicts available interpretations of Meulaboh-1 well and a depth-stretched seismic section covering the well location. Mismatch between water depths given in bathymetric (Fig. 2) and well data are due to limitation of bathymetry derived from satellite altimetry in shallow water. The lower portion of thegure illustrates the continuation of unconformities into the deeper basin part and the major units identied. SeeFig. 2for location of MCS prole.

4 K. Berglar et al. / Marine Geology 253 (2008) 1–13

prograding facies. It can be traced into the deep water basin parts where it correlates with young onlapping strata indicating the latest phase of subsidence.

5. Seismic stratigraphy

In the MCS data we identied three major unconformities in the Simeulue basin. By correlating these with industry wells (see Section 4) we dated them as Early Miocene, base Late Miocene and base Pleis-tocene age. Based on these unconformities we dened four major stratigraphic units, which we further divided into subunits. The subunits do not necessarily represent different ages of the strata, but account for diverse depositional conditions. The main reason for this approach was to achieve a clear differentiation between shallow water, platform and deep water environments.

5.1. Pre-Neogene Unit 1

The acoustic basement is represented by the pre-Neogene succes-sion. Its internal geometry is not clearly imaged in the MCS data, but along some proles reectivity can be observed in its uppermost parts.

Especially in undeformed areas a band of parallel low-frequency high-amplitude reectors is visible. Here, these reections are concordant to their upper boundary and overlie an angular unconformity above steeply dipping older strata (Fig. 4, km 20, km 90120;Fig. 5, km 210 220;Fig. 6, km 45 and km 65). Chaotic, hummocky and wavy reectors seem to be connected to later subsidence and deformation, so this internal conguration is presumably rather an expression of deforma-tion than of a sedimentary process.

5.2. Early and Mid Miocene Unit 2

The seismic facies of the Early and Mid Miocene succession reects differentiated basin evolution in the Simeulue area. The following subunits have been identied:

5.2.1. Unit 2a

Unit 2a consists of parallel to sub-parallel high-amplitude reectors of high continuity, concordantly overlying the basal Neogene uncon-formity. The thickness of this unit varies from 100 ms TWT (Fig. 4) to about 300 ms TWT (Fig. 7). Small buildup-like structures occasionally formed atop (Fig. 7, km 4550;Fig. 8, km 65) indicating that Unit 2a is restricted to areas under shallow water conditions. The eastern slope of the basin was located approximately 3040 km westwards during the time of deposition compared to its present day position (Figs. 6 to 8). At the edge of the former shallow water areas reef complexes terminate this unit (Fig. 6, km 45;Fig. 7, km 35) and lead over to deep water environments. Existing paleomorphology is leveled out by Unit 2a.

5.2.2. Unit 2b

The internal reection pattern of Unit 2b is composed of parallel to sub-parallel continuous reectors which are onlapping on their lower boundary. Occurrence of Unit 2b is restricted to the area west of the shallow depositional environment of Unit 2a indicating the location of a local depocenter which is narrowing northwards. Maximum thickness of up to 3 s TWT can be found in the central part of the basin (Fig. 5, km 90). The internal conguration exhibits several onlapping sequences showing phases of increased subsidence. Deformation of Unit 2b can be observed throughout the Simeulue basin and is caused by subsequent Fig. 4.MCS prole covering the southeastern part of the Simeulue basin. A large prograding body (Unit 2c) is covered by sediments derived from the opposite direction (Unit 3b), presumably from the rapidly emerging Banyak Island located at southern border of the Simeulue basin. SeeFig. 2for location of prole.

5 K. Berglar et al. / Marine Geology 253 (2008) 1–13

Fig. 5.MCS prole along the main axis of the Simeulue basin. Early and Mid Miocene deformation is reected in the half graben faulting of Unit 1. The southeastern section is located at the Late Miocene/Pliocene location of the inner shelf of the basin where a carbonate platform evolved (Unit 3a). Recent uplift and associated normal faulting and erosion might be caused by reactivation of the Early/Mid Miocene half grabens. SeeFig. 2for location of prole.

6K.Berglaretal./MarineGeology253(2008)113

subsidence and development of normal faults (e.g.Fig. 5, km 90100;

Fig. 6, km 2028;Fig. 7, km 25).

5.2.3. Unit 2c

Unit 2c overlies the shallow water sediments of Unit 2a. In the southern part of the Simeulue basin it is characterized by basinward prograding, downlapping reection patterns. This unit reects a large sedimentary body of up to 1 s TWT thickness with internal sigmoidal conguration in the shallow area of the basin exhibiting a southward direction of progradation (Figs. 4 and 7). In the southernmost part of the basin Unit 2c interngers with Unit 2b (Fig. 5, km 210). In the central (Fig. 6) and northern (Fig. 8) part of the basin the reectors of Unit 2c are sub-parallel and of good to moderate continuity.

5.3. Late Miocene and Pliocene Unit 3

The Late Miocene and the Pliocene are characterized by an ex-ceedingly expanded depocenter, the development of a large platform with patch reefs in the eastern and southern Simeulue basin during the Late Miocene and a high sedimentary input from the onshore regions.

The shelf of the basin evolved into its present form during this time with continuously prograding and aggrading sequences.

5.3.1. Unit 3a

Unit 3a is characterized by isolated mounded bodies with low-amplitude internal reectors. Below these pinnacles velocity pull-ups are present. Their occurrence is focused on the shallow areas of the eastern and southern Simeulue basin. This mound-like bodies inter-finger with the surrounding units differently. Well stratied (Fig. 5, km 130180), hummocky (Fig. 4, km 3050) or chaotic reection patterns (Fig. 7, km 4050) on the former inner shelf of the basin indicate locally diverse sedimentation environments. The maximum thickness of the patch reefs is about 0.5 s TWT (Fig. 7, km 38).

5.3.2. Unit 3b

Unit 3b is composed of downlapping strata with an internal pattern of parallel to sub-parallel reectors of good continuity. This unit is restricted to the southern part of the Simeulue basin where it reaches a thickness of up to 1 s TWT (Figs. 4 and 5). The northward direction of progradation suggests the Banyak Islands as the source area. With continued progradation the patch reefs of Unit 3a were buried (Fig. 5).

5.3.3. Unit 3c

Analog to Unit 2b this Late Miocene/Pliocene unit is an onlapping basinll facies restricted to the western basin part where strong Fig. 6.MCS profile across the central part of the Simeulue basin perpendicular to the main axis of the basin. Neogene sediment thickness of the Simeulue basin is maximal at its central-western border with 4 s TWT. Wrench faulting results in uplift of this area since the Late Miocene/Pliocene and the paleo-shelf break gets buried (Unit 3d) and subsidence expanded eastwards. SeeFig. 2for location of profile.

7 K. Berglar et al. / Marine Geology 253 (2008) 1–13

subsidence occurred resulting in a sedimentary package of up to 1.2 s TWT thickness. It consists of well stratied reectors with good con-tinuity in the southern basin part (Figs. 6 and 7). To the north, the reections are increasingly discontinuous (Fig. 8).

5.3.4. Unit 3d

Unit 3d is the upper succession of Unit 3c, but spreads over the entire basin with a maximum thickness of 0.4 s TWT. The patch reefs of Unit 3a were completely buried by this sedimentary unit. In the northern (Fig. 8) and central (Fig. 6) Simeulue basin the upper part of this unit is characterized by a chaotic reection pattern merged with stratied reectors in the direction of the inner shelf.

5.4. Pleistocene to recent Unit 4

The Pleistocene to recent sediments are characterized by a well stratied onlapping conguration. They reect the present-day depo-center which shows a northward narrowing triangular form. Recent sedimentation in the Simeulue basin is conned to a large area in the southern basin part (thickness of 0.6 s TWT,Fig. 7) and a narrow region in the north (thickness of 0.4s TWT,Fig. 8), whereas the central basin part is currently subject to inversion, indicated by erosional truncation of the uppermost layers (Fig. 5). The northern depocenter is

distinguished by a westward shift of sedimentation and a bypass of sediments on the slope of the basin (Fig. 8).

6. Discussion of subsidence history and tectonic structures

The subsidence evolution of the Simeulue forearc basin is documented by its sedimentary volume. Analogous to the sedimen-tary units we divide the subsidence evolution into four major stages, characterized by changing accumulation space, sedimentation rates, source areas and tectonic phases.

6.1. Pre-Neogene

The acoustic basement in the seismic sections is formed by pre-Neogene successions. It was penetrated by several wells in the Sumatran forearc basins andpresumably consists of several thousand feet of folded sedimentary rocks deposited at the continental margin of Sundaland(Rose, 1983). Where distinct, our seismic sections show mainly faulted and heavily deformed strata. According tode Smet and Barber (2005)the pre-Neogene unit was affected by a regional horst and graben stage which inuenced the entire Sunda mainland from the Late Eocene to Early Oligocene. In their interpretation the base-ment below the present-day forearc basin was part of a stable Fig. 7.MCS prole across the southern part of the Simeulue basin perpendicular to the main axis of the basin. The eastern basement block rotated after deposition of Unit 2a deforming the lower section of Unit 2b. Shortly after initiation of wrench faulting in Late Miocene/Pliocene time at the western border of the basin (upper part of Unit 3c) subsidence extended several tens of km eastward. SeeFig. 2for location of profile.

8 K. Berglar et al. / Marine Geology 253 (2008) 1–13

continental shelf reaching landwards to the volcanic arc. Other authors proposed that the acoustic basement is made up of sediments of a subduction system of late Mesozoic to early Tertiary age, evi-denced by mélange complexes exposed on the Banyak Islands (Karig et al., 1979) and the outer arc islands of Nias, Pagai and Sipora and by volcanics and intrusives intersecting these sediments (Simandjuntak and Barber, 1996). Paleogene subduction took place at the western Sumatran border according to several plate tectonic reconstructions (Hall, 1997; Longley, 1997; Morley, 2002). Thus these models imply an extensive oceanward migration of the subduction zone to its present location. A combination of both aspects was presented byKopp and Kukowski (2003), suggesting sediments of an older forearc basin overlying rocks originating from the arc framework corresponding to the stable continental shelf fromde Smet and Barber (2005). In this scenario the location of the forearc basin stayed more or less steady

continental shelf reaching landwards to the volcanic arc. Other authors proposed that the acoustic basement is made up of sediments of a subduction system of late Mesozoic to early Tertiary age, evi-denced by mélange complexes exposed on the Banyak Islands (Karig et al., 1979) and the outer arc islands of Nias, Pagai and Sipora and by volcanics and intrusives intersecting these sediments (Simandjuntak and Barber, 1996). Paleogene subduction took place at the western Sumatran border according to several plate tectonic reconstructions (Hall, 1997; Longley, 1997; Morley, 2002). Thus these models imply an extensive oceanward migration of the subduction zone to its present location. A combination of both aspects was presented byKopp and Kukowski (2003), suggesting sediments of an older forearc basin overlying rocks originating from the arc framework corresponding to the stable continental shelf fromde Smet and Barber (2005). In this scenario the location of the forearc basin stayed more or less steady